------------------------------------
Reading input data for simulation 1.
------------------------------------

	TITLE Example 11.--Transport and cation exchange.
	SOLUTION 0  CaCl2
	        units            mmol/kgw
	        temp             25.0
	        pH               7.0     charge
	        pe               12.5    O2(g)   -0.68
	        Ca               0.6
	        Cl               1.2
	SOLUTION 1-40  Initial solution for column
	        units            mmol/kgw
	        temp             25.0
	        pH               7.0     charge
	        pe               12.5    O2(g)   -0.68
	        Na               1.0
	        K                0.2
	        N(5)             1.2
	END
-----
TITLE
-----

 Example 11.--Transport and cation exchange.

-------------------------------------------
Beginning of initial solution calculations.
-------------------------------------------

Initial solution 0.	CaCl2

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                6.000e-04   6.000e-04
	Cl                1.200e-03   1.200e-03

----------------------------Description of solution----------------------------

                                       pH  =   6.995      Charge balance
                                       pe  =  13.632      Equilibrium with O2(g)
     Specific Conductance (µS/cm,  25°C)  = 152
                         Density (g/cm³)  =   0.99710
                               Volume (L)  =   1.00298
                        Viscosity (mPa s)  =   0.89067
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.800e-03
                       Mass of water (kg)  =   1.000e+00
                 Total alkalinity (eq/kg)  =   2.687e-21
                        Temperature (°C)  =  25.00
                  Electrical balance (eq)  =   4.260e-18
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   4
                                  Total H  = 1.110124e+02
                                  Total O  = 5.550675e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              1.057e-07   1.012e-07    -6.976    -6.995    -0.019      0.00
   OH-             1.049e-07   1.001e-07    -6.979    -7.000    -0.021     -4.10
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.07
Ca            6.000e-04
   Ca+2            6.000e-04   4.983e-04    -3.222    -3.302    -0.081    -18.10
   CaOH+           8.570e-10   8.175e-10    -9.067    -9.087    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.144e-03    -2.921    -2.941    -0.021     18.08
   HCl             3.981e-11   3.988e-11   -10.400   -10.399     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -44.404   -44.404     0.000     28.61
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     30.40

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(298 K,   1 atm)

  H2(g)           -41.30    -44.40   -3.10  H2
  H2O(g)           -1.50     -0.00    1.50  H2O
  O2(g)            -0.68     -3.57   -2.89  O2	 Pressure   0.2 atm, phi 1.000

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

Initial solution 1.	Initial solution for column

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 2.000e-04   2.000e-04
	N(5)              1.200e-03   1.200e-03
	Na                1.000e-03   1.000e-03

----------------------------Description of solution----------------------------

                                       pH  =   6.997      Charge balance
                                       pe  =  13.630      Equilibrium with O2(g)
     Specific Conductance (µS/cm,  25°C)  = 145
                         Density (g/cm³)  =   0.99711
                               Volume (L)  =   1.00301
                        Viscosity (mPa s)  =   0.89020
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.000e+00
                 Total alkalinity (eq/kg)  =  -9.247e-19
                        Temperature (°C)  =  25.00
                  Electrical balance (eq)  =  -9.243e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   3 (7 overall)
                                  Total H  = 1.110124e+02
                                  Total O  = 5.551035e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              1.045e-07   1.007e-07    -6.981    -6.997    -0.016      0.00
   OH-             1.045e-07   1.005e-07    -6.981    -6.998    -0.017     -4.10
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.07
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -44.404   -44.404     0.000     28.61
K             2.000e-04
   K+              2.000e-04   1.923e-04    -3.699    -3.716    -0.017      9.01
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     29.49
Na            1.000e-03
   Na+             1.000e-03   9.621e-04    -3.000    -3.017    -0.017     -1.47
   NaOH            9.666e-21   9.668e-21   -20.015   -20.015     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     30.40

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(298 K,   1 atm)

  H2(g)           -41.30    -44.40   -3.10  H2
  H2O(g)           -1.50     -0.00    1.50  H2O
  O2(g)            -0.68     -3.57   -2.89  O2	 Pressure   0.2 atm, phi 1.000

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

	EXCHANGE 1-40
	        equilibrate 1
	        X                0.0011
	END
