D:\jamie\Berple_X>build


NO is the default (<cr>) answer to all Y/N prompts


Enter name of computational option file to be created,
< 100 characters, left justified [default = in]:
in11.dat

Enter thermodynamic data file name, left justified, [default = hp98ver.dat]:
hp98ver.dat

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

Calculations with a saturated phase (Y/N)?
The phase is:   FLUID
Its compositional variable is: Y(CO2), X(O), etc.
y

Select the independent saturated phase components:
 H2O   CO2
Enter names, left justified, 1 per line, <cr> to finish:

For C-O-H fluids it is only necessary to select volatile species present in
the solids of interest. If the species listed here are H2O and CO2, then to
constrain O2 chemical potential to be consistent with C-O-H fluid speciation
treat O2 as a saturated component. Refer to the Perple_X Tutorial for details.

H2O
CO2


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   O2
Enter names, left justified, 1 per line, <cr> to finish:
O2
SIO2
AL2O3


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

Select thermodynamic components from the set:
 NA2O  MGO   K2O   CAO   TIO2  MNO   FEO
Enter names, left justified, 1 per line, <cr> to finish:
CAO
FEO


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
   9 - Max X(H2O) GCOH fluid Cesare & Connolly 1993
  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 - EoS Birch & Feeblebop (1993)
  15 - X(H2) low T H2-H2O hybrid-EoS
  16 - X(O) H-O HSMRK/MRK hybrid-EoS
  17 - X(O) 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
12

Choose a buffer:
  1 - Pyrite + Pyrrhotite
  2 - Pyrrhotite
  3 - f(S2)
2

Enter atomic Fe/S of pyrrhotite:
0.93

Compute f(H2) & f(O2) as the dependent fugacities
(do not unless you project through carbon) (Y/N)?
n

Reduce graphite activity (Y/N)?
n

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


Specify computational mode:
     1 - Unconstrained minimization [default]
     2 - Constrained minimization on a grid
     3 - Output pseudocompound data

Unconstrained optimization should be used for the calculation of composition,
mixed variable, and Schreinemakers diagrams, it may also be used for the
calculation of phase diagram sections for a fixed bulk composition. Gridded
minimization can be used to construct phase diagram sections for both fixed
and variable bulk composition. Gridded minimization is preferable for the
recovery of phase and bulk properties.

1

Specify number of independent potential variables:
     0 - Composition diagram [default]
     1 - Mixed-variable diagram
     2 - Sections and Schreinemakers-type diagrams
2

Select x-axis variable:
     1 - P(bars)
     2 - T(K)
     3 - X(O)
3

Enter minimum and maximum values, respectively, for: X(O)
0.2 1

Select y-axis variable:
     2 - T(K)
     3 - P(bars)
2

Enter minimum and maximum values, respectively, for: T(K)
753 873

Specify sectioning value for: P(bars)
6000

Constrain bulk composition (as in pseudosections, y/n)?
n

Do you want a print file (Y/N)?
y
Enter the print file name, < 100 characters, left justified [default = pr]:
print11
Long print file format (Y/N)?
n
Write full reaction equations (Y/N)?
n
Suppress console status messages (Y/N)?
n
Print dependent potentials for chemographies (Y/N)?
Answer no if you do not know what this means.
n

Do you want a plot file (Y/N)?
y
Enter the plot file name, < 100 characters, left justified [default = pl]:
plot11

Specify efficiency level [1-5, default = 3]:
  1 - gives lowest efficiency, highest reliability
  5 - gives highest efficiency, lowest reliability
High values increase probability that a curve
may be partially determined or skipped.


**warning ver013** phase iron     has null or negative composition and will be
 rejected from the composition space.
**warning ver013** phase gph      has null or negative composition and will be
 rejected from the composition space.
**warning ver013** phase diam     has null or negative composition and will be
 rejected from the composition space.
**warning ver013** phase CO       has null or negative composition and will be
 rejected from the composition space.
**warning ver013** phase CH4      has null or negative composition and will be
 rejected from the composition space.
**warning ver013** phase H2       has null or negative composition and will be
 rejected from the composition space.

Exclude phases (Y/N)?
n

Do you want to treat solution phases (Y/N)?
y
Enter solution model file name [default = solut.dat]
left justified, < 100 characters:
solut.dat

**warning ver113** F          is not a valid model because component H2O      or CO2
 is constrained or missing

... blah blah ...

**warning ver025** 1 endmembers for MnCtd      The solution will not be considered.

Select phases from the following list, enter 1 per line,
left justified, <cr> to finish

  aChl        Chl         sChl        T           hCrd        TrTs
  TrTsPg      GlTrTs      CzEpPs      EpCz        Opx(HP)     Cpx(l)
  Cpx(h)      GrAd(EW)    GrAd        GrPyAlSp(B  GrPyAlSp(G  GrPyAl(G)
  GtD         Gt(HP)      GrPyAl(B)   trgltsch    Mn-Opx(HP)  MnChl

Gt(HP)
EpCz

Calculate high variance phase fields (Y/N)?
y

Enter calculation title:
Test Problem 11

D:\jamie\Berple_X>vertex

Enter computational option file name (i.e. the file created
with BUILD), left justified:
in11.dat

Reading thermodynamic data from file: hp98ver.dat

Writing print output to file: print11

Writing plot output to file: plot11

Reading solution models from file: solut.dat

**warning ver013** phase iron     has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver013** phase gph      has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver013** phase diam     has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver013** phase CO       has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver013** phase CH4      has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver013** phase H2       has a negative composition and will be
 rejected from the saturated-component composition space.
**warning ver114** the following endmembers are missing for Gt(HP)
  py          spss
**warning ver501** Gt(HP)     will be recast with endmembers:
          111 - alm        211 - gr

 cycle            1           1           1
Initial number of divariant assemblages to be tested is:    1
Testing divariant assemblage     1,     0 assemblages remaining to be tested.
finished with equilibrium (    1) sid = fctd
finished with equilibrium (    2) fctd = fst
finished with equilibrium (    3) sid = fst
finished with equilibrium (    4) ma = an
finished with equilibrium (    4) ma = an
finished with equilibrium (    3) sid = fst
finished with equilibrium (    5) ank = ma sid
finished with equilibrium (    4) ma = an
finished with equilibrium (    6) ank = sid an
finished with equilibrium (    7) sid an = Gt(HP)(al80)
finished with equilibrium (    3) sid = fst
finished with equilibrium (    8) an fst = Gt(HP)(al80)
finished with equilibrium (    9) cc = ma
finished with equilibrium (   10) cc = an
finished with equilibrium (   11) sid an = Gt(HP)(al84)
finished with equilibrium (   12) sid Gt(HP)(al80) = Gt(HP)(al84)
finished with equilibrium (   13) Gt(HP)(al80) = an Gt(HP)(al84)
finished with equilibrium (    3) sid = fst
finished with equilibrium (   14) Gt(HP)(al80) fst = Gt(HP)(al84)
finished with equilibrium (   15) an fst = Gt(HP)(al84)
finished with equilibrium (   16) sid Gt(HP)(al84) = Gt(HP)(al87)
finished with equilibrium (    3) sid = fst
finished with equilibrium (   17) Gt(HP)(al84) fst = Gt(HP)(al87)
finished with equilibrium (   18) sid Gt(HP)(al87) = Gt(HP)(al91)
finished with equilibrium (    3) sid = fst
finished with equilibrium (   19) Gt(HP)(al87) fst = Gt(HP)(al91)
finished with equilibrium (   20) sid Gt(HP)(al91) = Gt(HP)(al94)
finished with equilibrium (    3) sid = fst
finished with equilibrium (   21) Gt(HP)(al91) fst = Gt(HP)(al94)
finished with equilibrium (   22) sid Gt(HP)(al94) = Gt(HP)(al98.0)
finished with equilibrium (    3) sid = fst
finished with equilibrium (   23) Gt(HP)(al94) fst = Gt(HP)(al98.0)
finished with equilibrium (   24) sid = Gt(HP)(alm)
finished with equilibrium (   25) fst = Gt(HP)(alm)
Metastable assemblage in FLIPIT
the assemblage is:     ma       fctd
v = 6000.00      836.449     0.200000      0.00000      0.00000
Testing divariant assemblage     2,     0 assemblages remaining to be tested.
Testing divariant assemblage     3,     1 assemblages remaining to be tested.
Testing divariant assemblage     4,     1 assemblages remaining to be tested.
Testing divariant assemblage     5,     1 assemblages remaining to be tested.
Testing divariant assemblage     6,     2 assemblages remaining to be tested.
Testing divariant assemblage     7,     2 assemblages remaining to be tested.
Testing divariant assemblage     8,     2 assemblages remaining to be tested.
Testing divariant assemblage     9,     3 assemblages remaining to be tested.
Testing divariant assemblage    10,     2 assemblages remaining to be tested.
Testing divariant assemblage    11,     2 assemblages remaining to be tested.
Testing divariant assemblage    12,     3 assemblages remaining to be tested.
Testing divariant assemblage    13,     3 assemblages remaining to be tested.
Testing divariant assemblage    14,     2 assemblages remaining to be tested.
Testing divariant assemblage    15,     3 assemblages remaining to be tested.
Testing divariant assemblage    16,     3 assemblages remaining to be tested.
Metastable assemblage in FLIPIT
the assemblage is:     al84     fst
v = 6000.00      873.000     0.439999      0.00000      0.00000
Testing divariant assemblage    17,     2 assemblages remaining to be tested.
Testing divariant assemblage    18,     3 assemblages remaining to be tested.
Testing divariant assemblage    19,     3 assemblages remaining to be tested.
Testing divariant assemblage    20,     2 assemblages remaining to be tested.
Testing divariant assemblage    21,     2 assemblages remaining to be tested.
Testing divariant assemblage    22,     2 assemblages remaining to be tested.
Testing divariant assemblage    23,     1 assemblages remaining to be tested.
Testing divariant assemblage    24,     1 assemblages remaining to be tested.
Testing divariant assemblage    25,     0 assemblages remaining to be tested.


D:\jamie\Berple_X>psvdraw

Enter the VERTEX plot file name:
plot11

PostScript will be written to file: plot11.ps


Modify the default plot (y/n)?


D:\jamie\Berple_X>