Miscibility.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from pycalphad import equilibrium\n",
    "from pycalphad import Database, Model\n",
    "import pycalphad.variables as v"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "ROSE_TDB = \"\"\"\n",
    "ELEMENT H                 TEST    0      0    0  !\n",
    "ELEMENT HE                TEST    0      0    0  !\n",
    "ELEMENT LI                TEST    0      0    0  !\n",
    "ELEMENT BE                TEST    0      0    0  !\n",
    "ELEMENT B                 TEST    0      0    0  !\n",
    "ELEMENT C                 TEST    0      0    0  !\n",
    "ELEMENT N                 TEST    0      0    0  !\n",
    "ELEMENT O                 TEST    0      0    0  !\n",
    "ELEMENT F                 TEST    0      0    0  !\n",
    "ELEMENT NE                TEST    0      0    0  !\n",
    "\n",
    "TYPE_DEFINITION % SEQ * !\n",
    "FUNCTION STR 300 200000; 6000 N !\n",
    "\n",
    " PHASE TEST % 1 1 !\n",
    " CONSTITUENT TEST : H,HE,LI,BE,B,C,N,O,F,NE : !\n",
    "\n",
    "PARAMETER G(TEST,H,HE,LI;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,LI;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,LI;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,BE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,BE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,BE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,HE,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,BE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,BE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,BE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,LI,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,BE,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,B,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,C,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,H,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,BE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,BE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,BE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,LI,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,BE,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,B,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,C,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,HE,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,B;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,B;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,B;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,BE,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,B,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,C,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,LI,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,C;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,C;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,C;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,B,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,C,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,BE,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,N;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,N;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,N;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,C,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,B,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,O;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,O;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,O;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,N,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,C,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,F;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,F;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,F;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,O,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,N,F,NE;2)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,O,F,NE;0)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,O,F,NE;1)     300  STR#;               6000 N !\n",
    "PARAMETER G(TEST,O,F,NE;2)     300  STR#;               6000 N !\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "dbf = Database(ROSE_TDB)\n",
    "my_phases_rose = ['TEST']\n",
    "comps = ['H', 'HE', 'LI', 'BE', 'B', 'C', 'N', 'O', 'F']\n",
    "conds = dict({v.T: 1000, v.P: 101325})\n",
    "for comp in comps[:-1]:\n",
    "    conds[v.X(comp)] = 1.0/float(len(comps))\n",
    "eqx = equilibrium(dbf, comps, my_phases_rose, conds, pbar=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:       (P: 1, T: 1, X_B: 1, X_BE: 1, X_C: 1, X_H: 1, X_HE: 1, X_LI: 1, X_N: 1, X_O: 1, component: 9, internal_dof: 9, vertex: 9)\n",
       "Coordinates:\n",
       "  * P             (P) float64 1.013e+05\n",
       "  * T             (T) float64 1e+03\n",
       "  * X_B           (X_B) float64 0.1111\n",
       "  * X_BE          (X_BE) float64 0.1111\n",
       "  * X_C           (X_C) float64 0.1111\n",
       "  * X_H           (X_H) float64 0.1111\n",
       "  * X_HE          (X_HE) float64 0.1111\n",
       "  * X_LI          (X_LI) float64 0.1111\n",
       "  * X_N           (X_N) float64 0.1111\n",
       "  * X_O           (X_O) float64 0.1111\n",
       "  * vertex        (vertex) int64 0 1 2 3 4 5 6 7 8\n",
       "  * component     (component) object 'B' 'BE' 'C' 'F' 'H' 'HE' 'LI' 'N' 'O'\n",
       "  * internal_dof  (internal_dof) int64 0 1 2 3 4 5 6 7 8\n",
       "Data variables:\n",
       "    MU            (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O, component) float64 -5.835e+03 ...\n",
       "    GM            (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O) float64 -5.835e+03 ...\n",
       "    NP            (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O, vertex) float64 0.2222 ...\n",
       "    X             (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O, vertex, component) float64 0.4956 ...\n",
       "    Y             (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O, vertex, internal_dof) float64 0.4956 ...\n",
       "    Phase         (P, T, X_B, X_BE, X_C, X_H, X_HE, X_LI, X_N, X_O, vertex) object 'TEST' ...\n",
       "Attributes:\n",
       "    solve_iterations: 15\n",
       "    hull_iterations: 1\n",
       "    engine: pycalphad 0.4.1+3.gfbd8ab2\n",
       "    created: 2016-08-09 19:02:13.963575"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "eqx"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'MU' (P: 1, T: 1, X_B: 1, X_BE: 1, X_C: 1, X_H: 1, X_HE: 1, X_LI: 1, X_N: 1, X_O: 1, component: 9)>\n",
       "array([[[[[[[[[[[-5835.09917627, -5835.09917627, -5835.09917627,\n",
       "                 -5835.09917627, -5835.09917627, -5835.09917627,\n",
       "                 -5835.09917627, -5835.09917627, -5835.09917627]]]]]]]]]]])\n",
       "Coordinates:\n",
       "  * P          (P) float64 1.013e+05\n",
       "  * T          (T) float64 1e+03\n",
       "  * X_B        (X_B) float64 0.1111\n",
       "  * X_BE       (X_BE) float64 0.1111\n",
       "  * X_C        (X_C) float64 0.1111\n",
       "  * X_H        (X_H) float64 0.1111\n",
       "  * X_HE       (X_HE) float64 0.1111\n",
       "  * X_LI       (X_LI) float64 0.1111\n",
       "  * X_N        (X_N) float64 0.1111\n",
       "  * X_O        (X_O) float64 0.1111\n",
       "  * component  (component) object 'B' 'BE' 'C' 'F' 'H' 'HE' 'LI' 'N' 'O'"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "eqx.MU"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'X' (P: 1, T: 1, X_B: 1, X_BE: 1, X_C: 1, X_H: 1, X_HE: 1, X_LI: 1, X_N: 1, X_O: 1, vertex: 9, component: 9)>\n",
       "array([[[[[[[[[[[[ 0.4956269 ,  0.4956269 ,  0.00124946,  0.00124946,\n",
       "                   0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.00124946],\n",
       "                 [ 0.00124946,  0.00124946,  0.4956269 ,  0.4956269 ,\n",
       "                   0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.00124946],\n",
       "                 [ 0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.4956269 ,  0.00124946,  0.00124946,  0.4956269 ,\n",
       "                   0.00124946],\n",
       "                 [ 0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.00124946,  0.4956269 ,  0.00124946,  0.4956269 ,\n",
       "                   0.00124946],\n",
       "                 [ 0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.4956269 ,  0.4956269 ,  0.00124946,  0.00124946,\n",
       "                   0.00124946],\n",
       "                 [ 0.00124946,  0.00124946,  0.00124946,  0.00124946,\n",
       "                   0.00124946,  0.00124946,  0.4956269 ,  0.00124946,\n",
       "                   0.4956269 ],\n",
       "                 [        nan,         nan,         nan,         nan,\n",
       "                          nan,         nan,         nan,         nan,\n",
       "                          nan],\n",
       "                 [        nan,         nan,         nan,         nan,\n",
       "                          nan,         nan,         nan,         nan,\n",
       "                          nan],\n",
       "                 [        nan,         nan,         nan,         nan,\n",
       "                          nan,         nan,         nan,         nan,\n",
       "                          nan]]]]]]]]]]]])\n",
       "Coordinates:\n",
       "  * P          (P) float64 1.013e+05\n",
       "  * T          (T) float64 1e+03\n",
       "  * X_B        (X_B) float64 0.1111\n",
       "  * X_BE       (X_BE) float64 0.1111\n",
       "  * X_C        (X_C) float64 0.1111\n",
       "  * X_H        (X_H) float64 0.1111\n",
       "  * X_HE       (X_HE) float64 0.1111\n",
       "  * X_LI       (X_LI) float64 0.1111\n",
       "  * X_N        (X_N) float64 0.1111\n",
       "  * X_O        (X_O) float64 0.1111\n",
       "  * vertex     (vertex) int64 0 1 2 3 4 5 6 7 8\n",
       "  * component  (component) object 'B' 'BE' 'C' 'F' 'H' 'HE' 'LI' 'N' 'O'"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "eqx.X"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'NP' (P: 1, T: 1, X_B: 1, X_BE: 1, X_C: 1, X_H: 1, X_HE: 1, X_LI: 1, X_N: 1, X_O: 1, vertex: 9)>\n",
       "array([[[[[[[[[[[ 0.22222222,  0.22222222,  0.11111111,  0.11111111,\n",
       "                  0.11111111,  0.22222222,         nan,         nan,\n",
       "                         nan]]]]]]]]]]])\n",
       "Coordinates:\n",
       "  * P        (P) float64 1.013e+05\n",
       "  * T        (T) float64 1e+03\n",
       "  * X_B      (X_B) float64 0.1111\n",
       "  * X_BE     (X_BE) float64 0.1111\n",
       "  * X_C      (X_C) float64 0.1111\n",
       "  * X_H      (X_H) float64 0.1111\n",
       "  * X_HE     (X_HE) float64 0.1111\n",
       "  * X_LI     (X_LI) float64 0.1111\n",
       "  * X_N      (X_N) float64 0.1111\n",
       "  * X_O      (X_O) float64 0.1111\n",
       "  * vertex   (vertex) int64 0 1 2 3 4 5 6 7 8"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "eqx.NP"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}