User expects solubility on the order of 10^-4 Fe in Al

User used point density of 2000. Here it can be seen that 500 (the default) and 10000 give the same answer, but 2000 does not.


In [9]:
from pycalphad import Database, equilibrium, variables as v
db_file = '/Users/brandon/Box Sync/databases/nims-databases/Al-Fe_sei.tdb'
dbf = Database(db_file)
comps = ['AL', 'FE', 'VA']
phases = ['FCC_A1', 'LIQUID']
conditions = {v.X('FE'): .0028156, v.T: 931.385693359375, v.P: 1.0e5}

In [13]:
eq = equilibrium(dbf, comps, phases, conditions, calc_opts={'pdens': 500 }, verbose=True)
eq.X


Calculation Backend: Compiled (autowrap)
Components: AL FE VA
Phases: FCC_A1 LIQUID [done]
('NEW_L_MULTIPLIERS', array([   -405.01747922,    1567.93152121,    7946.1014188 ,
        -37723.1363365 , -148381.49490613]))
('old_driving_force', -38031.778470333891)
(1.0, -38031.752189458384, 0.00080339025389464247)
(0.5, -38032.795578030498, 0.00017453753423712914)
(0.25, -38032.433091489511, 4.363438356069782e-05)
('alpha', 0.25)
('Phases', [CompositionSet(FCC_A1, [  1.00000000e+00   1.00000000e-12]), CompositionSet(LIQUID, [ 0.99649094  0.00350906])])
('step', array([  4.78188768e-12,  -4.78113553e-12,   1.67362776e-16,
         1.99580629e-03,  -1.99580629e-03,  -3.49808568e-01,
         3.49808568e-01]))
('Site fractions', array([  1.00000000e+00,   1.00000000e-12,   1.00000000e+00,
         9.96490936e-01,   3.50906446e-03]))
('Phase fractions', array([ 0.21005566,  0.78994434]))
('Chemical potentials', array([ -37723.1363365 , -148381.49490613]))
('Chem pot progress', array([  6.70741429e-07,  -1.03975750e+03]))
('Energy progress', 4.1738880748744123)
('Driving force', 2.2729193964405567)
('NEW_L_MULTIPLIERS', array([   -242.94450995,     940.50349659,    7783.51850924,
        -37723.13633664, -148147.21465611]))
('old_driving_force', -38032.422868815425)
(1.0, -38032.455768795407, 0.00032375989482125167)
(0.5, -38032.710183718191, 8.3213699967665811e-05)
('alpha', 0.5)
('Phases', [CompositionSet(FCC_A1, [  1.00000000e+00   1.00000000e-12]), CompositionSet(LIQUID, [ 0.99699955  0.00300045])])
('step', array([  2.82438561e-12,  -3.82475055e-12,  -1.80078222e-16,
         1.01722436e-03,  -1.01722436e-03,  -2.41427599e-01,
         2.41427599e-01]))
('Site fractions', array([  1.00000000e+00,   1.00000000e-12,   1.00000000e+00,
         9.96999548e-01,   3.00045228e-03]))
('Phase fractions', array([ 0.08934186,  0.91065814]))
('Chemical potentials', array([ -37723.13633664, -148147.21465611]))
('Chem pot progress', array([ -1.38184987e-07,   2.34280250e+02]))
('Energy progress', 4.083194631995866)
('Driving force', 1.336187833774602)
('NEW_L_MULTIPLIERS', array([    457.39189994,   -1770.68698082,    7081.01753868,
        -37723.13633957, -147659.57556242]))
('old_driving_force', -38032.669605465519)
(1.0, -38032.669441059246, 2.5380530732643791e-07)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  1.00000000e+00   1.82431704e-11]), CompositionSet(LIQUID, [ 0.99700793  0.00299207])])
('step', array([ -1.82427563e-11,   1.72431704e-11,   2.41937939e-16,
         8.38251673e-06,  -8.38251672e-06,  -3.02778704e-02,
         3.02778704e-02]))
('Site fractions', array([  1.00000000e+00,   1.82431704e-11,   1.00000000e+00,
         9.97007930e-01,   2.99206977e-03]))