-------------------------------------------------------
Beginning of initial exchange-composition calculations.
-------------------------------------------------------

Exchange 1.	

X                1.100e-03 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                5.507e-04   5.507e-04   5.006e-01    -0.017
	NaX               5.493e-04   5.493e-04   4.994e-01    -0.017

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 3.
------------------------------------

	SELECTED_OUTPUT
	        file            advect.sel
	        reset           false
	        totals          Na Cl K Ca
	USER_PUNCH
	  heading  Temperature Pressure Hyd_K
	  10 PUNCH TC, PRESSURE
	  20 PUNCH CALLBACK(cell_no, 0, "HYDRAULIC_K")
	END
------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 4.
------------------------------------

-------------------------------
End of Run after 13.67 Seconds.
-------------------------------

Number of threads:                                3
Number of MPI processes:                          1
MPI task number:                                  0
File prefix:                                      Species_c
Number of grid cells in the user's model:         40
Number of chemistry cells in the reaction module: 20
Number of components for transport:               9
       H2O        18.016
         H         1.008
         O        16.000
    Charge         1.000
        Ca        40.080
        Cl        35.453
         K        39.102
         N        14.007
        Na        22.990

Ca+2
    Charge: 2
    Dw:     7.93e-10
CaOH+
    Charge: 1
    Dw:     0
Cl-
    Charge: -1
    Dw:     2.03e-09
H+
    Charge: 1
    Dw:     9.31e-09
H2
    Charge: 0
    Dw:     5.13e-09
H2O
    Charge: 0
    Dw:     0
HCl
    Charge: 0
    Dw:     0
K+
    Charge: 1
    Dw:     1.96e-09
N2
    Charge: 0
    Dw:     1.96e-09
NH3
    Charge: 0
    Dw:     2.28e-09
NH4+
    Charge: 1
    Dw:     1.98e-09
NO2-
    Charge: -1
    Dw:     1.91e-09
NO3-
    Charge: -1
    Dw:     1.9e-09
Na+
    Charge: 1
    Dw:     1.33e-09
NaOH
    Charge: 0
    Dw:     0
O2
    Charge: 0
    Dw:     2.35e-09
OH-
    Charge: -1
    Dw:     5.27e-09

Time:                   10
Chemistry cell:         0
Grid cell(s) (0-based): 0 20 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      0
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 0.	Solution after simulation 1.
Using exchange 0.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	CaX2              1.100e-04   2.200e-04   1.000e+00    -0.080
	KX                5.579e-15   5.579e-15   2.536e-11    -0.020
	NaX               9.246e-22   9.246e-22   4.203e-18    -0.020

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                6.000e-04   1.196e-04
	Cl                1.200e-03   2.393e-04
	K                 2.946e-13   5.875e-14
	Na                2.522e-19   5.029e-20

----------------------------Description of solution----------------------------

                                       pH  =   7.079      Charge balance
                                       pe  =  13.969      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 136
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.19976
                        Viscosity (mPa s)  =   1.00227
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.800e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.294e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   8.803e-19
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   2
                                  Total H  = 2.213642e+01
                                  Total O  = 1.106831e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.713e-08   8.339e-08    -7.060    -7.079    -0.019      0.00
   OH-             8.613e-08   8.216e-08    -7.065    -7.085    -0.020     -4.31
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            6.000e-04
   Ca+2            6.000e-04   4.991e-04    -3.222    -3.302    -0.080    -18.16
   CaOH+           1.041e-09   9.933e-10    -8.983    -9.003    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.145e-03    -2.921    -2.941    -0.020     17.92
   HCl             3.389e-11   3.395e-11   -10.470   -10.469     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             2.946e-13
   K+              2.946e-13   2.811e-13   -12.531   -12.551    -0.020      8.85
Na            2.522e-19
   Na+             2.522e-19   2.407e-19   -18.598   -18.618    -0.020     -1.74
   NaOH            1.977e-36   1.978e-36   -35.704   -35.704     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite          -23.13    -21.56    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite         -16.37    -15.49    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.001 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         1
