This dialog illustrates how to set up a calculation in which oxygen chemical potential 
(fugacity) is buffered (QFM) as a function of pressure and temperature. See 

perplex.ethz.ch/perplex/faq/how_to_buffer_chemical_potentials.txt

for additional explanation. In the dialog below, after BUILD is run to set up the calculation,
MEEMUM is used to make a test calculation to verify that ferric iron (excess oxygen relative
to FEO is present in the stable minerals).

Ordinarily the user would first verify that the correct buffer is implemented as described 
in the above how_to_buffer... file. 


C:\jamie\Perple_X>build

Perple_X version 6.7.3, source updated Jun 13, 2016.


NO is the default (blank) answer to all Y/N prompts


Enter a name for this project (the name will be used as the
root for all output file names) [default = my_project]:
example29

The problem definition file will be named: example29.dat

Enter thermodynamic data file name [default = hp02ver.dat]:


Enter the computational option file name [default = perplex_option.dat]:
See: www.perplex.ethz.ch/perplex_options.html


Reading Perple_X options from: perplex_option.dat

The current data base components are:
 NA2O  MGO   AL2O3 SIO2  K2O   CAO   TIO2  MNO   FEO   NIO   ZRO2  CL2
 O2    H2O   CO2
Transform them (Y/N)?


Calculations with a saturated   FLUID  (Y/N)?


Calculations with saturated components (Y/N)?
y

**warning ver015** if you select > 1 saturated component, then the order you
enter the components determines the saturation heirarchy and may effect your
results (see Connolly 1990).


Select < 6 saturated components from the set:
 NA2O  MGO   AL2O3 SIO2  K2O   CAO   TIO2  MNO   FEO   NIO   ZRO2  CL2
 O2    H2O   CO2
Enter names, 1 per line, press <enter> to finish:
O2


Use chemical potentials, activities or fugacities as independent
variables (Y/N)?


Select thermodynamic components from the set:
 NA2O  MGO   AL2O3 SIO2  K2O   CAO   TIO2  MNO   FEO   NIO   ZRO2  CL2
 H2O   CO2
Enter names, 1 per line, press <enter> to finish:
SIO2
Al2O3

Al2O3 is invalid. Check spelling, upper/lower case match, and do not use leading blanks. Try again:

AL2O3
FEO
MGO
CAO
NA2O
K2O
H2O


**warning ver016** you are going to treat a saturated (fluid) phase component
as a thermodynamic component, this may not be what you want to do.


Select fluid equation of state:

   0 - X(CO2) Modified Redlich-Kwong (MRK/DeSantis/Holloway)
   1 - X(CO2) Kerrick & Jacobs 1981 (HSMRK)
   2 - X(CO2) Hybrid MRK/HSMRK
   3 - X(CO2) Saxena & Fei 1987 pseudo-virial expansion
   4 - Bottinga & Richet 1981 (CO2 RK)
   5 - X(CO2) Holland & Powell 1991, 1998 (CORK)
   6 - X(CO2) Hybrid Haar et al 1979/CORK (TRKMRK)
   7 - f(O2/CO2)-f(S2) Graphite buffered COHS MRK fluid
   8 - f(O2/CO2)-f(S2) Graphite buffered COHS hybrid-EoS fluid
  10 - X(O) GCOH-fluid hybrid-EoS Connolly & Cesare 1993
  11 - X(O) GCOH-fluid MRK Connolly & Cesare 1993
  12 - X(O)-f(S2) GCOHS-fluid hybrid-EoS Connolly & Cesare 1993
  13 - X(H2) H2-H2O hybrid-EoS
  14 - X(CO2) Pitzer & Sterner 1994; Holland & Powell mixing 2003
  15 - X(H2) low T H2-H2O hybrid-EoS
  16 - X(O) H-O HSMRK/MRK hybrid-EoS
  17 - X(O)-f(S2) H-O-S HSMRK/MRK hybrid-EoS
  18 - X(CO2) Delany/HSMRK/MRK hybrid-EoS, for P > 10 kb
  19 - X(O)-X(S) COHS hybrid-EoS Connolly & Cesare 1993
  20 - X(O)-X(C) COHS hybrid-EoS Connolly & Cesare 1993
  21 - X(CO2) Halbach & Chatterjee 1982, P > 10 kb, hybrid-Eos
  22 - X(CO2) DHCORK, hybrid-Eos
  23 - Toop-Samis Silicate Melt
  24 - f(O2/CO2)-N/C Graphite saturated COHN MRK fluid
  25 - H2O-CO2-NaCl Aranovich et al. 2010
  26 - O-Si Silicate fluid RK EoS
  27 - C-O-H hybrid EoS
5
Specify computational mode:

     1 - Convex-Hull minimization
     2 - Constrained minimization on a 2d grid [default]
     3 - Constrained minimization on a 1d grid
     4 - Output pseudocompound data
     5 - Phase fractionation calculations