('Phase fractions', array([ 0.05906399,  0.94093601]))
('Chemical potentials', array([ -37723.13633957, -147659.57556242]))
('Chem pot progress', array([ -2.93365883e-06,   4.87639094e+02]))
('Energy progress', -9.120150984956126)
('Driving force', 0.0039367861609207466)
('NEW_L_MULTIPLIERS', array([    472.53110841,   -1829.29492058,    7065.83163434,
        -37723.13637274, -147658.20186757]))
('old_driving_force', -38032.66944071059)
(1.0, -38032.669497073832, 1.2679815458889465e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  1.00000000e+00   3.03359433e-10]), CompositionSet(LIQUID, [ 0.99700144  0.00299856])])
('step', array([ -2.85115699e-10,   2.85116262e-10,  -1.71573871e-17,
        -6.48626365e-06,   6.48626364e-06,   1.95495797e-03,
        -1.95495797e-03]))
('Site fractions', array([  1.00000000e+00,   3.03359433e-10,   1.00000000e+00,
         9.97001444e-01,   2.99855603e-03]))
('Phase fractions', array([ 0.06101895,  0.93898105]))
('Chemical potentials', array([ -37723.13637274, -147658.20186757]))
('Chem pot progress', array([ -3.31664050e-05,   1.37369485e+00]))
('Energy progress', -0.026564505467831623)
('Driving force', 4.6069362724665552e-05)
('NEW_L_MULTIPLIERS', array([    472.66368287,   -1829.80809564,    7065.69894075,
        -37723.1367257 , -147658.06824193]))
('old_driving_force', -38032.669497072129)
(1.0, -38032.669515699556, 1.099120794378905e-12)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99999996e-01   4.04542765e-09]), CompositionSet(LIQUID, [ 0.99700138  0.00299862])])
('step', array([ -3.74206827e-09,   3.74206822e-09,  -2.93691911e-16,
        -6.79238063e-08,   6.79238081e-08,   1.71174709e-05,
        -1.71174709e-05]))
('Site fractions', array([  9.99999996e-01,   4.04542765e-09,   1.00000000e+00,
         9.97001376e-01,   2.99862395e-03]))
('Phase fractions', array([ 0.06103606,  0.93896394]))
('Chemical potentials', array([ -37723.1367257 , -147658.06824193]))
('Chem pot progress', array([-0.00035297,  0.13362564]))
('Energy progress', -0.0014124613444437273)
('Driving force', 3.1818635761737823e-06)
('NEW_L_MULTIPLIERS', array([    473.53413107,   -1833.17743017,    7064.82886815,
        -37723.13965201, -147657.09528569]))
('old_driving_force', -38032.669515699556)
(1.0, -38032.669667386377, 3.5612512938598684e-11)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99999956e-01   4.35273916e-08]), CompositionSet(LIQUID, [ 0.99700102  0.00299898])])
('step', array([ -3.94819641e-08,   3.94819639e-08,   4.92309392e-16,
        -3.56345511e-07,   3.56345511e-07,   1.12386471e-04,
        -1.12386471e-04]))
('Site fractions', array([  9.99999956e-01,   4.35273916e-08,   1.00000000e+00,
         9.97001020e-01,   2.99898030e-03]))
('Phase fractions', array([ 0.06114845,  0.93885155]))
('Chemical potentials', array([ -37723.13965201, -147657.09528569]))
('Chem pot progress', array([-0.00292631,  0.97295623]))
('Energy progress', -0.00014789276610827073)
('Driving force', 3.0107679776847363e-05)
('NEW_L_MULTIPLIERS', array([    478.90069935,   -1853.95095481,    7059.46487943,
        -37723.15765503, -147651.11029913]))
('old_driving_force', -38032.669667386166)
(1.0, -38032.670584166095, 1.2943292939174289e-09)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99999632e-01   3.67813173e-07]), CompositionSet(LIQUID, [ 0.99699883  0.00300117])])
('step', array([ -3.24285781e-07,   3.24285782e-07,  -6.07449609e-17,
        -2.19224959e-06,   2.19224959e-06,   6.92909202e-04,
        -6.92909202e-04]))
