------------------------------------
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 3.777 Seconds.
-------------------------------

SetValues variables and units:
  0                                      Concentrations   mol/L
  1                                         DensityUser   kg L-1
  2                                          FilePrefix   prefix
  3                                   NthSelectedOutput   id
  4                                      SaturationUser   unitless
  5                                                Time   s
  6                                            TimeStep   s
  7                                            Porosity   unitless
  8                                            Pressure   atm
  9                                    SelectedOutputOn   bool
  10                                         Temperature   C
GetValues variables and units:
  0                                      ComponentCount   count
  1                                          Components   names
  2                                      Concentrations   mol/L
  3                                   DensityCalculated   kg L-1
  4                                         ErrorString   error
  5                                          FilePrefix   prefix
  6                                                 Gfw   g mol-1
  7                                       GridCellCount   count
  8                                SaturationCalculated   unitless
  9                                      SelectedOutput   user specified
  10                           SelectedOutputColumnCount   count
  11                                 SelectedOutputCount   count
  12                              SelectedOutputHeadings   names
  13                              SelectedOutputRowCount   count
  14                                      SolutionVolume   L
  15                                                Time   s
  16                                            TimeStep   s
  17                     CurrentSelectedOutputUserNumber   id
  18                                            Porosity   unitless
  19                                            Pressure   atm
  20                                    SelectedOutputOn   bool
  21                                         Temperature   C
  22                                           Viscosity   mPa s
  23                exchange_X_species_log_molality_CaX2   log mol kgw-1
  24                  exchange_X_species_log_molality_KX   log mol kgw-1
  25                 exchange_X_species_log_molality_NaX   log mol kgw-1
  26                  exchange_X_species_log_molality_X-   log mol kgw-1
  27                           exchange_total_molality_X   mol kgw-1
  28                                 solution_alkalinity   eq kgw-1
  29                             solution_charge_balance   eq kgw-1
  30                             solution_ionic_strength   mol kgw-1
  31                                         solution_pe   -
  32                              solution_percent_error   -
  33                                         solution_ph   -
  34                       solution_specific_conductance   uS cm-1
  35                          solution_total_molality_Ca   mol kgw-1
  36                          solution_total_molality_Cl   mol kgw-1
  37                           solution_total_molality_H   mol kgw-1
  38                        solution_total_molality_H(0)   mol kgw-1
  39                        solution_total_molality_H(1)   mol kgw-1
  40                           solution_total_molality_K   mol kgw-1
  41                           solution_total_molality_N   mol kgw-1
  42                       solution_total_molality_N(-3)   mol kgw-1
  43                        solution_total_molality_N(0)   mol kgw-1
  44                        solution_total_molality_N(3)   mol kgw-1
  45                        solution_total_molality_N(5)   mol kgw-1
  46                          solution_total_molality_Na   mol kgw-1
  47                           solution_total_molality_O   mol kgw-1
  48                       solution_total_molality_O(-2)   mol kgw-1
  49                        solution_total_molality_O(0)   mol kgw-1
  50                                 solution_water_mass   kg
File prefix:                                      AdvectBMI_cpp
Number of grid cells in the user's model:         40
Number of components for transport:               8
         H    1.008
         O    16
    Charge    1
        Ca    40.08
        Cl    35.453
         K    39.102
         N    14.0067
        Na    22.9898

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.136e-21
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   8.144e-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 4.153 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.234e-12   1.234e-12   5.609e-09    -0.020
	NaX               1.560e-18   1.560e-18   7.089e-15    -0.020

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

	Elements           Molality       Moles

	Ca                6.000e-04   1.196e-04
	Cl                1.200e-03   2.393e-04
	K                 6.516e-11   1.299e-11
	Na                4.254e-16   8.482e-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.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)  =   1.259e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   5.003e-19
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   5
                                  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             6.516e-11
   K+              6.516e-11   6.217e-11   -10.186   -10.206    -0.020      8.85
Na            4.254e-16
   Na+             4.254e-16   4.061e-16   -15.371   -15.391    -0.020     -1.74
   NaOH            3.334e-33   3.336e-33   -32.477   -32.477     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.33    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite         -14.02    -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 4.158 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.614e-10   1.614e-10   7.337e-07    -0.020
	NaX               1.938e-15   1.938e-15   8.808e-12    -0.020

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

	Elements           Molality       Moles

	Ca                6.000e-04   1.196e-04
	Cl                1.200e-03   2.393e-04
	K                 8.524e-09   1.700e-09
	Na                5.285e-13   1.054e-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.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)  =   1.355e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   1.227e-18
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   6
                                  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             8.524e-09
   K+              8.524e-09   8.133e-09    -8.069    -8.090    -0.020      8.85