Use Convex-Hull minimization for Schreinemakers projections or phase diagrams
with > 2 independent variables. Use constrained minimization for phase diagrams
or phase diagram sections with < 3 independent variables.


The data base has P(bar)   and T(K)     as default independent potentials.
Make one dependent on the other, e.g., as along a geothermal gradient (y/n)?


Select x-axis variable:
     1 - P(bar)
     2 - T(K)
     3 - Composition X(C1)* (user defined)

*X(C1) can not be selected as the y-axis variable

2

Enter minimum and maximum values, respectively, for: T(K)
673 1473

Enter minimum and maximum values, respectively, for: P(bar)
3d4 8d4

For gridded minimization, grid resolution is determined by the number of levels
(grid_levels) and the resolution at the lowest level in the X- and Y-directions (x_nodes
and y_nodes) these parameters are currently set for the exploratory and autorefine cycles
as follows:

stage        grid_levels  xnodes  ynodes    effective resolution
exploratory       1          40      40        40 x  40 nodes
auto-refine       4          60      60       473 x 473 nodes

To change these options edit or create the file perplex_option.dat
See: www.perplex.ethz.ch/perplex_options.html#grid_parameters


All thermodynamic components must be constrained,
constrain saturated components also (Y/N)?


Specify component amounts by weight (Y/N)?


Enter molar  amounts of the components:
 SIO2  AL2O3 FEO   MGO   CAO   NA2O  K2O   H2O
for the bulk composition of interest:
45 15.26 9.84 6.54 12.66 2.03 0.55 2.63

Output a print file (Y/N)?
y

Summary of valid make definitions:

          NA2O  MGO   AL2O3 SIO2  K2O   CAO   FEO   O2    H2O
qfm       0.00  0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00
fbr       0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00  1.00
fchum     0.00  0.00  0.00  4.00  0.00  0.00  9.00  0.00  1.00
fphA      0.00  0.00  0.00  2.00  0.00  0.00  7.00  0.00  3.00
fper      0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00  0.00
fatg      0.00  0.00  0.00 34.00  0.00  0.00 48.00  0.00 31.00
atgts     0.00 44.00  4.00 30.00  0.00  0.00  0.00  0.00 31.00
sil8L     0.00  0.00  1.60  1.60  0.00  0.00  0.00  0.00  0.00
fo8L      0.00  4.00  0.00  2.00  0.00  0.00  0.00  0.00  0.00
fa8L      0.00  0.00  0.00  2.00  0.00  0.00  4.00  0.00  0.00
q8L       0.00  0.00  0.00  4.00  0.00  0.00  0.00  0.00  0.00
mrb       1.00  3.00  0.00  8.00  0.00  0.00  2.00  0.50  1.00
ts_dqf    0.00  3.00  2.00  6.00  0.00  2.00  0.00  0.00  1.00
parg_dqf  0.50  4.00  1.50  6.00  0.00  2.00  0.00  0.00  1.00
gl_dqf    1.00  3.00  1.00  8.00  0.00  0.00  0.00  0.00  1.00
mpa       0.50  6.00  1.50  6.00  0.00  0.00  0.00  0.00  1.00
ged_dqf   0.00  5.00  2.00  6.00  0.00  0.00  0.00  0.00  1.00
ogl_dqf   1.00  3.00  1.00  8.00  0.00  0.00  0.00  0.00  1.00
fanth_dq  0.00  0.00  0.00  8.00  0.00  0.00  7.00  0.00  1.00
omrb_dqf  1.00  3.00  0.00  8.00  0.00  0.00  2.00  0.50  1.00
mpa2      0.50  6.00  1.50  6.00  0.00  0.00  0.00  0.00  1.00
ged_dqf2  0.00  5.00  2.00  6.00  0.00  0.00  0.00  0.00  1.00
fbit      0.00  2.00  1.00  2.00  0.50  0.00  1.00  0.25  1.00
ocar      0.00  0.00  0.50  2.00  0.00  0.00  2.00  0.25  2.00
octd      0.00  0.00  0.75  1.00  0.00  0.00  1.50  0.12  1.00
fmu       0.00  0.00  1.00  3.00  0.50  0.00  1.00  0.25  1.00
ma_dqf    0.00  0.00  2.00  2.00  0.00  1.00  0.00  0.00  1.00
kho       0.00  3.00  0.00  3.00  0.00  0.00  2.00  0.50  0.00
mots      0.00  1.00  0.50  1.00  0.00  0.00  1.00  0.25  0.00
esn       0.00  0.00  0.50  1.00  0.00  1.00  1.00  0.25  0.00




**warning ver013** iron     has null or negative composition and will be
rejected from the composition space.

**warning ver013** H2       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** fesi     has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-diam  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-HCP   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-HCP   has null or negative composition and will be
rejected from the composition space.

Exclude pure and/or endmember phases (Y/N)?
y
Do you want to be prompted for phases (Y/N)?
n
Enter names, 1 per line, press <enter> to finish:
O2
h2oL


Include solution models (Y/N)?
y
Enter the solution model file name [default = solution_model.dat]:



**warning ver114** the following endmembers are missing for Chl(W)
  ff3cli    f3clin    mame      mafchl    f3mcli    mnchl

**warning ver114** the following endmembers are missing for Cpx(JH)
  cfs       crdi      cess      cenjh
**warning ver025** 0 endmembers for Melt(JH)   The solution will not be considered.

**warning ver114** the following endmembers are missing for Sp(JH)
  picr

**warning ver114** the following endmembers are missing for Grt(JH)
  fkno_d    knor      ckno_d

**warning ver114** the following endmembers are missing for Opx(JH)
  odijh     tsdi      oen       ofs       odif      cren      crfs      odicr


**warning ver114** the following endmembers are missing for Gt(W)
  spss      fmn_i     fkho_i    kho1      andr_i

**warning ver114** the following endmembers are missing for Mica(W)
  ma1_dqf

**warning ver114** the following endmembers are missing for Ctd(W)
  ctdo

**warning ver114** the following endmembers are missing for St(W)
  mstt      msto      mnst

**warning ver114** the following endmembers are missing for Opx(W)
  mnots     mots1     feots     modi      odi       fodi      mnts      mnopx


**warning ver114** the following endmembers are missing for Bi(W)
  mfbi      ffbi_i    fbi       mtbi      ftbi_i    tbi       mts       mnbi


**warning ver114** the following endmembers are missing for Sp(WPC)
  usp
**warning ver025** 1 endmembers for Ilm(W)     The solution will not be considered.

**warning ver114** the following endmembers are missing for Crd(W)
  mncrd     hmncrd_i

**warning ver114** the following endmembers are missing for Sa(WP)
  spr5      fsp5_d    ospr      fospr_d

**warning ver114** the following endmembers are missing for Bio(TCC)
  mfbit     mtbit     ftbit     tbit      mts       mnbi

**warning ver114** the following endmembers are missing for Bio(WPH)
  ffbi_i    fbi       ftbi_i    tbi1

**warning ver114** the following endmembers are missing for CrOpx(HP)
  crts      fcrts_d

**warning ver114** the following endmembers are missing for CrGt
  uv_d      fuv_d     knor

**warning ver114** the following endmembers are missing for Gt(MPF)
  uv_d      fuv_d     knor

**warning ver114** the following endmembers are missing for Chl(HP)
  mame_i    mafchl_i  mnchl

**warning ver114** the following endmembers are missing for melt(HP)
  zr8L

**warning ver114** the following endmembers are missing for pMELTS(G)
  tiGL

**warning ver114** the following endmembers are missing for mMELTS(G)
  crGL      tiGL
**warning ver025** 0 endmembers for MELTS(GS)  The solution will not be considered.

**warning ver114** the following endmembers are missing for Sapp(KWP)
  spr5      fsp5_i

**warning ver114** the following endmembers are missing for GaHcSp
  gah
**warning ver025** 1 endmembers for F          The solution will not be considered.
**warning ver025** 1 endmembers for F(salt)    The solution will not be considered.
**warning ver025** 0 endmembers for Scap       The solution will not be considered.

**warning ver114** the following endmembers are missing for St(HP)
  mnst

**warning ver114** the following endmembers are missing for Ctd(HP)
  mnctd

**warning ver114** the following endmembers are missing for hCrd
  mncrd     hmncrd_i
**warning ver025** 1 endmembers for A          The solution will not be considered.

**warning ver114** the following endmembers are missing for Amph(DHP)
  mfets     ffets_i

**warning ver114** the following endmembers are missing for Amph(DPW)
  mfets     ffets_i
**warning ver025** 1 endmembers for GCOHF      The solution will not be considered.
**warning ver025** 1 endmembers for COHF       The solution will not be considered.

**warning ver114** the following endmembers are missing for Mica(CHA1)
  tip       ntip      ctip      ftip      nftip     cftip

**warning ver114** the following endmembers are missing for Mica(CHA)
  tip

**warning ver114** the following endmembers are missing for Mica+(CHA)
  tip

**warning ver114** the following endmembers are missing for O(HP)
  teph
**warning ver025** 0 endmembers for Do(HP)     The solution will not be considered.
**warning ver025** 0 endmembers for M(HP)      The solution will not be considered.
**warning ver025** 0 endmembers for Do(AE)     The solution will not be considered.
**warning ver025** 0 endmembers for Cc(AE)     The solution will not be considered.
**warning ver025** 1 endmembers for Mt(W)      The solution will not be considered.
**warning ver025** 1 endmembers for IlHm(A)    The solution will not be considered.
**warning ver025** 0 endmembers for IlGkPy     The solution will not be considered.

**warning ver114** the following endmembers are missing for MtUl(A)
  usp

**warning ver114** the following endmembers are missing for Gt(B)
  spss

**warning ver114** the following endmembers are missing for Gt(GCT)
  spss

**warning ver114** the following endmembers are missing for Gt(HP)
  spss