('Site fractions', array([  9.99999632e-01,   3.67813173e-07,   1.00000000e+00,
         9.96998827e-01,   3.00117255e-03]))
('Phase fractions', array([ 0.06184136,  0.93815864]))
('Chemical potentials', array([ -37723.15765503, -147651.11029913]))
('Chem pot progress', array([-0.01800302,  5.98498656]))
('Energy progress', -0.0007784115296090022)
('Driving force', 0.00021433201618492603)
('NEW_L_MULTIPLIERS', array([    502.30168664,   -1944.53848252,    7036.07720966,
        -37723.23594463, -147625.10445358]))
('old_driving_force', -38032.670584132335)
(1.0, -38032.674536630439, 2.2701099111621659e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99997612e-01   2.38767030e-06]), CompositionSet(LIQUID, [ 0.99698929  0.00301071])])
('step', array([ -2.01985712e-06,   2.01985712e-06,   7.65625443e-16,
        -9.53295511e-06,   9.53295511e-06,   3.02153635e-03,
        -3.02153635e-03]))
('Site fractions', array([  9.99997612e-01,   2.38767030e-06,   1.00000000e+00,
         9.96989294e-01,   3.01070550e-03]))
('Phase fractions', array([ 0.0648629,  0.9351371]))
('Chemical potentials', array([ -37723.23594463, -147625.10445358]))
('Chem pot progress', array([ -0.07828959,  26.00584555]))
('Energy progress', -0.001599854207597673)
('Driving force', 0.0011089691979577765)
('NEW_L_MULTIPLIERS', array([    567.82755957,   -2198.2211098 ,    6970.60435865,
        -37723.45454269, -147552.75781008]))
('old_driving_force', -38032.674534983125)
(1.0, -38032.68578875719, 1.4812519200102514e-07)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99988506e-01   1.14936788e-05]), CompositionSet(LIQUID, [ 0.99696268  0.00303732])])
('step', array([ -9.10600854e-06,   9.10600854e-06,  -5.98929581e-16,
        -2.66123223e-05,   2.66123223e-05,   8.46124392e-03,
        -8.46124392e-03]))
('Site fractions', array([  9.99988506e-01,   1.14936788e-05,   1.00000000e+00,
         9.96962682e-01,   3.03731783e-03]))
('Phase fractions', array([ 0.07332414,  0.92667586]))
('Chemical potentials', array([ -37723.45454269, -147552.75781008]))
('Chem pot progress', array([ -0.21859807,  72.34664351]))
('Energy progress', 0.002521466572943609)
('Driving force', 0.0041402154747629538)
('NEW_L_MULTIPLIERS', array([    648.21334639,   -2509.55199414,    6890.34949876,
        -37723.71906602, -147466.25304497]))
('old_driving_force', -38032.685775904472)
(1.0, -38032.70136438754, 8.7100081875256252e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99964716e-01   3.52839033e-05]), CompositionSet(LIQUID, [ 0.9969305  0.0030695])])
('step', array([ -2.37902244e-05,   2.37902244e-05,   6.26605595e-16,
        -3.21792461e-05,   3.21792461e-05,   1.03826269e-02,
        -1.03826269e-02]))
('Site fractions', array([  9.99964716e-01,   3.52839033e-05,   1.00000000e+00,
         9.96930503e-01,   3.06949707e-03]))
('Phase fractions', array([ 0.08370677,  0.91629323]))
('Chemical potentials', array([ -37723.71906602, -147466.25304497]))
('Chem pot progress', array([ -0.26452332,  86.5047651 ]))
('Energy progress', -0.022285533108515665)
('Driving force', 0.0087802996786194853)
('NEW_L_MULTIPLIERS', array([    656.92363807,   -2543.53692444,    6881.76581108,
        -37723.73207313, -147462.48351413]))
('old_driving_force', -38032.701364058077)
(1.0, -38032.707794185248, 2.6258867524830976e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99939920e-01   6.00800007e-05]), CompositionSet(LIQUID, [ 0.99692893  0.00307107])])
('step', array([ -2.47960974e-05,   2.47960974e-05,  -1.79055572e-16,
        -1.57432878e-06,   1.57432878e-06,   1.13078669e-03,
        -1.13078669e-03]))
