Input file: ../examples/ex20a Output file: ex20a.out Database file: ../database/iso.dat ------------------ Reading data base. ------------------ SOLUTION_MASTER_SPECIES SOLUTION_SPECIES PHASES EXCHANGE_MASTER_SPECIES EXCHANGE_SPECIES SURFACE_MASTER_SPECIES SURFACE_SPECIES SOLUTION_MASTER_SPECIES ISOTOPES ISOTOPE_RATIOS ISOTOPE_ALPHAS NAMED_EXPRESSIONS CALCULATE_VALUES CALCULATE_VALUES SOLUTION_SPECIES PHASES ------------------------------------ Reading input data for simulation 1. ------------------------------------ DATABASE ../database/iso.dat TITLE Example 20a.--Calculate carbonate solid solution PRINT censor_species 1e-006 SOLUTION 1 # water to find composition of marine carbonate pH 8.2 Na 1 charge Ca 10 Calcite 0 C 2 [13C] 0 # permil [14C] 0 # pmc D 0 # permil [18O] 0 # permil END ----- TITLE ----- Example 20a.--Calculate carbonate solid solution ------------------------------------------- Beginning of initial solution calculations. ------------------------------------------- Initial solution 1. -----------------------------------Isotopes------------------------------------ Isotope Molality Moles Ratio Units H 1.10997e+02 1.10997e+02 D 1.72889e-02 1.72889e-02 0.00000e+00 permil T 0.00000e+00 0.00000e+00 0.00000e+00 TU C 1.97789e-03 1.97789e-03 [13C] 2.21132e-05 2.21132e-05 0.00000e+00 permil [14C] 0.00000e+00 0.00000e+00 0.00000e+00 pmc O 5.54011e+01 5.54011e+01 [18O] 1.11090e-01 1.11090e-01 0.00000e+00 permil -----------------------------Solution composition------------------------------ Elements Molality Moles C 2.000e-03 2.000e-03 Ca 3.089e-04 3.089e-04 Equilibrium with Calcite Na 1.380e-03 1.380e-03 Charge balance ----------------------------Description of solution---------------------------- pH = 8.200 pe = 4.000 Activity of water = 1.000 Ionic strength = 2.294e-03 Mass of water (kg) = 1.000e+00 Total alkalinity (eq/kg) = 1.998e-03 Total CO2 (mol/kg) = 2.000e-03 Temperature (deg C) = 25.00 Electrical balance (eq) = 1.660e-13 Percent error, 100*(Cat-|An|)/(Cat+|An|) = 0.00 Iterations = 8 Total H = 1.109971e+02 Total O = 5.540110e+01 ----------------------------Distribution of species---------------------------- Log Log Log Species Molality Activity Molality Activity Gamma OH- 1.673e-06 1.586e-06 -5.776 -5.800 -0.023 H3O+ 6.628e-09 6.310e-09 -8.179 -8.200 -0.021 H2O 5.551e+01 9.999e-01 1.744 -0.000 0.000 C(-4) 0.000e+00 CH4 0.000e+00 0.000e+00 -77.792 -77.792 0.000 C(4) 2.000e-03 HCO3- 1.944e-03 1.844e-03 -2.711 -2.734 -0.023 CO2 2.615e-05 2.616e-05 -4.583 -4.582 0.000 CO3-2 1.693e-05 1.371e-05 -4.771 -4.863 -0.092 CaHCO3+ 5.989e-06 5.687e-06 -5.223 -5.245 -0.022 CaCO3 5.562e-06 5.565e-06 -5.255 -5.255 0.000 NaHCO3 1.356e-06 1.356e-06 -5.868 -5.868 0.000 NaCO3- 3.520e-07 3.339e-07 -6.453 -6.476 -0.023 Ca 3.089e-04 Ca+2 2.974e-04 2.417e-04 -3.527 -3.617 -0.090 CaHCO3+ 5.989e-06 5.687e-06 -5.223 -5.245 -0.022 CaCO3 5.562e-06 5.565e-06 -5.255 -5.255 0.000 H(0) 5.634e-28 H2 2.817e-28 2.819e-28 -27.550 -27.550 0.000 Na 1.380e-03 Na+ 1.379e-03 1.308e-03 -2.861 -2.883 -0.023 NaHCO3 1.356e-06 1.356e-06 -5.868 -5.868 0.000 NaCO3- 3.520e-07 3.339e-07 -6.453 -6.476 -0.023 O(0) 1.049e-37 O2 5.243e-38 5.246e-38 -37.280 -37.280 0.000 ------------------------------Saturation indices------------------------------- Phase SI log IAP log K(298 K, 1 atm) Calcite 0.00* -8.48 -8.48 CaCO3 Delta_V 54.21 cm3/mol CH4(g) -74.93* -77.79 -2.86 CH4 Delta_V -180.69 cm3/mol CO2(g) -3.11 -4.58 -1.47 CO2 H2(g) -24.40* -27.55 -3.15 H2 Delta_V -36.14 cm3/mol H2O(g) -1.51* -0.00 1.51 H2O Delta_V 18.07 cm3/mol O2(g) -34.39* -37.28 -2.89 O2 Delta_V 108.41 cm3/mol * with Delta_V * (P - 1) / 2.3RT. ------------------ End of simulation. ------------------ ------------------------------------ Reading input data for simulation 2. ------------------------------------ SOLID_SOLUTION 1 No [14C] Calcite component Calcite 0 component CaCO2[18O](s) 0 component CaCO[18O]2(s) 0 component CaC[18O]3(s) 0 component Ca[13C]O3(s) 0 component Ca[13C]O2[18O](s) 0 component Ca[13C]O[18O]2(s) 0 component Ca[13C][18O]3(s) 0 END ------------------ End of simulation. ------------------ ------------------------------------ Reading input data for simulation 3. ------------------------------------ RUN_CELLS -cells 1 USER_PRINT start 10 PRINT pad("Component", 20), "Mole fraction" 20 t = LIST_S_S("Calcite", count, name$, moles) 30 for i = 1 to count 40 PRINT pad(name$(i),20), moles(i)/t 50 next i end END -------------------------- Beginning of run as cells. -------------------------- ----------------------------------------- Beginning of batch-reaction calculations. ----------------------------------------- Reaction step 1. Using solution 1. Using solid solution assemblage 1. No [14C] ----------------------------------User print----------------------------------- Component Mole fraction Calcite 9.8283e-01 Ca[13C]O3(s) 1.1011e-02 CaCO2[18O](s) 6.0825e-03 Ca[13C]O2[18O](s) 6.8147e-05 CaCO[18O]2(s) 1.2548e-05 Ca[13C]O[18O]2(s) 1.4058e-07 CaC[18O]3(s) 8.6284e-09 Ca[13C][18O]3(s) 9.6671e-11 --------------------------------Solid solutions-------------------------------- Solid solution Component Moles Delta moles Mole fract Calcite 2.54e-07 Calcite 2.49e-07 2.49e-07 9.83e-01 CaCO2[18O](s) 1.54e-09 1.54e-09 6.08e-03 CaCO[18O]2(s) 3.18e-12 3.18e-12 1.25e-05 CaC[18O]3(s) 2.19e-15 2.19e-15 8.63e-09 Ca[13C]O3(s) 2.79e-09 2.79e-09 1.10e-02 Ca[13C]O2[18O](s) 1.73e-11 1.73e-11 6.81e-05 Ca[13C]O[18O]2(s) 3.56e-14 3.56e-14 1.41e-07 Ca[13C][18O]3(s) 2.45e-17 2.45e-17 9.67e-11 --------------------------------Isotope Ratios--------------------------------- Isotope Ratio Ratio Input Units R(D) 1.55760e-04 3.4803e-13 permil R(18O) 2.00520e-03 -3.9446e-07 permil R(13C) 1.11802e-02 -0.00026756 permil R(D) H2O(l) 1.55760e-04 1.138e-05 permil R(18O) H2O(l) 2.00520e-03 -3.9135e-05 permil R(D) OH- 3.70888e-05 -761.88 permil R(18O) OH- 1.93086e-03 -37.072 permil R(D) H3O+ 1.62250e-04 41.665 permil R(18O) H3O+ 2.05157e-03 23.123 permil R(D) H2(aq) 1.55760e-04 1.138e-05 permil R(13C) CO2(aq) 1.10852e-02 -8.4965 permil R(18O) CO2(aq) 2.08959e-03 42.084 permil R(D) HCO3- 1.55760e-04 1.1382e-05 permil R(18O) HCO3- 2.00520e-03 -3.9135e-05 permil R(13C) HCO3- 1.11816e-02 0.12943 permil R(18O) CO3-2 2.00520e-03 -3.9137e-05 permil R(13C) CO3-2 1.11656e-02 -1.3058 permil R(D) CH4(aq) 1.55760e-04 1.138e-05 permil R(13C) CH4(aq) 1.10852e-02 -8.4965 permil R(18O) Calcite 2.06293e-03 28.79 permil R(13C) Calcite 1.12038e-02 2.1103 permil --------------------------------Isotope Alphas--------------------------------- 1000ln(Alpha) ---------------------- Isotope Ratio Solution alpha Solution 25.0 C Alpha D OH-/H2O(l) 0.23812 -1435 -1435 Alpha 18O OH-/H2O(l) 0.96293 -37.777 -37.777 Alpha D H3O+/H2O(l) 1.0417 40.82 40.82 Alpha 18O H3O+/H2O(l) 1.0231 22.86 22.86 Alpha D H2(aq)/H2O(l) 1 8.6597e-12 0 Alpha 18O CO2(aq)/H2O(l) 1.0421 41.223 41.223 Alpha D HCO3-/H2O(l) 1 2.2582e-09 0 Alpha 18O HCO3-/H2O(l) 1 1.9984e-12 0 Alpha 13C HCO3-/CO2(aq) 1.0087 8.6622 8.6622 Alpha 18O CO3-2/H2O(l) 1 -1.653e-09 0 Alpha 13C CO3-2/CO2(aq) 1.0073 7.2261 7.2261 Alpha D CH4(aq)/H2O(l) 1 -4.2188e-12 0 Alpha 13C CH4(aq)/CO2(aq) 1 1.5543e-12 0 Alpha 18O Calcite/H2O(l) 1.0288 28.383 28.383 Alpha 13C Calcite/CO2(aq) 1.0107 10.641 10.641 -----------------------------Solution composition------------------------------ Elements Molality Moles [13C] 2.218e-05 2.211e-05 [18O] 1.114e-01 1.111e-01 C 1.984e-03 1.978e-03 Ca 3.097e-04 3.087e-04 D 1.734e-02 1.729e-02 Na 1.385e-03 1.380e-03 ----------------------------Description of solution---------------------------- pH = 8.199 Charge balance pe = -2.783 Adjusted to redox equilibrium Activity of water = 0.998 Ionic strength = 2.300e-03 Mass of water (kg) = 9.968e-01 Total alkalinity (eq/kg) = 2.004e-03 Total CO2 (mol/kg) = 1.984e-03 Temperature (deg C) = 25.00 Electrical balance (eq) = 1.660e-13 Percent error, 100*(Cat-|An|)/(Cat+|An|) = 0.00 Iterations = 40 Total H = 1.109971e+02 Total O = 5.540110e+01 ----------------------------Distribution of species---------------------------- Log Log Log Species Molality Activity Molality Activity Gamma OH- 1.661e-06 1.575e-06 -5.780 -5.803 -0.023 H3O+ 6.648e-09 6.328e-09 -8.177 -8.199 -0.021 H2O 5.556e+01 9.977e-01 1.745 -0.001 0.000 C(-4) 3.063e-24 CH4 3.061e-24 3.063e-24 -23.514 -23.514 0.000 CH3D 1.907e-27 1.908e-27 -26.720 -26.719 0.000 C(4) 1.984e-03 HCO3- 1.916e-03 1.818e-03 -2.718 -2.741 -0.023 CO2 2.597e-05 2.598e-05 -4.586 -4.585 0.000 CO3-2 1.661e-05 1.344e-05 -4.780 -4.872 -0.092 CaHCO3+ 5.916e-06 5.618e-06 -5.228 -5.250 -0.023 CaCO3 5.466e-06 5.469e-06 -5.262 -5.262 0.000 HCO2[18O]- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 HCO[18O]O- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 HC[18O]O2- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 NaHCO3 1.341e-06 1.341e-06 -5.873 -5.872 0.000 NaCO3- 3.463e-07 3.285e-07 -6.460 -6.483 -0.023 DCO3- 2.985e-07 2.831e-07 -6.525 -6.548 -0.023 CO[18O] 1.085e-07 1.086e-07 -6.964 -6.964 0.000 CO2[18O]-2 9.991e-08 8.087e-08 -7.000 -7.092 -0.092 CaCO2[18O] 3.288e-08 3.290e-08 -7.483 -7.483 0.000 CaHCO2[18O]+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 CaHCO[18O]O+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 CaHC[18O]O2+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 HCO[18O]2- 7.705e-09 7.308e-09 -8.113 -8.136 -0.023 HC[18O]2O- 7.705e-09 7.308e-09 -8.113 -8.136 -0.023 HC[18O]O[18O]- 7.705e-09 7.308e-09 -8.113 -8.136 -0.023 NaHCO2[18O] 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaHCO[18O]O 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaHC[18O]O2 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaCO2[18O]- 2.083e-09 1.976e-09 -8.681 -8.704 -0.023 Ca 3.097e-04 Ca+2 2.981e-04 2.422e-04 -3.526 -3.616 -0.090 CaHCO3+ 5.916e-06 5.618e-06 -5.228 -5.250 -0.023 CaCO3 5.466e-06 5.469e-06 -5.262 -5.262 0.000 CaH[13C]O3+ 6.616e-08 6.281e-08 -7.179 -7.202 -0.023 Ca[13C]O3 6.104e-08 6.107e-08 -7.214 -7.214 0.000 CaCO2[18O] 3.288e-08 3.290e-08 -7.483 -7.483 0.000 CaHCO2[18O]+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 CaHCO[18O]O+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 CaHC[18O]O2+ 1.186e-08 1.126e-08 -7.926 -7.948 -0.023 CaDCO3+ 9.216e-10 8.750e-10 -9.035 -9.058 -0.023 Ca[13C]O2[18O] 3.672e-10 3.674e-10 -9.435 -9.435 0.000 D(0) 3.259e-18 HD 3.259e-18 3.261e-18 -17.487 -17.487 0.000 D2 2.538e-22 2.539e-22 -21.596 -21.595 0.000 D(1) 1.734e-02 HDO 1.731e-02 3.108e-04 -1.762 -3.507 -1.746 HD[18O] 3.470e-05 6.233e-07 -4.460 -6.205 -1.746 D2O 1.348e-06 2.421e-08 -5.870 -7.616 -1.746 DCO3- 2.985e-07 2.831e-07 -6.525 -6.548 -0.023 H(0) 2.093e-14 H2 1.046e-14 1.047e-14 -13.980 -13.980 0.000 HD 3.259e-18 3.261e-18 -17.487 -17.487 0.000 Na 1.385e-03 Na+ 1.383e-03 1.312e-03 -2.859 -2.882 -0.023 NaHCO3 1.341e-06 1.341e-06 -5.873 -5.872 0.000 NaCO3- 3.463e-07 3.285e-07 -6.460 -6.483 -0.023 NaH[13C]O3 1.499e-08 1.500e-08 -7.824 -7.824 0.000 Na[13C]O3- 3.867e-09 3.668e-09 -8.413 -8.436 -0.023 NaHCO2[18O] 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaHCO[18O]O 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaHC[18O]O2 2.688e-09 2.690e-09 -8.571 -8.570 0.000 NaCO2[18O]- 2.083e-09 1.976e-09 -8.681 -8.704 -0.023 O(0) 0.000e+00 O2 0.000e+00 0.000e+00 -64.422 -64.422 0.000 O[18O] 0.000e+00 0.000e+00 -66.819 -66.818 0.000 [13C](-4) 3.396e-26 [13C]H4 3.394e-26 3.395e-26 -25.469 -25.469 0.000 [13C]H3D 2.114e-29 2.115e-29 -28.675 -28.675 0.000 [13C](4) 2.218e-05 H[13C]O3- 2.143e-05 2.032e-05 -4.669 -4.692 -0.023 [13C]O2 2.879e-07 2.880e-07 -6.541 -6.541 0.000 [13C]O3-2 1.854e-07 1.501e-07 -6.732 -6.824 -0.092 CaH[13C]O3+ 6.616e-08 6.281e-08 -7.179 -7.202 -0.023 Ca[13C]O3 6.104e-08 6.107e-08 -7.214 -7.214 0.000 H[13C]O2[18O]- 4.296e-08 4.075e-08 -7.367 -7.390 -0.023 H[13C]O[18O]O- 4.296e-08 4.075e-08 -7.367 -7.390 -0.023 H[13C][18O]O2- 4.296e-08 4.075e-08 -7.367 -7.390 -0.023 NaH[13C]O3 1.499e-08 1.500e-08 -7.824 -7.824 0.000 Na[13C]O3- 3.867e-09 3.668e-09 -8.413 -8.436 -0.023 D[13C]O3- 3.337e-09 3.166e-09 -8.477 -8.500 -0.023 [13C]O[18O] 1.203e-09 1.204e-09 -8.920 -8.919 0.000 [13C]O2[18O]-2 1.116e-09 9.029e-10 -8.953 -9.044 -0.092 Ca[13C]O2[18O] 3.672e-10 3.674e-10 -9.435 -9.435 0.000 CaH[13C]O2[18O]+ 1.327e-10 1.260e-10 -9.877 -9.900 -0.023 CaH[13C]O[18O]O+ 1.327e-10 1.260e-10 -9.877 -9.900 -0.023 CaH[13C][18O]O2+ 1.327e-10 1.260e-10 -9.877 -9.900 -0.023 H[13C]O[18O]2- 8.615e-11 8.172e-11 -10.065 -10.088 -0.023 H[13C][18O]2O- 8.615e-11 8.172e-11 -10.065 -10.088 -0.023 H[13C][18O]O[18O]- 8.615e-11 8.172e-11 -10.065 -10.088 -0.023 NaH[13C]O2[18O] 3.006e-11 3.008e-11 -10.522 -10.522 0.000 NaH[13C]O[18O]O 3.006e-11 3.008e-11 -10.522 -10.522 0.000 NaH[13C][18O]O2 3.006e-11 3.008e-11 -10.522 -10.522 0.000 Na[13C]O2[18O]- 2.326e-11 2.207e-11 -10.633 -10.656 -0.023 [18O](-2) 1.114e-01 H2[18O] 1.114e-01 2.001e-03 -0.953 -2.699 -1.746 HD[18O] 3.470e-05 6.233e-07 -4.460 -6.205 -1.746 HCO2[18O]- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 HCO[18O]O- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 HC[18O]O2- 3.842e-06 3.645e-06 -5.415 -5.438 -0.023 [18O](0) 0.000e+00 O[18O] 0.000e+00 0.000e+00 -66.819 -66.818 0.000 [18O]2 0.000e+00 0.000e+00 -69.817 -69.817 0.000 ------------------------------Saturation indices------------------------------- Phase SI log IAP log K(298 K, 1 atm) [13C][18O]2(g) -10.43* -11.94 -1.50 [13C][18O]2 Delta_V -36.14 cm3/mol [13C]D4(g) -37.84* -40.70 -2.86 [13C]D4 Delta_V -216.82 cm3/mol [13C]H2D2(g) -29.45* -33.08 -3.64 [13C]H2D2 Delta_V -198.75 cm3/mol [13C]H3D(g) -25.81* -29.28 -3.46 [13C]H3D Delta_V -189.72 cm3/mol [13C]H4(g) -22.61* -25.47 -2.86 [13C]H4 Delta_V -180.69 cm3/mol [13C]HD3(g) -33.43* -36.89 -3.46 [13C]HD3 Delta_V -207.79 cm3/mol [13C]O2(g) -5.07 -6.54 -1.47 [13C]O2 [13C]O[18O](g) -7.45* -9.24 -1.79 [13C]O[18O] Delta_V -18.07 cm3/mol [18O]2(g) -67.53* -69.82 -2.29 [18O]2 Delta_V 72.27 cm3/mol C[18O]2(g) -8.48* -9.98 -1.50 C[18O]2 Delta_V -36.14 cm3/mol Ca[13C][18O]3(s) -10.01 -1.86 8.16 Ca[13C][18O]3 Ca[13C]O2[18O](s) -4.17* 3.54 7.71 Ca[13C]O2[18O] Delta_V 36.14 cm3/mol Ca[13C]O3(s) -1.96* 6.24 8.20 Ca[13C]O3 Delta_V 54.21 cm3/mol Ca[13C]O[18O]2(s) -6.85* 0.84 7.69 Ca[13C]O[18O]2 Delta_V 18.07 cm3/mol CaC[18O]3(s) -8.06 0.10 8.16 CaC[18O]3 CaCO2[18O](s) -2.22* 5.50 7.71 CaCO2[18O] Delta_V 36.14 cm3/mol CaCO[18O]2(s) -4.90* 2.80 7.70 CaCO[18O]2 Delta_V 18.07 cm3/mol Calcite -0.01* -8.49 -8.48 CaCO3 Delta_V 54.21 cm3/mol CD4(g) -35.88* -38.74 -2.86 CD4 Delta_V -216.82 cm3/mol CH2D2(g) -27.49* -31.13 -3.64 CH2D2 Delta_V -198.75 cm3/mol CH3D(g) -23.86* -27.32 -3.46 CH3D Delta_V -189.72 cm3/mol CH4(g) -20.65* -23.51 -2.86 CH4 Delta_V -180.69 cm3/mol CHD3(g) -31.47* -34.94 -3.46 CHD3 Delta_V -207.79 cm3/mol CO2(g) -3.12 -4.59 -1.47 CO2 CO[18O](g) -5.50* -7.28 -1.79 CO[18O] Delta_V -18.07 cm3/mol D2(g) -18.45* -21.60 -3.15 D2 Delta_V -54.21 cm3/mol D2[18O](g) -11.89* -10.31 1.58 D2[18O] Delta_V -18.07 cm3/mol D2O(g) -9.19 -7.62 1.58 D2O H2(g) -10.83* -13.98 -3.15 H2 Delta_V -36.14 cm3/mol H2[18O](g) -4.21 -2.70 1.51 H2[18O] H2O(g) -1.51* -0.00 1.51 H2O Delta_V 18.07 cm3/mol HD(g) -14.34* -17.79 -3.45 HD Delta_V -45.17 cm3/mol HD[18O](g) -7.75* -6.51 1.25 HD[18O] Delta_V -9.03 cm3/mol HDO(g) -5.05* -3.81 1.24 HDO Delta_V 9.03 cm3/mol O2(g) -61.53* -64.42 -2.89 O2 Delta_V 108.41 cm3/mol O[18O](g) -64.23* -67.12 -2.89 O[18O] Delta_V 90.34 cm3/mol * with Delta_V * (P - 1) / 2.3RT. ------------------ End of simulation. ------------------ ------------------------------------ Reading input data for simulation 4. ------------------------------------ ----------- End of run. -----------