**warning ver114** the following endmembers are missing for Gt(WPH)
  spss      fmn_i
**warning ver025** 0 endmembers for casmelt    The solution will not be considered.
**warning ver025** 0 endmembers for Toop-Melt  The solution will not be considered.
**warning ver025** 0 endmembers for oCcM(HP)   The solution will not be considered.
**warning ver025** 0 endmembers for Tour(V)    The solution will not be considered.
**warning ver025** 0 endmembers for DoDo       The solution will not be considered.
**warning ver025** 1 endmembers for Wus(fab)   The solution will not be considered.
**warning ver025** 1 endmembers for Aki(fab)   The solution will not be considered.
**warning ver025** 0 endmembers for Pv(fab)    The solution will not be considered.
**warning ver025** 0 endmembers for Ppv(og)    The solution will not be considered.
**warning ver025** 0 endmembers for Wad(stx)   The solution will not be considered.
**warning ver025** 0 endmembers for Ring(stx)  The solution will not be considered.

**warning ver114** the following endmembers are missing for Gt(stx)
  maj
**warning ver025** 0 endmembers for C2/c(stx)  The solution will not be considered.

**warning ver114** the following endmembers are missing for Cpx(stx)
  mdi

**warning ver114** the following endmembers are missing for Cpx(stx7)
  mdi

**warning ver114** the following endmembers are missing for Cpx(HP)
  ccrts

**warning ver114** the following endmembers are missing for CrSp
  mcrm      fcrm_d

**warning ver114** the following endmembers are missing for Eskol(C)
  esk

**warning ver114** the following endmembers are missing for Ca-Amph(D)
  mfets     ffets_i
**warning ver025** 1 endmembers for Wus(stx7)  The solution will not be considered.
**warning ver025** 1 endmembers for Aki(stx7)  The solution will not be considered.
**warning ver025** 0 endmembers for Pv(stx7)   The solution will not be considered.
**warning ver025** 0 endmembers for Wad(stx7)  The solution will not be considered.
**warning ver025** 0 endmembers for Ring(stx7) The solution will not be considered.

**warning ver114** the following endmembers are missing for Sp(WPC)
  usp
**warning ver025** 1 endmembers for Ilm(WPH)   The solution will not be considered.
**warning ver025** 0 endmembers for Wad(stx8)  The solution will not be considered.
**warning ver025** 0 endmembers for Ring(stx8) The solution will not be considered.

**warning ver114** the following endmembers are missing for Opx(stx8)
  odi

**warning ver114** the following endmembers are missing for Cpx(stx8)
  mdi       cts
**warning ver025** 1 endmembers for Aki(stx8)  The solution will not be considered.

**warning ver114** the following endmembers are missing for Gt(stx8)
  maj       jmaj
**warning ver025** 0 endmembers for Ppv(stx8)  The solution will not be considered.
**warning ver025** 0 endmembers for CF(stx8)   The solution will not be considered.

**warning ver114** the following endmembers are missing for Stlp(M)
  mnsp

**warning ver114** the following endmembers are missing for Carp(M)
  mnca
**warning ver025** 0 endmembers for Carb(M)    The solution will not be considered.

**warning ver114** the following endmembers are missing for Bio(HP)
  mtbi_i    ftbi_i    tbi       mnts_i    mnbi

**warning ver114** the following endmembers are missing for Ctd(SGH)
  mnctd

**warning ver114** the following endmembers are missing for Carp(SGH)
  mncar

**warning ver114** the following endmembers are missing for Gt(WPPH)
  spss      fmn_i
**warning ver025** 0 endmembers for FeSi(BCC)  The solution will not be considered.
**warning ver025** 0 endmembers for FeSi(fcc)  The solution will not be considered.
**warning ver025** 0 endmembers for oFeSi(bcc) The solution will not be considered.
**warning ver025** 0 endmembers for FeSi_liq   The solution will not be considered.
**warning ver025** 0 endmembers for ZrRu       The solution will not be considered.

**warning ver114** the following endmembers are missing for ZrGt(KP)
  zrg
**warning ver025** 0 endmembers for Si-Fluid   The solution will not be considered.

**warning ver114** the following endmembers are missing for Maj
  maj
**warning ver025** 0 endmembers for Wad        The solution will not be considered.
**warning ver025** 0 endmembers for Ring       The solution will not be considered.
**warning ver025** 1 endmembers for Aki        The solution will not be considered.
**warning ver025** 0 endmembers for Pv         The solution will not be considered.
**warning ver025** 0 endmembers for FeSiC_liq  The solution will not be considered.

**warning ver114** the following endmembers are missing for Sapp(TP)
  spr5      fsp5_d    ospr      fospr_d