('Site fractions', array([  9.99939920e-01,   6.00800007e-05,   1.00000000e+00,
         9.96928929e-01,   3.07107140e-03]))
('Phase fractions', array([ 0.08483755,  0.91516245]))
('Chemical potentials', array([ -37723.73207313, -147462.48351413]))
('Chem pot progress', array([-0.01300711,  3.76953085]))
('Energy progress', -0.018869996652938426)
('Driving force', 0.0047075011025299318)
('NEW_L_MULTIPLIERS', array([    646.24656678,   -2502.30988724,    6892.47538694,
        -37723.68742936, -147476.97659578]))
('old_driving_force', -38032.707793803507)
(1.0, -38032.708361645724, 2.0423099827623914e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99930507e-01   6.94927062e-05]), CompositionSet(LIQUID, [ 0.99693436  0.00306564])])
('step', array([ -9.41270557e-06,   9.41270557e-06,  -1.45314495e-16,
         5.42875068e-06,  -5.42875068e-06,  -1.37608463e-03,
         1.37608463e-03]))
('Site fractions', array([  9.99930507e-01,   6.94927062e-05,   1.00000000e+00,
         9.96934357e-01,   3.06564265e-03]))
('Phase fractions', array([ 0.08346147,  0.91653853]))
('Chemical potentials', array([ -37723.68742936, -147476.97659578]))
('Chem pot progress', array([  0.04464377, -14.49308166]))
('Energy progress', 0.004555657156743109)
('Driving force', 0.00042869464232353494)
('NEW_L_MULTIPLIERS', array([    644.78082645,   -2496.64371201,    6893.94315294,
        -37723.6818316 , -147478.80983012]))
('old_driving_force', -38032.708361683275)
(1.0, -38032.708366375038, 3.0574740005409606e-10)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929571e-01   7.04288567e-05]), CompositionSet(LIQUID, [ 0.99693504  0.00306496])])
('step', array([ -9.36150481e-07,   9.36150480e-07,   3.54018812e-16,
         6.81012090e-07,  -6.81012090e-07,  -1.89064111e-04,
         1.89064111e-04]))
('Site fractions', array([  9.99929571e-01,   7.04288567e-05,   1.00000000e+00,
         9.96935038e-01,   3.06496164e-03]))
('Phase fractions', array([ 0.0832724,  0.9167276]))
('Chemical potentials', array([ -37723.6818316 , -147478.80983012]))
('Chem pot progress', array([ 0.00559777, -1.83323434]))
('Energy progress', -0.0022126743278931826)
('Driving force', 3.6132623790763319e-06)
('NEW_L_MULTIPLIERS', array([    644.76695428,   -2496.59007178,    6893.95703319,
        -37723.68178008, -147478.82679363]))
('old_driving_force', -38032.708366375045)
(1.0, -38032.708366375424, 2.622251374373441e-14)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929563e-01   7.04372282e-05]), CompositionSet(LIQUID, [ 0.99693504  0.00306496])])
('step', array([ -8.37149423e-09,   8.37149540e-09,   1.57236134e-16,
         6.27309320e-09,  -6.27309335e-09,  -1.78971179e-06,
         1.78971179e-06]))
('Site fractions', array([  9.99929563e-01,   7.04372282e-05,   1.00000000e+00,
         9.96935045e-01,   3.06495537e-03]))
('Phase fractions', array([ 0.08327061,  0.91672939]))
('Chemical potentials', array([ -37723.68178008, -147478.82679363]))
('Chem pot progress', array([  5.15201245e-05,  -1.69635107e-02]))
('Energy progress', -3.355482476763427e-05)
('Driving force', 2.7648638933897018e-10)
('NEW_L_MULTIPLIERS', array([    644.76695302,   -2496.59006694,    6893.95703438,
        -37723.68178007, -147478.82679498]))
('old_driving_force', -38032.708366375424)
(1.0, -38032.708366375344, 1.3995559268800695e-15)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929563e-01   7.04372289e-05]), CompositionSet(LIQUID, [ 0.99693504  0.00306496])])
('step', array([ -6.69138407e-13,   6.67683197e-13,   1.09537994e-16,
         5.61952758e-13,  -5.61952758e-13,  -1.62212914e-10,
         1.62213912e-10]))