Grid cell(s) (0-based): 1 21 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      1
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 1.	Solution after simulation 1.
Using exchange 1.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	CaX2              1.100e-04   2.200e-04   1.000e+00    -0.080
	KX                1.224e-12   1.224e-12   5.566e-09    -0.020
	NaX               1.544e-18   1.544e-18   7.020e-15    -0.020

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                6.000e-04   1.198e-04
	Cl                1.200e-03   2.396e-04
	K                 6.466e-11   1.291e-11
	Na                4.213e-16   8.410e-17

----------------------------Description of solution----------------------------

                                       pH  =   7.079      Charge balance
                                       pe  =  13.969      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 136
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00227
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.800e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   2.526e-20
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   4.854e-19
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   5
                                  Total H  = 2.216319e+01
                                  Total O  = 1.108170e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.713e-08   8.339e-08    -7.060    -7.079    -0.019      0.00
   OH-             8.613e-08   8.216e-08    -7.065    -7.085    -0.020     -4.31
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            6.000e-04
   Ca+2            6.000e-04   4.991e-04    -3.222    -3.302    -0.080    -18.16
   CaOH+           1.041e-09   9.933e-10    -8.983    -9.003    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.145e-03    -2.921    -2.941    -0.020     17.92
   HCl             3.389e-11   3.395e-11   -10.470   -10.469     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             6.466e-11
   K+              6.466e-11   6.169e-11   -10.189   -10.210    -0.020      8.85
Na            4.213e-16
   Na+             4.213e-16   4.022e-16   -15.375   -15.396    -0.020     -1.74
   NaOH            3.302e-33   3.303e-33   -32.481   -32.481     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite          -19.90    -18.34    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite         -14.03    -13.15    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.005 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         2
Grid cell(s) (0-based): 2 22 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      2
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 2.	Solution after simulation 1.
Using exchange 2.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	CaX2              1.100e-04   2.200e-04   1.000e+00    -0.080
	KX                1.599e-10   1.599e-10   7.270e-07    -0.020
	NaX               1.914e-15   1.914e-15   8.702e-12    -0.020

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                6.000e-04   1.198e-04
	Cl                1.200e-03   2.396e-04
	K                 8.446e-09   1.686e-09
	Na                5.222e-13   1.042e-13

----------------------------Description of solution----------------------------

                                       pH  =   7.079      Charge balance
                                       pe  =  13.969      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 136
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00227
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.800e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   2.311e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   1.546e-18
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   6
                                  Total H  = 2.216319e+01
                                  Total O  = 1.108170e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.713e-08   8.339e-08    -7.060    -7.079    -0.019      0.00
   OH-             8.613e-08   8.216e-08    -7.065    -7.085    -0.020     -4.31
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            6.000e-04
   Ca+2            6.000e-04   4.991e-04    -3.222    -3.302    -0.080    -18.16
   CaOH+           1.041e-09   9.933e-10    -8.983    -9.003    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.145e-03    -2.921    -2.941    -0.020     17.92
   HCl             3.389e-11   3.395e-11   -10.470   -10.469     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             8.446e-09
   K+              8.446e-09   8.058e-09    -8.073    -8.094    -0.020      8.85
Na            5.222e-13
   Na+             5.222e-13   4.985e-13   -12.282   -12.302    -0.020     -1.74
   NaOH            4.093e-30   4.095e-30   -29.388   -29.388     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite          -16.81    -15.24    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite         -11.91    -11.03    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

-------------------------------
End of Run after 14.01 Seconds.
-------------------------------

Time:                   10
Chemistry cell:         3
Grid cell(s) (0-based): 3 23 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      3
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 3.	Solution after simulation 1.
Using exchange 3.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	CaX2              1.100e-04   2.200e-04   9.999e-01    -0.080
	KX                1.683e-08   1.683e-08   7.649e-05    -0.020
	NaX               2.243e-12   2.243e-12   1.019e-08    -0.020

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                5.996e-04   1.197e-04
	Cl                1.200e-03   2.396e-04
	K                 8.884e-07   1.774e-07
	Na                6.115e-10   1.221e-10

----------------------------Description of solution----------------------------

                                       pH  =   7.079      Charge balance
                                       pe  =  13.969      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 137
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00227
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.800e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   2.769e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   6.231e-18
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   6
                                  Total H  = 2.216319e+01
                                  Total O  = 1.108170e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.713e-08   8.339e-08    -7.060    -7.079    -0.019      0.00
   OH-             8.613e-08   8.216e-08    -7.065    -7.085    -0.020     -4.31
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            5.996e-04
   Ca+2            5.996e-04   4.988e-04    -3.222    -3.302    -0.080    -18.16
   CaOH+           1.040e-09   9.925e-10    -8.983    -9.003    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.145e-03    -2.921    -2.941    -0.020     17.92
   HCl             3.389e-11   3.395e-11   -10.470   -10.469     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             8.884e-07
   K+              8.884e-07   8.476e-07    -6.051    -6.072    -0.020      8.85
Na            6.115e-10
   Na+             6.115e-10   5.837e-10    -9.214    -9.234    -0.020     -1.74
   NaOH            4.793e-27   4.795e-27   -26.319   -26.319     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite          -13.74    -12.17    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -9.89     -9.01    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.015 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         4
Grid cell(s) (0-based): 4 24 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      4
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 4.	Solution after simulation 1.
Using exchange 4.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	CaX2              1.092e-04   2.184e-04   9.927e-01    -0.079
	KX                1.599e-06   1.599e-06   7.266e-03    -0.020
	NaX               2.604e-09   2.604e-09   1.184e-05    -0.020

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                5.588e-04   1.116e-04
	Cl                1.200e-03   2.396e-04
	K                 8.176e-05   1.632e-05
	Na                6.878e-07   1.373e-07

----------------------------Description of solution----------------------------

                                       pH  =   7.079      Charge balance
                                       pe  =  13.969      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 138
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00221
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.759e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   3.476e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   2.791e-16
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   7
                                  Total H  = 2.216319e+01
                                  Total O  = 1.108170e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.706e-08   8.335e-08    -7.060    -7.079    -0.019      0.00
   OH-             8.612e-08   8.220e-08    -7.065    -7.085    -0.020     -4.31
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            5.588e-04
   Ca+2            5.588e-04   4.657e-04    -3.253    -3.332    -0.079    -18.16
   CaOH+           9.711e-10   9.272e-10    -9.013    -9.033    -0.020     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.146e-03    -2.921    -2.941    -0.020     17.92
   HCl             3.390e-11   3.396e-11   -10.470   -10.469     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             8.176e-05
   K+              8.176e-05   7.805e-05    -4.087    -4.108    -0.020      8.85
Na            6.878e-07
   Na+             6.878e-07   6.570e-07    -6.163    -6.182    -0.020     -1.74
   NaOH            5.397e-24   5.399e-24   -23.268   -23.268     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite          -10.69     -9.12    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.92     -7.05    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.019 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         5
Grid cell(s) (0-based): 5 25 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      5
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 5.	Solution after simulation 1.
Using exchange 5.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.063e-04   1.063e-04   4.831e-01    -0.017
	CaX2              5.574e-05   1.115e-04   5.068e-01    -0.066
	NaX               2.219e-06   2.219e-06   1.009e-02    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                1.091e-05   2.179e-06
	Cl                1.200e-03   2.396e-04
	K                 1.063e-03   2.123e-04
	Na                1.147e-04   2.289e-05

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.966      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 154
                         Density (g/cm³)  =   0.99831
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00173
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.211e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   3.873e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   1.029e-13
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   8
                                  Total H  = 2.216295e+01
                                  Total O  = 1.108158e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   OH-             8.598e-08   8.269e-08    -7.066    -7.083    -0.017     -4.32
   H+              8.596e-08   8.286e-08    -7.066    -7.082    -0.016      0.00
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            1.091e-05
   Ca+2            1.091e-05   9.367e-06    -4.962    -5.028    -0.066    -18.19
   CaOH+           1.950e-11   1.876e-11   -10.710   -10.727    -0.017     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.154e-03    -2.921    -2.938    -0.017     17.91
   HCl             3.397e-11   3.401e-11   -10.469   -10.468     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.063e-03
   K+              1.063e-03   1.023e-03    -2.973    -2.990    -0.017      8.84
Na            1.147e-04
   Na+             1.147e-04   1.103e-04    -3.941    -3.957    -0.017     -1.75
   NaOH            9.120e-22   9.122e-22   -21.040   -21.040     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.573     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite           -8.46     -6.89    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -6.80     -5.93    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.022 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         6
Grid cell(s) (0-based): 6 26 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      6
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 6.	Solution after simulation 1.
Using exchange 6.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.816e-04   1.816e-04   8.256e-01    -0.017
	NaX               3.834e-05   3.834e-05   1.743e-01    -0.017
	CaX2              7.889e-09   1.578e-08   7.172e-05    -0.066

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                1.541e-10   3.077e-11
	Cl                1.200e-03   2.396e-04
	K                 5.741e-04   1.146e-04
	Na                6.259e-04   1.250e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.966      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 144
                         Density (g/cm³)  =   0.99830
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00179
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -8.230e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -7.166e-16
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   6
                                  Total H  = 2.216281e+01
                                  Total O  = 1.108151e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   OH-             8.597e-08   8.270e-08    -7.066    -7.082    -0.017     -4.32
   H+              8.594e-08   8.285e-08    -7.066    -7.082    -0.016      0.00
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            1.541e-10
   Ca+2            1.541e-10   1.324e-10    -9.812    -9.878    -0.066    -18.19
   CaOH+           2.755e-16   2.651e-16   -15.560   -15.577    -0.017     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.154e-03    -2.921    -2.938    -0.017     17.91
   HCl             3.397e-11   3.401e-11   -10.469   -10.468     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             5.741e-04
   K+              5.741e-04   5.523e-04    -3.241    -3.258    -0.017      8.84
Na            6.259e-04
   Na+             6.259e-04   6.024e-04    -3.203    -3.220    -0.017     -1.75
   NaOH            4.979e-21   4.980e-21   -20.303   -20.303     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.573     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite           -7.72     -6.16    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.07     -6.20    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.002 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         7
Grid cell(s) (0-based): 7 27 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      7
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 7.	Solution after simulation 1.
Using exchange 7.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.243e-04   1.243e-04   5.650e-01    -0.017
	NaX               9.569e-05   9.569e-05   4.350e-01    -0.017
	CaX2              1.240e-14   2.480e-14   1.127e-10    -0.066

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                9.129e-17   1.823e-17
	Cl                1.200e-03   2.396e-04
	K                 2.412e-04   4.815e-05
	Na                9.588e-04   1.914e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.966      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 137
                         Density (g/cm³)  =   0.99830
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00183
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =   5.075e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   7.587e-13
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   4
                                  Total H  = 2.216293e+01
                                  Total O  = 1.108157e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   OH-             8.597e-08   8.270e-08    -7.066    -7.082    -0.017     -4.32
   H+              8.594e-08   8.285e-08    -7.066    -7.082    -0.016      0.00
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            9.129e-17
   Ca+2            9.129e-17   7.840e-17   -16.040   -16.106    -0.066    -18.19
   CaOH+           1.632e-22   1.570e-22   -21.787   -21.804    -0.017     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.154e-03    -2.921    -2.938    -0.017     17.91
   HCl             3.397e-11   3.401e-11   -10.469   -10.468     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             2.412e-04
   K+              2.412e-04   2.320e-04    -3.618    -3.634    -0.017      8.84
Na            9.588e-04
   Na+             9.588e-04   9.228e-04    -3.018    -3.035    -0.017     -1.75
   NaOH            7.627e-21   7.629e-21   -20.118   -20.118     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.573     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite           -7.54     -5.97    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.45     -6.57    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.006 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         8
Grid cell(s) (0-based): 8 28 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      8
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 8.	Solution after simulation 1.
Using exchange 8.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.116e-04   1.116e-04   5.074e-01    -0.017
	NaX               1.084e-04   1.084e-04   4.926e-01    -0.017
	CaX2              2.913e-21   5.827e-21   2.648e-17    -0.066

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Ca                1.820e-23   3.634e-24
	Cl                1.200e-03   2.396e-04
	K                 1.995e-04   3.983e-05
	Na                1.000e-03   1.997e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.966      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 136
                         Density (g/cm³)  =   0.99830
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00184
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -5.853e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   1.476e-13
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   5
                                  Total H  = 2.216295e+01
                                  Total O  = 1.108158e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   OH-             8.597e-08   8.270e-08    -7.066    -7.082    -0.017     -4.32
   H+              8.594e-08   8.285e-08    -7.066    -7.082    -0.016      0.00
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Ca            1.820e-23
   Ca+2            1.820e-23   1.563e-23   -22.740   -22.806    -0.066    -18.19
   CaOH+           3.255e-29   3.131e-29   -28.488   -28.504    -0.017     (0)  
Cl            1.200e-03
   Cl-             1.200e-03   1.154e-03    -2.921    -2.938    -0.017     17.91
   HCl             3.397e-11   3.401e-11   -10.469   -10.468     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.995e-04
   K+              1.995e-04   1.919e-04    -3.700    -3.717    -0.017      8.84
Na            1.000e-03
   Na+             1.000e-03   9.629e-04    -3.000    -3.016    -0.017     -1.75
   NaOH            7.959e-21   7.961e-21   -20.099   -20.099     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.573     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite           -7.52     -5.95    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.53     -6.65    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

-------------------------------
End of Run after 14.01 Seconds.
-------------------------------

Time:                   10
Chemistry cell:         9
Grid cell(s) (0-based): 9 29 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      9
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 9.	Solution after simulation 1.
Using exchange 9.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.035e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.965e-01    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	Cl                1.200e-03   2.396e-04
	K                 1.970e-04   3.932e-05
	Na                1.003e-03   2.002e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.966      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 136
                         Density (g/cm³)  =   0.99830
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00184
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -4.369e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -1.850e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216295e+01
                                  Total O  = 1.108158e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   OH-             8.597e-08   8.270e-08    -7.066    -7.082    -0.017     -4.32
   H+              8.594e-08   8.285e-08    -7.066    -7.082    -0.016      0.00
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
Cl            1.200e-03
   Cl-             1.200e-03   1.154e-03    -2.921    -2.938    -0.017     17.91
   HCl             3.397e-11   3.401e-11   -10.469   -10.468     0.001     (0)  
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.970e-04
   K+              1.970e-04   1.895e-04    -3.706    -3.722    -0.017      8.84
Na            1.003e-03
   Na+             1.003e-03   9.653e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.979e-21   7.981e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.675e-04   2.676e-04    -3.573    -3.573     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  Halite           -7.52     -5.95    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.53     -6.66    0.87  KCl

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.016 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         10
Grid cell(s) (0-based): 10 30 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      10
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 10.	Solution after simulation 1.
Using exchange 10.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.021 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         11
Grid cell(s) (0-based): 11 31 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      11
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 11.	Solution after simulation 1.
Using exchange 11.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.002 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         12
Grid cell(s) (0-based): 12 32 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      12
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 12.	Solution after simulation 1.
Using exchange 12.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.007 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         13
Grid cell(s) (0-based): 13 33 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      13
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 13.	Solution after simulation 1.
Using exchange 13.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

-------------------------------
End of Run after 14.01 Seconds.
-------------------------------

Time:                   10
Chemistry cell:         14
Grid cell(s) (0-based): 14 34 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      14
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 14.	Solution after simulation 1.
Using exchange 14.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.014 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         15
Grid cell(s) (0-based): 15 35 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      15
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 15.	Solution after simulation 1.
Using exchange 15.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.018 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         16
Grid cell(s) (0-based): 16 36 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      16
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 16.	Solution after simulation 1.
Using exchange 16.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.022 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         17
Grid cell(s) (0-based): 17 37 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      17
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 17.	Solution after simulation 1.
Using exchange 17.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.026 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         18
Grid cell(s) (0-based): 18 38 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      18
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 18.	Solution after simulation 1.
Using exchange 18.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.029 Seconds.
--------------------------------

Time:                   10
Chemistry cell:         19
Grid cell(s) (0-based): 19 39 
------------------------------------
Reading input data for simulation 1.
------------------------------------

	RUN_CELLS
	  -start_time 777600
	  -time_step  86400
	  -cells      19
	END
--------------------------
Beginning of run as cells.
--------------------------

-----------------------------------------
Beginning of batch-reaction calculations.
-----------------------------------------

Reaction step 1.

Using solution 19.	Solution after simulation 1.
Using exchange 19.	Exchange assemblage after simulation 1.

-----------------------------Exchange composition------------------------------

X                2.200e-04 mol

	                               Equiv-    Equivalent      Log 
	Species             Moles      alents      Fraction     Gamma

	KX                1.108e-04   1.108e-04   5.034e-01    -0.017
	NaX               1.092e-04   1.092e-04   4.966e-01    -0.017
	NH4X              3.892e-65   3.892e-65   1.769e-61    -0.017

-----------------------------Solution composition------------------------------

	Elements           Molality       Moles

	K                 1.969e-04   3.931e-05
	N                 1.200e-03   2.396e-04
	Na                1.003e-03   2.003e-04

----------------------------Description of solution----------------------------

                                       pH  =   7.082      Charge balance
                                       pe  =  13.967      Adjusted to redox equilibrium
     Specific Conductance (µS/cm,  20°C)  = 131
                         Density (g/cm³)  =   0.99832
                               Volume (L)  =   0.20000
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.996e-01
                 Total alkalinity (eq/kg)  =  -1.012e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -4.165e-15
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.216266e+01
                                  Total O  = 1.108216e+01

----------------------------Distribution of species----------------------------

                                               Log       Log       Log    mole V
   Species          Molality    Activity  Molality  Activity     Gamma   cm³/mol

   H+              8.595e-08   8.287e-08    -7.066    -7.082    -0.016      0.00
   OH-             8.595e-08   8.268e-08    -7.066    -7.083    -0.017     -4.32
   H2O             5.551e+01   1.000e+00     1.744    -0.000     0.000     18.05
H(0)          0.000e+00
   H2              0.000e+00   0.000e+00   -45.225   -45.225     0.000     28.61
K             1.969e-04
   K+              1.969e-04   1.894e-04    -3.706    -3.723    -0.017      8.84
N(-3)         0.000e+00
   NH4+            0.000e+00   0.000e+00   -64.054   -64.071    -0.017     17.80
   NH3             0.000e+00   0.000e+00   -66.390   -66.390     0.000     24.19
N(0)          5.774e-20
   N2              2.887e-20   2.888e-20   -19.540   -19.539     0.000     29.29
N(3)          1.258e-16
   NO2-            1.258e-16   1.210e-16   -15.900   -15.917    -0.017     24.69
N(5)          1.200e-03
   NO3-            1.200e-03   1.154e-03    -2.921    -2.938    -0.017     28.98
Na            1.003e-03
   Na+             1.003e-03   9.654e-04    -2.999    -3.015    -0.017     -1.75
   NaOH            7.978e-21   7.980e-21   -20.098   -20.098     0.000     (0)  
O(0)          5.351e-04
   O2              2.676e-04   2.676e-04    -3.573    -3.572     0.000     29.98

------------------------------Saturation indices-------------------------------

  Phase               SI** log IAP   log K(293 K,   2 atm)

  H2(g)           -42.14    -45.22   -3.09  H2
  H2O(g)           -1.63     -0.00    1.63  H2O
  N2(g)           -16.39    -19.54   -3.14  N2
  NH3(g)          -68.29    -66.39    1.90  NH3
  O2(g)            -0.72     -3.57   -2.85  O2

**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
  For ideal gases, phi = 1.

------------------
End of simulation.
------------------

------------------------------------
Reading input data for simulation 2.
------------------------------------

--------------------------------
End of Run after 14.033 Seconds.
--------------------------------

