C:\jamie\perplex_f90>build NO is the default () answer to all Y/N prompts Enter name of computational option file to be created, < 100 characters, left justified [default = in]: in6.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)? 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, 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)? n Use chemical potentials, activities or fugacities as independent variables (Y/N)? y Chemical potentials are symbolized in Perple_X, by a lowercase u followed by the component name in parentheses, e.g., the potential of A is written u(A). To convert chemical potentials to activities or fugacities run the Perple_X program MU_2_F after running VERTEX (see on-line documentation for additional information). Select < 3 mobile components from the set: NA2O MGO AL2O3 SIO2 K2O CAO TIO2 MNO FEO O2 Enter names, left justified, 1 per line, to finish: O2 SIO2 Select thermodynamic components from the set: NA2O MGO AL2O3 K2O CAO TIO2 MNO FEO Enter names, left justified, 1 per line, to finish: AL2O3 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 5 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)? 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 - Y(CO2) 4 - u(O2) 5 - u(SIO2) 3 Enter minimum and maximum values, respectively, for: Y(CO2) 0 0.8 Select y-axis variable: 2 - T(K) 3 - P(bars) 4 - u(O2) 5 - u(SIO2) 5 Enter minimum and maximum values, respectively, for: u(SIO2) -928046 -895002 Specify sectioning value for: P(bars) 2000 Specify sectioning value for: u(O2) -291926 Specify sectioning value for: T(K) 863 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]: print6 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]: plot6 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. 3 **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 cons ... 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, 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 GrAd EpCz Calculate high variance phase fields (Y/N)? y Enter calculation title: Test Problem 6 C:\jamie\perplex_f90>vertex Enter computational option file name (i.e. the file created with BUILD), left justified: in6.dat Reading thermodynamic data from file: hp98ver.dat Writing print output to file: print6 Writing plot output to file: plot6 Reading solution models from file: solut.dat cycle 1 1 1 cycle 2 2 3 cycle 3 4 5 cycle 4 6 7 cycle 5 8 9 cycle 6 10 11 cycle 7 12 13 cycle 8 14 15 cycle 9 16 17 cycle 10 18 19 cycle 11 20 21 Initial number of divariant assemblages to be tested is: 21 Testing divariant assemblage 1, 20 assemblages remaining to be tested. finished with equilibrium ( 1) cor cc = GrAd(gr) finished with equilibrium ( 2) cor GrAd(gr) = an finished with equilibrium ( 3) cor cc = an finished with equilibrium ( 4) cc an = GrAd(gr) finished with equilibrium ( 5) cor cc GrAd(gr80) = GrAd(gr89) finished with equilibrium ( 6) cor GrAd(gr89) = an GrAd(gr80) finished with equilibrium ( 3) cor cc = an finished with equilibrium ( 7) cc an GrAd(gr80) = GrAd(gr89) finished with equilibrium ( 8) cc = wo finished with equilibrium ( 9) GrAd(gr) = an wo finished with equilibrium ( 8) cc = wo finished with equilibrium ( 10) cor cc GrAd(gr70) = GrAd(gr80) finished with equilibrium ( 3) cor cc = an finished with equilibrium ( 11) cor GrAd(gr80) = an GrAd(gr70) finished with equilibrium ( 12) cc an GrAd(gr70) = GrAd(gr80) finished with equilibrium ( 13) cor cc GrAd(gr60) = GrAd(gr70) finished with equilibrium ( 3) cor cc = an finished with equilibrium ( 14) cor GrAd(gr70) = an GrAd(gr60) finished with equilibrium ( 15) cc an GrAd(gr60) = GrAd(gr70) finished with equilibrium ( 16) cor cc hem = GrAd(gr60) finished with equilibrium ( 3) cor cc = an finished with equilibrium ( 17) cor GrAd(gr60) = an hem finished with equilibrium ( 18) cc an hem = GrAd(gr60) finished with equilibrium ( 19) cc an GrAd(gr50) = GrAd(gr60) finished with equilibrium ( 20) cc GrAd(gr60) hem = GrAd(gr50) finished with equilibrium ( 21) cc an hem = GrAd(gr50) finished with equilibrium ( 22) GrAd(gr60) hem = an GrAd(gr50) finished with equilibrium ( 23) cc an hem = GrAd(gr40) finished with equilibrium ( 24) cc an GrAd(gr40) = GrAd(gr50) finished with equilibrium ( 25) cc hem GrAd(gr50) = GrAd(gr40) finished with equilibrium ( 26) hem GrAd(gr50) = an GrAd(gr40) finished with equilibrium ( 27) cc an hem = GrAd(gr30) finished with equilibrium ( 28) cc an GrAd(gr30) = GrAd(gr40) finished with equilibrium ( 29) cc hem GrAd(gr40) = GrAd(gr30) finished with equilibrium ( 30) hem GrAd(gr40) = an GrAd(gr30) Testing divariant assemblage 2, 20 assemblages remaining to be tested. Testing divariant assemblage 3, 20 assemblages remaining to be tested. Testing divariant assemblage 4, 20 assemblages remaining to be tested. Testing divariant assemblage 5, 20 assemblages remaining to be tested. Testing divariant assemblage 6, 20 assemblages remaining to be tested. Testing divariant assemblage 7, 20 assemblages remaining to be tested. Testing divariant assemblage 8, 20 assemblages remaining to be tested. Testing divariant assemblage 9, 20 assemblages remaining to be tested. Testing divariant assemblage 10, 20 assemblages remaining to be tested. Testing divariant assemblage 11, 20 assemblages remaining to be tested. Testing divariant assemblage 12, 20 assemblages remaining to be tested. Testing divariant assemblage 13, 20 assemblages remaining to be tested. Testing divariant assemblage 14, 20 assemblages remaining to be tested. Testing divariant assemblage 15, 20 assemblages remaining to be tested. Testing divariant assemblage 16, 20 assemblages remaining to be tested. Testing divariant assemblage 17, 20 assemblages remaining to be tested. finished with equilibrium ( 31) hem GrAd(gr30) cc = GrAd(gr20) finished with equilibrium ( 32) hem GrAd(gr30) = GrAd(gr20) an Testing divariant assemblage 18, 22 assemblages remaining to be tested. Testing divariant assemblage 19, 22 assemblages remaining to be tested. finished with equilibrium ( 33) hem GrAd(gr20) cc = GrAd(gr10) Testing divariant assemblage 20, 22 assemblages remaining to be tested. Testing divariant assemblage 21, 22 assemblages remaining to be tested. finished with equilibrium ( 34) hem cc = GrAd(andr) Testing divariant assemblage 22, 22 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: cor gr89 cc v = 2000.00 863.000 0.625000E-02 -291930. -928050. Testing divariant assemblage 23, 21 assemblages remaining to be tested. Testing divariant assemblage 24, 20 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: cor gr80 cc v = 2000.00 863.000 0.750000E-02 -291930. -928050. Testing divariant assemblage 25, 19 assemblages remaining to be tested. Testing divariant assemblage 26, 18 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: cor gr70 cc v = 2000.00 863.000 0.750000E-02 -291930. -928050. Testing divariant assemblage 27, 17 assemblages remaining to be tested. Testing divariant assemblage 28, 16 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: cor gr60 cc v = 2000.00 863.000 0.750000E-02 -291930. -928050. Testing divariant assemblage 29, 15 assemblages remaining to be tested. Testing divariant assemblage 30, 14 assemblages remaining to be tested. Testing divariant assemblage 31, 15 assemblages remaining to be tested. Testing divariant assemblage 32, 14 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: gr60 hem cc v = 2000.00 863.000 0.750000E-02 -291930. -928050. Testing divariant assemblage 33, 13 assemblages remaining to be tested. Testing divariant assemblage 34, 12 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: hem gr50 cc v = 2000.00 863.000 0.750000E-02 -291930. -928050. Testing divariant assemblage 35, 11 assemblages remaining to be tested. Testing divariant assemblage 36, 10 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: hem gr40 cc v = 2000.00 863.000 0.625000E-02 -291930. -928050. Testing divariant assemblage 37, 9 assemblages remaining to be tested. Testing divariant assemblage 38, 8 assemblages remaining to be tested. Testing divariant assemblage 39, 7 assemblages remaining to be tested. Testing divariant assemblage 40, 7 assemblages remaining to be tested. Testing divariant assemblage 41, 6 assemblages remaining to be tested. Testing divariant assemblage 42, 5 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: hem gr20 cc v = 2000.00 863.000 0.500000E-02 -291930. -928050. Testing divariant assemblage 43, 4 assemblages remaining to be tested. Testing divariant assemblage 44, 3 assemblages remaining to be tested. Metastable assemblage in FLIPIT the assemblage is: hem gr10 cc v = 2000.00 863.000 0.500000E-02 -291930. -928050. Testing divariant assemblage 45, 2 assemblages remaining to be tested. Testing divariant assemblage 46, 2 assemblages remaining to be tested. finished with equilibrium ( 35) cor = and Testing divariant assemblage 47, 3 assemblages remaining to be tested. Testing divariant assemblage 48, 3 assemblages remaining to be tested. Testing divariant assemblage 49, 4 assemblages remaining to be tested. Testing divariant assemblage 50, 3 assemblages remaining to be tested. Testing divariant assemblage 51, 3 assemblages remaining to be tested. Testing divariant assemblage 52, 3 assemblages remaining to be tested. Testing divariant assemblage 53, 4 assemblages remaining to be tested. Testing divariant assemblage 54, 3 assemblages remaining to be tested. Testing divariant assemblage 55, 3 assemblages remaining to be tested. Testing divariant assemblage 56, 3 assemblages remaining to be tested. Testing divariant assemblage 57, 4 assemblages remaining to be tested. Testing divariant assemblage 58, 3 assemblages remaining to be tested. Testing divariant assemblage 59, 3 assemblages remaining to be tested. Testing divariant assemblage 60, 2 assemblages remaining to be tested. Testing divariant assemblage 61, 3 assemblages remaining to be tested. Testing divariant assemblage 62, 2 assemblages remaining to be tested. Testing divariant assemblage 63, 1 assemblages remaining to be tested. Testing divariant assemblage 64, 2 assemblages remaining to be tested. Testing divariant assemblage 65, 1 assemblages remaining to be tested. Testing divariant assemblage 66, 2 assemblages remaining to be tested. Testing divariant assemblage 67, 1 assemblages remaining to be tested. Testing divariant assemblage 68, 1 assemblages remaining to be tested. Testing divariant assemblage 69, 0 assemblages remaining to be tested. Testing divariant assemblage 70, 0 assemblages remaining to be tested. C:\jamie\perplex_f90>psvdraw Enter the VERTEX plot file name: plot6 PostScript will be written to file: plot6.ps Modify the default plot (y/n)? C:\jamie\perplex_f90>