('Site fractions', array([  9.99929563e-01,   7.04372289e-05,   1.00000000e+00,
         9.96935045e-01,   3.06495537e-03]))
('Phase fractions', array([ 0.08327061,  0.91672939]))
('Chemical potentials', array([ -37723.68178007, -147478.82679498]))
('Chem pot progress', array([  4.07453626e-09,  -1.34363654e-06]))
('Energy progress', -2.837623469531536e-09)
('Driving force', 7.2759576141834259e-11)
No progress
Out[13]:
<xarray.DataArray 'X' (P: 1, T: 1, X_FE: 1, vertex: 2, component: 2)>
array([[[[[  9.99929563e-01,   7.04372289e-05],
          [  9.96935045e-01,   3.06495537e-03]]]]])
Coordinates:
  * X_FE       (X_FE) float64 0.002816
  * component  (component) <U2 'AL' 'FE'
  * vertex     (vertex) int64 0 1
  * P          (P) float64 1e+05
  * T          (T) float64 931.4

In [16]:
eq = equilibrium(dbf, comps, phases, conditions, calc_opts={'pdens': 2000 }, verbose=True)
eq.X


Calculation Backend: Compiled (autowrap)
Components: AL FE VA
Phases: FCC_A1 LIQUID [done]
('NEW_L_MULTIPLIERS', array([    472.5425895 ,   -1829.33937232,    7065.81999672,
        -37723.13634075, -147658.20014453]))
('old_driving_force', -38032.669484192571)
(1.0, -38032.669495361733, 2.6818302023556839e-09)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  1.00000000e+00   2.66513395e-11]), CompositionSet(LIQUID, [ 0.99700143  0.00299857])])
('step', array([ -2.56502130e-11,   2.56513395e-11,  -1.44832580e-16,
         2.92931975e-06,  -2.92931976e-06,  -9.15504689e-04,
         9.15504689e-04]))
('Site fractions', array([  1.00000000e+00,   2.66513395e-11,   1.00000000e+00,
         9.97001429e-01,   2.99857143e-03]))
('Phase fractions', array([ 0.06102043,  0.93897957]))
('Chemical potentials', array([ -37723.13634075, -147658.20014453]))
('Chem pot progress', array([ -3.58490070e-06,  -6.61725708e-03]))
('Energy progress', 0.00028365789330564439)
('Driving force', 1.1037191143259406e-05)
No progress
Out[16]:
<xarray.DataArray 'X' (P: 1, T: 1, X_FE: 1, vertex: 2, component: 2)>
array([[[[[  1.00000000e+00,   2.66513395e-11],
          [  9.97001429e-01,   2.99857143e-03]]]]])
Coordinates:
  * X_FE       (X_FE) float64 0.002816
  * component  (component) <U2 'AL' 'FE'
  * vertex     (vertex) int64 0 1
  * P          (P) float64 1e+05
  * T          (T) float64 931.4

In [15]:
eq = equilibrium(dbf, comps, phases, conditions, calc_opts={'pdens': 10000 }, verbose=True)
eq.X


Calculation Backend: Compiled (autowrap)
Components: AL FE VA
Phases: FCC_A1 LIQUID [done]
('NEW_L_MULTIPLIERS', array([    649.21667531,   -2513.79796303,    6889.50458142,
        -37723.69860884, -147473.37738669]))
('old_driving_force', -38032.707168508248)
(1.0, -38032.708334909221, 9.5264208521683713e-08)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99931928e-01   6.80716004e-05]), CompositionSet(LIQUID, [ 0.99693306  0.00306694])])
('step', array([  1.85323370e-05,  -1.85323370e-05,   1.06637900e-16,
        -1.17765120e-05,   1.17765120e-05,   3.14311533e-03,
        -3.14311533e-03]))
('Site fractions', array([  9.99931928e-01,   6.80716004e-05,   1.00000000e+00,
         9.96933056e-01,   3.06694442e-03]))