**warning ver025** 0 endmembers for oCcM(EF)   The solution will not be considered.
**warning ver025** 0 endmembers for LIQ(NK)    The solution will not be considered.
**warning ver025** 0 endmembers for LIQ(EF)    The solution will not be considered.
**warning ver025** 0 endmembers for dis(EF)    The solution will not be considered.
**warning ver025** 1 endmembers for Si-vapor   The solution will not be considered.
**warning ver025** 1 endmembers for CCO-vapor  The solution will not be considered.
**warning ver025** 0 endmembers for ideal_gas  The solution will not be considered.
**warning ver025** 0 endmembers for FeSiC-BCC  The solution will not be considered.
**warning ver025** 0 endmembers for FeSiC-FCC  The solution will not be considered.
**warning ver025** 0 endmembers for FeCr(liq)  The solution will not be considered.
**warning ver025** 0 endmembers for FeCr-FCC   The solution will not be considered.
**warning ver025** 0 endmembers for FeCr-BCC   The solution will not be considered.
**warning ver025** 0 endmembers for FeCr-s     The solution will not be considered.
**warning ver025** 0 endmembers for FeCliq     The solution will not be considered.
**warning ver025** 0 endmembers for FeCfcc     The solution will not be considered.
**warning ver025** 0 endmembers for FeCbcc     The solution will not be considered.

Select models from the following list, enter 1 per line, press <enter> to finish

chlorite models:          Chl(W)     Chl(HP)    Chl(LWV)
epidote models:           Ep(HP11)   Ep(HP)
clinopyroxene models:     Cpx(JH)    Omph(GHP)  Cpx(l)     Cpx(h)     Cpx(stx)
                          Cpx(stx7)  Omph(HP)   Cpx(HP)    Cpx(m)     Cpx(stx8)
                          Omph(GHP2)
olivine models:           O(JH)      O(SG)      O(HP)      O(stx)     O(stx7)
                          Ol(m)      O(stx8)
spinel models:            Sp(JH)     Sp(WPC)    GaHcSp     MF         Sp(JR)
                          Sp(GS)     Sp(HP)     MtUl(A)    Sp(stx)    CrSp
                          Sp(stx7)   Sp(WPC)    Sp(stx8)
binary-feldspar models:   Pl(JH)     Pl(h)      Kf         San        San(TH)
                          Pl(stx8)
garnet models:            Grt(JH)    Gt(W)      CrGt       Gt(MPF)    Gt(B)
                          Gt(GCT)    Gt(HP)     Gt(EWHP)   Gt(WPH)    Gt(stx)
                          Gt(stx8)   Gt(WPPH)   ZrGt(KP)   Maj
orthopyroxene models:     Opx(JH)    Opx(W)     Opx(HP)    CrOpx(HP)  Opx(stx)
                          Opx(stx8)
liquid models:            melt(W)    melt(HP)   pMELTS(G)  mMELTS(G)
white-mica models:        Mica(W)    Pheng(HP)  MaPa       Mica(CF)   Mica(CHA1)
                          Mica(CHA)  Mica+(CHA) Mica(M)    Mica(SGH)
chloritoid models:        Ctd(W)     Ctd(HP)    Ctd(SGH)
staurolite models:        St(W)      St(HP)
biotite models:           Bi(W)      Bio(TCC)   Bio(WPH)   Bio(HP)
cordierite models:        Crd(W)     hCrd
sapphirine models:        Sa(WP)     Sapp(HP)   Sapp(KWP)  Sapp(TP)
orthoamphibole models:    oAmph(DP)  Anth       o-Amph     oAmph(DP2)
clinoamphibole models:    cAmph(DP)  Cumm       Gl         Tr         GlTrTsPg
                          Amph(DHP)  Amph(DPW)  Ca-Amph(D) Na-Amph(D) Act(M)
                          GlTrTsMr   cAmph(DP2)
osumilite models:         Osm(HP)
talc models:              T
carpholite models:        Carp       Carp(M)    Carp(SGH)
sudoite models:           Sud(Livi)  Sud        Sud(M)
ternary-feldspar models:  feldspar   feldspar_B Pl(I1,HP)  Fsp(C1)
feldspathoid models:      Neph(FB)
alphabet-phase models:    A-phase
clinohumite models:       Chum
serpentine models:        Atg(PN)
brucite models:           B
ilmenite models:          Eskol(C)
pumpellyite models:       Pu(M)
stilpnomelane models:     Stlp(M)
wuestite models:          Wus

For details on these models see:www.perplex.ethz.ch/perplex_solution_model_glossary.html
or read the commentary in the solution model file.

Chl(W)
Ep(HP11)
Omph(GHP2)
Sp(WPC)
Pl(h)
Gt(WPH)
Mica(SGH)
Opx(W)
Bio(TCC)
GlTrTsMr
T
Pu(M)
Stlp(M)


Enter calculation title:
buffered fo2 calculation

C:\jamie\Perple_X>meemum

Perple_X version 6.7.3, source updated Jun 13, 2016.

Enter the project name (the name assigned in BUILD) [default = my_project]:
example_29

Reading problem definition from file: example_29.dat
Reading thermodynamic data from file: hp02ver.dat
Reading solution models from file: solution_model.dat
Writing print output to file: example_29.prn
Writing plot output to file: example_29.plt
Writing bulk composition plot output to file: example_29.blk

