Microsoft Windows [Version 6.1.7601] Copyright (c) 2009 Microsoft Corporation. All rights reserved. C:\jamie\Perple_X>werami Perple_X version 6.8.4, source updated Sept 15, 2018. Copyright (C) 1986-2018 James A D Connolly . Enter the project name (the name assigned in BUILD) [default = my_project]: gloss1di Reading Perple_X options from: gloss1di_perplex_option.dat Writing Perple_X option summary to file: not requested Perple_X version 6.8.4, source updated Sept 15, 2018. Copyright (C) 1986-2018 James A D Connolly . Perple_X computational option settings for WERAMI: Keyword: Value: Permitted values [default]: Input/Output options: aqueous_output T [F] T aqeuous_species 20 [20] 0-100 aq_solvent_composition m [y] m: y => mol fraction, m => molality aq_solute_composition m y [m]: y => mol fraction, m => molality spreadsheet T [F] T logarithmic_p F [F] T bad_number NaN [0.0] composition_constant F [F] T composition_phase mol [mol] wt composition_system mol [wt] mol proportions vol [vol] wt mol interpolation on [on] off melt_is_fluid F [F] T solution_names mod [model] abbreviation full species_output T [T] F species_Gibbs_energies F [F] T seismic_output som [some] none all pause_on_error T [T] F poisson_test 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 [F] T hybrid_EoS_H2O 0 [7] 0-2, 4-8 hybrid_EoS_CO2 5 [7] 0-4 hybrid_EoS_CH4 4 [1] 0-1 fd_expansion_factor 2.0 >0 [2.] finite_difference_p 0.1D+05 >0 [1d4]; fraction = 0.1D-02 [1d-2] 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 To change these options see: www.perplex.ethz.ch/perplex_options.html **warning ver099** aq_output or aq_lagged_speciation is T, but dependent_potentials is off, dependent_potentials set = on (AQIDST) **warning ver099** aq_lagged_speciation is T, but refine_endmembers is F, refine_endmembers set = T (AQIDST) Select operational mode: 1 - properties at specified conditions 3 - properties along a 1d path 0 - EXIT 3 Select a property [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 - Composition (Mol, Mass, or Wt%) of the system 7 - Mode (Vol, Mol, or Wt proportion) of a phase 8 - Composition (Mol, Mass, or Wt%) of a solution phase 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) 36 - All phase &/or system properties 37 - Absolute amount (Vol, Mol, or Wt) of a phase 38 - Multiple property output 39 - Heat capacity ratio (Cp/Cv) 40 - Lagged or back-calculated aqueous solute chemistry 25 Output cumulative modes (y/n)? (see www.perplex.ethz.ch/perplex_options.html#cumulative_modes) y **warning ver210** the fancy_cumulative_modes option (perplex_option.dat) is T. This option may lead to inconsistent curves for coarsely sampled profiles. Include fluid in computation of aggregate (or modal) properties (y/n)? y **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Grt(JH) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Mica(CF) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Do(HP) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Omph(GHP) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: law is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: coe is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: pyr is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Grt(JH) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Mica(CF) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Do(HP) is missing or invalid and has been estimated from the default poisson ratio **warning ver178** at T(K)= 946.5 P(bar)= 0.4300E+05 the shear modulus of: Omph(GHP) is missing or invalid and has been estimated from the default poisson ratio **warning ver049** warning 178 will not be repeated for future instances of this problem. currently in routine: GETPHP Data ranges excluding values equal to bad_number ( NaN) specified in perplex_option.dat: Grt(JH) Mica(CF) Do(HP) Omph(GHP) COH-Fluid+ law coe pyr M(HP) ky min 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 max 56.1883 64.6615 10.3769 56.1883 9.27418 26.6859 100.000 10.3769 8.64021 8.64021 Output has been written to two files: plt format is in file: gloss1di_1.plt 1d tab format is in file: gloss1di_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: perplex.ethz.ch/perplex/faq/Perple_X_tab_file_format.txt Select operational mode: 1 - properties at specified conditions 3 - properties along a 1d path 0 - EXIT 3 Select a property [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 - Composition (Mol, Mass, or Wt%) of the system 7 - Mode (Vol, Mol, or Wt proportion) of a phase 8 - Composition (Mol, Mass, or Wt%) of a solution phase 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) 36 - All phase &/or system properties 37 - Absolute amount (Vol, Mol, or Wt) of a phase 38 - Multiple property output 39 - Heat capacity ratio (Cp/Cv) 40 - Lagged or back-calculated aqueous solute chemistry 40 This calculation was done with lagged speciation, but back-calculated speciation calculations are also enabled. Output the back-calculated results (y/n)? Data ranges excluding values equal to bad_number ( NaN) specified in perplex_option.dat: H2,mol_abs C,mol_abs Si,mol_abs Al,mol_abs Fe,mol_abs Mg,mol_abs Ca,mol_abs Na,mol_abs K,mol_abs O2,mol_abs S2,mol_abs m_{CO2} m_{H2O} m_{FeC4H6O4} m_{FeC2H3O2} m_{NaCO3-} m_{CO3-2} m_{BENZENE} m_{ACETATE} m_{ACETIC- A} m_{Al+3} m_{AlO2-} m_{Ca(HCO3)} m_{Ca(OH)+} m_{Ca+2} m_{CaCO3,aq} m_{CaSO4,aq} m_{CO,aq} m_{ETHANE,A} m_{ETHANOL,} m_{ETHYLENE} m_{Fe+2} m_{Fe+3} m_{FORMATE,} m_{FORMIC-A} m_{GLUTARIC} m_{GLYCOLIC} m_{H+} m_{H2,aq} m_{H2S,aq} m_{HAlO2,aq} m_{HCO3-} m_{HEXANE,A} m_{HO2-} m_{HS-} m_{HSiO3-} m_{HSO3-} m_{HSO4-} m_{HSO5-} m_{K+} m_{KOH,aq} m_{KSO4-} m_{LACTIC-A} m_{METHANE,} m_{METHANOL} m_{Mg(HCO3)} m_{Mg(HSiO3} m_{Mg+2} m_{MgCO3,aq} m_{MgOH+} m_{MgSO4,aq} m_{Na+} m_{NaHCO3,a} m_{NaHSiO3,} m_{NaOH,aq} m_{O2,aq} m_{OH-} m_{PROPANE,} m_{PROPANOA} m_{PROPANOI} m_{PROPANOL} m_{S2-2} m_{S2O3-2} m_{S2O4-2} m_{S2O5-2} m_{S2O6-2} m_{S2O8-2} m_{S3-2} m_{S3O6-2} m_{S4-2} m_{S4O6-2} m_{S5-2} m_{S5O6-2} m_{Si2O4,aq} m_{SiO2,aq} m_{SO2,aq} m_{SO3-2} m_{SO4-2} m_{TOLUENE,} pH-pH_0 pH error_pH permittivity I,m tot_solute_m min 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.633694E-001 0.00000 0.00000 0.00000 max 56.1883 64.6615 10.3769 56.1883 9.27418 26.6859 100.000 10.3769 8.64021 8.64021 0.105658E-003 1.14585 53.9019 0.193038E-009 0.198289E-007 0.305840E-003 0.423792E-002 0.272055E-021 0.211805E-004 0.245367E -004 0.910385E-010 0.116988E-002 1.27698 0.905453E-002 0.216780 0.166635E-001 0.124701 0.116823E-003 0.655435E-009 0.517089E-006 0.289508E-011 0.214286E-003 0.431223E-009 0.390134E-001 0.266316E-002 0.132049E-018 0.326274E-010 0.843846E-003 0.148798E-002 0.111875E -001 0.218879E-002 0.143291 0.134005E-006 0.689313E-018 0.697719E-003 2.77779 0.893982E-002 0.211103E-001 0.631614E-017 1.07943 0.187859E-001 0.136675E-003 0.791754E-014 0.144375E-004 0.272947E-005 0.148882E-002 0.502660E-001 0.190187E-001 0.326537E-009 0.878626E -002 0.198302E-004 1.60466 0.108881E-001 0.910996E-002 0.538470E-002 0.923642E-020 0.272921E-001 0.555845E-013 0.251615E-005 0.272378E-009 0.378046E-012 0.109077E-008 0.588972E-007 0.161294E-013 0.174227E-013 0.248968E-014 0.598167E-030 0.155324E-009 0.151974E-015 0.216711E -010 0.666743E-013 0.309643E-011 0.477627E-017 0.145200 0.664002 0.286892E-004 0.243869E-007 0.186555E-002 0.987577E-024 1.30604 4.16544 0.00000 31.2963 3.01219 6.86117 Output has been written to the 1d tab format file: gloss1di_2.tab 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: perplex.ethz.ch/perplex/faq/Perple_X_tab_file_format.txt Select operational mode: 1 - properties at specified conditions 3 - properties along a 1d path 0 - EXIT 3 Select a property [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 - Composition (Mol, Mass, or Wt%) of the system 7 - Mode (Vol, Mol, or Wt proportion) of a phase 8 - Composition (Mol, Mass, or Wt%) of a solution phase 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) 36 - All phase &/or system properties 37 - Absolute amount (Vol, Mol, or Wt) of a phase 38 - Multiple property output 39 - Heat capacity ratio (Cp/Cv) 40 - Lagged or back-calculated aqueous solute chemistry 36 In this mode you may tabulate: 1 - properties of the system 2 - properties of a phase 3 - properties of the system and its phases Output for option 1 & 2 can be plotted with PSPLOT, PYWERAMI or MatLab. Output for option 3 can only be plotted with PHEMGP. Select an option [default = 1]: 2 Enter solution or compound (left justified): COH-Fluid+ Include fluid in computation of aggregate (or modal) properties (y/n)? Data ranges excluding values equal to bad_number ( NaN) specified in perplex_option.dat: V,J/bar/mol H,J/mol Gruneisen_T Ks,bar Gs,bar v0,km/s vp,km/s vs,km/s vp/vs rho,kg/m3 G,J/mol cp,J/K/mol alpha,1/K beta,1/bar S,J/K/mol n,mol N,g Ks_{T},bar/K Gs_{T},bar/K Ks_{P} Gs_P v0_{T} vp_{T} vs_{T} v0_{P} vp_P vs_{P} cp/cv vol,% wt,% mol,% H2,mol_abs C,mol_abs Si,mol_abs Al,mol_abs Fe,mol_abs Mg,mol_abs Ca,mol_abs Na,mol_abs K,mol_abs O2,mol_abs S2,mol_abs mu[H2],J/mol mu[C],J/mol mu[Si],J/mol mu[Al],J/mol mu[Fe],J/mol mu[Mg],J/mol mu[Ca],J/mol mu[Na],J/mol mu[K],J/mol mu[O2],J/mol mu[S2],J/mol nom_ox min 0.00000 -184226. 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -313449. 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -113.936 0.00000 0.00000 0.00000 -0.933584E-003 -0.933584E-003 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -54127.9 -11025.0 -434361. -495054. -96793.9 -416745. -510142. -265016. -256410. -400406. -34732.7 -0.878266E-001 max 56.1883 64.6615 10.3769 197721. 9.27418 26.6859 100.000 10.3769 8.64021 1802.43 0.105658E-003 61.5599 53.9019 0.600944E-005 136.528 0.133438 24.3718 0.272055E-021 0.211805E-004 4.79990 0.910385E-010 0.116988E-002 1.27698 0.905453E-002 0.216780 0.166635E-001 0.124701 1.00809 0.655435E-009 0.517089E-006 0.289508E-011 80.3257 0.706762 3.31455 0.389985E-002 0.158873E-003 0.633503E-001 0.519752 0.786221 0.460642 44.8085 0.143291 0.134005E-006 10612.7 0.697719E-003 2.77779 0.893982E-002 0.211103E-001 0.631614E-017 1.07943 0.187859E-001 0.136675E-003 0.791754E-014 0.144375E-004 Output has been written to the 1d tab format file: gloss1di_3.tab 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: perplex.ethz.ch/perplex/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\Perple_X>