Na            5.285e-13
   Na+             5.285e-13   5.045e-13   -12.277   -12.297    -0.020     -1.74
   NaOH            4.143e-30   4.145e-30   -29.383   -29.383     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.80    -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 4.163 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.700e-08   1.700e-08   7.725e-05    -0.020
	NaX               2.274e-12   2.274e-12   1.034e-08    -0.020

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

	Elements           Molality       Moles

	Ca                5.996e-04   1.196e-04
	Cl                1.200e-03   2.393e-04
	K                 8.972e-07   1.789e-07
	Na                6.201e-10   1.236e-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.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)  =  -1.923e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   6.122e-18
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   6
                                  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            5.996e-04
   Ca+2            5.996e-04   4.987e-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.972e-07
   K+              8.972e-07   8.560e-07    -6.047    -6.068    -0.020      8.85
Na            6.201e-10
   Na+             6.201e-10   5.920e-10    -9.208    -9.228    -0.020     -1.74
   NaOH            4.861e-27   4.863e-27   -26.313   -26.313     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.88     -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 4.169 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.926e-01    -0.079
	KX                1.616e-06   1.616e-06   7.344e-03    -0.020
	NaX               2.646e-09   2.646e-09   1.203e-05    -0.020

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

	Elements           Molality       Moles

	Ca                5.583e-04   1.113e-04
	Cl                1.200e-03   2.393e-04
	K                 8.262e-05   1.647e-05
	Na                6.987e-07   1.393e-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.19976
                        Viscosity (mPa s)  =   1.00221
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.758e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -2.589e-21
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   2.787e-16
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   7
                                  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.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.583e-04
   Ca+2            5.583e-04   4.654e-04    -3.253    -3.332    -0.079    -18.16
   CaOH+           9.704e-10   9.265e-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.262e-05
   K+              8.262e-05   7.886e-05    -4.083    -4.103    -0.020      8.85
Na            6.987e-07
   Na+             6.987e-07   6.674e-07    -6.156    -6.176    -0.020     -1.74
   NaOH            5.482e-24   5.485e-24   -23.261   -23.261     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.68     -9.12    1.57  NaCl
  O2(g)            -0.72     -3.57   -2.85  O2
  Sylvite          -7.92     -7.04    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 4.173 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.072e-04   1.072e-04   4.871e-01    -0.017
	CaX2              5.530e-05   1.106e-04   5.027e-01    -0.066
	NaX               2.252e-06   2.252e-06   1.023e-02    -0.017

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

	Elements           Molality       Moles

	Ca                1.065e-05   2.124e-06
	Cl                1.200e-03   2.393e-04
	K                 1.063e-03   2.120e-04
	Na                1.154e-04   2.301e-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.19976
                        Viscosity (mPa s)  =   1.00173
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.211e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =   6.841e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   9.176e-14
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   8
                                  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

   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.065e-05
   Ca+2            1.065e-05   9.140e-06    -4.973    -5.039    -0.066    -18.19
   CaOH+           1.903e-11   1.831e-11   -10.721   -10.737    -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.154e-04
   Na+             1.154e-04   1.110e-04    -3.938    -3.955    -0.017     -1.75
   NaOH            9.177e-22   9.179e-22   -21.037   -21.037     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 4.178 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.815e-04   1.815e-04   8.250e-01    -0.017
	NaX               3.850e-05   3.850e-05   1.750e-01    -0.017
	CaX2              7.697e-09   1.539e-08   6.997e-05    -0.066

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

	Elements           Molality       Moles

	Ca                1.499e-10   2.988e-11
	Cl                1.200e-03   2.393e-04
	K                 5.727e-04   1.142e-04
	Na                6.273e-04   1.251e-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.19976
                        Viscosity (mPa s)  =   1.00179
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -5.540e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -7.026e-16
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   6
                                  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

   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.499e-10
   Ca+2            1.499e-10   1.287e-10    -9.824    -9.890    -0.066    -18.19
   CaOH+           2.679e-16   2.578e-16   -15.572   -15.589    -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.727e-04
   K+              5.727e-04   5.509e-04    -3.242    -3.259    -0.017      8.84
Na            6.273e-04
   Na+             6.273e-04   6.037e-04    -3.203    -3.219    -0.017     -1.75
   NaOH            4.990e-21   4.992e-21   -20.302   -20.302     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 4.155 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.242e-04   1.242e-04   5.647e-01    -0.017
	NaX               9.577e-05   9.577e-05   4.353e-01    -0.017
	CaX2              1.209e-14   2.418e-14   1.099e-10    -0.066

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

	Elements           Molality       Moles

	Ca                8.889e-17   1.773e-17
	Cl                1.200e-03   2.393e-04
	K                 2.409e-04   4.804e-05
	Na                9.591e-04   1.912e-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.19976
                        Viscosity (mPa s)  =   1.00183
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =   4.502e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   7.519e-13
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   4
                                  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

   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            8.889e-17
   Ca+2            8.889e-17   7.634e-17   -16.051   -16.117    -0.066    -18.19
   CaOH+           1.589e-22   1.529e-22   -21.799   -21.816    -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.409e-04
   K+              2.409e-04   2.318e-04    -3.618    -3.635    -0.017      8.84
Na            9.591e-04
   Na+             9.591e-04   9.230e-04    -3.018    -3.035    -0.017     -1.75
   NaOH            7.629e-21   7.632e-21   -20.118   -20.117     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 4.16 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.857e-21   5.713e-21   2.597e-17    -0.066

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

	Elements           Molality       Moles

	Ca                1.785e-23   3.559e-24
	Cl                1.200e-03   2.393e-04
	K                 1.995e-04   3.978e-05
	Na                1.000e-03   1.995e-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.19976
                        Viscosity (mPa s)  =   1.00184
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =   5.157e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =   1.445e-13
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =   0.00
                               Iterations  =   5
                                  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

   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.785e-23
   Ca+2            1.785e-23   1.533e-23   -22.748   -22.815    -0.066    -18.19
   CaOH+           3.191e-29   3.070e-29   -28.496   -28.513    -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 4.166 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.393e-04
	K                 1.970e-04   3.928e-05
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00184
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =   2.876e-19
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -1.628e-17
 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

   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 4.172 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.181 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.155 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.161 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.166 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.171 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.176 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.183 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.185 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.187 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.926e-05
	N                 1.200e-03   2.393e-04
	Na                1.003e-03   2.000e-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.19976
                        Viscosity (mPa s)  =   1.00175
                        Activity of water  =   1.000
                 Ionic strength (mol/kgw)  =   1.200e-03
                       Mass of water (kg)  =   1.994e-01
                 Total alkalinity (eq/kg)  =  -3.942e-18
                        Temperature (°C)  =  20.00
                           Pressure (atm)  =  2.00
                  Electrical balance (eq)  =  -6.789e-17
 Percent error, 100*(Cat-|An|)/(Cat+|An|)  =  -0.00
                               Iterations  =   2
                                  Total H  = 2.213587e+01
                                  Total O  = 1.106876e+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 4.189 Seconds.
-------------------------------

Current distribution of cells for workers
Worker      First cell        Last Cell
0           0                 5
1           6                 10
2           11                 19
Selected output sequence number: 0
Selected output user number:     1
Cell number 0
     Calculated Density: 0.998305
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 4.07186e-18
          3 Ca: 0.000598214
          4 Cl: 0.00119643
          5 K: 2.93754e-13
          6 N: 0
          7 Na: 2.51428e-19
     Selected output: 
          0 Na(mol/kgw): 2.52179e-19
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 2.94631e-13
          3 Ca(mol/kgw): 0.0006
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 1
     Calculated Density: 0.998305
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 2.50163e-18
          3 Ca: 0.000598214
          4 Cl: 0.00119643
          5 K: 6.4971e-11
          6 N: 0
          7 Na: 4.24098e-16
     Selected output: 
          0 Na(mol/kgw): 4.25364e-16
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 6.51649e-11
          3 Ca(mol/kgw): 0.0006
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 2
     Calculated Density: 0.998305
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 6.13391e-18
          3 Ca: 0.00059821
          4 Cl: 0.00119643
          5 K: 8.49904e-09
          6 N: 0
          7 Na: 5.26944e-13
     Selected output: 
          0 Na(mol/kgw): 5.28517e-13
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 8.52441e-09
          3 Ca(mol/kgw): 0.000599996
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 3
     Calculated Density: 0.998305
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 3.06082e-17
          3 Ca: 0.000597767
          4 Cl: 0.00119643
          5 K: 8.94576e-07
          6 N: 0
          7 Na: 6.18236e-10
     Selected output: 
          0 Na(mol/kgw): 6.20082e-10
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 8.97246e-07
          3 Ca(mol/kgw): 0.000599551
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 4
     Calculated Density: 0.998305
     Calculated Volume:  0.199759
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 1.39372e-15
          3 Ca: 0.000556681
          4 Cl: 0.00119643
          5 K: 8.23713e-05
          6 N: 0
          7 Na: 6.96665e-07
     Selected output: 
          0 Na(mol/kgw): 6.98744e-07
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 8.26171e-05
          3 Ca(mol/kgw): 0.000558342
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 5
     Calculated Density: 0.998306
     Calculated Volume:  0.199762
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 4.58801e-13
          3 Ca: 1.0618e-05
          4 Cl: 0.00119643
          5 K: 0.00106014
          6 N: 0
          7 Na: 0.00011505
     Selected output: 
          0 Na(mol/kgw): 0.000115394
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 0.00106331
          3 Ca(mol/kgw): 1.06497e-05
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 6
     Calculated Density: 0.998303
     Calculated Volume:  0.199761
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: -3.51279e-15
          3 Ca: 1.49424e-10
          4 Cl: 0.00119643
          5 K: 0.000570961
          6 N: 0
          7 Na: 0.000625467
     Selected output: 
          0 Na(mol/kgw): 0.000627334
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 0.000572666
          3 Ca(mol/kgw): 1.4987e-10
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 7
     Calculated Density: 0.998301
     Calculated Volume:  0.199761
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 3.75972e-12
          3 Ca: 8.86289e-17
          4 Cl: 0.00119643
          5 K: 0.000240192
          6 N: 0
          7 Na: 0.000956237
     Selected output: 
          0 Na(mol/kgw): 0.000959091
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 0.000240909
          3 Ca(mol/kgw): 8.88934e-17
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 8
     Calculated Density: 0.998301
     Calculated Volume:  0.199761
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: 7.2234e-13
          3 Ca: 1.77958e-23
          4 Cl: 0.00119643
          5 K: 0.000198915
          6 N: 0
          7 Na: 0.000997514
     Selected output: 
          0 Na(mol/kgw): 0.00100049
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 0.000199509
          3 Ca(mol/kgw): 1.78489e-23
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 9
     Calculated Density: 0.998301
     Calculated Volume:  0.199761
     Components: 
          0 H: 110.682
          1 O: 55.3416
          2 Charge: -8.142e-17
          3 Ca: 0
          4 Cl: 0.00119643
          5 K: 0.000196389
          6 N: 0
          7 Na: 0.00100004
     Selected output: 
          0 Na(mol/kgw): 0.00100303
          1 Cl(mol/kgw): 0.0012
          2 K(mol/kgw): 0.000196975
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 10
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 11
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 12
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 13
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 14
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 15
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 16
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 17
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 18
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
Cell number 19
     Calculated Density: 0.99832
     Calculated Volume:  0.199758
     Components: 
          0 H: 110.679
          1 O: 55.3438
          2 Charge: -3.39446e-16
          3 Ca: 0
          4 Cl: 0
          5 K: 0.000196302
          6 N: 0.0011964
          7 Na: 0.0010001
     Selected output: 
          0 Na(mol/kgw): 0.00100311
          1 Cl(mol/kgw): 0
          2 K(mol/kgw): 0.000196893
          3 Ca(mol/kgw): 0
          4 Temperature: 20
          5 Pressure: 2
          6 Hyd_K: 0
      pH      CaX2     KX       NaX
   7.07888  0.00011  0.00000  0.00000
   7.07888  0.00011  0.00000  0.00000
   7.07888  0.00011  0.00000  0.00000
   7.07888  0.00011  0.00000  0.00000
   7.07907  0.00011  0.00000  0.00000
   7.08166  0.00006  0.00011  0.00000
   7.08171  0.00000  0.00018  0.00004
   7.08171  0.00000  0.00012  0.00010
   7.08171  0.00000  0.00011  0.00011
   7.08171  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
   7.08162  0.00000  0.00011  0.00011
