PyPI - pyEQL - Versions diffs - 1.4.0rc9__cp310-cp310-win_amd64.whl - Mend

pyEQL 1.4.0rc9__cp310-cp310-win_amd64.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (519) hide show

pyEQL/phreeqc/database/phreeqc_rates.dat ADDED Viewed

@@ -0,0 +1,3158 @@
+# phreeqc_rates.dat for calculating temperature and pressure dependence of reactions, and the specific conductance and viscosity of the solution. Augmented with kinetic rates for minerals from compilations. Based on:
+#   diffusion coefficients and molal volumina of aqueous species, solubility and volume of minerals, and critical temperatures and pressures of gases in Peng-Robinson's EOS.
+# Details are given at the end of this file.
+SOLUTION_MASTER_SPECIES
+#
+#element  species	 alk	gfw_formula	element_gfw
+#
+H           H+      -1  H               1.008
+H(0)        H2      0   H
+H(1)        H+      -1  H
+E           e-      1   0               0
+O           H2O     0   O               16
+O(0)        O2      0   O
+O(-2)       H2O     0   0
+Ca          Ca+2    0   Ca              40.08
+Mg          Mg+2    0   Mg              24.312
+Na          Na+     0   Na              22.9898
+K           K+      0   K               39.102
+Fe          Fe+2    0   Fe              55.847
+Fe(+2)      Fe+2    0   Fe
+Fe(+3)      Fe+3    -2  Fe
+Mn          Mn+2    0   Mn              54.938
+Mn(+2)      Mn+2    0   Mn
+Mn(+3)      Mn+3    0   Mn
+Al          Al+3    0   Al              26.9815
+Ba          Ba+2    0   Ba              137.34
+Sr          Sr+2    0   Sr              87.62
+Si          H4SiO4  0   SiO2            28.0843
+Cl          Cl-     0   Cl              35.453
+C           CO3-2   2   HCO3            12.0111
+C(+4)       CO3-2   2   HCO3
+C(-4)       CH4     0   CH4
+Alkalinity  CO3-2   1   Ca0.5(CO3)0.5   50.05
+S           SO4-2   0   SO4             32.064
+S(6)        SO4-2   0   SO4
+S(-2)       HS-     1   S
+N           NO3-    0   N               14.0067
+N(+5)       NO3-    0   N
+N(+3)       NO2-    0   N
+N(0)        N2      0   N
+N(-3)       NH4+    0   N               14.0067
+#Amm         AmmH+   0   AmmH           17.031
+B           H3BO3   0   B               10.81
+P           PO4-3   2   P               30.9738
+F           F-      0   F               18.9984
+Li          Li+     0   Li              6.939
+Br          Br-     0   Br              79.904
+Zn          Zn+2    0   Zn              65.37
+Cd          Cd+2    0   Cd              112.4
+Pb          Pb+2    0   Pb              207.19
+Cu          Cu+2    0   Cu              63.546
+Cu(+2)      Cu+2    0   Cu
+Cu(+1)      Cu+1    0   Cu
+#   redox-uncoupled gases
+Hdg         Hdg     0   Hdg             2.016   #   H2  gas
+Oxg         Oxg     0   Oxg             32      #   O2  gas
+Mtg         Mtg     0   Mtg             16.032  #   CH4 gas
+Sg          H2Sg    0   H2Sg            32.064  #   H2S gas
+Ntg         Ntg     0   Ntg             28.0134 #   N2  gas
+SOLUTION_SPECIES
+H+ = H+
+	-gamma 9 0
+	-viscosity 9.35e-2 -8.31e-2 2.487e-2 4.49e-4 2.01e-2 1.57 # for viscosity parameters see ref. 4
+	-dw 9.31e-9 838 6.96 -2.285 0.206 24.01 0
+#	   Dw(25 C) dw_T a    a2    visc  a3    a_v_dif
+# Dw(TK) = 9.31e-9 * exp(838 / TK - 838 / 298.15) * viscos_0_25 / viscos_0_tc
+# a = DH ion size, a2 = exponent, visc = viscosity exponent, a3(H+) = 24.01 = new dw calculation from A.D. 2024, a_v_dif = exponent in (viscos_0_tc / viscos)^a_v_dif for tracer diffusion.
+# For SC, Dw(TK) *=  (viscos_0_tc / viscos)^visc (visc = 0.206 for H+)
+# a3 > 5 or a3 = 0 or not defined ? ka = DH_B * a * (1 + (vm - v0))^a2 * mu^0.5, in Onsager-Falkenhagen eqn. (For H+, the reference ion, vm = v0 = 0, a *= (1 + mu)^a2.)
+# a3 = -10 ? ka = DH_B * a * mu^a2  (Define a3 = -10, not used in this database.) (a3 = 24.01 for H+, a flag.)
+# -3 < a3 < 4 ? ka = DH_B * a2 * mu^0.5 / (1 + mu^a3), Appelo, 2017: Dw(I) = Dw(TK) * exp(-a * DH_A * z * sqrt_mu / (1 + ka)) (Sr+2 in this database)
+# If a_v_dif <> 0, Dw(TK) *= (viscos_0_tc / viscos)^a_v_dif in TRANSPORT.
+e- = e-
+H2O = H2O
+	-dw 2.299e-9 -254
+# H2O + 0.01e- = H2O-0.01; -log_k -9 # aids convergence
+Li+ = Li+
+	-gamma 6 0 # The apparent volume parameters are defined in ref. 1 & 2
+	-Vm -0.419 -0.069 13.16 -2.78 0.416 0 0.296 -12.4 -2.74e-3 1.26 # ref. 2 and Ellis, 1968, J. Chem. Soc. A, 1138
+	-viscosity 0.162 -2.45e-2 3.73e-2 9.7e-4 8.1e-4 2.087 # < 10 M LiCl
+	-dw 1.03e-9 -14 4.03 0.8341 1.679
+Na+ = Na+
+	-gamma 4 0.075
+	-gamma 4.08 0.082 # halite solubility
+	-Vm 2.28 -4.38 -4.1 -0.586 0.09 4 0.3 52 -3.33e-3 0.566
+	# -Vm   2.28  -4.38  -4.1  -0.586  0.09  4  0.3  52  -3.33e-3  0.45 # for densities (rho) when I > 3.
+	-viscosity 0.1387 -8.66e-2 1.25e-2 1.45e-2 7.5e-3 1.062
+	-dw 1.33e-9 75 3.627 0 0.7037
+K+ = K+
+	-gamma 3.5 0.015
+	-Vm 3.322 -1.473 6.534 -2.712 9.06e-2 3.5 0 29.7 0 1
+	-viscosity 0.116 -0.191 1.52e-2 1.4e-2 2.59e-2 0.9028
+	-dw 1.96e-9 254 3.484 0 0.1964
+Mg+2 = Mg+2
+	-gamma 5.5 0.2
+	-Vm -1.41 -8.6 11.13 -2.39 1.332 5.5 1.29 -32.9 -5.86e-3 1
+	-viscosity 0.426 0 0 1.66e-3 4.32e-3 2.461
+	-dw 0.705e-9 -4 5.569 0 1.047
+Ca+2 = Ca+2
+	-gamma 5 0.165
+	-Vm -0.3456 -7.252 6.149 -2.479 1.239 5 1.6 -57.1 -6.12e-3 1
+	-viscosity 0.359 -0.158 4.2e-2 1.5e-3 8.04e-3 2.3 # ref. 4, CaCl2 < 6 M
+	-dw 0.792e-9 34 5.411 0 1.046
+Sr+2 = Sr+2
+	-gamma 5.26 0.121
+	-Vm -1.57e-2 -10.15 10.18 -2.36 0.86 5.26 0.859 -27 -4.1e-3 1.97
+	-viscosity 0.472 -0.252 5.51e-3 3.67e-3 0 1.876
+	-dw 0.794e-9 149 0.805 1.961 1e-9 0.7876
+Ba+2 = Ba+2
+	-gamma 5 0
+	-gamma 4 0.153 # Barite solubility
+	-Vm 2.063 -10.06 1.9534 -2.36 0.4218 5 1.58 -12.03 -8.35e-3 1
+	-viscosity 0.338 -0.227 1.39e-2 3.07e-2 0 0.768
+	-dw 0.848e-9 174 10.53 0 3
+Fe+2 = Fe+2
+	-gamma 6 0
+	-Vm -0.3255 -9.687 1.536 -2.379 0.3033 6 -4.21e-2 39.7 0 1
+	-dw 0.719e-9
+Mn+2 = Mn+2
+	-gamma 6 0
+	-Vm -1.1 -8.03 4.08 -2.45 1.4 6 8.07 0 -1.51e-2 0.118
+	-dw 0.688e-9
+Al+3 = Al+3
+	-gamma 9 0
+	-Vm -2.28 -17.1 10.9 -2.07 2.87 9 0 0 5.5e-3 1 # ref. 2 and Barta and Hepler, 1986, Can. J.C. 64, 353
+	-dw 0.559e-9
+H4SiO4 = H4SiO4
+	-Vm 10.5 1.7 20 -2.7 0.1291 # supcrt  2*H2O in a1
+	-dw 1.1e-9
+Cl- = Cl-
+	-gamma 3.5 0.015
+	-gamma 3.63 0.017 # cf. pitzer.dat
+	-Vm 4.465 4.801 4.325 -2.847 1.748 0 -0.331 20.16 0 1
+	-viscosity 0 0 0 0 0 0 1 # the reference solute
+	-dw 2.033e-9 216 3.16 0.2071 0.7432
+CO3-2 = CO3-2
+	-gamma 5.4 0
+	-Vm 6.09 -2.78 -0.405 -5.3 5.02 0 0.169 101 -1.38e-2 0.9316
+	-viscosity -0.5 0.6521 5.44e-3 1.06e-3 -2.18e-2 1.208 -2.147
+	-dw 0.955e-9 -103 2.246 7.13e-2 0.3686
+SO4-2 = SO4-2
+	-gamma 5 -0.04
+	-Vm -7.77 43.17 176 -51.45 3.794 0 42.99 -541 -0.145 0.45 # with analytical_expressions for log K of NaSO4-, KSO4- & MgSO4, 0 - 200 oC
+	-viscosity -0.3 0.501 2.57e-3 0.195 3.14e-2 2.015 0.605
+	-dw 1.07e-9 -114 17 6.02e-2 4.94e-2
+NO3- = NO3-
+	-gamma 3 0
+	-Vm 6.32 6.78 0 -3.06 0.346 0 0.93 0 -0.012 1
+	-viscosity 8.37e-2 -0.458 1.54e-2 0.34 1.79e-2 5.02e-2 0.7381
+	-dw 1.9e-9 104 1.11
+# AmmH+ = AmmH+
+	# -gamma 2.50
+	# -Vm 5.35 2.345 3.72 -2.88 1.55 2.5 -4.54 217 2.344e-2 0.569
+	# -viscosity 9.9e-2 -0.159 1.36e-2 6.51e-3 3.21e-2 0.972
+	# -dw 1.98e-9 203 1.47 2.644 6.81e-2
+H3BO3 = H3BO3
+	-Vm 7.0643 8.8547 3.5844 -3.1451 -0.2 # supcrt
+	-dw 1.1e-9
+PO4-3 = PO4-3
+	-gamma 4 0
+	-Vm 1.24 -9.07 9.31 -2.4 5.61 0 0 0 -1.41e-2 1
+	-dw 0.612e-9
+F- = F-
+	-gamma 3.5 0
+	-Vm 0.928 1.36 6.27 -2.84 1.84 0 0 -0.318 0 1
+	-viscosity 0 2.85e-2 1.35e-2 6.11e-2 4.38e-3 1.384 0.586
+	-dw 1.46e-9 -36 4.352
+Br- = Br-
+	-gamma 3 0
+	-Vm 6.72 2.85 4.21 -3.14 1.38 0 -9.56e-2 7.08 -1.56e-3 1
+	-viscosity -6.98e-2 -0.141 1.78e-2 0.159 7.76e-3 6.25e-2 0.859
+	-dw 2.09e-9 208 3.5 0 0.5737
+Zn+2 = Zn+2
+	-gamma 5 0
+	-Vm -1.96 -10.4 14.3 -2.35 1.46 5 -1.43 24 1.67e-2 1.11
+	-dw 0.715e-9
+Cd+2 = Cd+2
+	-Vm 1.63 -10.7 1.01 -2.34 1.47 5 0 0 0 1
+	-dw 0.717e-9
+Pb+2 = Pb+2
+	-Vm -0.0051 -7.7939 8.8134 -2.4568 1.0788 4.5 # supcrt
+	-dw 0.945e-9
+Cu+2 = Cu+2
+	-gamma 6 0
+	-Vm -1.13 -10.5 7.29 -2.35 1.61 6 9.78e-2 0 3.42e-3 1
+	-dw 0.733e-9
+# redox-uncoupled gases
+Hdg = Hdg # H2
+	-Vm 6.52 0.78 0.12 # supcrt
+	-dw 5.13e-9
+Oxg = Oxg # O2
+	-Vm 5.7889 6.3536 3.2528 -3.0417 -0.3943 # supcrt
+	-dw 2.35e-9
+Mtg = Mtg # CH4
+	-Vm 9.01 -1.11 0 -1.85 -1.5 # Hnedkovsky et al., 1996, JCT 28, 125
+	-dw 1.85e-9
+Ntg = Ntg # N2
+	-Vm 7 # Pray et al., 1952, IEC 44, 1146
+	-dw 1.96e-9 -90 # Cadogan et al. 2014, JCED 59, 519
+H2Sg = H2Sg # H2S
+	-Vm 1.39 28.3 0 -7.22 -0.59 # Hnedkovsky et al., 1996, JCT 28, 125
+	-dw 2.1e-9
+# aqueous species
+H2O = OH- + H+
+	-analytic 293.29227 0.1360833 -10576.913 -123.73158 0 -6.996455e-5
+	-gamma 3.5 0
+	-Vm -9.66 28.5 80 -22.9 1.89 0 1.09 0 0 1
+	-viscosity -2.26e-2 0.106 2.184e-2  -3.2e-3 0 0.4082 -1.634 # < 5 M Li,Na,KOH
+	-dw 5.27e-9 478 0.8695
+2 H2O = O2 + 4 H+ + 4 e-
+	-log_k -86.08
+	-delta_h 134.79 kcal
+	-Vm 5.7889 6.3536 3.2528 -3.0417 -0.3943 # supcrt
+	-dw 2.35e-9
+2 H+ + 2 e- = H2
+	-log_k -3.15
+	-delta_h -1.759 kcal
+	-Vm 6.52 0.78 0.12 # supcrt
+	-dw 5.13e-9
+H+ + Cl- = HCl
+	-log_k -0.5
+	-analytical_expression 0.334 -2.684e-3 1.015 # from Pitzer.dat, up to 15 M HCl, 0 - 50�C
+	-gamma 0 0.4256
+	-viscosity 0.921 -0.765 8.32e-3 8.25e-4 2.53e-3 4.223
+CO3-2 + H+ = HCO3-
+	-log_k 10.329;  -delta_h -3.561 kcal
+	-analytic 107.8871 0.03252849 -5151.79 -38.92561 563713.9
+	-gamma 5.4 0
+	-Vm 10.26 -2.92 -12.58 -0.241 2.23 0 -5.49 320 2.83e-2 1.144
+	-viscosity -0.6 1.366 -1.216e-2 0e-2 3.139e-2 -1.135 1.253
+	-dw 1.18e-9 -190 11.386
+CO3-2 + 2 H+ = CO2 + H2O
+	-log_k 16.681
+	-delta_h -5.738 kcal
+	-analytic 464.1965 0.09344813 -26986.16 -165.75951 2248628.9
+	-Vm 7.29 0.92 2.07 -1.23 -1.6 # McBride et al. 2015, JCED 60, 171
+	-gamma 0 0.066 # Rumpf et al. 1994, J. Sol. Chem. 23, 431
+	-viscosity 6.8e-3 9.03e-2 3.27e-2 0 0 0 0.18
+	-dw 1.92e-9 -120 # TK dependence from Cadogan et al. 2014, , JCED 59, 519
+2 CO2 = (CO2)2 # activity correction for CO2 solubility at high P, T
+	-log_k -1.8
+	-analytical_expression 8.68 -0.0103 -2190
+	-Vm 14.58 1.84 4.14 -2.46 -3.2
+	-viscosity 1.36e-2 0.1806 3.27e-2 0 0 0 0.36
+	-dw 1.92e-9 -120 # TK dependence from Cadogan et al. 2014, , JCED 59, 519
+CO3-2 + 10 H+ + 8 e- = CH4 + 3 H2O
+	-log_k 41.071
+	-delta_h -61.039 kcal
+	-Vm 9.01 -1.11 0 -1.85 -1.5 # Hnedkovsky et al., 1996, JCT 28, 125
+	-dw 1.85e-9
+SO4-2 + H+ = HSO4-
+	-log_k 1.988;  -delta_h 3.85 kcal
+	-analytic -56.889 0.006473 2307.9 19.8858
+	-Vm 8.2 9.259 2.1108 -3.1618 1.1748 0 -0.3 15 0 1
+	-viscosity 0.5 -6.97e-2 6.07e-2 1e-5 -0.1333 0.4865 0.7987
+	-dw 1.22e-9 1000 15 2.861
+HS- = S-2 + H+
+	-log_k -12.918
+	-delta_h 12.1 kcal
+	-gamma 5 0
+	-dw 0.731e-9
+SO4-2 + 9 H+ + 8 e- = HS- + 4 H2O
+	-log_k 33.65
+	-delta_h -60.14 kcal
+	-gamma 3.5 0
+	-Vm 5.0119 4.9799 3.4765 -2.9849 1.441 # supcrt
+	-dw 1.73e-9
+HS- + H+ = H2S
+	-log_k 6.994; -delta_h -5.3 kcal
+	-analytical -11.17 0.02386 3279
+	-Vm 1.39 28.3 0 -7.22 -0.59 # Hnedkovsky et al., 1996, JCT 28, 125
+	-dw 2.1e-9
+2 H2S = (H2S)2 # activity correction for H2S solubility at high P, T
+	-analytical_expression 10.227 -0.01384 -2200
+	-Vm 36.41 -71.95 0 0 2.58
+	-dw 2.1e-9
+H2Sg = HSg- + H+
+	-log_k -6.994; -delta_h 5.3 kcal
+	-analytical_expression 11.17 -0.02386 -3279
+	-gamma 3.5 0
+	-Vm 5.0119 4.9799 3.4765 -2.9849 1.441 # supcrt
+	-dw 1.73e-9
+2 H2Sg = (H2Sg)2 # activity correction for H2S solubility at high P, T
+	-analytical_expression 10.227 -0.01384 -2200
+	-Vm 36.41 -71.95 0 0 2.58
+	-dw 2.1e-9
+NO3- + 2 H+ + 2 e- = NO2- + H2O
+	-log_k 28.57
+	-delta_h -43.76 kcal
+	-gamma 3 0
+	-Vm 5.5864 5.859 3.4472 -3.0212 1.1847 # supcrt
+	-dw 1.91e-9
+2 NO3- + 12 H+ + 10 e- = N2 + 6 H2O
+	-log_k 207.08
+	-delta_h -312.13 kcal
+	-Vm 7 # Pray et al., 1952, IEC 44 1146
+	-dw 1.96e-9 -90 # Cadogan et al. 2014, JCED 59, 519
+NO3- + 10 H+ + 8 e- = NH4+ + 3 H2O
+	-log_k 119.077
+	-delta_h -187.055 kcal
+	-gamma 2.5 0
+	-Vm 5.35 2.345 3.72 -2.88 1.55 2.5 -4.54 217 2.344e-2 0.569
+	-viscosity 9.9e-2 -0.159 1.36e-2 6.51e-3 3.21e-2 0.972
+	-dw 1.98e-9 203 1.47 2.644 6.81e-2
+#AmmH+ = Amm + H+
+NH4+ = NH3 + H+
+	-log_k -9.252
+	-delta_h 12.48 kcal
+	-analytic 0.6322 -0.001225 -2835.76
+	-Vm 6.69 2.8 3.58 -2.88 1.43
+	-viscosity 0.08 0 0 7.82e-3 -0.134 -0.986
+	-dw 2.28e-9
+#AmmH+ + SO4-2 = AmmHSO4-
+NH4+ + SO4-2 = NH4SO4-
+	-gamma 2.08 -0.0416
+	-log_k	1.211; -delta_h 8.56 kJ
+	-Vm -8.78 0 -36.09 0 -8.60 0 87.62 0 -0.3123 0.1172
+	-viscosity 0 0.116 -8.6e-3 0.159 -9.3e-3 0.522 0.627
+	-dw 0.9e-9 100 2.1 2 0
+H3BO3 = H2BO3- + H+
+	-log_k -9.24
+	-delta_h 3.224 kcal
+H3BO3 + F- = BF(OH)3-
+	-log_k -0.4
+	-delta_h 1.85 kcal
+H3BO3 + 2 F- + H+ = BF2(OH)2- + H2O
+	-log_k 7.63
+	-delta_h 1.618 kcal
+H3BO3 + 2 H+ + 3 F- = BF3OH- + 2 H2O
+	-log_k 13.67
+	-delta_h -1.614 kcal
+H3BO3 + 3 H+ + 4 F- = BF4- + 3 H2O
+	-log_k 20.28
+	-delta_h -1.846 kcal
+PO4-3 + H+ = HPO4-2
+	-log_k 12.346
+	-delta_h -3.53 kcal
+	-gamma 5 0
+	-dw 0.69e-9
+	-Vm 3.52 1.09 8.39 -2.82 3.34 0 0 0 0 1
+PO4-3 + 2 H+ = H2PO4-
+	-log_k 19.553
+	-delta_h -4.52 kcal
+	-gamma 5.4 0
+	-Vm 5.58 8.06 12.2 -3.11 1.3 0 0 0 1.62e-2 1
+	-dw 0.846e-9
+PO4-3 + 3 H+ = H3PO4
+	log_k 21.721 # log_k and delta_h from minteq.v4.dat, NIST46.3
+	delta_h -10.1 kJ
+	-Vm 7.47 12.4 6.29 -3.29 0
+H+ + F- = HF
+	-log_k 3.18
+	-delta_h 3.18 kcal
+	-analytic -2.033 0.012645 429.01
+	-Vm 3.4753 .7042 5.4732 -2.8081 -.0007 # supcrt
+H+ + 2 F- = HF2-
+	-log_k 3.76
+	-delta_h 4.55 kcal
+	-Vm 5.2263 4.9797 3.7928 -2.9849 1.2934 # supcrt
+Ca+2 + H2O = CaOH+ + H+
+	-log_k -12.78
+Ca+2 + CO3-2 = CaCO3
+	-log_k 3.224; -delta_h 3.545 kcal
+	-analytic -1228.732 -0.29944 35512.75 485.818
+	-dw 4.46e-10 # complexes: calc'd with the Pikal formula
+	-Vm -.243 -8.3748 9.0417 -2.4328 -.03 # supcrt
+Ca+2 + CO3-2 + H+ = CaHCO3+
+	-log_k 10.91; -delta_h 4.38 kcal
+	-analytic -6.009 3.377e-2 2044
+	-gamma 6 0
+	-Vm 30.19 .01 5.75 -2.78 .308 5.4
+	-dw 5.06e-10
+Ca+2 + SO4-2 = CaSO4
+	-log_k 2.25
+	-delta_h 1.325 kcal
+	-dw 4.71e-10
+	-Vm 2.791 -.9666 6.13 -2.739 -.001 # supcrt
+Ca+2 + HSO4- = CaHSO4+
+	-log_k 1.08
+Ca+2 + PO4-3 = CaPO4-
+	-log_k 6.459
+	-delta_h 3.1 kcal
+	-gamma 5.4 0
+Ca+2 + HPO4-2 = CaHPO4
+	-log_k 2.739
+	-delta_h 3.3 kcal
+Ca+2 + H2PO4- = CaH2PO4+
+	-log_k 1.408
+	-delta_h 3.4 kcal
+	-gamma 5.4 0
+# Ca+2 + F- = CaF+
+	# -log_k	0.94
+	# -delta_h 4.120	kcal
+	# -gamma	5.5  0.0
+	# -Vm   .9846  -5.3773  7.8635  -2.5567  .6911 5.5 # supcrt
+Mg+2 + H2O = MgOH+ + H+
+	-log_k -11.44
+	-delta_h 15.952 kcal
+	-gamma 6.5 0
+Mg+2 + CO3-2 = MgCO3
+	-log_k 2.98
+	-delta_h 2.713 kcal
+	-analytic 0.991 0.00667
+	-Vm -0.5837 -9.2067 9.3687 -2.3984 -.03 # supcrt
+	-dw 4.21e-10
+Mg+2 + H+ + CO3-2 = MgHCO3+
+	-log_k 11.399
+	-delta_h -2.771 kcal
+	-analytic 48.6721 0.03252849 -2614.335 -18.00263 563713.9
+	-gamma 4 0
+	-Vm 2.7171 -1.1469 6.2008 -2.7316 .5985 4 # supcrt
+	-dw 4.78e-10
+Mg+2 + SO4-2 = MgSO4
+	-gamma 0 0.2
+	-log_k 2.42; -delta_h 19 kJ
+	-analytical_expression 0 9.64e-3 -136 # mean salt gamma from Pitzer.dat and epsomite/hexahydrite/kieserite solubilities, 0 - 200 oC
+	-Vm 8.65 -10.21 29.58 -18.6 1.061
+	-viscosity 0.318 -5.4e-4 -3.42e-2 0.708 3.7e-3 0.696
+	-dw 4.45e-10
+SO4-2 + MgSO4 = Mg(SO4)2-2
+	-gamma 7 0.047
+	-log_k 0.52; -delta_h -13.6 kJ
+	-analytical_expression 0 -1.51e-3 0 0 8.604e4 # mean salt gamma from Pitzer.dat and epsomite/hexahydrite/kieserite solubilities, 0 - 200 oC
+	-Vm -8.14 -62.2 -15.96 3.29 -3.01 0 150 0 0.153 3.79e-2
+	-viscosity -0.169 5e-4 -5.69e-2 0.11 2.03e-3 2.027 -1e-3
+	-dw 0.845e-9 -200 8 0 0.965
+Mg+2 + PO4-3 = MgPO4-
+	-log_k 6.589
+	-delta_h 3.1 kcal
+	-gamma 5.4 0
+Mg+2 + HPO4-2 = MgHPO4
+	-log_k 2.87
+	-delta_h 3.3 kcal
+Mg+2 + H2PO4- = MgH2PO4+
+	-log_k 1.513
+	-delta_h 3.4 kcal
+	-gamma 5.4 0
+Mg+2 + F- = MgF+
+	-log_k 1.82
+	-delta_h 3.2 kcal
+	-gamma 4.5 0
+	-Vm .6494 -6.1958 8.1852 -2.5229 .9706 4.5 # supcrt
+# Na+ + OH- = NaOH
+	# -log_k -14.7 # remove this complex
+Na+ + HCO3- = NaHCO3
+	-log_k -0.06; -delta_h 21 kJ
+	-gamma 0 0.2
+	-Vm 7.95 0 0 0 0.609
+	-viscosity -4e-2 -2.717 1.67e-5
+	-dw 6.73e-10
+Na+ + SO4-2 = NaSO4-
+	-gamma 5.5 0
+	-log_k 0.6; -delta_h -14.4 kJ
+	-analytical_expression 255.903 0.10057 0 -1.11138e2 -8.5983e5 # mirabilite/thenardite solubilities, 0 - 200 oC
+	-Vm 1.99 -10.78 21.88 -12.7 1.601 5 32.38 501 1.565e-2 0.2325
+	-viscosity 0.2 -5.93e-2 -4e-4 8.46e-3 1.78e-3 2.308 -0.208
+	-dw 1.13e-9 -23 8.5 0.392 0.521
+Na+ + HPO4-2 = NaHPO4-
+	-log_k 0.29
+	-gamma 5.4 0
+	-Vm 5.2 8.1 13 -3 0.9 0 0 1.62e-2 1
+Na+ + F- = NaF
+	-log_k -0.24
+	-Vm 2.7483 -1.0708 6.1709 -2.7347 -.03 # supcrt
+K+ + HCO3- = KHCO3
+	-log_k -0.35;  -delta_h 12 kJ
+	-gamma 0 9.4e-3
+	-Vm 9.48 0 0 0 -0.542
+	-viscosity 0.7 -1.289 9e-2
+K+ + SO4-2 = KSO4-
+	-gamma 5.4 0.19
+	-log_k 0.6; -delta_h -10.4 kJ
+	-analytical_expression -3.0246 9.986e-3 0 0 1.093e5 # arcanite solubility, 0 - 200 oC
+	-Vm 13.48 -18.03 61.74 -19.6 2.046 5.4 -17.32 0 0.1522 1.919
+	-viscosity -1 1.06 1e-4 -0.464 3.78e-2 0.539 -0.69
+	-dw 0.9e-9 63 8.48 0 1.8
+K+ + HPO4-2 = KHPO4-
+	-log_k 0.29
+	-gamma 5.4 0
+	-Vm 5.4 8.1 19 -3.1 0.7 0 0 0 1.62e-2 1
+Fe+2 + H2O = FeOH+ + H+
+	-log_k -9.5
+	-delta_h 13.2 kcal
+	-gamma 5 0
+Fe+2 + 3 H2O = Fe(OH)3- + 3 H+
+	-log_k -31
+	-delta_h 30.3 kcal
+	-gamma 5 0
+Fe+2 + Cl- = FeCl+
+	-log_k 0.14
+Fe+2 + CO3-2 = FeCO3
+	-log_k 4.38
+Fe+2 + HCO3- = FeHCO3+
+	-log_k 2
+Fe+2 + SO4-2 = FeSO4
+	-log_k 2.25
+	-delta_h 3.23 kcal
+	-Vm -13 0 123
+Fe+2 + HSO4- = FeHSO4+
+	-log_k 1.08
+Fe+2 + 2 HS- = Fe(HS)2
+	-log_k 8.95
+Fe+2 + 3 HS- = Fe(HS)3-
+	-log_k 10.987
+Fe+2 + HPO4-2 = FeHPO4
+	-log_k 3.6
+Fe+2 + H2PO4- = FeH2PO4+
+	-log_k 2.7
+	-gamma 5.4 0
+Fe+2 + F- = FeF+
+	-log_k 1
+Fe+2 = Fe+3 + e-
+	-log_k -13.02
+	-delta_h 9.68 kcal
+	-gamma 9 0
+Fe+3 + H2O = FeOH+2 + H+
+	-log_k -2.19
+	-delta_h 10.4 kcal
+	-gamma 5 0
+Fe+3 + 2 H2O = Fe(OH)2+ + 2 H+
+	-log_k -5.67
+	-delta_h 17.1 kcal
+	-gamma 5.4 0
+Fe+3 + 3 H2O = Fe(OH)3 + 3 H+
+	-log_k -12.56
+	-delta_h 24.8 kcal
+Fe+3 + 4 H2O = Fe(OH)4- + 4 H+
+	-log_k -21.6
+	-delta_h 31.9 kcal
+	-gamma 5.4 0
+Fe+2 + 2 H2O = Fe(OH)2 + 2 H+
+	-log_k -20.57
+	-delta_h 28.565 kcal
+2 Fe+3 + 2 H2O = Fe2(OH)2+4 + 2 H+
+	-log_k -2.95
+	-delta_h 13.5 kcal
+3 Fe+3 + 4 H2O = Fe3(OH)4+5 + 4 H+
+	-log_k -6.3
+	-delta_h 14.3 kcal
+Fe+3 + Cl- = FeCl+2
+	-log_k 1.48
+	-delta_h 5.6 kcal
+	-gamma 5 0
+Fe+3 + 2 Cl- = FeCl2+
+	-log_k 2.13
+	-gamma 5 0
+Fe+3 + 3 Cl- = FeCl3
+	-log_k 1.13
+Fe+3 + SO4-2 = FeSO4+
+	-log_k 4.04
+	-delta_h 3.91 kcal
+	-gamma 5 0
+Fe+3 + HSO4- = FeHSO4+2
+	-log_k 2.48
+Fe+3 + 2 SO4-2 = Fe(SO4)2-
+	-log_k 5.38
+	-delta_h 4.6 kcal
+Fe+3 + HPO4-2 = FeHPO4+
+	-log_k 5.43
+	-delta_h 5.76 kcal
+	-gamma 5 0
+Fe+3 + H2PO4- = FeH2PO4+2
+	-log_k 5.43
+	-gamma 5.4 0
+Fe+3 + F- = FeF+2
+	-log_k 6.2
+	-delta_h 2.7 kcal
+	-gamma 5 0
+Fe+3 + 2 F- = FeF2+
+	-log_k 10.8
+	-delta_h 4.8 kcal
+	-gamma 5 0
+Fe+3 + 3 F- = FeF3
+	-log_k 14
+	-delta_h 5.4 kcal
+Mn+2 + H2O = MnOH+ + H+
+	-log_k -10.59
+	-delta_h 14.4 kcal
+	-gamma 5 0
+Mn+2 + 3 H2O = Mn(OH)3- + 3 H+
+	-log_k -34.8
+	-gamma 5 0
+Mn+2 + Cl- = MnCl+
+	-log_k 0.61
+	-gamma 5 0
+	-Vm 7.25 -1.08 -25.8 -2.73 3.99 5 0 0 0 1
+Mn+2 + 2 Cl- = MnCl2
+	-log_k 0.25
+	-Vm 1e-5 0 144
+Mn+2 + 3 Cl- = MnCl3-
+	-log_k -0.31
+	-gamma 5 0
+	-Vm 11.8 0 0 0 2.4 0 0 0 3.6e-2 1
+Mn+2 + CO3-2 = MnCO3
+	-log_k 4.9
+Mn+2 + HCO3- = MnHCO3+
+	-log_k 1.95
+	-gamma 5 0
+Mn+2 + SO4-2 = MnSO4
+	-log_k 2.25
+	-delta_h 3.37 kcal
+	-Vm -1.31 -1.83 62.3 -2.7
+Mn+2 + 2 NO3- = Mn(NO3)2
+	-log_k 0.6
+	-delta_h -0.396 kcal
+	-Vm 6.16 0 29.4 0 0.9
+Mn+2 + F- = MnF+
+	-log_k 0.84
+	-gamma 5 0
+Mn+2 = Mn+3 + e-
+	-log_k -25.51
+	-delta_h 25.8 kcal
+	-gamma 9 0
+Al+3 + H2O = AlOH+2 + H+
+	-log_k -5
+	-delta_h 11.49 kcal
+	-analytic -38.253 0 -656.27 14.327
+	-gamma 5.4 0
+	-Vm -1.46 -11.4 10.2 -2.31 1.67 5.4 0 0 0 1 # Barta and Hepler, 1986, Can. J. Chem. 64, 353
+Al+3 + 2 H2O = Al(OH)2+ + 2 H+
+	-log_k -10.1
+	-delta_h 26.9 kcal
+	-gamma 5.4 0
+	-analytic 88.5 0 -9391.6 -27.121
+Al+3 + 3 H2O = Al(OH)3 + 3 H+
+	-log_k -16.9
+	-delta_h 39.89 kcal
+	-analytic 226.374 0 -18247.8 -73.597
+Al+3 + 4 H2O = Al(OH)4- + 4 H+
+	-log_k -22.7
+	-delta_h 42.3 kcal
+	-analytic 51.578 0 -11168.9 -14.865
+	-gamma 4.5 0
+	-dw 1.04e-9 # Mackin & Aller, 1983, GCA 47, 959
+Al+3 + SO4-2 = AlSO4+
+	-log_k 3.5
+	-delta_h 2.29 kcal
+	-gamma 4.5 0
+Al+3 + 2 SO4-2 = Al(SO4)2-
+	-log_k 5
+	-delta_h 3.11 kcal
+	-gamma 4.5 0
+Al+3 + HSO4- = AlHSO4+2
+	-log_k 0.46
+Al+3 + F- = AlF+2
+	-log_k 7
+	-delta_h 1.06 kcal
+	-gamma 5.4 0
+Al+3 + 2 F- = AlF2+
+	-log_k 12.7
+	-delta_h 1.98 kcal
+	-gamma 5.4 0
+Al+3 + 3 F- = AlF3
+	-log_k 16.8
+	-delta_h 2.16 kcal
+Al+3 + 4 F- = AlF4-
+	-log_k 19.4
+	-delta_h 2.2 kcal
+	-gamma 4.5 0
+# Al+3 + 5 F- = AlF5-2
+	# log_k 20.6
+	# delta_h 1.840 kcal
+# Al+3 + 6 F- = AlF6-3
+	# log_k 20.6
+	# delta_h -1.670	kcal
+H4SiO4 = H3SiO4- + H+
+	-log_k -9.83
+	-delta_h 6.12 kcal
+	-analytic -302.3724 -0.050698 15669.69 108.18466 -1119669
+	-gamma 4 0
+	-Vm 7.94 1.0881 5.3224 -2.824 1.4767 # supcrt  H2O in a1
+H4SiO4 = H2SiO4-2 + 2 H+
+	-log_k -23
+	-delta_h 17.6 kcal
+	-analytic -294.0184 -0.07265 11204.49 108.18466 -1119669
+	-gamma 5.4 0
+H4SiO4 + 4 H+ + 6 F- = SiF6-2 + 4 H2O
+	-log_k 30.18
+	-delta_h -16.26 kcal
+	-gamma 5 0
+	-Vm 8.5311 13.0492 .6211 -3.3185 2.7716 # supcrt
+Ba+2 + H2O = BaOH+ + H+
+	-log_k -13.47
+	-gamma 5 0
+Ba+2 + CO3-2 = BaCO3
+	-log_k 2.71
+	-delta_h 3.55 kcal
+	-analytic 0.113 0.008721
+	-Vm .2907 -7.0717 8.5295 -2.4867 -.03 # supcrt
+Ba+2 + HCO3- = BaHCO3+
+	-log_k 0.982
+	-delta_h 5.56 kcal
+	-analytic -3.0938 0.013669
+Ba+2 + SO4-2 = BaSO4
+	-log_k 2.7
+Sr+2 + H2O = SrOH+ + H+
+	-log_k -13.29
+	-gamma 5 0
+Sr+2 + CO3-2 + H+ = SrHCO3+
+	-log_k 11.509
+	-delta_h 2.489 kcal
+	-analytic 104.6391 0.04739549 -5151.79 -38.92561 563713.9
+	-gamma 5.4 0
+Sr+2 + CO3-2 = SrCO3
+	-log_k 2.81
+	-delta_h 5.22 kcal
+	-analytic -1.019 0.012826
+	-Vm -.1787 -8.2177 8.9799 -2.4393 -.03 # supcrt
+Sr+2 + SO4-2 = SrSO4
+	-log_k 2.29
+	-delta_h 2.08 kcal
+	-Vm 6.791 -.9666 6.13 -2.739 -.001 # celestite solubility
+Li+ + SO4-2 = LiSO4-
+	-log_k 0.64
+	-gamma 5 0
+Cu+2 + e- = Cu+
+	-log_k 2.72
+	-delta_h 1.65 kcal
+	-gamma 2.5 0
+Cu+ + 2 Cl- = CuCl2-
+	-log_k 5.5
+	-delta_h -0.42 kcal
+	-gamma 4 0
+Cu+ + 3 Cl- = CuCl3-2
+	-log_k 5.7
+	-delta_h 0.26 kcal
+	-gamma 5 0
+Cu+2 + CO3-2 = CuCO3
+	-log_k 6.73
+Cu+2 + 2 CO3-2 = Cu(CO3)2-2
+	-log_k 9.83
+Cu+2 + HCO3- = CuHCO3+
+	-log_k 2.7
+Cu+2 + Cl- = CuCl+
+	-log_k 0.43
+	-delta_h 8.65 kcal
+	-gamma 4 0
+	-Vm -4.19 0 30.4 0 0 4 0 0 1.94e-2 1
+Cu+2 + 2 Cl- = CuCl2
+	-log_k 0.16
+	-delta_h 10.56 kcal
+	-Vm 26.8 0 -136
+Cu+2 + 3 Cl- = CuCl3-
+	-log_k -2.29
+	-delta_h 13.69 kcal
+	-gamma 4 0
+Cu+2 + 4 Cl- = CuCl4-2
+	-log_k -4.59
+	-delta_h 17.78 kcal
+	-gamma 5 0
+Cu+2 + F- = CuF+
+	-log_k 1.26
+	-delta_h 1.62 kcal
+Cu+2 + H2O = CuOH+ + H+
+	-log_k -8
+	-gamma 4 0
+Cu+2 + 2 H2O = Cu(OH)2 + 2 H+
+	-log_k -13.68
+Cu+2 + 3 H2O = Cu(OH)3- + 3 H+
+	-log_k -26.9
+Cu+2 + 4 H2O = Cu(OH)4-2 + 4 H+
+	-log_k -39.6
+2 Cu+2 + 2 H2O = Cu2(OH)2+2 + 2 H+
+	-log_k -10.359
+	-delta_h 17.539 kcal
+	-analytical 2.497 0 -3833
+Cu+2 + SO4-2 = CuSO4
+	-log_k 2.31
+	-delta_h 1.22 kcal
+	-Vm 5.21 0 -14.6
+Cu+2 + 3 HS- = Cu(HS)3-
+	-log_k 25.9
+Zn+2 + H2O = ZnOH+ + H+
+	-log_k -8.96
+	-delta_h 13.4 kcal
+Zn+2 + 2 H2O = Zn(OH)2 + 2 H+
+	-log_k -16.9
+Zn+2 + 3 H2O = Zn(OH)3- + 3 H+
+	-log_k -28.4
+Zn+2 + 4 H2O = Zn(OH)4-2 + 4 H+
+	-log_k -41.2
+Zn+2 + Cl- = ZnCl+
+	-log_k 0.43
+	-delta_h 7.79 kcal
+	-gamma 4 0
+	-Vm 14.8 -3.91 -105.7 -2.62 0.203 4 0 0 -5.05e-2 1
+Zn+2 + 2 Cl- = ZnCl2
+	-log_k 0.45
+	-delta_h 8.5 kcal
+	-Vm -10.1 4.57 241 -2.97 -1e-3
+Zn+2 + 3 Cl- = ZnCl3-
+	-log_k 0.5
+	-delta_h 9.56 kcal
+	-gamma 4 0
+	-Vm 0.772 15.5 -0.349 -3.42 1.25 0 -7.77 0 0 1
+Zn+2 + 4 Cl- = ZnCl4-2
+	-log_k 0.2
+	-delta_h 10.96 kcal
+	-gamma 5 0
+	-Vm 28.42 28 -5.26 -3.94 2.67 0 0 0 4.62e-2 1
+Zn+2 + H2O + Cl- = ZnOHCl + H+
+	-log_k -7.48
+Zn+2 + 2 HS- = Zn(HS)2
+	-log_k 14.94
+Zn+2 + 3 HS- = Zn(HS)3-
+	-log_k 16.1
+Zn+2 + CO3-2 = ZnCO3
+	-log_k 5.3
+Zn+2 + 2 CO3-2 = Zn(CO3)2-2
+	-log_k 9.63
+Zn+2 + HCO3- = ZnHCO3+
+	-log_k 2.1
+Zn+2 + SO4-2 = ZnSO4
+	-log_k 2.37
+	-delta_h 1.36 kcal
+	-Vm 2.51 0 18.8
+Zn+2 + 2 SO4-2 = Zn(SO4)2-2
+	-log_k 3.28
+	-Vm 10.9 0 -98.7 0 0 0 24 0 -0.236 1
+Zn+2 + Br- = ZnBr+
+	-log_k -0.58
+Zn+2 + 2 Br- = ZnBr2
+	-log_k -0.98
+Zn+2 + F- = ZnF+
+	-log_k 1.15
+	-delta_h 2.22 kcal
+Cd+2 + H2O = CdOH+ + H+
+	-log_k -10.08
+	-delta_h 13.1 kcal
+Cd+2 + 2 H2O = Cd(OH)2 + 2 H+
+	-log_k -20.35
+Cd+2 + 3 H2O = Cd(OH)3- + 3 H+
+	-log_k -33.3
+Cd+2 + 4 H2O = Cd(OH)4-2 + 4 H+
+	-log_k -47.35
+2 Cd+2 + H2O = Cd2OH+3 + H+
+	-log_k -9.39
+	-delta_h 10.9 kcal
+Cd+2 + H2O + Cl- = CdOHCl + H+
+	-log_k -7.404
+	-delta_h 4.355 kcal
+Cd+2 + NO3- = CdNO3+
+	-log_k 0.4
+	-delta_h -5.2 kcal
+	-Vm 5.95 0 -1.11 0 2.67 7 0 0 1.53e-2 1
+Cd+2 + Cl- = CdCl+
+	-log_k 1.98
+	-delta_h 0.59 kcal
+	-Vm 5.69 0 -30.2 0 0 6 0 0 0.112 1
+Cd+2 + 2 Cl- = CdCl2
+	-log_k 2.6
+	-delta_h 1.24 kcal
+	-Vm 5.53
+Cd+2 + 3 Cl- = CdCl3-
+	-log_k 2.4
+	-delta_h 3.9 kcal
+	-Vm 4.6 0 83.9 0 0 0 0 0 0 1
+Cd+2 + CO3-2 = CdCO3
+	-log_k 2.9
+Cd+2 + 2 CO3-2 = Cd(CO3)2-2
+	-log_k 6.4
+Cd+2 + HCO3- = CdHCO3+
+	-log_k 1.5
+Cd+2 + SO4-2 = CdSO4
+	-log_k 2.46
+	-delta_h 1.08 kcal
+	-Vm 10.4 0 57.9
+Cd+2 + 2 SO4-2 = Cd(SO4)2-2
+	-log_k 3.5
+	-Vm -6.29 0 -93 0 9.5 7 0 0 0 1
+Cd+2 + Br- = CdBr+
+	-log_k 2.17
+	-delta_h -0.81 kcal
+Cd+2 + 2 Br- = CdBr2
+	-log_k 2.9
+Cd+2 + F- = CdF+
+	-log_k 1.1
+Cd+2 + 2 F- = CdF2
+	-log_k 1.5
+Cd+2 + HS- = CdHS+
+	-log_k 10.17
+Cd+2 + 2 HS- = Cd(HS)2
+	-log_k 16.53
+Cd+2 + 3 HS- = Cd(HS)3-
+	-log_k 18.71
+Cd+2 + 4 HS- = Cd(HS)4-2
+	-log_k 20.9
+Pb+2 + H2O = PbOH+ + H+
+	-log_k -7.71
+Pb+2 + 2 H2O = Pb(OH)2 + 2 H+
+	-log_k -17.12
+Pb+2 + 3 H2O = Pb(OH)3- + 3 H+
+	-log_k -28.06
+Pb+2 + 4 H2O = Pb(OH)4-2 + 4 H+
+	-log_k -39.7
+2 Pb+2 + H2O = Pb2OH+3 + H+
+	-log_k -6.36
+Pb+2 + Cl- = PbCl+
+	-log_k 1.6
+	-delta_h 4.38 kcal
+	-Vm 2.8934 -.7165 6.0316 -2.7494 .1281 6 # supcrt
+Pb+2 + 2 Cl- = PbCl2
+	-log_k 1.8
+	-delta_h 1.08 kcal
+	-Vm 6.5402 8.1879 2.5318 -3.1175 -.03 # supcrt
+Pb+2 + 3 Cl- = PbCl3-
+	-log_k 1.7
+	-delta_h 2.17 kcal
+	-Vm 11.0396 19.1743 -1.7863 -3.5717 .7356 # supcrt
+Pb+2 + 4 Cl- = PbCl4-2
+	-log_k 1.38
+	-delta_h 3.53 kcal
+	-Vm 16.415 32.2997 -6.9452 -4.1143 2.3118 # supcrt
+Pb+2 + CO3-2 = PbCO3
+	-log_k 7.24
+Pb+2 + 2 CO3-2 = Pb(CO3)2-2
+	-log_k 10.64
+Pb+2 + HCO3- = PbHCO3+
+	-log_k 2.9
+Pb+2 + SO4-2 = PbSO4
+	-log_k 2.75
+Pb+2 + 2 SO4-2 = Pb(SO4)2-2
+	-log_k 3.47
+Pb+2 + 2 HS- = Pb(HS)2
+	-log_k 15.27
+Pb+2 + 3 HS- = Pb(HS)3-
+	-log_k 16.57
+3 Pb+2 + 4 H2O = Pb3(OH)4+2 + 4 H+
+	-log_k -23.88
+	-delta_h 26.5 kcal
+Pb+2 + NO3- = PbNO3+
+	-log_k 1.17
+Pb+2 + Br- = PbBr+
+	-log_k 1.77
+	-delta_h 2.88 kcal
+Pb+2 + 2 Br- = PbBr2
+	-log_k 1.44
+Pb+2 + F- = PbF+
+	-log_k 1.25
+Pb+2 + 2 F- = PbF2
+	-log_k 2.56
+Pb+2 + 3 F- = PbF3-
+	-log_k 3.42
+Pb+2 + 4 F- = PbF4-2
+	-log_k 3.1
+PHASES
+Calcite
+	CaCO3 = CO3-2 + Ca+2
+	-log_k -8.48
+	-delta_h -2.297 kcal
+	-analytic 17.118 -0.046528 -3496 # 0 - 250�C, Ellis, 1959, Plummer and Busenberg, 1982
+	-Vm 36.9 cm3/mol # MW (100.09 g/mol) / rho (2.71 g/cm3)
+Aragonite
+	CaCO3 = CO3-2 + Ca+2
+	-log_k -8.336
+	-delta_h -2.589 kcal
+	-analytic -171.9773 -0.077993 2903.293 71.595
+	-Vm 34.04
+Dolomite
+	CaMg(CO3)2 = Ca+2 + Mg+2 + 2 CO3-2
+	-log_k -17.09
+	-delta_h -9.436 kcal
+	-analytic 31.283 -0.0898 -6438 # 25�C: Hemingway and Robie, 1994; 50�175�C: B�n�zeth et al., 2018, GCA 224, 262-275
+	-Vm 64.5
+Siderite
+	FeCO3 = Fe+2 + CO3-2
+	-log_k -10.89
+	-delta_h -2.48 kcal
+	-Vm 29.2
+Rhodochrosite
+	MnCO3 = Mn+2 + CO3-2
+	-log_k -11.13
+	-delta_h -1.43 kcal
+	-Vm 31.1
+Strontianite
+	SrCO3 = Sr+2 + CO3-2
+	-log_k -9.271
+	-delta_h -0.4 kcal
+	-analytic 155.0305 0 -7239.594 -56.58638
+	-Vm 39.69
+Witherite
+	BaCO3 = Ba+2 + CO3-2
+	-log_k -8.562
+	-delta_h 0.703 kcal
+	-analytic 607.642 0.121098 -20011.25 -236.4948
+	-Vm 46
+Gypsum
+	CaSO4:2H2O = Ca+2 + SO4-2 + 2 H2O
+	-log_k -4.58
+	-delta_h -0.109 kcal
+	-analytic 68.2401 0 -3221.51 -25.0627
+	-analytical_expression 93.7 5.99E-3 -4e3 -35.019 # better fits the appendix data of Appelo, 2015, AG 55, 62
+	-Vm 73.9 # 172.18 / 2.33 (Vm H2O = 13.9 cm3/mol)
+Anhydrite
+	CaSO4 = Ca+2 + SO4-2
+	-log_k -4.36
+	-delta_h -1.71 kcal
+	-analytic 84.9 0 -3135.12 -31.79 # 50 - 160oC, 1 - 1e3 atm, anhydrite dissolution, Blount and Dickson, 1973, Am. Mineral. 58, 323
+	-Vm 46.1 # 136.14 / 2.95
+Celestite
+	SrSO4 = Sr+2 + SO4-2
+	-log_k -6.63
+	-delta_h -4.037 kcal
+#	  -analytic	  -14805.9622	-2.4660924	 756968.533	 5436.3588	  -40553604.0
+	-analytic -7.14 6.11e-3 75 0 0 -1.79e-5 # Howell et al., 1992, JCED 37, 464
+	-Vm 46.4
+Barite
+	BaSO4 = Ba+2 + SO4-2
+	-log_k -9.97
+	-delta_h 6.35 kcal
+	-analytical_expression -282.43 -8.972e-2 5822 113.08 # Blount 1977; Templeton, 1960
+	-Vm 52.9
+Arcanite
+	K2SO4 = SO4-2 + 2 K+
+	log_k -1.776; -delta_h 5 kcal
+	-analytical_expression 674.142 0.30423 -18037 -280.236 0 -1.44055e-4 # ref. 3
+	# Note, the Linke and Seidell data may give subsaturation in other xpt's, SI = -0.06
+	-Vm 65.5
+Mirabilite
+	Na2SO4:10H2O = SO4-2 + 2 Na+ + 10 H2O
+	-analytical_expression -301.9326 -0.16232 0 141.078 # ref. 3
+	Vm 216
+Thenardite
+	Na2SO4 = 2 Na+ + SO4-2
+	-analytical_expression 57.185 8.6024e-2 0 -30.8341 0 -7.6905e-5 # ref. 3
+	-Vm 52.9
+Epsomite
+	MgSO4:7H2O = Mg+2 + SO4-2 + 7 H2O
+	log_k -1.74; -delta_h 10.57 kJ
+	-analytical_expression -3.59 6.21e-3
+	Vm 147
+Hexahydrite
+	MgSO4:6H2O = Mg+2 + SO4-2 + 6 H2O
+	log_k -1.57; -delta_h 2.35 kJ
+	-analytical_expression -1.978 1.38e-3
+	Vm 132
+Kieserite
+	MgSO4:H2O = Mg+2 + SO4-2 + H2O
+	log_k -1.16; -delta_h 9.22 kJ
+	-analytical_expression 29.485 -5.07e-2 0 -2.662 -7.95e5
+	Vm 53.8
+Hydroxyapatite
+	Ca5(PO4)3OH + 4 H+ = H2O + 3 HPO4-2 + 5 Ca+2
+	-log_k -3.421
+	-delta_h -36.155 kcal
+	-Vm 128.9
+Fluorite
+	CaF2 = Ca+2 + 2 F-
+	-log_k -10.6
+	-delta_h 4.69 kcal
+	-analytic 66.348 0 -4298.2 -25.271
+	-Vm 15.7
+SiO2(a)
+	SiO2 + 2 H2O = H4SiO4
+	-log_k -2.71
+	-delta_h 3.34 kcal
+	-analytic -0.26 0 -731
+Chalcedony
+	SiO2 + 2 H2O = H4SiO4
+	-log_k -3.55
+	-delta_h 4.72 kcal
+	-analytic -0.09 0 -1032
+	-Vm 23.1
+Quartz
+	SiO2 + 2 H2O = H4SiO4
+	-log_k -3.98
+	-delta_h 5.99 kcal
+	-analytic 0.41 0 -1309
+	-Vm 22.67
+Gibbsite
+	Al(OH)3 + 3 H+ = Al+3 + 3 H2O
+	-log_k 8.11
+	-delta_h -22.8 kcal
+	-Vm 32.22
+Al(OH)3(a)
+	Al(OH)3 + 3 H+ = Al+3 + 3 H2O
+	-log_k 10.8
+	-delta_h -26.5 kcal
+Kaolinite
+	Al2Si2O5(OH)4 + 6 H+ = H2O + 2 H4SiO4 + 2 Al+3
+	-log_k 7.435
+	-delta_h -35.3 kcal
+	-Vm 99.35
+Albite
+	NaAlSi3O8 + 8 H2O = Na+ + Al(OH)4- + 3 H4SiO4
+	-log_k -18.002
+	-delta_h 25.896 kcal
+	-Vm 101.31
+Anorthite
+	CaAl2Si2O8 + 8 H2O = Ca+2 + 2 Al(OH)4- + 2 H4SiO4
+	-log_k -19.714
+	-delta_h 11.58 kcal
+	-Vm 105.05
+K-feldspar
+	KAlSi3O8 + 8 H2O = K+ + Al(OH)4- + 3 H4SiO4
+	-log_k -20.573
+	-delta_h 30.82 kcal
+	-Vm 108.15
+K-mica
+	KAl3Si3O10(OH)2 + 10 H+ = K+ + 3 Al+3 + 3 H4SiO4
+	-log_k 12.703
+	-delta_h -59.376 kcal
+Chlorite(14A)
+	Mg5Al2Si3O10(OH)8 + 16 H+ = 5 Mg+2 + 2 Al+3 + 3 H4SiO4 + 6 H2O
+	-log_k 68.38
+	-delta_h -151.494 kcal
+Ca-Montmorillonite
+	Ca0.165Al2.33Si3.67O10(OH)2 + 12 H2O = 0.165 Ca+2 + 2.33 Al(OH)4- + 3.67 H4SiO4 + 2 H+
+	-log_k -45.027
+	-delta_h 58.373 kcal
+	-Vm 156.16
+Talc
+	Mg3Si4O10(OH)2 + 4 H2O + 6 H+ = 3 Mg+2 + 4 H4SiO4
+	-log_k 21.399
+	-delta_h -46.352 kcal
+	-Vm 68.34
+Illite
+	K0.6Mg0.25Al2.3Si3.5O10(OH)2 + 11.2 H2O = 0.6 K+ + 0.25 Mg+2 + 2.3 Al(OH)4- + 3.5 H4SiO4 + 1.2 H+
+	-log_k -40.267
+	-delta_h 54.684 kcal
+	-Vm 141.48
+Chrysotile
+	Mg3Si2O5(OH)4 + 6 H+ = H2O + 2 H4SiO4 + 3 Mg+2
+	-log_k 32.2
+	-delta_h -46.8 kcal
+	-analytic 13.248 0 10217.1 -6.1894
+	-Vm 106.5808 # 277.11/2.60
+Sepiolite
+	Mg2Si3O7.5OH:3H2O + 4 H+ + 0.5 H2O = 2 Mg+2 + 3 H4SiO4
+	-log_k 15.76
+	-delta_h -10.7 kcal
+	-Vm 143.765
+Sepiolite(d)
+	Mg2Si3O7.5OH:3H2O + 4 H+ + 0.5 H2O = 2 Mg+2 + 3 H4SiO4
+	-log_k 18.66
+Hematite
+	Fe2O3 + 6 H+ = 2 Fe+3 + 3 H2O
+	-log_k -4.008
+	-delta_h -30.845 kcal
+	-Vm 30.39
+Goethite
+	FeOOH + 3 H+ = Fe+3 + 2 H2O
+	-log_k -1
+	-delta_h -14.48 kcal
+	-Vm 20.84
+Fe(OH)3(a)
+	Fe(OH)3 + 3 H+ = Fe+3 + 3 H2O
+	-log_k 4.891
+Pyrite
+	FeS2 + 2 H+ + 2 e- = Fe+2 + 2 HS-
+	-log_k -18.479
+	-delta_h 11.3 kcal
+	-Vm 23.48
+FeS(ppt)
+	FeS + H+ = Fe+2 + HS-
+	-log_k -3.915
+Mackinawite
+	FeS + H+ = Fe+2 + HS-
+	-log_k -4.648
+	-Vm 20.45
+Sulfur
+	S + 2 H+ + 2 e- = H2S
+	-log_k 4.882
+	-delta_h -9.5 kcal
+Vivianite
+	Fe3(PO4)2:8H2O = 3 Fe+2 + 2 PO4-3 + 8 H2O
+	-log_k -36
+Pyrolusite # H2O added for surface calc's
+	MnO2:H2O + 4 H+ + 2 e- = Mn+2 + 3 H2O
+	-log_k 41.38
+	-delta_h -65.11 kcal
+Hausmannite
+	Mn3O4 + 8 H+ + 2 e- = 3 Mn+2 + 4 H2O
+	-log_k 61.03
+	-delta_h -100.64 kcal
+Manganite
+	MnOOH + 3 H+ + e- = Mn+2 + 2 H2O
+	-log_k 25.34
+Pyrochroite
+	Mn(OH)2 + 2 H+ = Mn+2 + 2 H2O
+	-log_k 15.2
+Halite
+	NaCl = Cl- + Na+
+	log_k 1.57
+	-delta_h 1.37
+	#-analytic -713.4616   -.1201241   37302.21    262.4583    -2106915.
+	-Vm 27.1
+Sylvite
+	KCl = K+ + Cl-
+	log_k 0.9
+	-delta_h 8.5
+	# -analytic	3.984	0.0   -919.55
+	Vm 37.5
+# Gases...
+CO2(g)
+	CO2 = CO2
+	-log_k -1.468
+	-delta_h -4.776 kcal
+	-analytic 10.5624 -2.3547e-2 -3972.8 0 5.8746e5 1.9194e-5
+	-T_c 304.2 # critical T, K
+	-P_c 72.86 # critical P, atm
+	-Omega 0.225 # acentric factor
+H2O(g)
+	H2O = H2O
+	-log_k 1.506; delta_h -44.03 kJ
+	-T_c 647.3; -P_c 217.6; -Omega 0.344
+	-analytic -16.5066 -2.0013E-3 2710.7 3.7646 0 2.24E-6
+O2(g)
+	O2 = O2
+	-log_k -2.8983
+	-analytic -7.5001 7.8981e-3 0 0 2.0027e5
+	-T_c 154.6; -P_c 49.8; -Omega 0.021
+H2(g)
+	H2 = H2
+	-log_k -3.105
+	-delta_h -4.184 kJ
+	-analytic -9.3114 4.6473e-3 -49.335 1.4341 1.2815e5
+	-T_c 33.2; -P_c 12.8; -Omega -0.225
+N2(g)
+	N2 = N2
+	-log_k -3.1864
+	-analytic -58.453 1.818e-3 3199 17.909 -27460
+	-T_c 126.2; -P_c 33.5; -Omega 0.039
+H2S(g)
+	H2S = H+ + HS-
+	log_k -7.93
+	-delta_h 9.1
+	-analytic -45.07 -0.02418 0 17.9205 # H2S solubilities, 0 - 300�C, 1 - 987 atm, Jiang et al., 2020, CG 555, 119816
+	-T_c 373.2; -P_c 88.2; -Omega 0.1
+CH4(g)
+	CH4 = CH4
+	-log_k -2.8
+	-analytic 10.44 -7.65e-3 -6669 0 1.014e6 # CH4 solubilities 25 - 100�C
+	-T_c 190.6; -P_c 45.4; -Omega 0.008
+#Amm(g)
+#	Amm = Amm
+NH3(g)
+	NH3 = NH3
+	-log_k 1.7966
+	-analytic -18.758 3.367e-4 2.5113e3 4.8619 39.192
+	-T_c 405.6; -P_c 111.3; -Omega 0.25
+# redox-uncoupled gases
+Oxg(g)
+	Oxg = Oxg
+	-analytic -7.5001 7.8981e-3 0 0 2.0027e5
+	-T_c 154.6; -P_c 49.8; -Omega 0.021
+Hdg(g)
+	Hdg = Hdg
+	-analytic -9.3114 4.6473e-3 -49.335 1.4341 1.2815e5
+	-T_c 33.2; -P_c 12.8; -Omega -0.225
+Ntg(g)
+	Ntg = Ntg
+	-analytic -58.453 1.818e-3 3199 17.909 -27460
+	T_c 126.2; -P_c 33.5; -Omega 0.039
+Mtg(g)
+	Mtg = Mtg
+	-log_k -2.8
+	-analytic 10.44 -7.65e-3 -6669 0 1.014e6 # CH4 solubilities 25 - 100�C
+	-T_c 190.6; -P_c 45.4; -Omega 0.008
+H2Sg(g)
+	H2Sg = H+ + HSg-
+	log_k -7.93
+	-delta_h 9.1
+	-analytic -45.07 -0.02418 0 17.9205 # H2S solubilities, 0 - 300�C, 1 - 987 atm, Jiang et al., 2020, CG 555, 119816
+	-T_c 373.2; -P_c 88.2; -Omega 0.1
+Melanterite
+	FeSO4:7H2O = 7 H2O + Fe+2 + SO4-2
+	-log_k -2.209
+	-delta_h 4.91 kcal
+	-analytic 1.447 -0.004153 0 0 -214949
+Alunite
+	KAl3(SO4)2(OH)6 + 6 H+ = K+ + 3 Al+3 + 2 SO4-2 + 6 H2O
+	-log_k -1.4
+	-delta_h -50.25 kcal
+Jarosite-K
+	KFe3(SO4)2(OH)6 + 6 H+ = 3 Fe+3 + 6 H2O + K+ + 2 SO4-2
+	-log_k -9.21
+	-delta_h -31.28 kcal
+Zn(OH)2(e)
+	Zn(OH)2 + 2 H+ = Zn+2 + 2 H2O
+	-log_k 11.5
+Smithsonite
+	ZnCO3 = Zn+2 + CO3-2
+	-log_k -10
+	-delta_h -4.36 kcal
+Sphalerite
+	ZnS + H+ = Zn+2 + HS-
+	-log_k -11.618
+	-delta_h 8.25 kcal
+Willemite 289
+	Zn2SiO4 + 4 H+ = 2 Zn+2 + H4SiO4
+	-log_k 15.33
+	-delta_h -33.37 kcal
+Cd(OH)2
+	Cd(OH)2 + 2 H+ = Cd+2 + 2 H2O
+	-log_k 13.65
+Otavite 315
+	CdCO3 = Cd+2 + CO3-2
+	-log_k -12.1
+	-delta_h -0.019 kcal
+CdSiO3 328
+	CdSiO3 + H2O + 2 H+ = Cd+2 + H4SiO4
+	-log_k 9.06
+	-delta_h -16.63 kcal
+CdSO4 329
+	CdSO4 = Cd+2 + SO4-2
+	-log_k -0.1
+	-delta_h -14.74 kcal
+Cerussite 365
+	PbCO3 = Pb+2 + CO3-2
+	-log_k -13.13
+	-delta_h 4.86 kcal
+Anglesite 384
+	PbSO4 = Pb+2 + SO4-2
+	-log_k -7.79
+	-delta_h 2.15 kcal
+Pb(OH)2 389
+	Pb(OH)2 + 2 H+ = Pb+2 + 2 H2O
+	-log_k 8.15
+	-delta_h -13.99 kcal
+GAS_BINARY_PARAMETERS
+H2O(g)  CO2(g)           0.19
+H2O(g)  H2S(g)           0.19
+H2O(g)  H2Sg(g)          0.19
+H2O(g)  CH4(g)           0.49
+H2O(g)  Mtg(g)           0.49
+H2O(g)  Methane(g)       0.49
+H2O(g)  N2(g)            0.49
+H2O(g)  Ntg(g)           0.49
+H2O(g)  Ethane(g)        0.49
+H2O(g)  Propane(g)       0.55
+EXCHANGE_MASTER_SPECIES
+	X X-
+EXCHANGE_SPECIES
+	X- = X-
+	-log_k 0
+	Na+ + X- = NaX
+	-log_k 0
+	-gamma 4.08 0.082
+	K+ + X- = KX
+	-log_k 0.7
+	-gamma 3.5 0.015
+	-delta_h -4.3 # Jardine & Sparks, 1984
+	Li+ + X- = LiX
+	-log_k -0.08
+	-gamma 6 0
+	-delta_h 1.4 # Merriam & Thomas, 1956
+# !!!!!
+#	  H+ + X- = HX
+#	  -log_k  1.0
+#	  -gamma  9.0	0
+#	AmmH+ + X- = AmmHX
+	NH4+ + X- = NH4X
+	-log_k 0.6
+	-gamma 2.5 0
+	-delta_h -2.4 # Laudelout et al., 1968
+	Ca+2 + 2 X- = CaX2
+	-log_k 0.8
+	-gamma 5 0.165
+	-delta_h 7.2 # Van Bladel & Gheyl, 1980
+	Mg+2 + 2 X- = MgX2
+	-log_k 0.6
+	-gamma 5.5 0.2
+	-delta_h 7.4 # Laudelout et al., 1968
+	Sr+2 + 2 X- = SrX2
+	-log_k 0.91
+	-gamma 5.26 0.121
+	-delta_h 5.5 # Laudelout et al., 1968
+	Ba+2 + 2 X- = BaX2
+	-log_k 0.91
+	-gamma 4 0.153
+	-delta_h 4.5 # Laudelout et al., 1968
+	Mn+2 + 2 X- = MnX2
+	-log_k 0.52
+	-gamma 6 0
+	Fe+2 + 2 X- = FeX2
+	-log_k 0.44
+	-gamma 6 0
+	Cu+2 + 2 X- = CuX2
+	-log_k 0.6
+	-gamma 6 0
+	Zn+2 + 2 X- = ZnX2
+	-log_k 0.8
+	-gamma 5 0
+	Cd+2 + 2 X- = CdX2
+	-log_k 0.8
+	-gamma 0 0
+	Pb+2 + 2 X- = PbX2
+	-log_k 1.05
+	-gamma 0 0
+	Al+3 + 3 X- = AlX3
+	-log_k 0.41
+	-gamma 9 0
+	AlOH+2 + 2 X- = AlOHX2
+	-log_k 0.89
+	-gamma 0 0
+SURFACE_MASTER_SPECIES
+	Hfo_s Hfo_sOH
+	Hfo_w Hfo_wOH
+SURFACE_SPECIES
+# All surface data from
+# Dzombak and Morel, 1990
+#
+#
+# Acid-base data from table 5.7
+#
+# strong binding site--Hfo_s,
+	Hfo_sOH = Hfo_sOH
+	-log_k 0
+	Hfo_sOH + H+ = Hfo_sOH2+
+	-log_k 7.29 # = pKa1,int
+	Hfo_sOH = Hfo_sO- + H+
+	-log_k -8.93 # = -pKa2,int
+# weak binding site--Hfo_w
+	Hfo_wOH = Hfo_wOH
+	-log_k 0
+	Hfo_wOH + H+ = Hfo_wOH2+
+	-log_k 7.29 # = pKa1,int
+	Hfo_wOH = Hfo_wO- + H+
+	-log_k -8.93 # = -pKa2,int
+###############################################
+# CATIONS #
+###############################################
+#
+# Cations from table 10.1 or 10.5
+#
+# Calcium
+	Hfo_sOH + Ca+2 = Hfo_sOHCa+2
+	-log_k 4.97
+	Hfo_wOH + Ca+2 = Hfo_wOCa+ + H+
+	-log_k -5.85
+# Strontium
+	Hfo_sOH + Sr+2 = Hfo_sOHSr+2
+	-log_k 5.01
+	Hfo_wOH + Sr+2 = Hfo_wOSr+ + H+
+	-log_k -6.58
+	Hfo_wOH + Sr+2 + H2O = Hfo_wOSrOH + 2 H+
+	-log_k -17.6
+# Barium
+	Hfo_sOH + Ba+2 = Hfo_sOHBa+2
+	-log_k 5.46
+	Hfo_wOH + Ba+2 = Hfo_wOBa+ + H+
+	-log_k -7.2 # table 10.5
+#
+# Cations from table 10.2
+#
+# Cadmium
+	Hfo_sOH + Cd+2 = Hfo_sOCd+ + H+
+	-log_k 0.47
+	Hfo_wOH + Cd+2 = Hfo_wOCd+ + H+
+	-log_k -2.91
+# Zinc
+	Hfo_sOH + Zn+2 = Hfo_sOZn+ + H+
+	-log_k 0.99
+	Hfo_wOH + Zn+2 = Hfo_wOZn+ + H+
+	-log_k -1.99
+# Copper
+	Hfo_sOH + Cu+2 = Hfo_sOCu+ + H+
+	-log_k 2.89
+	Hfo_wOH + Cu+2 = Hfo_wOCu+ + H+
+	-log_k 0.6 # table 10.5
+# Lead
+	Hfo_sOH + Pb+2 = Hfo_sOPb+ + H+
+	-log_k 4.65
+	Hfo_wOH + Pb+2 = Hfo_wOPb+ + H+
+	-log_k 0.3 # table 10.5
+#
+# Derived constants table 10.5
+#
+# Magnesium
+	Hfo_wOH + Mg+2 = Hfo_wOMg+ + H+
+	-log_k -4.6
+# Manganese
+	Hfo_sOH + Mn+2 = Hfo_sOMn+ + H+
+	-log_k -0.4 # table 10.5
+	Hfo_wOH + Mn+2 = Hfo_wOMn+ + H+
+	-log_k -3.5 # table 10.5
+# Iron, strong site: Appelo, Van der Weiden, Tournassat & Charlet, EST 36, 3096
+	Hfo_sOH + Fe+2 = Hfo_sOFe+ + H+
+	-log_k -0.95
+# Iron, weak site: Liger et al., GCA 63, 2939, re-optimized for D&M
+	Hfo_wOH + Fe+2 = Hfo_wOFe+ + H+
+	-log_k -2.98
+	Hfo_wOH + Fe+2 + H2O = Hfo_wOFeOH + 2 H+
+	-log_k -11.55
+###############################################
+# ANIONS #
+###############################################
+#
+# Anions from table 10.6
+#
+# Phosphate
+	Hfo_wOH + PO4-3 + 3 H+ = Hfo_wH2PO4 + H2O
+	-log_k 31.29
+	Hfo_wOH + PO4-3 + 2 H+ = Hfo_wHPO4- + H2O
+	-log_k 25.39
+	Hfo_wOH + PO4-3 + H+ = Hfo_wPO4-2 + H2O
+	-log_k 17.72
+#
+# Anions from table 10.7
+#
+# Borate
+	Hfo_wOH + H3BO3 = Hfo_wH2BO3 + H2O
+	-log_k 0.62
+#
+# Anions from table 10.8
+#
+# Sulfate
+	Hfo_wOH + SO4-2 + H+ = Hfo_wSO4- + H2O
+	-log_k 7.78
+	Hfo_wOH + SO4-2 = Hfo_wOHSO4-2
+	-log_k 0.79
+#
+# Derived constants table 10.10
+#
+	Hfo_wOH + F- + H+ = Hfo_wF + H2O
+	-log_k 8.7
+	Hfo_wOH + F- = Hfo_wOHF-
+	-log_k 1.6
+#
+# Carbonate: Van Geen et al., 1994 reoptimized for D&M model
+#
+	Hfo_wOH + CO3-2 + H+ = Hfo_wCO3- + H2O
+	-log_k 12.56
+	Hfo_wOH + CO3-2 + 2 H+ = Hfo_wHCO3 + H2O
+	-log_k 20.62
+#
+# Silicate: Swedlund, P.J. and Webster, J.G., 1999. Water Research 33, 3413-3422.
+#
+	Hfo_wOH + H4SiO4 = Hfo_wH3SiO4 + H2O	; log_K 4.28
+	Hfo_wOH + H4SiO4 = Hfo_wH2SiO4- + H+ + H2O; log_K -3.22
+	Hfo_wOH + H4SiO4 = Hfo_wHSiO4-2 + 2 H+ + H2O; log_K -11.69
+MEAN_GAMMAS
+CaCl2    Ca+2  1  Cl-    2
+CaSO4    Ca+2  1  SO4-2  1
+CaCO3    Ca+2  1  CO3-2  1
+Ca(OH)2  Ca+2  1  OH-    2
+MgCl2    Mg+2  1  Cl-    2
+MgSO4    Mg+2  1  SO4-2  1
+MgCO3    Mg+2  1  CO3-2  1
+Mg(OH)2  Mg+2  1  OH-    2
+NaCl     Na+   1  Cl-    1
+Na2SO4   Na+   2  SO4-2  1
+NaHCO3   Na+   1  HCO3-  1
+Na2CO3   Na+   2  CO3-2  1
+NaOH     Na+   1  OH-    1
+KCl      K+    1  Cl-    1
+K2SO4    K+    2  SO4-2  1
+HCO3     K+    1  HCO3-  1
+K2CO3    K+    2  CO3-2  1
+KOH      K+    1  OH-    1
+HCl      H+    1  Cl-    1
+H2SO4    H+    2  SO4-2  1
+HBr      H+    1  Br-    1
+RATES
+###########
+#Quartz
+###########
+#
+#######
+# Example of quartz kinetic rates block:
+#	  KINETICS
+#	  Quartz
+#		  -m0  158.8	 # 90 % Qu
+#		  -parms 0.146  1.5
+#		  -step 3.1536e8 in 10
+#		  -tol 1e-12
+Quartz
+  -start
+1 REM Specific rate k from Rimstidt and Barnes, 1980, GCA 44,1683
+2 REM k = 10^-13.7 mol/m2/s (25 C), Ea = 90 kJ/mol
+3 REM sp. rate * parm(2) due to salts (Dove and Rimstidt, MSA Rev. 29, 259)
+4 REM PARM(1) = Specific area of Quartz, m^2/mol Quartz
+5 REM PARM(2) = salt correction: (1 + 1.5 * c_Na (mM)), < 35
+10 dif_temp = 1/TK - 1/298
+20 pk_w = 13.7 + 4700.4 * dif_temp
+40 moles = PARM(1) * M0 * PARM(2) * (M/M0)^0.67 * 10^-pk_w * (1 - SR("Quartz"))
+#			  Integrate...
+50 SAVE moles * TIME
+  -end
+###########
+#K-feldspar
+###########
+#
+# Sverdrup and Warfvinge, 1995, Estimating field weathering rates
+# using laboratory kinetics: Reviews in mineralogy and geochemistry,
+# vol. 31, p. 485-541.
+#
+# As described in:
+# Appelo and Postma, 2005, Geochemistry, groundwater
+# and pollution, 2nd Edition: A.A. Balkema Publishers,
+# p. 162-163 and 395-399.
+#
+# Assume soil is 10% K-feldspar by mass in 1 mm spheres (radius 0.05 mm)
+# Assume density of rock and Kspar is 2600 kg/m^3 = 2.6 kg/L
+# GFW Kspar 0.278 kg/mol
+#
+# Moles of Kspar per liter pore space calculation:
+#   Mass of rock per liter pore space = 0.7*2.6/0.3	  = 6.07	kg rock/L pore space
+#   Mass of Kspar per liter pore space 6.07x0.1	   = 0.607    kg Kspar/L pore space
+#   Moles of Kspar per liter pore space 0.607/0.278	  = 2.18	mol Kspar/L pore space
+#
+# Specific area calculation:
+#   Volume of sphere 4/3 x pi x r^3			  = 5.24e-13 m^3 Kspar/sphere
+#   Mass of sphere 2600 x 5.24e-13			= 1.36e-9  kg Kspar/sphere
+#   Moles of Kspar in sphere 1.36e-9/0.278		= 4.90e-9  mol Kspar/sphere
+#   Surface area of one sphere 4 x pi x r^2		  = 3.14e-8  m^2/sphere
+#   Specific area of K-feldspar in sphere 3.14e-8/4.90e-9 = 6.41 m^2/mol Kspar
+#
+#
+# Example of KINETICS data block for K-feldspar rate:
+#	  KINETICS 1
+#	  K-feldspar
+#		  -m0 2.18	    # 10% Kspar, 0.1 mm cubes
+#		  -m  2.18	    # Moles per L pore space
+#		  -parms 6.41  0.1    # m^2/mol Kspar, fraction adjusts lab rate to field rate
+#		  -time 1.5 year in 40
+K-feldspar
+ -start
+1 REM Sverdrup and Warfvinge, 1995, mol m^-2 s^-1
+2 REM PARM(1) = Specific area of Kspar m^2/mol Kspar
+3 REM PARM(2) = Adjusts lab rate to field rate
+4 REM temp corr: from A&P, p. 162: E (kJ/mol) / R / 2.303 = H in H*(1/T-1/281)
+5 REM K-Feldspar parameters
+10 DATA 11.7, 0.5, 4e-6, 0.4, 500e-6, 0.15, 14.5, 0.14, 0.15, 13.1, 0.3
+20 RESTORE 10
+30 READ pK_H, n_H, lim_Al, x_Al, lim_BC, x_BC, pK_H2O, z_Al, z_BC, pK_OH, o_OH
+40 DATA 3500, 2000, 2500, 2000
+50 RESTORE 40
+60 READ e_H, e_H2O, e_OH, e_CO2
+70 pk_CO2 = 13
+80 n_CO2 = 0.6
+100 REM Generic rate follows
+110 dif_temp = 1/TK - 1/281
+120 BC = ACT("Na+") + ACT("K+") + ACT("Mg+2") + ACT("Ca+2")
+130 REM rate by H+
+140 pk_H = pk_H + e_H * dif_temp
+150 rate_H = 10^-pk_H * ACT("H+")^n_H / ((1 + ACT("Al+3") / lim_Al)^x_Al * (1 + BC / lim_BC)^x_BC)
+160 REM rate by hydrolysis
+170 pk_H2O = pk_H2O + e_H2O * dif_temp
+180 rate_H2O = 10^-pk_H2O / ((1 + ACT("Al+3") / lim_Al)^z_Al * (1 + BC / lim_BC)^z_BC)
+190 REM rate by OH-
+200 pk_OH = pk_OH + e_OH * dif_temp
+210 rate_OH = 10^-pk_OH * ACT("OH-")^o_OH
+220 REM rate by CO2
+230 pk_CO2 = pk_CO2 + e_CO2 * dif_temp
+240 rate_CO2 = 10^-pk_CO2 * (SR("CO2(g)"))^n_CO2
+250 rate = rate_H + rate_H2O + rate_OH + rate_CO2
+260 area = PARM(1) * M0 *(M/M0)^0.67
+270 rate = PARM(2) * area * rate * (1-SR("K-feldspar"))
+280 moles = rate * TIME
+290 SAVE moles
+ -end
+###########
+#Albite
+###########
+#
+# Sverdrup and Warfvinge, 1995, Estimating field weathering rates
+# using laboratory kinetics: Reviews in mineralogy and geochemistry,
+# vol. 31, p. 485-541.
+#
+# As described in:
+# Appelo and Postma, 2005, Geochemistry, groundwater
+# and pollution, 2nd Edition: A.A. Balkema Publishers,
+# p. 162-163 and 395-399.
+#
+# Example of KINETICS data block for Albite rate:
+#	  KINETICS 1
+#	  Albite
+#		  -m0 0.46	    # 2% Albite, 0.1 mm cubes
+#		  -m  0.46	    # Moles per L pore space
+#		  -parms 6.04  0.1    # m^2/mol Albite, fraction adjusts lab rate to field rate
+#		  -time 1.5 year in 40
+#
+# Assume soil is 2% Albite by mass in 1 mm spheres (radius 0.05 mm)
+# Assume density of rock and Albite is 2600 kg/m^3 = 2.6 kg/L
+# GFW Albite 0.262 kg/mol
+#
+# Moles of Albite per liter pore space calculation:
+#   Mass of rock per liter pore space = 0.7*2.6/0.3	  = 6.07	kg rock/L pore space
+#   Mass of Albite per liter pore space 6.07x0.02	 = 0.121    kg Albite/L pore space
+#   Moles of Albite per liter pore space 0.607/0.262	 = 0.46	mol Albite/L pore space
+#
+# Specific area calculation:
+#   Volume of sphere 4/3 x pi x r^3			  = 5.24e-13 m^3 Albite/sphere
+#   Mass of sphere 2600 x 5.24e-13			= 1.36e-9  kg Albite/sphere
+#   Moles of Albite in sphere 1.36e-9/0.262		  = 5.20e-9  mol Albite/sphere
+#   Surface area of one sphere 4 x pi x r^2		  = 3.14e-8  m^2/sphere
+#   Specific area of Albite in sphere 3.14e-8/5.20e-9	= 6.04 m^2/mol Albite
+Albite
+ -start
+1 REM Sverdrup and Warfvinge, 1995, mol m^-2 s^-1
+2 REM PARM(1) = Specific area of Albite m^2/mol Albite
+3 REM PARM(2) = Adjusts lab rate to field rate
+4 REM temp corr: from A&P, p. 162 E (kJ/mol) / R / 2.303 = H in H*(1/T-1/281)
+5 REM Albite parameters
+10 DATA 11.5, 0.5, 4e-6, 0.4, 500e-6, 0.2, 13.7, 0.14, 0.15, 11.8, 0.3
+20 RESTORE 10
+30 READ pK_H, n_H, lim_Al, x_Al, lim_BC, x_BC, pK_H2O, z_Al, z_BC, pK_OH, o_OH
+40 DATA 3500, 2000, 2500, 2000
+50 RESTORE 40
+60 READ e_H, e_H2O, e_OH, e_CO2
+70 pk_CO2 = 13
+80 n_CO2 = 0.6
+100 REM Generic rate follows
+110 dif_temp = 1/TK - 1/281
+120 BC = ACT("Na+") + ACT("K+") + ACT("Mg+2") + ACT("Ca+2")
+130 REM rate by H+
+140 pk_H = pk_H + e_H * dif_temp
+150 rate_H = 10^-pk_H * ACT("H+")^n_H / ((1 + ACT("Al+3") / lim_Al)^x_Al * (1 + BC / lim_BC)^x_BC)
+160 REM rate by hydrolysis
+170 pk_H2O = pk_H2O + e_H2O * dif_temp
+180 rate_H2O = 10^-pk_H2O / ((1 + ACT("Al+3") / lim_Al)^z_Al * (1 + BC / lim_BC)^z_BC)
+190 REM rate by OH-
+200 pk_OH = pk_OH + e_OH * dif_temp
+210 rate_OH = 10^-pk_OH * ACT("OH-")^o_OH
+220 REM rate by CO2
+230 pk_CO2 = pk_CO2 + e_CO2 * dif_temp
+240 rate_CO2 = 10^-pk_CO2 * (SR("CO2(g)"))^n_CO2
+250 rate = rate_H + rate_H2O + rate_OH + rate_CO2
+260 area = PARM(1) * M0 *(M/M0)^0.67
+270 rate = PARM(2) * area * rate * (1-SR("Albite"))
+280 moles = rate * TIME
+290 SAVE moles
+ -end
+########
+#Calcite
+########
+# Example of KINETICS data block for calcite rate,
+#   in mmol/cm2/s, Plummer et al., 1978, AJS 278, 179; Appelo et al., AG 13, 257
+# KINETICS 1
+# Calcite
+# -tol   1e-8
+# -m0    3.e-3
+# -m	3.e-3
+# -parms 1.67e5   0.6  # cm^2/mol calcite, exp factor
+# -time  1 day
+Calcite
+   -start
+1 REM PARM(1) = specific surface area of calcite, cm^2/mol calcite
+2 REM PARM(2) = exponent for M/M0
+10 si_cc = SI("Calcite")
+20 IF (M <= 0 and si_cc < 0) THEN GOTO 200
+30 k1 = 10^(0.198 - 444 / TK )
+40 k2 = 10^(2.84 - 2177 /TK )
+50 IF TC <= 25 THEN k3 = 10^(-5.86 - 317 / TK)
+60 IF TC > 25 THEN k3 = 10^(-1.1 - 1737 / TK )
+80 IF M0 > 0 THEN area = PARM(1)*M0*(M/M0)^PARM(2) ELSE area = PARM(1)*M
+110 rate = area * (k1 * ACT("H+") + k2 * ACT("CO2") + k3 * ACT("H2O"))
+120 rate = rate * (1 - 10^(2/3*si_cc))
+130 moles = rate * 0.001 * TIME # convert from mmol to mol
+200 SAVE moles
+   -end
+#######
+#Pyrite
+#######
+#
+# Williamson, M.A. and Rimstidt, J.D., 1994,
+# Geochimica et Cosmochimica Acta, v. 58, p. 5443-5454,
+# rate equation is mol m^-2 s^-1.
+#
+# Example of KINETICS data block for pyrite rate:
+#	  KINETICS 1
+#	  Pyrite
+#		  -tol    1e-8
+#		  -m0	5.e-4
+#		  -m	 5.e-4
+#		  -parms  0.3	0.67	.5	 -0.11
+#		  -time 1 day in 10
+Pyrite
+  -start
+1 REM Williamson and Rimstidt, 1994
+2 REM PARM(1) = log10(specific area), log10(m^2 per mole pyrite)
+3 REM PARM(2) = exp for (M/M0)
+4 REM PARM(3) = exp for O2
+5 REM PARM(4) = exp for H+
+10 REM Dissolution in presence of DO
+20 if (M <= 0) THEN GOTO 200
+30 if (SI("Pyrite") >= 0) THEN GOTO 200
+40 log_rate = -8.19 + PARM(3)*LM("O2") + PARM(4)*LM("H+")
+50 log_area = PARM(1) + LOG10(M0) + PARM(2)*LOG10(M/M0)
+60 moles = 10^(log_area + log_rate) * TIME
+200 SAVE moles
+  -end
+##########
+#Organic_C
+##########
+#
+# Example of KINETICS data block for SOC (sediment organic carbon):
+#	  KINETICS 1
+#	  Organic_C
+#		  -formula C
+#		  -tol    1e-8
+#		  -m	 5e-3   # SOC in mol
+#		  -time 30 year in 15
+Organic_C
+ -start
+1 REM Additive Monod kinetics for SOC (sediment organic carbon)
+2 REM Electron acceptors: O2, NO3, and SO4
+10 if (M <= 0) THEN GOTO 200
+20 mO2 = MOL("O2")
+30 mNO3 = TOT("N(5)")
+40 mSO4 = TOT("S(6)")
+50 k_O2 = 1.57e-9 # 1/sec
+60 k_NO3 = 1.67e-11 # 1/sec
+70 k_SO4 = 1.e-13 # 1/sec
+80 rate = k_O2 * mO2/(2.94e-4 + mO2)
+90 rate = rate + k_NO3 * mNO3/(1.55e-4 + mNO3)
+100 rate = rate + k_SO4 * mSO4/(1.e-4 + mSO4)
+110 moles = rate * M * (M/M0) * TIME
+200 SAVE moles
+ -end
+###########
+#Pyrolusite
+###########
+#
+# Postma, D. and Appelo, C.A.J., 2000, GCA, vol. 64, pp. 1237-1247.
+# Rate equation given as mol L^-1 s^-1
+#
+# Example of KINETICS data block for Pyrolusite
+#	  KINETICS 1-12
+#	  Pyrolusite
+#		  -tol    1.e-7
+#		  -m0	0.1
+#		  -m	 0.1
+#		  -time 0.5 day in 10
+Pyrolusite
+  -start
+10 if (M <= 0) THEN GOTO 200
+20 sr_pl = SR("Pyrolusite")
+30 if (sr_pl > 1) THEN GOTO 100
+40 REM sr_pl <= 1, undersaturated
+50 Fe_t = TOT("Fe(2)")
+60 if Fe_t < 1e-8 then goto 200
+70 moles = 6.98e-5 * Fe_t * (M/M0)^0.67 * TIME * (1 - sr_pl)
+80 GOTO 200
+100 REM sr_pl > 1, supersaturated
+110 moles = 2e-3 * 6.98e-5 * (1 - sr_pl) * TIME
+200 SAVE moles * SOLN_VOL
+  -end
+#
+#  Additional definition of PHASES, RATE parameters, and RATES examples
+#
+#  RATE_PARAMETERS_PK has parameters from Palandri and Kharaka (2004). A compilation of rate parameters of water-mineral interaction kinetics for application to geochemical modeling. USGS Open-File Report 2004-1068.
+#
+#  RATE_PARAMETERS_SVD has two examples from Sverdrup, Oelkers, Lampa, Belyazid, Kurz, and Akselsson (2019). Reviews and Syntheses: weathering of silicate minerals in soils and watersheds: parameterization of the weathering kinetics module in the PROFILE and ForSAFE models. Biogeosciences Discuss. 1-58.
+#
+#  RATE_PARAMETERS_HERMANSKA has parameters from Hermanska, Voigt, Marieni, Declercq, and Oelkers (2022, 2023). A comprehensive and internally consistent mineral dissolution rate database: Part I: Primary silicate minerals and glasses. Chemical Geology, 597, p.120807, Part II: Secondary silicate minerals. Chemical Geology, p.121632.
+#
+#  Example RATES definitions and input files with KINETICS show the application in
+#  Albite_PK # Palandri and Kharaka
+#  Albite_Svd # Sverdrup
+#  Albite_Hermanska
+#  Calcite_PK # Palandri and Kharaka
+#  Calcite # Plummer, Wigley, Parkhurst 1978, AJS 278, 179-216.
+#  Quartz_PK # Palandri and Kharaka
+#  Quartz_Svd # Sverdrup
+#  Quartz_Hermanska #
+#  Quartz_Rimstidt_Barnes
+#  Montmorillonite # Na, K, Mg, Ca exchange, Hermanska rate for the TOT layer
+#
+PHASES # defined for formulas and affinities of kinetic (mostly) dissolving minerals
+#        Unless noted otherwise, data from ThermoddemV1.10_15Dec2020.dat,
+#           tidied with lsp.exe from https://phreeplot.org/lsp/lsp.html
+Actinolite # Hornblende, Ferroactinolite
+Ca2(Mg2.25Fe2.5Al0.25)(Si7.75Al0.25)O22(OH)2 + 15 H+ + 7 H2O = 0.5 Al+3 + 2 Ca+2 + 2.5 Fe+2 + 2.25 Mg+2 + 7.75 H4SiO4
+	log_k 7.128
+	delta_h -181.662 #kJ/mol #19bla/lac
+	-analytic -5.0954182E3 -6.949504E-1 3.0825312E5 1.8133351E3 -1.8767155E7
+Almandine # (alpha)
+Fe3Al2Si3O12 + 12 H+ = 2 Al+3 + 3 Fe+2 + 3 H4SiO4
+	log_k 42.18
+	delta_h -458.683 #kJ/mol #95rob/hem
+	-analytic -3.0848427E3 -4.4981168E-1 1.9672956E5 1.0990475E3 -1.0509115E7
+Analcime
+Na0.99Al0.99Si2.01O6:H2O + 3.96 H+ + 1.04 H2O = 0.99 Al+3 + 0.99 Na+ + 2.01 H4SiO4
+	log_k 6.654
+	delta_h -98 #kJ/mol #04neu/hov
+	-analytic -1.3403358E3 -1.8135021E-1 8.3684586E4 4.7527556E2 -4.9476886E6
+Andalusite
+Al2SiO5 + 6 H+ = 2 Al+3 + H4SiO4 + H2O
+	log_k 16.206
+	delta_h -244.61 #kJ/mol #Internal calculation
+	-analytic -1.339469E3 -2.048042E-1 8.5279067E4 4.7661954E2 -4.3249835E6
+Andesine # defined for elemental release
+Na0.6Ca0.4Si2.6Al1.4O8 + 8 H2O = 0.6 Na+ + 0.4 Ca+2 + 2.6 H4SiO4 + 1.4 Al(OH)4-
+Andradite
+Ca3Fe2Si3O12 + 12 H+ = 3 Ca+2 + 2 Fe+3 + 3 H4SiO4
+	log_k 33.787
+	delta_h -327.864 #kJ/mol #Internal calculation
+	-analytic -2.9077837E3 -4.2372897E-1 1.7981493E5 1.040602E3 -9.7870213E6
+Anglesite
+PbSO4 = Pb+2 + SO4-2
+	log_k -7.848
+	delta_h 11.55 #kJ/mol #89cox/wag
+	-analytic -1.6531905E3 -2.6395706E-1 9.1051907E4 5.9877724E2 -5.5987833E6
+Annite
+KFe3(AlSi3)O10(OH)2 + 10 H+ = Al+3 + 3 Fe+2 + K+ + 3 H4SiO4
+	log_k 32.771
+	delta_h -306.153 #kJ/mol #92cir/nav
+	-analytic -2.6382558E3 -3.7460641E-1 1.6621477E5 9.4111433E2 -9.2002058E6
+Anorthite
+Ca(Al2Si2)O8 + 8 H+ = 2 Al+3 + Ca+2 + 2 H4SiO4
+	log_k 24.235
+	delta_h -303.522 #kJ/mol #95rob/hem
+	-analytic -1.9788284E3 -2.9190197E-1 1.2612201E5 7.0425974E2 -6.7173266E6
+Anthophyllite
+Mg7Si8O22(OH)2 + 14 H+ + 8 H2O = 7 Mg+2 + 8 H4SiO4
+	log_k 73.783
+	delta_h -583.247 #kJ/mol #95rob/hem
+	-analytic -5.2321622E3 -7.0079895E-1 3.3845592E5 1.8579984E3 -1.9360477E7
+Antigorite
+Mg48Si34O85(OH)62 + 96 H+ = 48 Mg+2 + 34 H4SiO4 + 11 H2O
+	log_k 500.08
+	delta_h -3743.421 #kJ/mol #98hol/pow
+	-analytic -2.9383249E4 -4.0195982 1.8738549E6 1.0481455E4 -1.0123582E8
+# As2S3  # Orpiment # no As in phreeqc.dat
+# As2S3 + 6H2O = 2H2AsO3- + 3HS- + 5H+
+# log_k -65.110
+# delta_h 334.975 #kJ/mol        #Internal calculation
+# -analytic -2.5599772E+3 -4.2267991E-1 1.1988784E+5 9.3328822E+2 -8.0517057E+6
+Augite # Pyroxene(CaAl)
+CaAl(AlSi)O6 + 8 H+ = 2 Al+3 + Ca+2 + H4SiO4 + 2 H2O
+	log_k 36.234
+	delta_h -370.792 #kJ/mol #Internal calculation
+	-analytic -1.5908243E3 -2.4603865E-1 1.0453251E5 5.681931E2 -4.9909659E6
+Biotite # defined for elemental release
+KFe3(AlSi3)O10(OH)2 + 10 H+ = Al+3 + K+ + 3 Fe+2 + 3 H4SiO4
+Bronzite # defined for elemental release
+Mg0.8Fe0.2SiO3 + 2 H+ + H2O = 0.8 Mg+2 + 0.2 Fe+2 + H4SiO4
+Brucite
+Mg(OH)2 + 2 H+ = Mg+2 + 2 H2O
+	log_k 17.112
+	delta_h -114.518 #kJ/mol #08bla
+	-analytic -3.5641635E2 -5.3167189E-2 2.4317829E4 1.2873122E2 -9.5286882E5
+Bytownite # defined for elemental release
+Na0.2Ca0.8Si2.2Al1.8O8 + 8 H2O = 0.2 Na+ + 0.8 Ca+2 + 2.2 H4SiO4 + 1.8 Al(OH)4-
+Chabazite
+Ca(Al2Si4)O12:6H2O + 8 H+ = 2 Al+3 + Ca+2 + 4 H4SiO4 + 2 H2O
+	log_k 11.541
+	delta_h -200.464 #kJ/mol #08bla
+	-analytic -2.5875779E3 -3.5298441E-1 1.6180839E5 9.1700928E2 -9.5494778E6
+Chamosite(Daphnite)
+Fe5Al(AlSi3)O10(OH)8 + 16 H+ = 2 Al+3 + 5 Fe+2 + 3 H4SiO4 + 6 H2O
+	log_k 47.603
+	delta_h -497.518 #kJ/mol #01vid/par
+	-analytic -3.7422355E3 -5.4789298E-1 2.3185338E5 1.338448E3 -1.2120616E7
+Chrysotile
+Mg3Si2O5(OH)4 + 6 H+ = 3 Mg+2 + 2 H4SiO4 + H2O
+	log_k 33.182
+	delta_h -244.552 #kJ/mol #04eva
+	-analytic -1.8039877E3 -2.4743291E-1 1.1552931E5 6.4375706E2 -6.1763163E6
+Clinochlore
+Mg5Al(AlSi3)O10(OH)8 + 16 H+ = 2 Al+3 + 5 Mg+2 + 3 H4SiO4 + 6 H2O
+	log_k 61.706
+	delta_h -593.773 #kJ/mol #05vid/par
+	-analytic -3.933293E3 -5.6860144E-1 2.4698841E5 1.4055516E3 -1.2607E7
+Clinoptilolite(Ca)
+Ca0.55(Si4.9Al1.1)O12:3.9H2O + 4.4 H+ + 3.7 H2O = 1.1 Al+3 + 0.55 Ca+2 + 4.9 H4SiO4
+	log_k -2.085
+	delta_h -58.407 #kJ/mol #09bla
+	-analytic -2.3815518E3 -3.0085981E-1 1.4942318E5 8.390927E2 -9.6254008E6
+Clinoptilolite(K)
+K1.1(Si4.9Al1.1)O12:2.7H2O + 4.4 H+ + 4.9 H2O = 1.1 Al+3 + 1.1 K+ + 4.9 H4SiO4
+	log_k -1.142
+	delta_h -49.035 #kJ/mol #09bla
+	-analytic -2.3148616E3 -2.905299E-1 1.4612903E5 8.1530832E2 -9.5298429E6
+Clinoptilolite(Na)
+Na1.1(Si4.9Al1.1)O12:3.5H2O + 4.4 H+ + 4.1 H2O = 1.1 Al+3 + 1.1 Na+ + 4.9 H4SiO4
+	log_k -0.113
+	delta_h -50.769 #kJ/mol #09bla
+	-analytic -2.3846087E3 -2.9645291E-1 1.4988094E5 8.401942E2 -9.6738611E6
+Cordierite
+Mg2Al3(AlSi5)O18 + 16 H+ + 2 H2O = 4 Al+3 + 2 Mg+2 + 5 H4SiO4
+	log_k 49.433
+	delta_h -648.745 #kJ/mol #95rob/hem
+	-analytic -4.3696636E3 -6.2958321E-1 2.8022776E5 1.5507866E3 -1.5147654E7
+Cristobalite # (alpha)
+SiO2 + 2 H2O = H4SiO4
+	log_k -3.158
+	delta_h 18.829 #kJ/mol #04fab/sax
+	-analytic -3.544017E2 -4.1702635E-2 2.2114271E4 1.2427357E2 -1.6001472E6
+# Cristobalite(beta)
+# SiO2 + 2H2O = 1H4SiO4
+# log_k -3.096
+# #delta_h   0     #kJ/mol
+# -analytic -3.6088361E+2 -4.1957223E-2 2.2873339E+4 1.2628239E+2 -1.6799304E+6
+Dawsonite
+NaAlCO3(OH)2 + 3 H+ = Al+3 + HCO3- + Na+ + 2 H2O
+	log_k 4.327
+	delta_h -76.33 #kJ/mol #76fer/stu
+	-analytic -1.21599E3 -1.9110794E-1 6.8919359E4 4.3970018E2 -3.7220307E6
+Diaspore
+AlO(OH) + 3 H+ = Al+3 + 2 H2O
+	log_k 6.866
+	delta_h -108.76 #kJ/mol #95rob/hem
+	-analytic -4.8201662E2 -7.7930965E-2 2.9964822E4 1.7237439E2 -1.3257386E6
+Diopside
+CaMg(SiO3)2 + 4 H+ + 2 H2O = Ca+2 + Mg+2 + 2 H4SiO4
+	log_k 21.743
+	delta_h -153.574 #kJ/mol #Internal calculation
+	-analytic -1.332806E3 -1.8198553E-1 8.603858E4 4.749095E2 -4.8802351E6
+Dolomite(disordered)
+CaMg(CO3)2 + 2 H+ = 2 HCO3- + Ca+2 + Mg+2
+	log_k 4.299
+	delta_h -73.162 #kJ/mol #78hel/del,92ajoh
+	-analytic -1.7814432E3 -2.8852695E-1 9.9263747E4 6.4714027E2 -5.5533944E6
+Edenite # (alpha)
+Na(Ca2Mg5)(AlSi7)O22(OH)2 + 18 H+ + 4 H2O = Al+3 + 2 Ca+2 + 5 Mg+2 + Na+ + 7 H4SiO4
+	log_k 81.946
+	delta_h -679.296 #kJ/mol #97got
+	-analytic -5.4623009E3 -7.5241996E-1 3.5051336E5 1.9444511E3 -1.942E7
+Enstatite # (alpha)
+MgSiO3 + 2 H+ + H2O = Mg+2 + H4SiO4
+	log_k 11.844
+	delta_h -93.265 #kJ/mol #78hel/del
+	-analytic -7.0139177E2 -9.4618096E-2 4.5846726E4 2.4912172E2 -2.5565294E6
+Epidote
+Ca2FeAl2Si3O12(OH) + 13 H+ = 2 Al+3 + 2 Ca+2 + Fe+3 + 3 H4SiO4 + H2O
+	log_k 32.23
+	delta_h -411.613 #kJ/mol #04got
+	-analytic -3.1567388E3 -4.6487997E-1 1.9676775E5 1.1260692E3 -1.0558252E7
+Fayalite
+Fe2SiO4 + 4 H+ = 2 Fe+2 + H4SiO4
+	log_k 19.03
+	delta_h -157.157 #kJ/mol #Internal calculation
+	-analytic -1.0258478E3 -1.4618015E-1 6.6129821E4 3.6618221E2 -3.5053712E6
+Ferroactinolite # = Ferrotremolite
+(Ca2Fe5)Si8O22(OH)2 + 14 H+ + 8 H2O = 2 Ca+2 + 5 Fe+2 + 8 H4SiO4
+	log_k 53.699
+	delta_h -412.225 #kJ/mol #Internal calculation
+	-analytic -4.942592E3 -6.6976495E-1 3.1400258E5 1.7585882E3 -1.8552107E7
+Fluorapatite # (Natur)
+Ca5(PO4)3F + 6 H+ = 5 Ca+2 + F- + 3 H2PO4-
+	log_k -0.91
+	delta_h -115.601 #kJ/mol #Internal calculation
+	-analytic -3.7675938E3 -6.2227437E-1 2.0719593E5 1.369906E3 -1.1775417E7
+Forsterite
+Mg2SiO4 + 4 H+ = 2 Mg+2 + H4SiO4
+	log_k 28.609
+	delta_h -217.115 #kJ/mol #Internal calculation
+	-analytic -1.0983766E3 -1.5385695E-1 7.321503E4 3.91599E2 -3.7061609E6
+Glauconite
+(K0.75Mg0.25Fe1.5Al0.25)(Al0.25Si3.75)O10(OH)2 + 7 H+ + 3 H2O = 0.5 Al+3 + 1.25 Fe+3 + 0.75 K+ + 0.25 Mg+2 + 3.75 H4SiO4 + 0.25 Fe+2
+	log_k 1.873
+	delta_h -120.903 #kJ/mol #15bla/vie
+	-analytic -2.3976207E3 -3.2091227E-1 1.4807364E5 8.4865741E2 -9.0151175E6
+Glaucophane
+Na2(Mg3Al2)Si8O22(OH)2 + 14 H+ + 8 H2O = 2 Al+3 + 3 Mg+2 + 2 Na+ + 8 H4SiO4
+	log_k 37.026
+	delta_h -378.727 #kJ/mol #95rob/hem
+	-analytic -5.095188E3 -6.8518568E-1 3.2040873E5 1.8087612E3 -1.9006796E7
+Grossular
+Ca3Al2Si3O12 + 12 H+ = 2 Al+3 + 3 Ca+2 + 3 H4SiO4
+	log_k 49.372
+	delta_h -442.383 #kJ/mol #95rob/hem
+	-analytic -2.9566754E3 -4.3410622E-1 1.8868769E5 1.057027E3 -1.0038715E7
+# Hornblende # see Actinolite, Edenite, Pargasite, Ferroactinolite
+Heulandite(Ca)
+Ca1.07Al2.14Si6.86O18:6.17H2O + 8.56 H+ + 3.27 H2O = 2.14 Al+3 + 1.07 Ca+2 + 6.86 H4SiO4
+	log_k 2.457
+	delta_h -139.108 #kJ/mol #09bla
+	-analytic -3.7607701E3 -5.0483789E-1 2.3083824E5 1.3337643E3 -1.4294418E7
+# Ilmenite # Ti not in phreeqc.dat
+# FeTiO3 + 2H+ + 1H2O = 1Fe+2 + 1Ti(OH)4
+# log_k 1.816
+# delta_h -87.445 #kJ/mol        #Internal calculation
+# -analytic -7.7719505E+2 -8.1479565E-2 4.34898E+4 2.7302259E+2 -1.612373E+6
+Heulandite(Na)
+Na2.14Al2.14Si6.86O18:6.17H2O + 8.56 H+ + 3.27 H2O = 2.14 Al+3 + 2.14 Na+ + 6.86 H4SiO4
+	log_k 2.797
+	delta_h -126.775 #kJ/mol #09bla
+	-analytic -3.7890714E3 -4.9720069E-1 2.3269508E5 1.3423841E3 -1.4400431E7
+Jadeite
+NaAl(SiO3)2 + 4 H+ + 2 H2O = Al+3 + Na+ + 2 H4SiO4
+	log_k 7.561
+	delta_h -95.502 #kJ/mol #95rob/hem
+	-analytic -1.3237509E3 -1.8118316E-1 8.2628986E4 4.7016122E2 -4.9060741E6
+Kyanite
+Al2SiO5 + 6 H+ = 2 Al+3 + H4SiO4 + H2O
+	log_k 15.936
+	delta_h -240.322 #kJ/mol #Internal calculation
+	-analytic -1.3447799E3 -2.0581745E-1 8.5324148E4 4.7877192E2 -4.3369481E6
+Labradorite # defined for elemental release
+Na0.4Ca0.6Si2.4Al1.6O8 + 8 H2O = 0.4 Na+ + 0.6 Ca+2 + 2.4 H4SiO4 + 1.6 Al(OH)4-
+Larnite(alpha)
+Ca2SiO4 + 4 H+ = 2 Ca+2 + H4SiO4
+	log_k 39.044
+	delta_h -238.161 #kJ/mol #95rob/hem
+	-analytic -8.9908942E2 -1.301379E-1 6.3335055E4 3.2296168E2 -3.0793446E6
+Larnite(beta)
+Ca2SiO4 + 4 H+ = 2 Ca+2 + H4SiO4
+	log_k 39.322
+#delta_h 0 #kJ/mol
+	-analytic -9.0365527E2 -1.3027777E-1 6.4015139E4 3.243254E2 -3.1477489E6
+Larnite(gamma)
+Ca2SiO4 + 4 H+ = 2 Ca+2 + H4SiO4
+	log_k 41.444
+#delta_h 0 #kJ/mol
+	-analytic -8.7896206E2 -1.2907359E-1 6.3430487E4 3.1585123E2 -3.1477489E6
+Laumontite
+Ca(Al2Si4)O12:4H2O + 8 H+ = 2 Al+3 + Ca+2 + 4 H4SiO4
+	log_k 11.695
+	delta_h -204.244 #kJ/mol #96kis/nav
+	-analytic -2.6447429E3 -3.6684244E-1 1.6419074E5 9.3900001E2 -9.6343473E6
+Leonhardtite
+MgSO4:4H2O = Mg+2 + SO4-2 + 4 H2O
+	log_k -0.886
+	delta_h -24.03 #kJ/mol #74nau/ryz
+	-analytic -1.8009396E3 -2.6450971E-1 9.9216758E4 6.5010323E2 -5.5554353E6
+Leucite # minteq.dat
+KAlSi2O6 + 2 H2O + 4 H+ = 2 H4SiO4 + Al+3 + K+
+	log_k 6.423
+	delta_h -22.085 kcal
+Lizardite
+Mg3Si2O5(OH)4 + 6 H+ = 3 Mg+2 + 2 H4SiO4 + H2O
+	log_k 33.093
+	delta_h -242.552 #kJ/mol #04eva
+	-analytic -1.8045338E3 -2.475614E-1 1.1546724E5 6.4405193E2 -6.1786442E6
+Magnetite
+Fe3O4 + 8 H+ = 2 Fe+3 + Fe+2 + 4 H2O
+	log_k 10.362
+	delta_h -215.594 #kJ/mol #90hem
+	-analytic -1.3520774E3 -2.1498134E-1 8.0017747E4 4.8502632E2 -3.7344997E6
+Microcline
+K(AlSi3)O8 + 4 H+ + 4 H2O = Al+3 + K+ + 3 H4SiO4
+	log_k 0.015
+	delta_h -49.203 #kJ/mol #95rob/hem
+	-analytic -1.6018728E3 -2.1339241E-1 9.9207574E4 5.6723025E2 -6.2943433E6
+Montmorillonite(HcCa)
+Ca0.3Mg0.6Al1.4Si4O10(OH)2 + 6 H+ + 4 H2O = 1.4 Al+3 + 0.3 Ca+2 + 0.6 Mg+2 + 4 H4SiO4
+	log_k 6.903
+	delta_h -154.564 #kJ/mol #15bla/vie
+	-analytic -2.3616529E3 -3.1379357E-1 1.4899818E5 8.3431323E2 -9.0744862E6
+Montmorillonite(HcK)
+K0.6Mg0.6Al1.4Si4O10(OH)2 + 6 H+ + 4 H2O = 1.4 Al+3 + 0.6 K+ + 0.6 Mg+2 + 4 H4SiO4
+	log_k 4.449
+	delta_h -119.628 #kJ/mol #15bla/vie
+	-analytic -2.3324885E3 -3.0832834E-1 1.4605682E5 8.2462838E2 -9.022722E6
+Montmorillonite(HcMg)
+Mg0.3Mg0.6Al1.4Si4O10(OH)2 + 6 H+ + 4 H2O = 1.4 Al+3 + 0.9 Mg+2 + 4 H4SiO4
+	log_k 5.996
+	delta_h -156.964 #kJ/mol #15bla/vie
+	-analytic -2.3909331E3 -3.1726069E-1 1.5070041E5 8.4429278E2 -9.163021E6
+Montmorillonite(HcNa)
+Na0.6Mg0.6Al1.4Si4O10(OH)2 + 6 H+ + 4 H2O = 1.4 Al+3 + 0.6 Mg+2 + 0.6 Na+ + 4 H4SiO4
+	log_k 5.472
+	delta_h -135.658 #kJ/mol #15bla/vie
+	-analytic -2.3671642E3 -3.1193536E-1 1.486659E5 8.3634354E2 -9.1085654E6
+Montmorillonite(MgCa)
+Ca0.17Mg0.34Al1.66Si4O10(OH)2 + 6 H+ + 4 H2O = 1.66 Al+3 + 0.17 Ca+2 + 0.34 Mg+2 + 4 H4SiO4
+	log_k 4.222
+	delta_h -146.668 #kJ/mol #15bla/vie
+	-analytic -2.3648299E3 -3.1580182E-1 1.4861699E5 8.3532612E2 -9.0862785E6
+Montmorillonite(MgK)
+K0.34Mg0.34Al1.66Si4O10(OH)2 + 6 H+ + 4 H2O = 1.66 Al+3 + 0.34 K+ + 0.34 Mg+2 + 4 H4SiO4
+	log_k 2.83
+	delta_h -126.865 #kJ/mol #15bla/vie
+	-analytic -2.3483045E3 -3.1270489E-1 1.4694997E5 8.2983827E2 -9.056946E6
+Montmorillonite(MgMg)
+Mg0.17Mg0.34Al1.66Si4O10(OH)2 + 6 H+ + 4 H2O = 1.66 Al+3 + 0.51 Mg+2 + 4 H4SiO4
+	log_k 3.708
+	delta_h -148.028 #kJ/mol #15bla/vie
+	-analytic -2.3814282E3 -3.1776702E-1 1.4958186E5 8.4098328E2 -9.1364559E6
+Montmorillonite(MgNa)
+Na0.34Mg0.34Al1.66Si4O10(OH)2 + 6 H+ + 4 H2O = 1.66 Al+3 + 0.34 Mg+2 + 0.34 Na+ + 4 H4SiO4
+	log_k 3.411
+	delta_h -135.953 #kJ/mol #15bla/vie
+	-analytic -2.3679565E3 -3.1474933E-1 1.4842879E5 8.3647775E2 -9.1055977E6
+MordeniteB # (Ca)
+Ca0.515Al1.03Si4.97O12:3.1H2O + 4.12 H+ + 4.78 H2O = 1.03 Al+3 + 0.515 Ca+2 + 4.97 H4SiO4
+	log_k -2.898
+	delta_h -56.278 #kJ/mol #09bla
+	-analytic -2.3577543E3 -2.9682032E-1 1.4847577E5 8.2993876E2 -9.6241393E6
+MordeniteJ
+Ca0.289Na0.362Al0.94Si5.06O12:3.468H2O + 3.76 H+ + 4.772 H2O = 0.94 Al+3 + 0.289 Ca+2 + 0.362 Na+ + 5.06 H4SiO4
+	log_k -4.16
+	delta_h -29.442 #kJ/mol #92joh/tas
+	-analytic -2.3112502E3 -2.9430315E-1 1.4403365E5 8.1541676E2 -9.418252E6
+Muscovite # (ordered)
+KAl2(AlSi3)O10(OH)2 + 10 H+ = 3 Al+3 + K+ + 3 H4SiO4
+	log_k 11.353
+	delta_h -253.923 #kJ/mol #06bla/pia
+	-analytic -2.5862792E3 -3.7607072E-1 1.5907562E5 9.2024545E2 -8.9668534E6
+Natrolite
+Na2(Al2Si3)O10:2H2O + 8 H+ = 2 Al+3 + 2 Na+ + 3 H4SiO4
+	log_k 19.326
+	delta_h -215.463 #kJ/mol #83joh/flo
+	-analytic -2.303612E3 -3.1993458E-1 1.4352482E5 8.1980235E2 -8.1431211E6
+Nepheline
+Na(AlSi)O4 + 4 H+ = Al+3 + Na+ + H4SiO4
+	log_k 14.077
+	delta_h -144.506 #kJ/mol #Internal calculation
+	-analytic -9.7409139E2 -1.3955693E-1 6.2423687E4 3.467383E2 -3.3400695E6
+Oligoclase # defined for elemental release
+Na0.8Ca0.2Si2.8Al1.2O8 + 8 H2O = 0.8 Na+ + 0.2 Ca+2 + 2.8 H4SiO4 + 1.2 Al(OH)4-
+Palygorskite # defined for elemental release
+Mg2Al2Si8O20(OH)2:8H2O + 10 H+ + 2 H2O = 2 Mg+2 + 2 Al+3 + 8 H4SiO4
+Paragonite
+NaAl2(AlSi3)O10(OH)2 + 10 H+ = 3 Al+3 + Na+ + 3 H4SiO4
+	log_k 16.804
+	delta_h -294.623 #kJ/mol #96rou/hov
+	-analytic -2.6452559E3 -3.8247258E-1 1.64246E5 9.4070011E2 -9.1107641E6
+Pargasite # Hornblende
+Na(Ca2Mg4Al)(Al2Si6)O22(OH)2 + 22 H+ = 3 Al+3 + 2 Ca+2 + 4 Mg+2 + Na+ + 6 H4SiO4
+	log_k 104.557
+	delta_h -940.614 #kJ/mol #Internal calculation
+	-analytic -5.7962939E3 -8.2700886E-1 3.7555969E5 2.0652064E3 -1.9772394E7
+Phlogopite
+KMg3(AlSi3)O10(OH)2 + 10 H+ = Al+3 + K+ + 3 Mg+2 + 3 H4SiO4
+	log_k 41.098
+	delta_h -353.123 #kJ/mol #92cir/nav
+	-analytic -2.7194067E3 -3.8106546E-1 1.7318081E5 9.69566E2 -9.4102646E6
+Prehnite
+Ca2Al2Si3O10(OH)2 + 10 H+ = 2 Al+3 + 2 Ca+2 + 3 H4SiO4
+	log_k 32.596
+	delta_h -339.617 #kJ/mol #98cha/kru
+	-analytic -2.6255465E3 -3.8041883E-1 1.6586587E5 9.3642007E2 -9.0549681E6
+Pyrophyllite
+Al2Si4O10(OH)2 + 6 H+ + 4 H2O = 2 Al+3 + 4 H4SiO4
+	log_k -0.418
+	delta_h -128.924 #kJ/mol #95rob/hem
+	-analytic -2.3595061E3 -3.237303E-1 1.4585394E5 8.3524091E2 -8.9193526E6
+Pyrrhotite(Hx) # Pyrrhotite
+FeS + H+ = Fe+2 + HS-
+	log_k -3.679
+	delta_h -10.009 #kJ/mol #05wal/pel
+	-analytic -1.1321823E3 -1.8235764E-1 6.1304821E4 4.1103628E2 -3.5403537E6
+Pyrrhotite(Mc) # Pyrrhotite
+FeS + H+ = Fe+2 + HS-
+	log_k -3.679
+	delta_h -10.009 #kJ/mol #05wal/pel
+	-analytic -1.1321823E3 -1.8235764E-1 6.1304821E4 4.1103628E2 -3.5403537E6
+Rhyolite # a mixture of minerals, defined for elemental release...
+Na0.078K0.046Al0.26Si1.23O2.912 + 3.048 H2O = 0.136 H+ + 0.078 Na+ + 0.046 K+ + 0.26 Al(OH)4- + 1.23 H4SiO4
+Riebeckite
+Na2(Fe3Fe2)Si8O22(OH)2 + 14 H+ + 8 H2O = 3 Fe+2 + 2 Na+ + 8 H4SiO4 + 2 Fe+3
+	log_k 9.199
+	delta_h -197.377 #kJ/mol #98hol/pow
+	-analytic -5.0079102E3 -6.7170777E-1 3.0608951E5 1.7785742E3 -1.8686839E7
+Saponite(Ca)
+Ca0.17Mg3Al0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + 0.17 Ca+2 + 3 Mg+2 + 3.66 H4SiO4
+	log_k 29.355
+	delta_h -262.766 #kJ/mol #15bla/vie
+	-analytic -2.5667428E3 -3.4039957E-1 1.6475488E5 9.099285E2 -9.472597E6
+Saponite(FeCa)
+Ca0.17Mg2FeAl0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + 0.17 Ca+2 + Fe+2 + 2 Mg+2 + 3.66 H4SiO4
+	log_k 26.569
+	delta_h -250.636 #kJ/mol #15bla/vie
+	-analytic -2.5356344E3 -3.373844E-1 1.6236385E5 8.9871835E2 -9.386812E6
+Saponite(FeK)
+K0.34Mg2FeAl0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + Fe+2 + 0.34 K+ + 2 Mg+2 + 3.66 H4SiO4
+	log_k 25.398
+	delta_h -232.093 #kJ/mol #15bla/vie
+	-analytic -2.515955E3 -3.3384661E-1 1.6058454E5 8.9209651E2 -9.3470003E6
+Saponite(FeMg)
+Mg0.17Mg2FeAl0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + Fe+2 + 2.17 Mg+2 + 3.66 H4SiO4
+	log_k 26.022
+	delta_h -251.806 #kJ/mol #15bla/vie
+	-analytic -2.5507675E3 -3.3914471E-1 1.6323608E5 9.0384868E2 -9.4321235E6
+Saponite(FeNa)
+Na0.34Mg2FeAl0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + Fe+2 + 2 Mg+2 + 0.34 Na+ + 3.66 H4SiO4
+	log_k 25.896
+	delta_h -240.711 #kJ/mol #15bla/vie
+	-analytic -2.5368817E3 -3.3606965E-1 1.6211086E5 8.9919435E2 -9.3999007E6
+Saponite(K)
+K0.33Mg3Al0.33Si3.67O10(OH)2 + 7.32 H+ + 2.68 H2O = 0.33 Al+3 + 0.33 K+ + 3 Mg+2 + 3.67 H4SiO4
+	log_k 27.43
+	delta_h -239.483 #kJ/mol #15bla/vie
+	-analytic -2.544416E3 -3.3629993E-1 1.6263915E5 9.0231366E2 -9.4312976E6
+Saponite(Mg)
+Mg0.17Mg3Al0.34Si3.66O10(OH)2 + 7.36 H+ + 2.64 H2O = 0.34 Al+3 + 3.17 Mg+2 + 3.66 H4SiO4
+	log_k 28.81
+	delta_h -263.946 #kJ/mol #15bla/vie
+	-analytic -2.5818719E3 -3.4215988E-1 1.6562747E5 9.1505763E2 -9.5179085E6
+Saponite(Na)
+Na0.33Mg3Al0.33Si3.67O10(OH)2 + 7.32 H+ + 2.68 H2O = 0.33 Al+3 + 3 Mg+2 + 0.33 Na+ + 3.67 H4SiO4
+	log_k 27.971
+	delta_h -248.219 #kJ/mol #15bla/vie
+	-analytic -2.5647603E3 -3.3846001E-1 1.6414122E5 9.0921188E2 -9.482682E6
+Saponite(SapCa)
+(Na0.394K0.021Ca0.038)(Si3.569Al0.397)(Mg2.949Fe0.055)O10(OH)2 + 7.724 H+ + 2.276 H2O = 0.397 Al+3 + 0.038 Ca+2 + 0.034 Fe+3 + 0.021 K+ + 2.949 Mg+2 + 0.394 Na+ + 3.569 H4SiO4 + 0.021 Fe+2
+	log_k 31.473
+	delta_h -277.172 #kJ/mol #13gai/bla
+	-analytic -2.5790231E3 -3.508959E-1 1.6429225E5 9.168404E2 -9.2969386E6
+Scolecite
+CaAl2Si3O10:3H2O + 8 H+ = 2 Al+3 + Ca+2 + 3 H4SiO4 + H2O
+	log_k 16.647
+	delta_h -233.213 #kJ/mol #83joh/flo
+	-analytic -2.3692738E3 -3.4026162E-1 1.4623007E5 8.4431312E2 -8.2035956E6
+Smectite # (MX80)
+Na0.409K0.024Ca0.009(Si3.738Al0.262)(Al1.598Mg0.214Fe0.208)O10(OH)2 + 7.048 H+ + 2.952 H2O = 1.86 Al+3 + 0.009 Ca+2 + 0.173 Fe+3 + 0.024 K+ + 0.214 Mg+2 + 0.409 Na+ + 3.738 H4SiO4 + 0.035 Fe+2
+	log_k 5.278
+	delta_h -175.308 #kJ/mol #12gai/bla
+	-analytic -2.4267042E3 -3.3712249E-1 1.5038583E5 8.6021197E2 -8.9284687E6
+Smectite(MX80:3.989H2O)
+Na0.409K0.024Ca0.009(Si3.738Al0.262)(Al1.598Mg0.214Fe0.208)O10(OH)2:3.989H2O + 7.048 H+ = 1.86 Al+3 + 0.009 Ca+2 + 0.173 Fe+3 + 0.024 K+ + 0.214 Mg+2 + 0.409 Na+ + 3.738 H4SiO4 + 0.035 Fe+2 + 1.037 H2O
+	log_k 1.774
+	delta_h -148.524 #kJ/mol #12gai/bla
+	-analytic -2.3838609E3 -3.2232449E-1 1.4844358E5 8.4261556E2 -8.9910004E6
+Smectite(MX80:5.189H2O)
+Na0.409K0.024Ca0.009(Si3.738Al0.262)(Al1.598Mg0.214Fe0.208)O10(OH)2:5.189H2O + 7.048 H+ = 1.86 Al+3 + 0.009 Ca+2 + 0.173 Fe+3 + 0.024 K+ + 0.214 Mg+2 + 0.409 Na+ + 3.738 H4SiO4 + 0.035 Fe+2 + 2.237 H2O
+	log_k 1.435
+	delta_h -140.43 #kJ/mol #12gai/bla
+	-analytic -2.3706061E3 -3.2008903E-1 1.4737914E5 8.3812012E2 -8.9524821E6
+Spodumene # from core10.dat
+LiAlSi2O6 + 4 H+ + 2 H2O = Al+3 + Li+ + 2 H4SiO4
+	log_k 6.9972
+	delta_h -89.1817
+	-analytic -9.8111 2.1191e-3 9.692e3 -3.0484 -7.8822e5
+	-Vm 58.37
+Staurolite
+Fe2Al9Si4O23(OH) + 31 H+ = 9 Al+3 + 2 Fe+2 + 4 H4SiO4 + 8 H2O
+	log_k 216.34
+	delta_h -1956.484 #kJ/mol #87woo/gar
+	-analytic -6.5297334E3 -1.0061427 4.5225123E5 2.3281295E3 -2.0588442E7
+Stilbite
+NaCa2(Al5Si13)O36:16H2O + 20 H+ = 5 Al+3 + 2 Ca+2 + Na+ + 13 H4SiO4
+	log_k 23.044
+	delta_h -403.823 #kJ/mol #01fri/neu
+	-analytic -7.4700792E3 -1.0099722 4.6170528E5 2.6510812E3 -2.7934606E7
+Thomsonite # defined for elemental release
+Na0.5CaAl2.5Si2.5O10:3H2O + 10 H+ = 2.5 Al+3 + 0.5 Na+ + Ca+2 + 2.5 H4SiO4 + 3 H2O
+Tourmaline # defined for elemental release
+NaFe1.5Mg1.5Al6B3Si6O27(OH)4 + 26 H2O + H+ = Na+ + 1.5 Fe+2 + 1.5 Mg+2 + 6 Al(OH)4- + 3 H3BO3 + 6 H4SiO4
+Tremolite
+(Ca2Mg5)Si8O22(OH)2 + 14 H+ + 8 H2O = 2 Ca+2 + 5 Mg+2 + 8 H4SiO4
+	log_k 67.281
+	delta_h -502.247 #kJ/mol #95rob/hem
+	-analytic -5.0977019E3 -6.8545317E-1 3.2680746E5 1.8129659E3 -1.8919407E7
+# Uraninite
+#        UO2 + 4 H+ = U+4 + 2 H2O
+#        log_k          -3.490
+#        delta_h        -18.630 kcal
+Wollastonite
+CaSiO3 + 2 H+ + H2O = Ca+2 + H4SiO4
+	log_k 14.047
+	delta_h -85.986 #kJ/mol #78hel/del,92ajoh
+	-analytic -6.3184784E2 -8.6944016E-2 4.1722732E4 2.2563038E2 -2.3494013E6
+Zoisite
+Ca2Al3Si3O12(OH) + 13 H+ = 3 Al+3 + 2 Ca+2 + 3 H4SiO4 + H2O
+	log_k 43.848
+	delta_h -485.113 #kJ/mol #01sme/fra
+	-analytic -3.1722373E3 -4.6912132E-1 2.0150433E5 1.1315082E3 -1.0643978E7
+RATE_PARAMETERS_PK
+#               Acid                    Neutral         Base
+#               log K   E       n(H+)   log K   E       log K   E       n(OH-)
+#               ================================================================
+Quartz          -30     0       0       -13.4   90.9    -30     0       0       # Table 4
+#
+SiO2(a)         -30     0       0       -12.31  76      -30     0       0       # Table 6
+Cristobalite    -30     0       0       -12.31  65      -30     0       0
+#
+Albite          -10.16  65      0.317   -12.56  65      -15.6   66.5    -0.471  # Table 1
+Oligoclase      -9.67   65      0.457   -11.84  69.8    -30     0       0       # Table 13
+Andesine        -8.88   53.5    0.541   -11.47  57.4    -30     0       0
+Labradorite     -7.87   42.1    0.626   -10.91  45.2    -30     0       0
+Bytownite       -5.85   29.3    1.018   -9.82   31.5    -30     0       0
+Anorthite       -3.5    16.6    1.411   -9.12   17.8    -30     0       0
+#
+K-feldspar      -10.06  51.7    0.5     -12.41  38      -21.2   94.1    -0.823  # Table 15
+#
+Nepheline       -2.73   62.9    1.13    -8.56   65.4    -10.76  37.8    -0.2    # Table 18
+Leucite         -6      132.2   0.7     -9.2    75.5    -10.66  56.6    -0.2
+#
+Forsterite      -6.85   67.2    0.47    -10.64  79      -30     0       0       # Table 23
+Fayalite        -4.8    94.4    0       -12.8   94.4    -30     0       0
+Almandine       -5.2    94.4    1       -10.7   103.8   -13.71  37.8    -0.35
+Grossular       -5.1    85      1       -10.7   103.8   -30     0       0
+Andradite       -5.2    94.4    1       -10.7   103.8   -30     0       0
+Kyanite         -10.17  -53.9   1.268   -17.44  53.9    -30     0       0
+Staurolite      -6.9    18.9    1       -12.2   56.6    -14.9   47.2    -0.3
+Epidote         -10.6   71.1    0.338   -11.99  70.7    -17.33  79.1    -0.556
+Zoisite         -7.5    66.1    0.5     -11.2   66.1    -30     0       0
+#
+Cordierite      -3.8    113.3   1       -11.2   28.3    -30     0       0       # Table 25
+Tourmaline      -6.5    75.5    1       -11.2   85      -30     0       0
+#
+augite          -6.82   78      0.7     -11.97  78      -30     0       0       # Table 26
+bronzite        -8.3    47.2    0.65    -11.7   66.1    -30     0       0
+diopside        -6.36   96.1    0.71    -11.11  40.6    -30     0       0
+enstatite       -9.02   80      0.6     -12.72  80      -30     0       0
+jadeite         -6      132.2   0.7     -9.5    94.4    -30     0       0
+spodumene       -4.6    94.4    0.7     -9.3    66.1    -30     0       0
+wollastonite    -5.37   54.7    0.4     -8.88   54.7    -30     0       0
+#
+anthophyllite   -11.94  51      0.44    -14.24  51      -30     0       0       # Table 27
+glaucophane     -5.6    85      0.7     -10.1   94.4    -30     0       0
+hornblende      -7      75.5    0.6     -10.3   94.4    -30     0       0
+riebeckite      -7.7    56.6    0.7     -12.2   47.2    -30     0       0
+tremolite       -8.4    18.9    0.7     -10.6   94.4    -30     0       0
+#
+biotite         -9.84   22      0.525   -12.55  22      -30     0       0       # Table 28
+glauconite      -4.8    85      0.7     -9.1    85      -30     0       0
+muscovite       -11.85  22      0.37    -13.55  22      -14.55  22      -0.22
+muscovite       -30     0       0       -13     22      -30     0       0
+paragonite      -30     0       0       -13     22      -30     0       0
+phlogopite      -30     0       0       -12.4   29      -30     0       0
+pyrophyllite    -30     0       0       -12.4   29      -30     0       0
+#
+kaolinite       -11.31  65.9    0.777   -13.18  22.2    -17.05  17.9    -0.472  # Table 29
+montmorillonite -12.71  48      0.22    -14.41  48      -14.41  48      -0.13   # Montmorillonite, K0.318(Si3.975Al0.025)(Al1.509Fe0.205Mg0.283)(OH)2.
+smectite        -10.98  23.6    0.34    -12.78  35      -16.52  58.9    -0.4    # Smectite, K0.04Ca0.5(Al2.8Fe0.53Mg0.7)(Si7.65Al0.35)O20(OH)4.
+#
+lizardite       -5.7    75.5    0.8     -12.4   56.6    -30     0       0       # Table 30
+chrysotile      -30     0       0       -12     73.5    -13.58  73.5    -0.23
+chlorite(14A)   -11.11  88      0.5     -12.52  88      -30     0       0
+talc            -30     0       0       -12     42      -30     0       0
+prehnite        -10.66  80.5    0.256   -13.16  93.4    -14.86  93.4    -0.2
+#
+goethite        -30     0       0       -7.94   86.5    -30     0       0       # Table 31
+hematite        -9.39   66.2    1       -14.6   66.2    -30     0       0
+magnetite       -8.59   18.6    0.279   -10.78  18.6    -30     0       0
+ilmenite        -8.35   37.9    0.421   -11.16  37.9    -30     0       0
+uraninite       -30     0       0       -7.98   32      -30     0       0
+#
+brucite         -4.73   59      0.5     -8.24   42      -30     0       0       # Table 32
+gibbsite        -7.65   47.5    0.992   -11.5   61.2    -16.65  80.1    -0.784
+diaspore        -30     0       0       -13.33  47.5    -23.6   47.5    -1.503
+#
+anglesite       -5.58   31.3    0.298   -6.5    31.3    -30     0       0       # Table 34
+anhydrite       -30     0       0       -3.19   14.3    -30     0       0
+gypsum          -30     0       0       -2.79   0       -30     0       0
+barite          -6.9    30.8    0.22    -7.9    30.8    -30     0       0
+celestite       -5.66   23.8    0.109   -30     0       -30     0       0
+#
+hydroxyapatite  -4.29   250     0.171   -6      250     -30     0       0       # Table 36
+fluorapatite    -3.73   250     0.613   -8      250     -30     0       0
+#
+halite          -30     0       0       -0.21   7.4     -30     0       0       # Table 37
+fluorite        -7.14   73      1       -13.79  73      -30     0       0
+#
+#               Acid                    Neutral         P_CO2
+#               log K   E       n(H+)   log K   E       log K   E       n(P_CO2)        Table
+#               ================================================================================
+calcite         -0.3    14.4    1       -5.81   23.5    -3.48   35.4    1               33      # specify Table number for P_CO2^n(P_CO2)
+dawsonite       -30     0       0       -7      62.8    -30     0       0               33
+dolomite(d)     -3.19   36.1    0.5     -7.53   52.2    -5.11   34.8    0.5             33
+dolomite        -3.76   56.7    0.5     -8.6    95.3    -5.37   45.7    0.5             33
+magnesite       -6.38   14.4    1       -9.34   23.5    -5.22   62.8    1               33
+#
+#               Acid and Fe+3                           Neutral and O2                  Base
+#               log K   E       n(H+)   n(Fe+3) log K   E       n(O2)   log K   E       n(OH-)  Table
+#               =========================================================================================
+pyrite          -7.52   56.9    -0.5    0.5     -4.55   56.9    0.5     -30     0       0       35      # specify Table number for Fe+3 and O2
+pyrrhotite(Mc)  -8.04   50.8    -0.597  0.355   -30     0       0       -30     0       0       35
+pyrrhotite(Hx)  -6.79   63      -0.09   0.356   -30     0       0       -30     0       0       35
+As2S3(a)        -30     0       0       0       -9.83   8.7     0.18    -17.39  8.7     -1.208  35
+RATE_PARAMETERS_SVD
+#               Table 4: E's                                    Table 3: H+-reaction                            H2O-reaction                                            CO2-reaction    Organic_acids           OH--reaction                                                            Table 5
+#               H+      H2O     CO2     Organic acids   OH-     pkH     nH      yAl     CAl     xBC     CBC     pkH2O   yAl     CAl     xBC     CBC     zSi     CSi     pkCO2   nCO2    pkOrg   nOrg    COrg    pkOH-   wOH-    yAl     CAl     xBC     CBC     zSi     CSi     #       Num     Mineral Formula
+#               =================================================================================================================================================================================================================================================================================================
+Albite          3350    2500    1680    1200            3100    14.6    0.5     0.4     0.4     0.4     0.5     16.8    0.15    4       0.15    200     3       900     16.05   0.6     14.7    0.5     5       15.4    0.3     0.1     12      0.5     5       3       900     #       1.6     Albite  NaAlSi3O8
+Quartz          3890    0       2200    2000            3320    18.4    0.3     0.3     5       0       500     17.8    0       5       0       5000    4       900     18      0.5     16.3    0.5     5       14.1    0.3     0.4     200     0       5000    1       900     #       8.3     Quartz  SiO2
+RATE_PARAMETERS_HERMANSKA
+#               Acid mechanism                          Neutral mechanism       Basic mechanism
+#               logk25  Aa              Eaa     n(H+)   logk25  Ab              Eab     logk25  Ac              Eac     n(OH)   #       Formula
+#               ================================================================================================================================
+# Amphiboles
+Anthophyllite   -12.4   5.70E-04        52      0.4     -13.7   5.00E-06        48      0       0               0       0
+Ferroactinolite -11.3   3.00E-03        50      0.2     -13.1   2.00E-05        48      0       0               0       0
+Riebeckite      -11.3   3.00E-03        50      0.2     -13.1   2.00E-05        48      0       0               0       0
+Tremolite       -11.3   3.00E-03        50      0.2     -13.1   2.00E-05        48      0       0               0       0
+Glaucophane     -6.1    2.20E+02        50      0.7     0       0               0       -12.8   1.00E-04        48      -0.1    #       Na0.14K0.09Ca2Fe1.78Mg2Al2Si7O22(OH)2
+Hornblende      -10.7   5.00E-03        50      0.2     0       0               0       -13.4   2.10E-05        48      -0.1    #       Ca2Mg4Al0.75Fe0.25(Si7AlO22)(OH)2
+# Feldspars
+Albite          -10.32  0.7             58      0.3     -11.19  0.21            60      -13.58  1.50E-05        50      -0.3
+Andesine        -7.99   147             58      0.7     -11.23  0.19            60      -13.58  1.50E-05        50      -0.3
+Anorthite       -5.17   9.80E+04        58      1.2     -11.34  0.15            60      -13.58  1.50E-05        50      -0.3
+Bytownite       -5.88   1.90E+04        58      1.1     -11.28  0.17            60      -13.58  1.50E-05        50      -0.3
+K-feldspar      -10.36  5.00E-02        51.7    0.5     -12.48  1.10E-02        60      -20.78  1.20E-10        62      -0.8    #       or Microcline
+Labradorite     -6.39   5.90E+03        58      1       -11.28  0.17            60      -13.58  1.50E-05        50      -0.3
+Oligoclase      -9.33   6.8             58      0.4     -11.21  0.2             60      -13.58  1.50E-05        50      -0.3
+# Glass
+Rhyolite        -9.1    1.60E-03        36      0.5     0       0               0       -16.27  7.00E-08        52      -0.6
+# Mica                                                                                                                          #       Also valid for
+Annite          -9.42   5.90E-07        18.2    0.5     -12.2   5.00E-09        22      -13.9   4.00E-10        25.5    -0.2    #       Biotite, Phlogopite
+Muscovite       -11.1   1.26E-04        41.3    0.4     -12.1   6.31E-06        39      -14.5   3.16E-05        57      -0.2
+# Olivines
+Fayalite        -6.26   1.20E+06        70.4    0.4     0       0               0       -7.39   1.91E+03        60.9    0.2
+Forsterite      -7.16   1.48E+05        70.4    0.4     0       0               0       -8.33   2.20E+02        60.9    0.2
+Larnite         -3.61   5.25E+08        70.4    0.4     0       0               0       -4.75   8.25E+05        60.9    0.2
+# Pyroxenes
+Augite          -8.2    1.52E+06        81.8    0.7     -12.8   350             83      0       0               0       0
+Bronzite        -9.8    9.50E-04        38.5    0.6     -11.7   7.60E-01        66.1    0       0               0       0
+Diopside        -9.8    8.55E-05        32.7    0.3     -11.01  4.30E-05        43.9    0       0               0       0
+Enstatite       -8.3    0.574           46.1    0.5     -11.9   6.30E+03        89.5    0       0               0       0
+# SiO2  polymorphs
+Quartz          -11.4   4.03E-04        45.6    0.3     0       0               0       -15     0.105           80      -0.4    #       Cristobalite
+SiO2(a)         -10.6   4.56E-04        41.6    0.3     0       0               0       -14.2   3.53E-02        73      -0.4
+# 2023, Table 1                                                                                                                         Also valid for
+Almandine       -5.21   2.00E+05        60      1       -11.2   2.31E-04        43.2    -14.6   6.00E-08        42.3    -0.4    #       Grossular
+Analcime        -3.3    5.00E+07        63      1       -11.3   1.00E-01        58.5    -14.3   7.50E-05        58      -0.4    #       Nepheline
+Andalusite      -10.57  3.90E-01        60      0.15    -12.61  8.00E-03        43.2    -22.82  8.80E-15        42.3    -1.2
+Andradite       -5.1    2.60E+05        60      1       -11.1   3.20E-04        43.2    0       0               0       0
+Antigorite      -10.3   2.80E-06        27      0.25    -12.4   2.00E-08        27      0       0               0       0       #       Chrysotile, Lizardite
+Chabazite       -6.56   2.21E-01        33.7    0.82    -11.55  1.56E-04        44.2    -12.05  4.94E-05        44.2    -0.2    #       Laumontite, Leonhardite
+Clinochlore     -9.08   1.50E-04        30      0.74    -13     4.70E-11        15      -14.3   2.00E-12        15      -0.2    #       Chamosite, Daphnite
+Clinoptilolite  -7.51   2.48E-02        33.7    0.82    -12.6   1.39E-05        44.2    -13.2   3.50E-06        44.2    -0.2    #       Heulandite, Mordenite, Stilbite
+Epidote         -10.47  1.09            60      0.3     -11.9   5.13E-05        43.2    -16.3   1.40E-09        42.3    -0.4    #       Zoisite
+Glauconite      -11.68  9.55E-07        32.3    0.37    -13.53  1.10E-07        37.5    0       0               0       0
+Illite          -11.9   7.30E-04        50      0.55    -14.68  3.84E-03        70      -20.19  6.00E-08        74      -0.6
+Jadeite         -6.68   25              46.1    0.5     -10.26  2.70E+05        89.5    0       0               0       0
+Kaolinite       -12.3   2.85            73      0.45    -14.1   4.15E-03        67      -21.3   2.40E-11        61      -0.76
+Kyanite         -11.1   1.15E-01        60      0.15    -13.5   1.00E-03        43.2    -21.6   1.50E-13        50      -1
+Mesolite        -5.61   1.97E+00        33.7    0.82    -10.7   1.11E-03        44.2    -11     5.54E-04        44.2    -0.2    #       Natrolite, Scolecite, Thomsonite
+Montmorillonite -11.7   1.66E-03        50.8    0.55    -14.3   9.00E-10        30      -17.2   1.50E-09        48      -0.3    #       Saponite, Smectite
+Paragonite      -11.9   7.30E-04        50      0.55    -14.68  3.84E-03        70      -20.19  6.00E-08        74      -0.6
+Prehnite        -10.4   1.30E+03        77      0.35    -14     1               80      -12.8   15              80      -0.075
+Pyrophyllite    -8.6    1.60E+04        73      0.7     -12.6   1.50E-01        67      -18.4   2.00E-08        61      -0.7
+Sepiolite       -11     5.89E-03        50.2    0.25    -13.2   8.00E-07        40.7    0       0               0       0       #       Palygorskite
+Spodumene       -5.38   4.90E+02        46.1    0.5     -8.95   5.40E+06        89.5    0       0               0       0
+Talc            -11.1   4.42E-03        50.2    0.36    -12.9   1.56E-06        40.7    0       0               0       0
+Wollastonite    -6.97   700             56      0.4     0       0               0       -7.81   200             52      0.15
+# # Example input files with RATES for KINETICS calculations
+# #
+# # compare Albite kinetics using rates from the compilations
+ # # for the PARMS, see https://www.hydrochemistry.eu/exmpls/kin_silicates.html
+# # =========================================================
+#
+# RATES
+# Albite_PK # Palandri and Kharaka, 2004
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Albite") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_PK("Albite")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# KINETICS 1
+# Albite_PK
+# -formula NaAlSi3O8; -parms 0 1 1 0.67
+# -m0 1; -time 1 # default
+# END
+# SOLUTION 1
+# PHASES
+  # Fix_pH; H+ = H+
+  # LiBr; LiBr = Li+ + Br-; -log_k -20 # (very) unsoluble phase with base cation and acid anion, permits to use HBr or LiOH as reactant
+# SELECTED_OUTPUT 1
+        # -file kinetic_rates_pH.inc
+        # -reset false
+# USER_PUNCH 1    # write out the pH's to equilibrate...
+        # 10 FOR i = 0 to 14 STEP 0.5
+        # 20      punch EOL$ + 'USE solution 1'
+        # 30      punch EOL$ + 'EQUILIBRIUM_PHASES 1'
+        # 40      punch EOL$ + ' LiBr'
+        # 50      punch EOL$ + ' Fix_pH ' + TRIM(STR$(-i)) + ' LiOH 10' # ...or HBr as reactant
+        # 60      punch EOL$ + 'USE kinetics 1'
+        # 70      punch EOL$ + 'END'
+        # 80 NEXT i
+# END
+# PRINT; -reset false
+# SELECTED_OUTPUT 1; -active false
+# USER_GRAPH 1; -headings pH Palandri
+# -axis_titles  pH "log10(initial rate / (mol / m2 / s))"
+# -axis_scale x_axis 0 14
+# 10 graph_x -la("H+")
+# 20 graph_sy log10(tot("Al"))
+# INCLUDE$ kinetic_rates_pH.inc
+# END
+# RATES
+# Albite_Svd # Sverdrup, 2019
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Albite") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_SVD("Albite")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# KINETICS 1
+# Albite_Svd
+# -formula NaAlSi3O8; -parms 0 1 20 0.67 # roughness = 20
+# USER_GRAPH 1; -headings pH Sverdup*20
+# INCLUDE$ kinetic_rates_pH.inc
+# END
+# KINETICS 1
+# Albite # from Sverdrup and Warfvinge, 1995
+# -formula NaAlSi3O8; -parms 1 20 # roughness = 20
+# USER_GRAPH 1; -headings pH Sverdup`95*20
+# INCLUDE$ kinetic_rates_pH.inc
+# END
+# RATES
+# Albite_Hermanska # 2022
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Albite") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_HERMANSKA("Albite")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# KINETICS 1
+# Albite_Hermanska
+# -formula NaAlSi3O8; -parms 0 1 1 0.67
+# USER_GRAPH 1; -headings pH Hermanska
+# INCLUDE$ kinetic_rates_pH.inc
+# END
+# USE solution 1
+# REACTION_TEMPERATURE 1; 25 25 in 21
+# USER_GRAPH 1; -headings Albite_data
+# 10 data 1.1, 2.05, 2.45, 2.9, 3, 3.5, 4.1, 5.1, 5.35, 5.47, 5.63, 5.63, 5.73, 7.73, 9.95, 9.95, 9.95, 10.6, 11.2, 11.55, 12.3
+# 20 data -10.25, -10.55, -10.82, -11.25, -11.1, -11.4, -11.47, -11.82, -11.75, -11.65, -11.83, -11.92, -11.92, -11.83, -10.97, -11.05, -11.13, -10.95, -10.55, -10.6, -10.38 # Chou, L., Wollast, R., 1985. Steady-state kinetics and dissolution mechanisms of albite. Am. J. Sci. 285, 963�993.
+# 30 restore 10 : dim ph(21) : for i = 1 to step_no : read ph(i) : next i
+# 40 restore 20 : dim lk(21) : for i = 1 to step_no : read lk(i) : next i
+# 50 i = step_no : plot_xy ph(i), lk(i), line_width = 0, color = Black, y_axis = 2, symbol_size = 10, symbol = Circle
+# END
+# # compare rates for calcite dissolution
+# #   of Palandri and Kharaka, 2004 and Plummer, Wigley and Parkhurst, 1978
+# #   at different initial CO2 concentrations.
+# # =====================================
+# USER_GRAPH 1; -active false
+# RATES
+# Calcite_PK # Palandri and Kharaka, 2004
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("calcite") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_PK("calcite")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# SOLUTION 1
+# pH 7 charge; C(4) 1 CO2(g) -2.5
+# KINETICS 1
+# calcite_PK
+# -formula CaCO3; -parms 0 1e-2 1 0.67
+# -time 0.1 10*1 hour
+# INCREMENTAL_REACTIONS true
+# USER_GRAPH 2; -headings h Palandri_SI(CO2_g).=.-2.5
+# -axis_titles "time / hours" "Calcite dissolved / (mmol/kgw)"
+# 10 graph_x total_time / 3600 : graph_sy tot("Ca") * 1e3
+# END
+# USE solution 1
+# KINETICS 1
+# Calcite
+# -parms  1e2  0.67  # cm^2/mol calcite, exp factor
+# -time  0.1 10*1 hour
+# USER_GRAPH 2; -headings h Plummer.Wigley.Parkhurst
+# END
+# SOLUTION 1
+# pH 7 charge; C(4) 1 CO2(g) -1.5
+# KINETICS 1
+# calcite_PK
+# -formula CaCO3
+# -parms  0  1e-2  1  0.67
+# -time 0.1 10*1 hour
+# USER_GRAPH 2; -headings h Palandri_SI(CO2_g).=.-1.5
+# END
+# USE solution 1
+# KINETICS 1
+# Calcite
+# -parms  1e2  0.67
+# -time  0.1 10*1 hour
+# USER_GRAPH 2; -headings h Plummer.Wigley.Parkhurst
+# END
+# # compare rates for quartz dissolution
+# #   and the effect of NaCl
+# # =====================================
+# USER_GRAPH 2; -active false
+# RATES
+# Quartz_PK # Palandri and Kharaka, 2004
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Quartz") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_PK("Quartz")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# SOLUTION 1
+# pH 7 charge
+# KINETICS 1
+# Quartz_PK
+# -formula SiO2
+# -parms  0  6  1  0.67
+# -time 0.1 10*1 year
+# INCREMENTAL_REACTIONS true
+# USER_GRAPH 3; -headings h Palandri
+# -axis_titles "time / years" "Quartz dissolved / (mmol/kgw)"
+# 10 graph_x total_time / 3.15e7 : graph_sy tot("Si") * 1e3
+# END
+# RATES
+# Quartz_Hermanska #
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Quartz") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_HERMANSKA("Quartz")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# USE solution 1
+# KINETICS 1
+# Quartz_Hermanska
+# -formula SiO2
+# -parms  0  6  1  0.67
+# -time 0.1 10*1 year
+# USER_GRAPH 3
+# -headings H Hermanska
+# END
+# RATES
+# Quartz_Svd # Sverdrup, 2019
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Quartz") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = RATE_SVD("Quartz")
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# USE solution 1
+# KINETICS 1
+# Quartz_Svd
+# -formula SiO2
+# -parms  0  6  1  0.67
+# -time 0.1 10*1 year
+# USER_GRAPH 3
+# -headings H Sverdup
+# END
+# RATES
+# Quartz_Rimstidt
+# #1 rem Specific rate k = 10^-13.34 mol/m2/s (25 C), Ea = 74 kJ/mol, Rimstidt, 2015, GCA 167, 195
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Quartz") : if affinity < parm(1) then SAVE 0 : END
+# 20 rate = 10^-(13.3 + 4700 * (1 / 298 - 1 / TK)) * (1 + 1500*tot("Na")) # salt correction, Dove and Rimstidt, 1994, MSA Rev. 29, 259
+# 20 rate = 10^-(13.3 + 4700 * (1 / 298 - 1 / TK)) + 11.2e3 * act("Na+")^0.33 * act("OH-")^0.44 * exp(-71.6/(8.314e-3 * TK)) # salt correction,  Rimstidt, 2015, GCA 167, 195
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# USE solution 1
+# KINETICS 1
+# Quartz_Rimstidt
+# -formula SiO2
+# -parms  0  6  1  0.67
+# -time 0.1 10*1 year
+# USER_GRAPH 3
+# -headings H Rimstidt
+# END
+# SOLUTION 1
+# pH 7 charge; Na 1; Cl 1
+# KINETICS 1
+# Quartz_Rimstidt
+# -formula SiO2
+# -parms  0  6  1  0.67
+# -time 0.1 10*1 year
+# USER_GRAPH 3
+# -headings H Rimstidt_1.mM.NaCl
+# END
+# # Example input file for calculating kinetic dissolution of Montmorillonite,
+# # a solid solution with exchangeable cations reacting fast;
+# # their ratios are related to the changing solution composition,
+# # and their amounts are connected to the kinetic reacting TOT layer.
+# #
+# # The affinity is related to a solid solution member, given by the fraction of the
+# # exchangeable cation (here Na+ or Ca+2). For the Gapon exchange formula,
+# # the exchange species and their log_k`s are from the solid solution members in ThermoddemV1
+# # For the Gaines Thomas formula, the Mg+2 and Ca+2 species are redefined.
+# # It also shows how the default X exchanger can be invkoed.
+# # # ==============================================================
+# USER_GRAPH 3; -active false
+# RATES
+# Montmorillonite
+# 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# # Gapon and Gaines-Tomas exchange formulas
+# 7  f_Na = (mol("Na0.34X_montm_mg") / tot("X_montm_mg"))
+# 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Montmorillonite(MgNa)") / f_Na
+# 20 rate = RATE_HERMANSKA("Montmorillonite") / f_Na
+# # #  Gapon, with Ca as exchange species...
+# # 7  f_Ca = (mol("Ca0.17X_montm_mg") / tot("X_montm_mg"))
+# # #  use SR("Montmorillonite(Mgca)")
+# # 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Montmorillonite(MgCa)") / f_Ca
+# # 20 rate = RATE_HERMANSKA("Montmorillonite") / f_Ca
+# # #  Gaines-Thomas exchange formula, with Ca as exchange species, uncomment the Gaines-Thomas EXCHANGE_SPECIES
+# # 7  f_Ca = (mol("Ca0.34X_montm_mg2") / 2 / tot("X_montm_mg")) : ex = 0.5
+# # 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Montmorillonite(MgCa)") / f_Ca^ex
+# # 20 rate = RATE_HERMANSKA("Montmorillonite") / f_Ca^ex
+# 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# 40 SAVE area * rate * affinity * TIME
+# -end
+# EXCHANGE_MASTER_SPECIES
+# X_montm_mg   X_montm_mg-0.34
+# EXCHANGE_SPECIES
+# # The Gapon formulation is easiest, with constants from Montmorillonite(Mg..) in PHASES
+              # X_montm_mg-0.34 =  X_montm_mg-0.34
+# 0.34 Na+  +   X_montm_mg-0.34 = Na0.34X_montm_mg; log_k -3.411 # 0 #
+# 0.34 K+   +   X_montm_mg-0.34 =  K0.34X_montm_mg; log_k -2.83 # 0.581 #
+# 0.17 Mg+2 +   X_montm_mg-0.34 = Mg0.17X_montm_mg; log_k -3.708 # -0.297 #
+# 0.17 Ca+2 +   X_montm_mg-0.34 = Ca0.17X_montm_mg; log_k -4.222 # -0.811 #
+# # # The divalent cations have rather low log_k, cf. A&P, p.254, log_k Ca0.5X ~ log_k KX
+# # # uncomment the following lines to see the effect...
+# # 0.17 Mg+2 +   X_montm_mg-0.34 = Mg0.17X_montm_mg; log_k -2.86
+# # 0.17 Ca+2 +   X_montm_mg-0.34 = Ca0.17X_montm_mg; log_k -2.83
+# # # also adapt the log_k`s of the solids...
+# # PHASES
+# # Montmorillonite(MgMg)
+# # Mg0.17Mg0.34Al1.66Si4O10(OH)2 + 6H+ + 4H2O = 1.660Al+3 + 0.510Mg+2 + 4H4SiO4
+     # # log_k     2.86
+# # Montmorillonite(MgCa)
+# # Ca0.17Mg0.34Al1.66Si4O10(OH)2 + 6H+ + 4H2O = 1.660Al+3 + 0.170Ca+2 + 0.340Mg+2 + 4H4SiO4
+     # # log_k     2.83
+# # # The divalent cations can be defined with the Gaines-Thomas convention...
+# # EXCHANGE_SPECIES
+# # # undefine the previous set...
+# # 0.17 Mg+2 +   X_montm_mg-0.34 = Mg0.17X_montm_mg; log_k -3.708e10
+# # 0.17 Ca+2 +   X_montm_mg-0.34 = Ca0.17X_montm_mg; log_k -4.222e10
+# # # write the Gaines-Thomas formulas...
+# # 0.34 Mg+2 + 2  X_montm_mg-0.34 = Mg0.34X_montm_mg2 ; log_k -7.416 # -0.297 #
+# # 0.34 Ca+2 + 2  X_montm_mg-0.34 = Ca0.34X_montm_mg2 ; log_k -8.444 # -0.811 #
+# # # The default exchanger X can be used, uncomment the following lines
+# # # redefine f_Na in the rate...
+# # RATES
+# # Montmorillonite
+# # 5  REM PARMS: 1 affinity, 2 m^2/mol, 3 roughness, 4 exponent
+# # 7  f_Na = (mol("NaX") / tot("X")) # when running with the default X exchange
+# # 10 if parm(1) = 1 then affinity = 1 else affinity = 1 - SR("Montmorillonite(MgNa)") / f_Na
+# # 20 rate = RATE_HERMANSKA("Montmorillonite") / f_Na
+# # 30 IF M > 0 THEN area = M * parm(2) * parm(3) * (M/M0)^parm(4) ELSE area = 0
+# # 40 SAVE area * rate * affinity * TIME
+# # -end
+# # # adapt log_k`s of the solids with default exchanger X:
+# # PHASES
+# # Montmorillonite(MgK)
+# # K0.34Mg0.34Al1.66Si4O10(OH)2 + 6H+ + 4H2O = 1.660Al+3 + 0.340K+ + 0.340Mg+2 + 4H4SiO4
+     # # log_k   2.6 # 3.41 - 0.7 * 0.34 = 3.17 expected, but is fraction-dependent, A&P, problems p. 305
+# # Montmorillonite(MgMg)
+# # Mg0.34(Mg0.34Al1.66Si4O10(OH)2)2 + 12 H+ + 8 H2O = 3.32 Al+3 + 1.02 Mg+2 + 8 H4SiO4
+     # # log_k   6.27 # 3.41 * 2 - 0.6 * 0.34 = 6.62
+# # Montmorillonite(MgCa)
+# # Ca0.34(Mg0.34Al1.66Si4O10(OH)2)2 + 12 H+ + 8 H2O = 3.32 Al+3 + 0.68 Mg+2 + 8 H4SiO4 + 0.34 Ca+2
+     # # log_k   6.2 #  3.41 * 2 - 0.8 * 0.34 = 6.55
+# # # in EXCHANGE 1, comment X_montm_mg, uncomment X...
+# END
+# SOLUTION 1
+# pH 7 charge
+# Na 1e-5
+# K  1e-5
+# Mg 1e-5
+# Ca 1e-5
+# END
+# # Give the solution composition for calculating the ininitial exchanger
+# SOLUTION 99
+# pH 7 charge
+# EQUILIBRIUM_PHASES 1
+# # solid solution of the end-members, SI = log10(fraction = 0.25)
+# Montmorillonite(MgNa)  -0.602  1e-2
+# Montmorillonite(MgCa)  -0.602  1e-2
+# Montmorillonite(MgK)   -0.602  1e-2
+# Montmorillonite(MgMg)  -0.602  1e-2
+# Kaolinite 0 0
+# SAVE solution 99
+# END
+# # # with Gapon only, initial exchanger can be defined explicitly
+# EXCHANGE 1
+# Na0.34X_montm_mg  1e-2
+# Ca0.17X_montm_mg  1e-2
+# K0.34X_montm_mg   1e-2
+# Mg0.17X_montm_mg  1e-2
+# END
+# USE solution 1
+# EQUILIBRIUM_PHASES 1
+# Kaolinite 0 0
+# # USE EXCHANGE 1 # with Gapon only, uncomment in KINETICS: # X_montm_mg -1
+# EXCHANGE 1
+# X_montm_mg Montmorillonite kin 1; -equil 99 # comment in KINETICS: # X_montm_mg -1
+# # X Montmorillonite kin 0.34; -equil 99 # default exchanger X, comment in KINETICS: # X_montm_mg -1
+# KINETICS 1
+# Montmorillonite
+# -formula Mg0.34Al1.66Si4O10(OH)2 1 # X_montm_mg -1
+# -m  4e-2
+# -parms 0  2.5e5 1 0.67
+# -step 30 100 1e3 1e4 2e4 2e4 3e4 3e4 3e4 3e4 1e5 1e5 1e5 3e5 6e5 1e6 3e6
+# -cvode true
+# INCREMENTAL_REACTIONS true
+# USER_GRAPH 4
+    # -headings               time Na K Mg Ca mm_diss
+    # -axis_titles            "Time / days" "Molality" "Montmorillonite dissolved / (mmol/kgw)"
+    # -axis_scale x_axis      auto auto auto auto log
+    # -axis_scale y_axis      auto auto auto auto log
+# 1  t = TOTAL_TIME / (3600 * 24) : put(t, 1)
+# 10 GRAPH_X t
+# 12 mg = tot("Mg") : if mg < 1e-24 then mg = 1e-24
+# 14 ca = tot("Ca") : if ca < 1e-24 then ca = 1e-24
+# 20 GRAPH_Y TOT("Na"), TOT("K"), mg, ca
+# 30 GRAPH_SY (4e-2 - kin("Montmorillonite")) * 1e3
+# END
+# USE solution 99; REACTION
+# USER_GRAPH 4; -connect_simulations false; -headings Solution_99
+# 1 t = get(1)
+# 10 plot_xy t, tot("Na"), symbol = Circle , symbol_size = 15, color = Red
+# 20 plot_xy t, tot("K"),  symbol = Circle , symbol_size = 15, color = Green
+# 30 plot_xy t, tot("Mg"), symbol = Circle , symbol_size = 15, color = Blue
+# 40 plot_xy t, tot("Ca"), symbol = Circle , symbol_size = 15, color = Orange
+# =============================================================================================
+#(a) means amorphous. (d) means disordered, or less crystalline.
+#(14A) refers to 14 angstrom spacing of clay planes. FeS(ppt),
+#precipitated, indicates an initial precipitate that is less crystalline.
+#Zn(OH)2(e) indicates a specific crystal form, epsilon.
+# =============================================================================================
+# For the reaction aA + bB = cC + dD,
+#   with delta_v = c*Vm(C) + d*Vm(D) - a*Vm(A) - b*Vm(B),
+# PHREEQC adds the pressure term to log_k: -= delta_v * (P - 1) / (2.3RT).
+#   Vm(A) is volume of A, cm3/mol, P is pressure, atm, R is the gas constant, T is Kelvin.
+# Gas-pressures and fugacity coefficients are calculated with Peng-Robinson's EOS.
+#   These binary interaction coefficients from Soreide and Whitson, 1992, FPE 77, 217 are
+#    hard-coded in calc_PR():
+#    kij    CH4    CO2    H2S    N2
+#    H2O    0.49   0.19   0.19   0.49
+#    but are overwritten by the data block GAS_BINARY_PARAMETERS of this file.
+# =============================================================================================
+# The molar volumes of solids are entered with
+#			 -Vm vm cm3/mol
+#    vm is the molar volume, cm3/mol (default), but dm3/mol and m3/mol are permitted.
+# Data for minerals' vm (= MW (g/mol) / rho (g/cm3)) are defined using rho from
+#   Deer, Howie and Zussman, The rock-forming minerals, Longman.
+#			   --------------------
+# Temperature- and pressure-dependent volumina of aqueous species are calculated with a Redlich-
+#   type equation (cf. Redlich and Meyer, Chem. Rev. 64, 221), from parameters entered with
+#			-Vm a1 a2 a3 a4 W a0 i1 i2 i3 i4
+# The volume (cm3/mol) is
+#    Vm(T, pb, I) = 41.84 * (a1 * 0.1 + a2 * 100 / (2600 + pb)  + a3 / (T - 228) +
+#			    a4 * 1e4 / (2600 + pb) / (T - 228) - W * QBrn)
+#		   + z^2 / 2 * Av * f(I^0.5)
+#		   + (i1 + i2 / (T - 228) + i3 * (T - 228)) * I^i4
+#   Volumina at I = 0 are obtained using supcrt92 formulas (Johnson et al., 1992, CG 18, 899).
+#   41.84 transforms cal/bar/mol into cm3/mol.
+#   pb is pressure in bar.
+#   W * QBrn is the energy of solvation, calculated from W and the pressure dependence of the Born equation,
+#	W is fitted on measured solution densities.
+#   z is charge of the solute species.
+#   Av is the Debye-H�ckel limiting slope (DH_AV in PHREEQC basic).
+#   a0 is the ion-size parameter in the extended Debye-H�ckel equation:
+#	f(I^0.5) = I^0.5 / (1 + a0 * DH_B * I^0.5),
+#	a0 = -gamma	  x for cations, = 0 for anions.
+# For details, consult ref. 1.
+# =============================================================================================
+# The viscosity is calculated with a (modified) Jones-Dole equation:
+#    viscos / viscos_0 = 1 + A * Sum(0.5 z_i m_i) + fan * Sum(B_i m_i + D_i m_i n_i)
+# Parameters are for calculating the B and D terms:
+#    -viscosity 9.35e-2 -8.31e-2 2.487e-2 4.49e-4 2.01e-2 1.570    0
+#    #	   b0	  b1	  b2	 d1	  d2	 d3	tan
+#    z_i is absolute charge number, m_i is molality of i
+#    B_i = b0 + b1 exp(-b2 * tc)
+#    fan = (2 - tan V_i / V_Cl-), corrects for the volume of anions
+#    D_i = d1 * exp(-d2 tc)
+#    n_i = (I^d3 * (1 + fI) + ((z_i^2 + z_i) / 2 � m_i)^d3) / (2 + fI), fI is an ionic strength term.
+# For details, consult ref. 4.
+#
+# ref. 1: Appelo, Parkhurst and Post, 2014. Geochim. Cosmochim. Acta 125, 49�67.
+# ref. 2: Procedures from ref. 1 using data compiled by Lalibert�, 2009, J. Chem. Eng. Data 54, 1725.
+# ref. 3: Appelo, 2017, Cem. Concr. Res. 101, 102-113.
+# ref. 4: Appelo and Parkhurst in prep., for details see subroutine viscosity in transport.cpp
+#
+# =============================================================================================
+# It remains the responsibility of the user to check the calculated results, for example with
+#   measured solubilities as a function of (P, T).