('Phase fractions', array([ 0.08384473,  0.91615527]))
('Chemical potentials', array([ -37723.69860884, -147473.37738669]))
('Chem pot progress', array([-0.03059329,  9.89886713]))
('Energy progress', 0.0092888152139494196)
('Driving force', 0.0014694959245389327)
('NEW_L_MULTIPLIERS', array([    644.85384905,   -2496.92608661,    6893.8700623 ,
        -37723.6821    , -147478.72098726]))
('old_driving_force', -38032.708335419848)
(1.0, -38032.7083663593, 2.3982501540231738e-09)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929618e-01   7.03820462e-05]), CompositionSet(LIQUID, [ 0.99693501  0.00306499])])
('step', array([ -2.31044584e-06,   2.31044584e-06,  -2.89679208e-16,
         1.95003683e-06,  -1.95003683e-06,  -5.62905503e-04,
         5.62905503e-04]))
('Site fractions', array([  9.99929618e-01,   7.03820462e-05,   1.00000000e+00,
         9.96935006e-01,   3.06499438e-03]))
('Phase fractions', array([ 0.08328183,  0.91671817]))
('Chemical potentials', array([ -37723.6821    , -147478.72098726]))
('Chem pot progress', array([ 0.01650884, -5.34360057]))
('Energy progress', -0.010223446377494838)
('Driving force', 2.1127787476871163e-05)
('NEW_L_MULTIPLIERS', array([    644.76700161,   -2496.59025489,    6893.95698643,
        -37723.68178025, -147478.82673662]))
('old_driving_force', -38032.708366359555)
(1.0, -38032.708366375336, 1.0549690461492123e-12)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929563e-01   7.04372002e-05]), CompositionSet(LIQUID, [ 0.99693504  0.00306496])])
('step', array([ -5.51540171e-08,   5.51540166e-08,   4.71625210e-16,
         3.89973901e-08,  -3.89973899e-08,  -1.12049689e-05,
         1.12049689e-05]))
('Site fractions', array([  9.99929563e-01,   7.04372002e-05,   1.00000000e+00,
         9.96935045e-01,   3.06495539e-03]))
('Phase fractions', array([ 0.08327062,  0.91672938]))
('Chemical potentials', array([ -37723.68178025, -147478.82673662]))
('Chem pot progress', array([ 0.00031975, -0.10574936]))
('Energy progress', -0.00026312009867979214)
('Driving force', 1.2536474969238043e-08)
('NEW_L_MULTIPLIERS', array([    644.76695304,   -2496.59006701,    6893.95703436,
        -37723.68178007, -147478.82679496]))
('old_driving_force', -38032.708366375336)
(1.0, -38032.708366375402, 1.1748414741052926e-15)
('alpha', 1.0)
('Phases', [CompositionSet(FCC_A1, [  9.99929563e-01   7.04372289e-05]), CompositionSet(LIQUID, [ 0.99693504  0.00306496])])
('step', array([ -2.86612571e-11,   2.86616223e-11,  -3.94094232e-16,
         2.19335299e-11,  -2.19325624e-11,  -6.26962270e-09,
         6.26962313e-09]))
('Site fractions', array([  9.99929563e-01,   7.04372289e-05,   1.00000000e+00,
         9.96935045e-01,   3.06495537e-03]))
('Phase fractions', array([ 0.08327061,  0.91672939]))
('Chemical potentials', array([ -37723.68178007, -147478.82679496]))
('Chem pot progress', array([  1.77278707e-07,  -5.83306246e-05]))
('Energy progress', -1.1589872883632779e-07)
('Driving force', -7.2759576141834259e-11)
No progress
Out[15]:
<xarray.DataArray 'X' (P: 1, T: 1, X_FE: 1, vertex: 2, component: 2)>
array([[[[[  9.99929563e-01,   7.04372289e-05],
          [  9.96935045e-01,   3.06495537e-03]]]]])
Coordinates:
  * X_FE       (X_FE) float64 0.002816
  * component  (component) <U2 'AL' 'FE'
  * vertex     (vertex) int64 0 1
  * P          (P) float64 1e+05
  * T          (T) float64 931.4

In [ ]: