C:\jamie\Berple_X\NewVersion>build Perple_X version 6.6.5.9, compiled 6/20/2011. 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]: in24 The problem definition file will be named: in24.dat Enter thermodynamic data file name [default = hp02ver.dat]: sfo05ver.dat Enter the computational option file name [default = perplex_option.dat]: See: www.perplex.ethz.ch/perplex_options.html Reading computational options from: perplex_option.dat The current data base components are: MGO AL2O3 SIO2 CAO FEO Transform them (Y/N)? Calculations with saturated components (Y/N)? Use chemical potentials, activities or fugacities as independent variables (Y/N)? Select thermodynamic components from the set: MGO AL2O3 SIO2 CAO FEO Enter names, 1 per line, press to finish: MGO AL2O3 SIO2 CAO FEO Specify computational mode: 1 - Unconstrained 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 unconstrained 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. 3 Enter path coordinates from a file (Y/N)? y In this mode VERTEX/WERAMI read path coordinatesfrom a file the file must be formatted so that the coordinates of each point are on a separate line, the coordinates are, in order: P(bar) T(K) Enter the coordinate file name [default = coor.dat]: aerotherm.dat Specify component amounts by weight (Y/N)? y Enter weight amounts of the components: MGO AL2O3 SIO2 CAO FEO for the bulk composition of interest: 31.5 2.31 42.71 2.49 18,71 Output a print file (Y/N)? Exclude endmember phases (Y/N)? y Do you want to be prompted for phases (Y/N)? n Enter names, 1 per line, press to finish: stv Include solution phases (Y/N)? y Enter the solution model file name [default = solution_model.dat]: **warning ver025** 0 endmembers for Bio(TCC) The solution will not be considered. **warning ver114** the following endmembers are missing for Opx(HP) mgts fets_d **warning ver114** the following endmembers are missing for CrOpx(HP) crts fcrts_d mgts fets_d **warning ver025** 0 endmembers for Si-vapor The solution will not be considered. **warning ver114** the following endmembers are missing for Gt(MPF) skiag kho andr **warning ver025** 0 endmembers for oAmph(DP) The solution will not be considered. **warning ver114** the following endmembers are missing for Omph(GHP) jd acm **warning ver025** 0 endmembers for cAmph(DP) The solution will not be considered. **warning ver025** 0 endmembers for Chl(HP) The solution will not be considered. **warning ver025** 0 endmembers for melt(HP) The solution will not be considered. **warning ver025** 0 endmembers for pMELTS(G) The solution will not be considered. **warning ver025** 0 endmembers for MELTS(GS) The solution will not be considered. **warning ver025** 0 endmembers for Ep(HP) The solution will not be considered. **warning ver025** 0 endmembers for Pheng(HP) The solution will not be considered. **warning ver025** 0 endmembers for Sapp(HP) The solution will not be considered. **warning ver025** 0 endmembers for Sapp(KWP) The solution will not be considered. **warning ver025** 0 endmembers for Osm(HP) The solution will not be considered. **warning ver114** the following endmembers are missing for GaHcSp gah **warning ver025** 0 endmembers for F The solution will not be considered. **warning ver025** 0 endmembers for F(salt) The solution will not be considered. **warning ver025** 0 endmembers for T The solution will not be considered. **warning ver025** 0 endmembers for Scap The solution will not be considered. **warning ver025** 0 endmembers for St(HP) The solution will not be considered. **warning ver025** 0 endmembers for Ctd(HP) The solution will not be considered. **warning ver025** 0 endmembers for Carp The solution will not be considered. **warning ver025** 0 endmembers for hCrd The solution will not be considered. **warning ver025** 0 endmembers for Sud(Livi) The solution will not be considered. **warning ver025** 0 endmembers for Sud The solution will not be considered. **warning ver025** 0 endmembers for Cumm The solution will not be considered. **warning ver025** 0 endmembers for Anth The solution will not be considered. **warning ver025** 0 endmembers for A The solution will not be considered. **warning ver025** 0 endmembers for Gl The solution will not be considered. **warning ver025** 0 endmembers for Tr The solution will not be considered. **warning ver025** 1 endmembers for GlTrTsPg The solution will not be considered. **warning ver025** 1 endmembers for Amph(DHP) The solution will not be considered. **warning ver025** 1 endmembers for Amph(DPW) The solution will not be considered. **warning ver025** 1 endmembers for feldspar The solution will not be considered. **warning ver025** 1 endmembers for feldspar_B The solution will not be considered. **warning ver025** 1 endmembers for Pl(h) The solution will not be considered. **warning ver025** 0 endmembers for Kf The solution will not be considered. **warning ver025** 0 endmembers for San The solution will not be considered. **warning ver025** 0 endmembers for San(TH) The solution will not be considered. **warning ver025** 0 endmembers for GCOHF The solution will not be considered. **warning ver025** 0 endmembers for MaPa The solution will not be considered. **warning ver025** 0 endmembers for Mica(CF) The solution will not be considered. **warning ver025** 0 endmembers for Mica(CHA1) The solution will not be considered. **warning ver025** 0 endmembers for Mica(CHA) The solution will not be considered. **warning ver114** the following endmembers are missing for O(HP) teph **warning ver114** the following endmembers are missing for Cpx(l) cats jd **warning ver114** the following endmembers are missing for Cpx(h) cats jd **warning ver114** the following endmembers are missing for Mont mont **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** 0 endmembers for MF The solution will not be considered. **warning ver025** 0 endmembers for Mt(W) The solution will not be considered. **warning ver025** 0 endmembers for IlHm(A) The solution will not be considered. **warning ver025** 0 endmembers for IlGkPy The solution will not be considered. **warning ver025** 0 endmembers for MtUl(A) The solution will not be considered. **warning ver025** 0 endmembers for Neph(FB) The solution will not be considered. **warning ver114** the following endmembers are missing for GrPyAlSp(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(EWHP) andr FA_d MA_d **warning ver114** the following endmembers are missing for Gt(WPH) spss fmn_i fkho_i kho_i andr **warning ver025** 0 endmembers for casmelt The solution will not be considered. **warning ver025** 0 endmembers for A-phase The solution will not be considered. **warning ver025** 0 endmembers for Chum The solution will not be considered. **warning ver025** 0 endmembers for Atg The solution will not be considered. **warning ver025** 0 endmembers for B The solution will not be considered. **warning ver025** 1 endmembers for P The solution will not be considered. **warning ver025** 0 endmembers for Toop-Melt The solution will not be considered. **warning ver025** 1 endmembers for OrFsp(C1) The solution will not be considered. **warning ver025** 1 endmembers for AbFsp(C1) The solution will not be considered. **warning ver025** 1 endmembers for Pl(I1,HP) The solution will not be considered. **warning ver025** 1 endmembers for Fsp(C1) 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** 0 endmembers for o-Amph The solution will not be considered. **warning ver114** the following endmembers are missing for Omph(HP) jd cats acm **warning ver114** the following endmembers are missing for CrGt uv_d fuv_d knor **warning ver114** the following endmembers are missing for Cpx(HP) ccrts cats jd acm **warning ver114** the following endmembers are missing for CrSp mcrm fcrm_d **warning ver025** 1 endmembers for Eskol(C) The solution will not be considered. **warning ver025** 1 endmembers for Ca-Amph(D) The solution will not be considered. **warning ver025** 0 endmembers for Na-Amph(D) The solution will not be considered. **warning ver114** the following endmembers are missing for Sp(WPC) usp mt **warning ver025** 0 endmembers for Ilm(WPH) The solution will not be considered. **warning ver114** the following endmembers are missing for Cpx(m) cats **warning ver025** 1 endmembers for Ol(m) The solution will not be considered. **warning ver025** 1 endmembers for Pl(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) jd cts **warning ver114** the following endmembers are missing for Gt(stx8) jmaj **warning ver025** 0 endmembers for CF(stx8) The solution will not be considered. **warning ver025** 0 endmembers for Pu(M) The solution will not be considered. **warning ver025** 0 endmembers for Act(M) The solution will not be considered. **warning ver025** 0 endmembers for Stlp(M) The solution will not be considered. **warning ver025** 0 endmembers for Mica(M) The solution will not be considered. **warning ver025** 0 endmembers for Carp(M) The solution will not be considered. **warning ver025** 0 endmembers for Sud(M) The solution will not be considered. **warning ver025** 1 endmembers for GlTrTsMr The solution will not be considered. **warning ver025** 0 endmembers for Carb(M) The solution will not be considered. **warning ver025** 0 endmembers for TiBio(HP) The solution will not be considered. Select phases from the following list, enter 1 per line, press to finish Opx(HP) CrOpx(HP) Gt(MPF) Omph(GHP) O(SG) GaHcSp O(HP) Cpx(l) Cpx(h) Mont Sp(JR) Sp(GS) Sp(HP) GrPyAlSp(B Gt(GCT) Gt(HP) Gt(EWHP) Gt(WPH) Wus(fab) Aki(fab) Pv(fab) Ppv(og) O(stx) Wad(stx) Ring(stx) Sp(stx) Gt(stx) C2/c(stx) Opx(stx) Cpx(stx) Cpx(stx7) Omph(HP) CrGt Cpx(HP) CrSp Wus(stx7) Aki(stx7) Pv(stx7) O(stx7) Wad(stx7) Ring(stx7) Sp(stx7) Sp(WPC) Cpx(m) Sp(stx8) O(stx8) Wad(stx8) Ring(stx8) Opx(stx8) Cpx(stx8) Aki(stx8) Gt(stx8) Ppv(stx8) For details on these models see:www.perplex.ethz.ch/perplex_solution_model_glossary.html or read the commentary in the solution model file. Wus(fab) Pv(fab) O(stx) Wad(stx) Ring(stx) C2/c(stx) Opx(stx) Cpx(stx) Sp(stx) Gt(stx) Aki(fab) Enter calculation title: Martian SNC model C:\jamie\Berple_X\NewVersion>vertex Perple_X version 6.6.5.9, compiled 6/20/2011. Enter the project name (the name assigned in BUILD) [default = my_project]: in24.dat file/project names cannot include . characters, try again in24 Reading problem definition from file: in24.dat Reading thermodynamic data from file: sfo05ver.dat Reading solution models from file: solution_model.dat Writing print output to file: none requested Writing plot output to file: in24.plt Writing bulk composition plot output to file: in24.blk Reading computational options from: perplex_option.dat Writing auto refine summary to file: not requested Writing pseudocompound glossary to file: not requested Writing computational option summary to file: not requested Perple_X computational option settings for VERTEX: Keyword: Value: Permitted values [default]: Auto-refine options: auto_refine aut off [manual] auto auto_refine_factor_I 2.0 >=1 [3] auto_refine_slop 0.010 [0.01] Free energy minimization options: zero_mode 0.1E-05 0->1 [1e-6]; < 0 => off zero_bulk 0.1E-05 0->1 [1e-6]; < 0 => off iteration limit 4 >0 [4]; iteration keyword value 1 refinement factor 3 >2*r [12]; iteration keyword value 2 refinement points 4 3->7 [4]; iteration keyword value 3 reach_factor 0.9 >0.5 [0.9] Fractionation path options: 1d_path 40 /150 [20/150], >0, <2048 Solution subdivision options: initial_resolution 0.07 0->1 [0.1], 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.1E-02 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 [T] F site_check T [T] F Input/Output options: dependent_potentials on off [on] logarithmic_p F [F] T bad_number NaN [0.0] Information file output options: option_list_files F [F] T; echo computational options pseudocompound_file F [F] T; echo static pseudocompound compositions auto_refine_file F [F] T; echo auto-refine compositions Worst case (Cartesian) compositional resolution (mol %): Exploratory stage: 0.098 Auto-refine stage: 0.049 To change these options see: www.perplex.ethz.ch/perplex_options.html Writing data for auto-refinement to file: in24.arf 14 pseudocompounds generated for: Wus(fab) 133 pseudocompounds generated for: Aki(fab) 133 pseudocompounds generated for: Pv(fab) 14 pseudocompounds generated for: O(stx) 14 pseudocompounds generated for: Wad(stx) 14 pseudocompounds generated for: Ring(stx) 14 pseudocompounds generated for: Sp(stx) 812 pseudocompounds generated for: Gt(stx) 14 pseudocompounds generated for: C2/c(stx) 133 pseudocompounds generated for: Opx(stx) 133 pseudocompounds generated for: Cpx(stx) Total number of pseudocompounds: 1428 ** Starting exploratory computational stage ** Fractionate phases (y/n)? Endmember compositional ranges for model: Wus(fab) Endmember Minimum Maximum per 0.62562 0.63427 Endmember compositional ranges for model: Aki(fab) Endmember Minimum Maximum cor 0.00372 0.00744 aki 0.79358 0.80102 Endmember compositional ranges for model: Pv(fab) Endmember Minimum Maximum aperov 0.00521 0.01638 perov 0.80996 0.81591 Endmember compositional ranges for model: O(stx) Endmember Minimum Maximum fo 0.73781 0.79953 Endmember compositional ranges for model: Wad(stx) Endmember Minimum Maximum wad 0.73104 0.82484 Endmember compositional ranges for model: Ring(stx) Endmember Minimum Maximum ring 0.46304 0.72465 Endmember compositional ranges for model: Sp(stx) Endmember Minimum Maximum sp 0.31416 0.34691 Endmember compositional ranges for model: Gt(stx) Endmember Minimum Maximum gr 0.05360 0.15115 alm 0.10869 0.30540 maj 0.02159 0.74319 Endmember compositional ranges for model: C2/c(stx) Endmember Minimum Maximum c2/c 0.74891 0.75784 Endmember compositional ranges for model: Opx(stx) Endmember Minimum Maximum en 0.68936 0.76678 fs 0.23322 0.28587 Endmember compositional ranges for model: Cpx(stx) Endmember Minimum Maximum di 0.68489 0.89184 hed 0.07593 0.17271 Reading data for auto-refinement from file: in24.arf 2 pseudocompounds generated for: Wus(fab) 4 pseudocompounds generated for: Aki(fab) 4 pseudocompounds generated for: Pv(fab) 3 pseudocompounds generated for: O(stx) 4 pseudocompounds generated for: Wad(stx) 9 pseudocompounds generated for: Ring(stx) 2 pseudocompounds generated for: Sp(stx) 592 pseudocompounds generated for: Gt(stx) 2 pseudocompounds generated for: C2/c(stx) 9 pseudocompounds generated for: Opx(stx) 29 pseudocompounds generated for: Cpx(stx) Total number of pseudocompounds: 660 ** Starting auto_refine computational stage ** The fractionated bulk composition will be written to file: fractionated_bulk.dat WARNING: The compositions of the following solutions reached internal limits that were automatically relaxed. Aki(fab) Pv(fab) O(stx) Wad(stx) Ring(stx) Gt(stx) Opx(stx) Cpx(stx) Endmember compositional ranges for model: Wus(fab) Endmember Minimum Maximum per 0.62678 0.63319 Endmember compositional ranges for model: Aki(fab) Endmember Minimum Maximum cor 0.00372 0.00707 aki 0.79469 0.80065 Endmember compositional ranges for model: Pv(fab) Endmember Minimum Maximum aperov 0.00521 0.01563 perov 0.80996 0.81666 Endmember compositional ranges for model: O(stx) Endmember Minimum Maximum fo 0.73781 0.79916 Endmember compositional ranges for model: Wad(stx) Endmember Minimum Maximum wad 0.73179 0.82484 Endmember compositional ranges for model: Ring(stx) Endmember Minimum Maximum ring 0.46304 0.72441 Endmember compositional ranges for model: Sp(stx) Endmember Minimum Maximum sp 0.31639 0.34617 Endmember compositional ranges for model: Gt(stx) Endmember Minimum Maximum gr 0.05360 0.15113 alm 0.10943 0.30540 maj 0.02159 0.74304 Endmember compositional ranges for model: C2/c(stx) Endmember Minimum Maximum c2/c 0.74966 0.75747 Endmember compositional ranges for model: Opx(stx) Endmember Minimum Maximum en 0.68936 0.76641 fs 0.23322 0.28587 Endmember compositional ranges for model: Cpx(stx) Endmember Minimum Maximum di 0.68489 0.89147 hed 0.07668 0.17271 C:\jamie\Berple_X\NewVersion>pssect Perple_X version 6.6.5.9, compiled 6/20/2011. Enter the project name (the name assigned in BUILD) [default = my_project]: in24 Reading computational options from: perplex_option.dat Perple_X plot options are currently set as: Keyword: Value: Permitted values [default]: axis_label_scale 1.20 [1.2] (rel) bounding_box : 0 [0] x-min (pts) 0 [0] y-min (pts) 800 [800] x-length (pts) 800 [800] y-length (pts) field_fill T [T] F field_label T [T] F field_label_scale 0.75 [0.72] (rel) font Helvetica grid F [F] T half_ticks T [T] F line_width 1.00 0-99 [1.] (pts) picture_transformation : 0.180 [0.18] x-scale (rel) 0.180 [0.18] y-scale (rel) 130. [0.18] x-translation (pts) 220. [0.18] y-translation (pts) 0.00 [0.0] rotation (deg) plot_aspect_ratio 1.000 [1.0] x_axis_length/y_axis_length replicate_label 0.250 0->1 [0.025] splines T [T] F tenth_ticks F [F] T text_scale 1.000 [1.] (rel) To change these options edit or create the plot option file See: www.perplex.ethz.ch/perplex_plot_options.html PostScript will be written to file: in24.ps Modify the default plot (y/n)? C:\jamie\Berple_X\NewVersion>werami Perple_X version 6.6.5.9, compiled 6/20/2011. Enter the project name (the name assigned in BUILD) [default = my_project]: in24 Reading computational options from: perplex_option.dat Writing computational option summary to file: not requested Perple_X computational option settings for WERAMI: Keyword: Value: Permitted values [default]: Input/Output options: spreadsheet T [F] T logarithmic_p F [F] T bad_number NaN [0.0] compositions mol wt [mol] proportions vol wt [vol] mol interpolation on off [on ] extrapolation off on [off] melt_is_fluid F [F] T Information file output options: option_list_files F [F] T; echo computational options Thermodynamic options: approx_alpha T [T] F Anderson-Gruneisen F [T] F Seismic velocity options: bounds VRH HS [VRH] vrh/hs_weighting 0.5 0->1 [0.5] explicit_bulk_modulus T [F] T poisson_ratio on off [on ] all; Poisson ratio = 0.35 To change these options see: www.perplex.ethz.ch/perplex_options.html Select operational mode: 1 - properties at specified conditions 3 - properties along a 1d path 0 - EXIT 3 The independent variable for this calculation is: node # Choose the independent variable for data plots: 1 - node # [default] 2 - P(bar) 3 - T(K) 2 Enter the node # coordinate at the beginning of the profile: 1 Enter the node # coordinate at the end of1 the profile: 95 How many points along the profile? 95 Select properties [enter 0 to finish]: 1 - Specific Enthalpy (J/m3) 2 - Density (kg/m3) 3 - Specific heat capacity (J/K/m3) 4 - Expansivity (1/K, for volume) 5 - Compressibility (1/bar, for volume) 6 - Weight (%) of a component 7 - Mode (Vol, Mol, or Wt proportion) of a phase 8 - Composition (Mol or Wt) of a solution 9 - Grueneisen thermal ratio 10 - Adiabatic bulk modulus (bar) 11 - Adiabatic shear modulus (bar) 12 - Sound velocity (km/s) 13 - P-wave velocity (Vp, km/s) 14 - S-wave velocity (Vs, km/s) 15 - Vp/Vs 16 - Specific entropy (J/K/m3) 17 - Entropy (J/K/kg) 18 - Enthalpy (J/kg) 19 - Heat Capacity (J/K/kg) 20 - Specific mass of a phase (kg/m3-system) 21 - Poisson ratio 22 - Molar Volume (J/bar) 23 - Dependent potentials (J/mol, bar, K) 24 - Assemblage Index 25 - Modes of all phases 26 - Sound velocity T derivative (km/s/K) 27 - P-wave velocity T derivative (km/s/K) 28 - S-wave velocity T derivative (km/s/K) 29 - Adiabatic bulk modulus T derivative (bar/K) 30 - Shear modulus T derivative (bar/K) 31 - Sound velocity P derivative (km/s/bar) 32 - P-wave velocity P derivative (km/s/bar) 33 - S-wave velocity P derivative (km/s/bar) 34 - Adiabatic bulk modulus P derivative (unitless) 35 - Shear modulus P derivative (unitless) 37 - Absolute amount (Vol, Mol, or Wt) of a phase 38 - Multiple property output (PHEMGP format) 25 Output cumulative modes (y/n)? (see www.perplex.ethz.ch/perplex_options.html#cumulative_modes) y Data ranges excluding values equal to bad_number ( NaN ) specified in perplex_option.dat: O(stx) Sp(stx) Opx(stx) Cpx(stx) Gt(stx) C2/c(stx) Ring(stx) Wad(stx) ca-pv Aki(fab) Pv(fab) Wus(fab) min 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 max 60.76207 60.76207 89.04490 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 Output has been written to two files: plt format is in file: in24_1.plt 1d tab format is in file: in24_1.tab plt format files can be plotted with: PSVDRAW 1d tab format files can be processed with: PSTABLE - a Perple_X plotting program PERPLE_X_PLOT - a Matlab plotting script spread-sheet programs, e.g., EXCEL for details on tab format refer to: www.perplex.ethz.ch/faq/perple_x_tab_file_format.txt Select operational mode: 1 - properties at specified conditions 3 - properties along a 1d path 0 - EXIT 0 C:\jamie\Berple_X\NewVersion>psvdraw Perple_X version 6.6.5.9, compiled 6/20/2011. Enter the project or plot file name [i.e., without the .plt suffix]: in24_1 Perple_X plot options are currently set as: Keyword: Value: Permitted values [default]: axis_label_scale 1.20 [1.2] (rel) bounding_box : 0 [0] x-min (pts) 0 [0] y-min (pts) 800 [800] x-length (pts) 800 [800] y-length (pts) field_fill T [T] F field_label T [T] F field_label_scale 0.75 [0.72] (rel) font Helvetica grid F [F] T half_ticks T [T] F line_width 1.00 0-99 [1.] (pts) picture_transformation : 0.180 [0.18] x-scale (rel) 0.180 [0.18] y-scale (rel) 130. [0.18] x-translation (pts) 220. [0.18] y-translation (pts) 0.00 [0.0] rotation (deg) plot_aspect_ratio 1.000 [1.0] x_axis_length/y_axis_length replicate_label 0.250 0->1 [0.025] splines T [T] F tenth_ticks F [F] T text_scale 1.000 [1.] (rel) To change these options edit or create the plot option file See: www.perplex.ethz.ch/perplex_plot_options.html PostScript will be written to file: in24_1.ps Modify the default plot (y/n)? C:\jamie\Berple_X\NewVersion>