Reading Perple_X options from: perplex_option.dat
Writing pseudocompound glossary to file: not requested
Writing Perple_X option summary to file: example_29_MEEMUM_options.txt

Perple_X computational option settings for MEEMUM:

    Keyword:               Value:     Permitted values [default]:

  Free energy minimization options:

    final_resolution:
      exploratory stage    0.10E-01   [1e-2], target value, see actual values below
      auto-refine stage    0.25E-03   [2.5e-4], target value, see actual values below
    resolution factor       3         >2 [3]; iteration keyword value 1
    refinement points       4         1->7 [4]; iteration keyword value 2
    solvus_tolerance_II    0.25       0->1 [0.25]
    global_reach_increment  0         >= 0 [0]
    reach_increment_switch on         [on] off all
    zero_mode              0.1E-05    0->1 [1e-6]; < 0 => off
    zero_bulk              0.1E-05    0->1 [1e-6]; < 0 => off

  Solution subdivision options:

    initial_resolution     0.067      0->1 [0.067], 0 => off
    stretch_factor         0.016      >0 [0.0164]
    subdivision_override   off        [off] lin str
    hard_limits            off        [off] on

  Thermodynamic options:

    solvus_tolerance       aut        [aut] or 0->1; aut = automatic, 0 => p=c pseudocompounds, 1 => homogenize
    speciation_tolerance   0.2E-04    0->1 [1e-3]; order-disorder speciation precision
    T_stop (K)                0.0     [0]
    T_melt (K)              873.0     [873]
    order_check            on         off [on]
    approx_alpha           T          [T] F
    Anderson-Gruneisen     F          [F] T
    site_check             T          [T] F
    speciation_max_it       100       [100]
    fd_expansion_factor    2.0        >0 [2.]
    finite_difference_p    0.1D+05    >0 [1d4]; fraction = 0.1D-01 [1d-2]

  Input/Output options:

    logarithmic_p          F          [F] T
    bad_number             NaN        [0.0]
    composition_constant   F          [F] T
    composition_phase      mol        [mol] wt
    composition_system     wt         [wt] mol
    proportions            vol        [vol] wt mol
    melt_is_fluid          F          [F] T
    solution_names         abb        [mod] abb ful
    species_output         T          [T] F
    seismic_output         som        [some] none all
    pause_on_error         T          [T] F
    poisson_test           F          [F] T

  Seismic velocity options:

    bounds                 VRH        [VRH] HS
    vrh/hs_weighting       0.5        [0.5] 0->1
    explicit_bulk_modulus  T          [F] T
    poisson_ratio          on         [on] all off; Poisson ratio = 0.35

  Information file output options:

    option_list_files      T          [F] T; echo computational options
    pseudocompound_file    F          [F] T; echo static pseudocompound compositions

  Worst case (Cartesian) compositional resolution (mol):

    Exploratory stage:    0.372E-02
    Auto-refine stage:    0.207E-03

  Adapative minimization will be done with:

    3 iterations in the exploratory stage
    5 iterations in the auto-refine stage

To change these options see: www.perplex.ethz.ch/perplex_options.html

Context specific options are echoed in: example_29_MEEMUM_options.txt


Summary of valid make definitions:

          NA2O  MGO   AL2O3 SIO2  K2O   CAO   FEO   O2    H2O
qfm       0.00  0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00
fbr       0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00  1.00
fchum     0.00  0.00  0.00  4.00  0.00  0.00  9.00  0.00  1.00
fphA      0.00  0.00  0.00  2.00  0.00  0.00  7.00  0.00  3.00
fper      0.00  0.00  0.00  0.00  0.00  0.00  1.00  0.00  0.00
fatg      0.00  0.00  0.00 34.00  0.00  0.00 48.00  0.00 31.00
atgts     0.00 44.00  4.00 30.00  0.00  0.00  0.00  0.00 31.00
sil8L     0.00  0.00  1.60  1.60  0.00  0.00  0.00  0.00  0.00
fo8L      0.00  4.00  0.00  2.00  0.00  0.00  0.00  0.00  0.00
fa8L      0.00  0.00  0.00  2.00  0.00  0.00  4.00  0.00  0.00
q8L       0.00  0.00  0.00  4.00  0.00  0.00  0.00  0.00  0.00
mrb       1.00  3.00  0.00  8.00  0.00  0.00  2.00  0.50  1.00
ts_dqf    0.00  3.00  2.00  6.00  0.00  2.00  0.00  0.00  1.00
parg_dqf  0.50  4.00  1.50  6.00  0.00  2.00  0.00  0.00  1.00
gl_dqf    1.00  3.00  1.00  8.00  0.00  0.00  0.00  0.00  1.00
mpa       0.50  6.00  1.50  6.00  0.00  0.00  0.00  0.00  1.00
ged_dqf   0.00  5.00  2.00  6.00  0.00  0.00  0.00  0.00  1.00
ogl_dqf   1.00  3.00  1.00  8.00  0.00  0.00  0.00  0.00  1.00
fanth_dq  0.00  0.00  0.00  8.00  0.00  0.00  7.00  0.00  1.00
omrb_dqf  1.00  3.00  0.00  8.00  0.00  0.00  2.00  0.50  1.00
mpa2      0.50  6.00  1.50  6.00  0.00  0.00  0.00  0.00  1.00
ged_dqf2  0.00  5.00  2.00  6.00  0.00  0.00  0.00  0.00  1.00
fbit      0.00  2.00  1.00  2.00  0.50  0.00  1.00  0.25  1.00
ocar      0.00  0.00  0.50  2.00  0.00  0.00  2.00  0.25  2.00
octd      0.00  0.00  0.75  1.00  0.00  0.00  1.50  0.12  1.00
fmu       0.00  0.00  1.00  3.00  0.50  0.00  1.00  0.25  1.00
ma_dqf    0.00  0.00  2.00  2.00  0.00  1.00  0.00  0.00  1.00
kho       0.00  3.00  0.00  3.00  0.00  0.00  2.00  0.50  0.00
mots      0.00  1.00  0.50  1.00  0.00  0.00  1.00  0.25  0.00
esn       0.00  0.00  0.50  1.00  0.00  1.00  1.00  0.25  0.00



**warning ver013** iron     has null or negative composition and will be
rejected from the composition space.

**warning ver013** H2       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** fesi     has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-diam  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-HCP   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-HCP   has null or negative composition and will be
rejected from the composition space.

**warning ver013** iron     has null or negative composition and will be
rejected from the composition space.

**warning ver013** H2       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si       has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-BCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-FCC   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si_LIQ   has null or negative composition and will be
rejected from the composition space.

**warning ver013** fesi     has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-diam  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-CBCC  has null or negative composition and will be
rejected from the composition space.

**warning ver013** Fe-HCP   has null or negative composition and will be
rejected from the composition space.

**warning ver013** Si-HCP   has null or negative composition and will be
rejected from the composition space.

**warning ver114** the following endmembers are missing for Chl(W)
  ff3cli    f3clin    mame      mafchl    f3mcli    mnchl

**warning ver050** reformulating reciprocal solution: Chl(W)     because of missing endmembers.
(reformulation can be controlled explicitly by excluding additional endmembers).

    2176 pseudocompounds generated for: Chl(W)
      16 pseudocompounds generated for: Ep(HP11)

**warning ver114** the following endmembers are missing for Opx(W)
  mnots     mots1     feots     modi      odi       fodi      mnts      mnopx


**warning ver050** reformulating reciprocal solution: Opx(W)     because of missing endmembers.
(reformulation can be controlled explicitly by excluding additional endmembers).

      64 pseudocompounds generated for: Opx(W)

**warning ver114** the following endmembers are missing for Sp(WPC)
  usp
      59 pseudocompounds generated for: Sp(WPC)

**warning ver114** the following endmembers are missing for Bio(TCC)
  mfbit     mtbit     ftbit     tbit      mts       mnbi

**warning ver050** reformulating reciprocal solution: Bio(TCC)   because of missing endmembers.
(reformulation can be controlled explicitly by excluding additional endmembers).

    3088 pseudocompounds generated for: Bio(TCC)
     145 pseudocompounds generated for: T
      16 pseudocompounds generated for: Pl(h)
         a reach_increment of  3 is specified for Pl(h)

**warning ver114** the following endmembers are missing for Gt(WPH)
  spss      fmn_i

**warning ver050** reformulating reciprocal solution: Gt(WPH)    because of missing endmembers.
(reformulation can be controlled explicitly by excluding additional endmembers).

    2176 pseudocompounds generated for: Gt(WPH)
     136 pseudocompounds generated for: Pu(M)

**warning ver114** the following endmembers are missing for Stlp(M)
  mnsp
      16 pseudocompounds generated for: Stlp(M)
    4422 pseudocompounds generated for: GlTrTsMr
   16800 pseudocompounds generated for: Mica(SGH)
     816 pseudocompounds generated for: Omph(GHP2)

Total number of pseudocompounds:   29930

Interactively enter bulk compositions (y/n)?
If you answer no, MEEMUM uses the bulk composition specified in the input file.



Enter (zeroes to quit) T(K)     P(bar)
750 3e4

**warning ver046** temperature ( 750.000     K) is out of range for endmember/phase: nasGL
nasGL    will be destabilized at this condition. In some cases this problem can be corrected by
setting Anderson_Gruneisen to TRUE in the Perple_X option file.

----------------------------------------

Stable phases at:
                             T(K)     =  750.000
                             P(bar)   =  30000.0


Phase Compositions (molar  proportions):
                   wt %      vol %     mol %     mol        H2O      AL2O3    FEO      MGO      CAO      NA2O     K2O      SIO2     O2
 Gt                54.14     50.19     35.63     7.21      0.00000  0.99256  1.29239  0.48644  1.23606  0.00000  0.00000  3.00000  0.00372
 Mica               7.15      8.84      5.47     1.11      1.00000  1.19370  0.03408  0.27358  0.00000  0.00279  0.49721  3.30494  0.00068
 Cpx                1.32      1.41      1.95    0.394      0.00000  0.40851  0.06572  0.11725  0.13400  0.43300  0.00000  2.00000  0.01224
 Cpx               22.70     24.47     32.87     6.65      0.00000  0.26335  0.06944  0.40386  0.44167  0.27917  0.00000  2.00000  0.00791
 ky                10.80     10.60     20.33     4.11      0.00000  1.00000  0.00000  0.00000  0.00000  0.00000  0.00000  1.00000  0.00000
 law                3.88      4.48      3.77    0.762      2.00000  1.00000  0.00000  0.00000  1.00000  0.00000  0.00000  2.00000  0.00000

Phase speciation (molar proportions):

 Gt                alm: 0.42583, py: 0.16215, gr: 0.40458, andr: 0.00744
 Mica              mu: 0.68675, pa: 0.00558, ma_dqf: 0.00000, cel: 0.27358, fmu: 0.00273, fcel: 0.03136
 Cpx               di: 0.08605, jd: 0.78778, acm: 0.04288, hed: 0.00870, om: 0.04630, cfm: 0.01609, jac: 0.01219
 Cpx               di: -.02176, jd: 0.12465, acm: 0.02158, hed: 0.00417, om: 0.78398, cfm: 0.06727, jac: 0.02011

Molar Properties and Density:
                    N(g)        H(J)      S(J/K)      V(J/bar)      Cp(J/K)      Alpha(1/K)  Beta(1/bar)    Cp/Cv     Density(kg/m3)
 Gt                463.34   -0.50124E+07   694.74       11.808       473.91      0.22533E-04  0.58920E-06   1.0164       3924.0
 Mica              398.80   -0.48625E+07   668.16       13.556       478.55      0.29535E-04  0.15931E-05   1.0118       2942.0
 Cpx               206.01   -0.25292E+07   325.55       6.0655       229.79      0.26453E-04  0.77189E-06   1.0183       3396.4
 Cpx               210.61   -0.25800E+07   328.91       6.2413       241.50      0.28466E-04  0.78040E-06   1.0205       3374.4
 ky                162.05   -0.22135E+07   227.15       4.3726       185.89      0.22380E-04  0.62414E-06   1.0144       3706.0
 law               314.24   -0.39821E+07   549.36       9.9765       424.52      0.31889E-04  0.93874E-06   1.0195       3149.8
 System           6167.90   -0.71792E+08   9414.1       169.57       6729.1      0.25063E-04  0.74671E-06   1.0162       3637.4

Seismic Properties:
                 Gruneisen_T       Ks(bar)       Mu(bar)    V0(km/s)     Vp(km/s)     Vs(km/s)     Poisson ratio
 Gt                0.96849      0.17250E+07  0.99312E+06   6.6303       8.8151       5.0308      0.25849
 Mica              0.53134      0.63511E+06  0.36395E+06   4.6463       6.1711       3.5172      0.25943
 Cpx               0.92112      0.13192E+07  0.80065E+06   6.2322       8.3828       4.8553      0.24760
 Cpx               0.96207      0.13077E+07  0.75467E+06   6.2252       8.2809       4.7291      0.25801
 ky                0.85556      0.16252E+07  0.10192E+07   6.6223       8.9735       5.2443      0.24065
 law               0.81387      0.10860E+07  0.50772E+06   5.8719       7.4814       4.0149      0.29776
 System            0.90869      0.14215E+07  0.81856E+06   6.2514       8.3118       4.7439      0.25844

Bulk Composition:

              mol       mol %      wt %
    H2O       2.630       2.78      0.77
    AL2O3    15.260      16.13     25.23
    FEO       9.840      10.40     11.46
    MGO       6.540       6.91      4.27
    CAO      12.660      13.38     11.51
    NA2O      2.030       2.15      2.04
    K2O       0.550       0.58      0.84
    SIO2     45.000      47.57     43.84
    O2        0.085       0.09      0.04

Other Bulk Properties:

 Enthalpy (J/kg) = -.116396E+08
 Specific Enthalpy (J/m3) = -.423377E+11
 Entropy (J/K/kg) =  1526.31
 Specific Entropy (J/K/m3) = 0.555178E+07
 Heat Capacity (J/K/kg) =  1090.99
 Specific Heat Capacity (J/K/m3) = 0.396835E+07


Chemical Potentials (J/mol):
      H2O           AL2O3         FEO           MGO           CAO           NA2O          K2O           SIO2
    -242979.     -0.155614E+07  -273314.      -596020.      -696612.      -683958.      -709863.      -827721.

Variance (c-p+2) =  4

----------------------------------------


Enter (zeroes to quit) T(K)     P(bar)
0 0

C:\jamie\Perple_X>