Unknown Cantera EOS to VCSnonideal: 45012
vcs_Cantera_convert: Species Type -1 not known 
vcs_Cantera_convert: Species Type 8 not known 
vcs_Cantera_convert: Species Type 64 not known 
vcs_Cantera_convert: Species Type 8 not known 
vcs_Cantera_convert: Species Type 8 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 8 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 4 not known 
vcs_Cantera_convert: Species Type 8 not known 

================================================================================
================ Cantera_to_vprob: START OF PROBLEM STATEMENT ====================
================================================================================
             Phase IDs of species
            species     phaseID        phaseName    Initial_Estimated_kMols
          H2O(L)          0   NaCl_electrolyte             2         
             Cl-          0   NaCl_electrolyte             0         
              H+          0   NaCl_electrolyte             0         
             Na+          0   NaCl_electrolyte             0         
             OH-          0   NaCl_electrolyte             0         
              O2          1                air             0         
              H2          1                air             0         
             CO2          1                air             0         
             H2O          1                air             0         
            NaCl          1                air             0         
              N2          1                air             4         
              OH          1                air             0         
         NaCl(S)          2            NaCl(S)             5         

--------------------------------------------------------------------------------
             Information about phases
  PhaseName    PhaseNum SingSpec GasPhase EqnState NumSpec  TMolesInert       Tmoles(kmol)
NaCl_electrolyte     0     0        0 UnkType: -1             5     0.000000e+00     2.000000e+00
             air     1     0        1 Ideal Gas               7     0.000000e+00     4.000000e+00
         NaCl(S)     2     1        0 Stoich Sub              1     0.000000e+00     5.000000e+00

================================================================================
================ Cantera_to_vprob: END OF PROBLEM STATEMENT ====================
================================================================================


================================================================================
==================== Cantera_to_vprob: START OF PROBLEM STATEMENT ====================
================================================================================

             Phase IDs of species
            species     phaseID        phaseName    Initial_Estimated_kMols
          H2O(L)          0   NaCl_electrolyte             2         
             Cl-          0   NaCl_electrolyte             0         
              H+          0   NaCl_electrolyte             0         
             Na+          0   NaCl_electrolyte             0         
             OH-          0   NaCl_electrolyte             0         
              O2          1                air             0         
              H2          1                air             0         
             CO2          1                air             0         
             H2O          1                air             0         
            NaCl          1                air             0         
              N2          1                air             4         
              OH          1                air             0         
         NaCl(S)          2            NaCl(S)             5         

--------------------------------------------------------------------------------
             Information about phases
  PhaseName    PhaseNum SingSpec GasPhase EqnState NumSpec  TMolesInert       Tmoles(kmol)
NaCl_electrolyte     0     0        0 UnkType: -1             5     0.000000e+00     2.000000e+00
             air     1     0        1 Ideal Gas               7     0.000000e+00     4.000000e+00
         NaCl(S)     2     1        0 Stoich Sub              1     0.000000e+00     5.000000e+00

================================================================================
==================== Cantera_to_vprob: END OF PROBLEM STATEMENT ====================
================================================================================


================================================================================
==================== VCS_PROB: PROBLEM STATEMENT ===============================
================================================================================

	Solve a constant T, P problem:
		T    = 298.15 K
		Pres = 1 atm

             Phase IDs of species
            species     phaseID        phaseName    Initial_Estimated_Moles   Species_Type
          H2O(L)          0   NaCl_electrolyte                N/A                 Mol_Num
             Cl-          0   NaCl_electrolyte                N/A                 Mol_Num
              H+          0   NaCl_electrolyte                N/A                 Mol_Num
             Na+          0   NaCl_electrolyte                N/A                 Mol_Num
             OH-          0   NaCl_electrolyte                N/A                 Mol_Num
              O2          1                air                N/A                 Mol_Num
              H2          1                air                N/A                 Mol_Num
             CO2          1                air                N/A                 Mol_Num
             H2O          1                air                N/A                 Mol_Num
            NaCl          1                air                N/A                 Mol_Num
              N2          1                air                N/A                 Mol_Num
              OH          1                air                N/A                 Mol_Num
         NaCl(S)          2            NaCl(S)                N/A                 Mol_Num

--------------------------------------------------------------------------------
             Information about phases
  PhaseName    PhaseNum SingSpec  GasPhase    EqnState    NumSpec  TMolesInert      TKmoles
NaCl_electrolyte     0     0        0 UnkType: -1             5     0.000000e+00    N/A
             air     1     0        1 Ideal Gas               7     0.000000e+00    N/A
         NaCl(S)     2     1        0 Stoich Sub              1     0.000000e+00    N/A

Elemental Abundances:             Target_kmol    ElemType ElActive
                          O   2.000000000000E+00    0         1
                          H   4.000000000000E+00    0         1
                          C   0.000000000000E+00    0         1
                          Fe  0.000000000000E+00    0         1
                          Si  0.000000000000E+00    0         1
                          N   8.000000000000E+00    0         1
                          Na  5.000000000000E+00    0         1
                          Cl  5.000000000000E+00    0         1
                          E   0.000000000000E+00    1         0
                          cn  0.000000000000E+00    2         1
                          Ca  0.000000000000E+00    0         1

Chemical Potentials:  (J/kmol)
             Species       (phase)        SS0ChemPot       StarChemPot
          H2O(L) NaCl_electrolyte    -3.06686e+08       -3.06686e+08
             Cl-                     -1.83974e+08       -1.83974e+08
              H+                                0                  0
             Na+                     -2.57752e+08       -2.57752e+08
             OH-                     -2.26784e+08       -2.26784e+08
              O2              air     -6.1165e+07        -6.1165e+07
              H2                     -3.89624e+07       -3.89624e+07
             CO2                     -4.57249e+08       -4.57249e+08
             H2O                     -2.98124e+08       -2.98124e+08
            NaCl                     -2.49904e+08       -2.49904e+08
              N2                     -5.71282e+07       -5.71282e+07
              OH                     -2.26793e+08       -2.26793e+08
         NaCl(S)          NaCl(S)     -4.3262e+08        -4.3262e+08

================================================================================
==================== VCS_PROB: END OF PROBLEM STATEMENT ========================
================================================================================

   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to calculate the number of components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C       Fe       Si        N       Na       Cl        E cn_NaCl_       Ca
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---         Cl- |         0        0        0        0        0        0        0        1        1       -1        0
   ---          H+ |         0        1        0        0        0        0        0        0       -1        1        0
   ---         Na+ |         0        0        0        0        0        0        1        0       -1        1        0
   ---         OH- |         1        1        0        0        0        0        0        0        1       -1        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        0        0        1        1        0        0        0
   ---          N2 |         0        0        0        0        0        2        0        0        0        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0
   ---     NaCl(S) |         0        0        0        0        0        0        1        1        0        0        0

   ---   NaCl(S)     (        5) replaces H2O(L)      (        2) as component   0
   ---   N2          (        4) replaces Cl-         (        0) as component   1
   ---   H2O(L)      (        2) replaces H+          (        0) as component   2
   ---   H+          (        0) replaces Na+         (        0) as component   3
   ---   Na+         (        0) replaces OH-         (        0) as component   4
   ---   OH          (        0) replaces O2          (        0) as component   5
   ---   CO2         (        0) replaces H2          (        0) as component   6
   ---   Total number of components found =   7 (ne = 11)
    -----------------------------------------------------------------------------
   --- Subroutine elem_rearrange() called to check stoich. coefficent matrix
   ---    and to rearrange the element ordering once
   ---   N (        8) replaces Fe(        0) as element   3
   ---   Na(        5) replaces Si(        0) as element   4
   ---   Cl(        5) replaces Fe(        0) as element   5
   ---   cn(        0) replaces Si(        0) as element   6
   --- vcs_inest: Go find an initial estimate for the equilibrium problem
   --- call setInitialMoles
 --- Summary of current progress:
 ---                   Name           Moles  -       SSGibbs 
 -------------------------------------------------------------------------------------
 ---                   NaCl(S)            5  -       -174.5
 ---                        N2            4  -       -23.05
 ---                    H2O(L)            2  -       -123.7
 ---                        H+            0  -            0
 ---                       Na+            0  -         -104
 ---                        OH            0  -       -91.49
 ---                       CO2            0  -       -184.5
 ---                        H2            0  -       -15.72
 ---                       H2O            0  -       -120.3
 ---                      NaCl            0  -       -100.8
 ---                       Cl-            0  -       -74.21
 ---                        O2            0  -       -24.67
 ---                       OH-            0  -       -91.48
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---  Total sum to be minimized = -1212.2
 --- seMolesLinProg  Mole numbers failing element abundances
 --- seMolesLinProg Call vcs_elcorr to attempt fix
   --- vcsc_elcorr: Element abundances correction routine (m_numComponents != m_numElemConstraints)
   ---    Elem_Abund:  Correct             Initial                Final
   ---       O    2.000000000000E+00   0.000000000000E+00   2.000000000000E+00
   ---       H    4.000000000000E+00   0.000000000000E+00   4.000000000000E+00
   ---       C    0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       N    8.000000000000E+00   0.000000000000E+00   8.000000000000E+00
   ---       Na   5.000000000000E+00   0.000000000000E+00   5.000000000000E+00
   ---       Cl   5.000000000000E+00   0.000000000000E+00   5.000000000000E+00
   ---       cn   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Fe   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       E    0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Si   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Ca   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---            Diff_Norm:           3.490246149173E+00   0.000000000000E+00
   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to reevaluate the components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C        N       Na       Cl cn_NaCl_       Fe        E       Si       Ca
   ---     NaCl(S) |         0        0        0        0        1        1        0        0        0        0        0
   ---          N2 |         0        0        0        2        0        0        0        0        0        0        0
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---          H+ |         0        1        0        0        0        0        1        0       -1        0        0
   ---         Na+ |         0        0        0        0        1        0        1        0       -1        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        1        1        0        0        0        0        0
   ---         Cl- |         0        0        0        0        0        1       -1        0        1        0        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---         OH- |         1        1        0        0        0        0       -1        0        1        0        0

   ---   NaCl(S)     (        5) remains                          as component   0
   ---   N2          (        4) remains                          as component   1
   ---   H2O(L)      (        2) remains                          as component   2
   ---   O2          (        0) replaces H+          (        0) as component   3
   ---   H+          (        0) replaces Na+         (        0) as component   4
   ---   Na+         (        0) replaces OH          (        0) as component   5
   ---   CO2         (        0) remains                          as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:         5         4         2         0         0         0         0
   ---   NonComponent|   Moles  |       NaCl(S)   N2        H2O(L)    O2        H+        Na+       CO2       
   ---   7 H2        |         0|      -0.00      -0.00      -1.00       0.50      -0.00       0.00       0.00
   ---   8 H2O       |         0|      -0.00      -0.00      -1.00       0.00      -0.00       0.00       0.00
   ---   9 NaCl      |         0|      -1.00      -0.00       0.00       0.00      -0.00       0.00       0.00
   ---  10 Cl-       |         0|      -1.00      -0.00       0.00       0.00      -0.00       1.00       0.00
   ---  11 OH        |         0|      -0.00      -0.00      -0.50      -0.25      -0.00       0.00       0.00
   ---  12 OH-       |         0|      -0.00      -0.00      -1.00       0.00       1.00       0.00       0.00
   -----------------------------------------------------------------------------
iteration 0
 --- Summary of current progress:
 ---                   Name           Moles  -       SSGibbs 
 -------------------------------------------------------------------------------------
 ---                   NaCl(S)            0  -       -174.5
 ---                        N2            4  -       -23.05
 ---                    H2O(L)            2  -       -123.7
 ---                        O2            0  -       -24.67
 ---                        H+            0  -            0
 ---                       Na+            5  -         -104
 ---                       CO2            0  -       -184.5
 ---                        H2            0  -       -15.72
 ---                       H2O            0  -       -120.3
 ---                      NaCl            0  -       -100.8
 ---                       Cl-            5  -       -74.21
 ---                        OH            0  -       -91.49
 ---                       OH-            0  -       -91.48
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---  Total sum to be minimized = -1230.56
   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to reevaluate the components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C        N       Na       Cl cn_NaCl_       Fe        E       Si       Ca
   ---     NaCl(S) |         0        0        0        0        1        1        0        0        0        0        0
   ---          N2 |         0        0        0        2        0        0        0        0        0        0        0
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---          H+ |         0        1        0        0        0        0        1        0       -1        0        0
   ---         Na+ |         0        0        0        0        1        0        1        0       -1        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        1        1        0        0        0        0        0
   ---         Cl- |         0        0        0        0        0        1       -1        0        1        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0
   ---         OH- |         1        1        0        0        0        0       -1        0        1        0        0

   ---   Na+         (        5) replaces NaCl(S)     (        0) as component   0
   ---   Cl-         (        5) replaces N2          (        4) as component   1
   ---   N2          (        4) replaces H2O(L)      (        2) as component   2
   ---   H2O(L)      (        2) replaces O2          (        0) as component   3
   ---   OH          (        0) replaces H+          (        0) as component   4
   ---   H+          (        0) replaces NaCl(S)     (        0) as component   5
   ---   CO2         (        0) remains                          as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:         5         5         4         2         0         0         0
   ---   NonComponent|   Moles  |       Na+       Cl-       N2        H2O(L)    OH        H+        CO2       
   ---   7 H2        |         0|      -0.00      -0.00      -0.00      -2.00       2.00      -0.00      -0.00
   ---   8 H2O       |         0|      -0.00      -0.00      -0.00      -1.00      -0.00      -0.00      -0.00
   ---   9 NaCl      |         0|      -1.00      -1.00      -0.00       0.00      -0.00      -0.00      -0.00
   ---  10 O2        |         0|      -0.00      -0.00      -0.00       2.00      -4.00      -0.00      -0.00
   ---  11 NaCl(S)   |         0|      -1.00      -1.00      -0.00       0.00      -0.00      -0.00      -0.00
   ---  12 OH-       |         0|      -0.00      -0.00      -0.00      -1.00      -0.00       1.00      -0.00
   -----------------------------------------------------------------------------
iteration 1
 --- Summary of current progress:
 ---                   Name           Moles  -       SSGibbs 
 -------------------------------------------------------------------------------------
 ---                       Na+            5  -         -104
 ---                       Cl-            5  -       -74.21
 ---                        N2            4  -       -23.05
 ---                    H2O(L)            2  -       -123.7
 ---                        OH            0  -       -91.49
 ---                        H+            0  -            0
 ---                       CO2            0  -       -184.5
 ---                        H2            0  -       -15.72
 ---                       H2O            0  -       -120.3
 ---                      NaCl            0  -       -100.8
 ---                        O2            0  -       -24.67
 ---                   NaCl(S)            0  -       -174.5
 ---                       OH-            0  -       -91.48
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---                                      0  -            0
 ---  Total sum to be minimized = -1230.56
   --- vcs_inest:  Mole Numbers returned from linear programming (vcs_inest initial guess):
   --- vcs_inest:      SPECIES          MOLE_NUMBER      -SS_ChemPotential
   --- vcs_inest:      Na+                        5           104
   --- vcs_inest:      Cl-                        5          74.2
   --- vcs_inest:      N2                         4            23
   --- vcs_inest:      H2O(L)                     2           124
   --- vcs_inest:      OH                         0          91.5
   --- vcs_inest:      H+                         0            -0
   --- vcs_inest:      CO2                        0           184
   --- vcs_inest:      H2                         0          15.7
   --- vcs_inest:      H2O                        0           120
   --- vcs_inest:      NaCl                       0           101
   --- vcs_inest:      O2                         0          24.7
   --- vcs_inest:      NaCl(S)                    0           175
   --- vcs_inest:      OH-                        0          91.5
   --- vcs_inest:  Element Abundance Agreement returned from linear programming (vcs_inest initial guess):
   --- vcs_inest:      Element           Goal         Actual
   --- vcs_inest:         O                    2            2
   --- vcs_inest:         H                    4            4
   --- vcs_inest:         C                    0            0
   --- vcs_inest:         N                    8            8
   --- vcs_inest:         Na                   5            5
   --- vcs_inest:         Cl                   5            5
   --- vcs_inest:         cn_NaCl_e            0            0
   --- vcs_inest:         Fe                   0            0
   --- vcs_inest:         Si                   0            0
   --- vcs_inest:         Ca                   0            0

   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to reevaluate the components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C        N       Na       Cl cn_NaCl_       Fe        E       Si       Ca
   ---         Na+ |         0        0        0        0        1        0        1        0       -1        0        0
   ---         Cl- |         0        0        0        0        0        1       -1        0        1        0        0
   ---          N2 |         0        0        0        2        0        0        0        0        0        0        0
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0
   ---          H+ |         0        1        0        0        0        0        1        0       -1        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        1        1        0        0        0        0        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---     NaCl(S) |         0        0        0        0        1        1        0        0        0        0        0
   ---         OH- |         1        1        0        0        0        0       -1        0        1        0        0

   ---   Na+         (        5) remains                          as component   0
   ---   Cl-         (        5) remains                          as component   1
   ---   N2          (        4) remains                          as component   2
   ---   H2O(L)      (        2) remains                          as component   3
   ---   OH-         (    1e-30) replaces OH          (    1e-30) as component   4
   ---   CO2         (    1e-30) replaces H+          (    1e-30) as component   5
   ---   H2          (    1e-30) replaces H+          (    1e-30) as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:         5         5         4         2     1e-30     1e-30     1e-30
   ---   NonComponent|   Moles  |       Na+       Cl-       N2        H2O(L)    OH-       CO2       H2        
   ---   7 H+        |     1e-30|      -0.00      -0.00      -0.00      -1.00       1.00       0.00       0.00
   ---   8 H2O       |     1e-30|      -0.00      -0.00      -0.00      -1.00      -0.00       0.00       0.00
   ---   9 NaCl      |     1e-30|      -1.00      -1.00      -0.00       0.00      -0.00       0.00       0.00
   ---  10 O2        |     1e-30|      -0.00      -0.00      -0.00      -2.00      -0.00       0.00       2.00
   ---  11 NaCl(S)   |         0|      -1.00      -1.00      -0.00       0.00      -0.00       0.00       0.00
   ---  12 OH        |     1e-30|      -0.00      -0.00      -0.00      -1.00      -0.00       0.00       0.50
   -----------------------------------------------------------------------------
   --- Subroutine vcs_deltag called for all noncomponents
   --- vcs_inest: Na+         fe* =         -104.85 ff =         -103.98
   --- vcs_inest: Cl-         fe* =          -75.09 ff =         -74.214
   --- vcs_inest: N2          fe* =         -23.045 ff =         -23.045
   --- vcs_inest: H2O(L)      fe* =         -125.51 ff =         -123.72
   --- vcs_inest: OH-         fe* =         -163.05 ff =         -91.483
   --- vcs_inest: CO2         fe* =         -254.92 ff =         -184.45
   --- vcs_inest: H2          fe* =         -86.181 ff =         -15.717
   --- vcs_inest: H+          fe* =               0 ff =               0 dg* =         -37.539
   --- vcs_inest: H2O         fe* =         -120.26 ff =         -120.26 dg* =          5.2455
   --- vcs_inest: NaCl        fe* =         -100.81 ff =         -100.81 dg* =          79.131
   --- vcs_inest: O2          fe* =         -24.674 ff =         -24.674 dg* =          53.979
   --- vcs_inest: NaCl(S)     fe* =         -174.52 ff =         -174.52 dg* =          5.4241
   --- vcs_inest: OH          fe* =         -91.487 ff =         -91.487 dg* =         -9.0703
   --- vcs_inest: direction (Na+         ) = -8.28614e-26
   --- vcs_inest: direction (Cl-         ) = -8.28614e-26
   --- vcs_inest: direction (N2          ) = 0
   --- vcs_inest: direction (H2O(L)      ) = -2.81169e+17
   --- vcs_inest: direction (OH-         ) = 2.81169e+17
   --- vcs_inest: direction (CO2         ) = 0
   --- vcs_inest: direction (H2          ) = 43466.7
   --- vcs_inest: direction (H+          ) = 2.81169e+17
   --- vcs_inest: direction (H2O         ) = 0.0527108
   --- vcs_inest: direction (NaCl        ) = 8.28614e-26
   --- vcs_inest: direction (O2          ) = 3.60859e-23
   --- vcs_inest: direction (NaCl(S)     ) = 0 (ssPhase doesn't exist -> stability not checked)
   --- vcs_inest: direction (OH          ) = 86933.4
   --- Subroutine vcs_dfe called for all species
   --- Subroutine vcs_dfe called for all species
   --- vcs_inest:      Final Mole Numbers produced by inest:
   --- vcs_inest:      SPECIES      MOLE_NUMBER
   --- vcs_inest:      Na+          5
   --- vcs_inest:      Cl-          5
   --- vcs_inest:      N2           4
   --- vcs_inest:      H2O(L)       1.87965
   --- vcs_inest:      OH-          0.120348
   --- vcs_inest:      CO2          1e-30
   --- vcs_inest:      H2           1.8605e-14
   --- vcs_inest:      H+           0.120348
   --- vcs_inest:      H2O          2.25617e-20
   --- vcs_inest:      NaCl         3.5467e-44
   --- vcs_inest:      O2           1.54458e-41
   --- vcs_inest:      NaCl(S)      0
   --- vcs_inest:      OH           3.721e-14
   --- vcs_inest: Initial guess failed element abundances
   --- vcs_inest: Call vcs_elcorr to attempt fix
   --- vcsc_elcorr: Element abundances correction routine (m_numComponents != m_numElemConstraints)
   ---    Elem_Abund:  Correct             Initial                Final
   ---       O    2.000000000000E+00   2.000000000000E+00   2.000000000000E+00
   ---       H    4.000000000000E+00   4.000000000000E+00   4.000000000000E+00
   ---       C    0.000000000000E+00   1.000000000000E-30   0.000000000000E+00
   ---       N    8.000000000000E+00   8.000000000000E+00   8.000000000000E+00
   ---       Na   5.000000000000E+00   5.000000000000E+00   5.000000000000E+00
   ---       Cl   5.000000000000E+00   5.000000000000E+00   5.000000000000E+00
   ---       cn   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Fe   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       E    0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Si   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       Ca   0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---            Diff_Norm:           3.015113445778E-31   0.000000000000E+00
   --- vcs_inest: Initial guess now satisfies element abundances
   --- vcs_inest: Total Dimensionless Gibbs Free Energy =  -1.2474176E+03
VCS CALCULATION METHOD

 MultiPhase Object


   13 SPECIES
   11 ELEMENTS
    7 COMPONENTS
    3 PHASES
 PRESSURE                101325 Pa 
 TEMPERATURE            298.150 K
 PHASE1 INERTS            0.000
 PHASE2 INERTS            0.000

 ELEMENTAL ABUNDANCES             CORRECT          FROM ESTIMATE           Type

                          O   2.000000000000E+00  2.000000000000E+00       0
                          H   4.000000000000E+00  4.000000000000E+00       0
                          C   0.000000000000E+00  0.000000000000E+00       0
                          N   8.000000000000E+00  8.000000000000E+00       0
                          Na  5.000000000000E+00  5.000000000000E+00       0
                          Cl  5.000000000000E+00  5.000000000000E+00       0
                          cn  0.000000000000E+00  0.000000000000E+00       2
                          Fe  0.000000000000E+00  0.000000000000E+00       0
                          E   0.000000000000E+00  0.000000000000E+00       1
                          Si  0.000000000000E+00  0.000000000000E+00       0
                          Ca  0.000000000000E+00  0.000000000000E+00       0

 MODIFIED LINEAR PROGRAMMING ESTIMATE OF EQUILIBRIUM - forced
 Stan. Chem. Pot. in J/kmol

 SPECIES       FORMULA VECTOR                                STAN_CHEM_POT   EQUILIBRIUM_EST.  Species_Type

               O  H  C  N  Na Cl cn Fe E  Si Ca SI(I)
 Na+           0  0  0  0  1  0  1  0 -1  0  0  0              -2.57752E+08   5.00000E+00       Mol_Num
 Cl-           0  0  0  0  0  1 -1  0  1  0  0  0              -1.83974E+08   5.00000E+00       Mol_Num
 N2            0  0  0  2  0  0  0  0  0  0  0  1              -5.71282E+07   4.00000E+00       Mol_Num
 H2O(L)        1  2  0  0  0  0  0  0 -0  0  0  0              -3.06686E+08   1.87965E+00       Mol_Num
 OH-           1  1  0  0  0  0 -1  0  1  0  0  0              -2.26784E+08   1.20348E-01       Mol_Num
 CO2           2  0  1  0  0  0  0  0  0  0  0  1              -4.57249E+08   0.00000E+00       Mol_Num
 H2            0  2  0  0  0  0  0  0  0  0  0  1              -3.89624E+07   1.86050E-14       Mol_Num
 H+            0  1  0  0  0  0  1  0 -1  0  0  0               0.00000E+00   1.20348E-01       Mol_Num
 H2O           1  2  0  0  0  0  0  0  0  0  0  1              -2.98124E+08   2.25617E-20       Mol_Num
 NaCl          0  0  0  0  1  1  0  0  0  0  0  1              -2.49904E+08   3.54670E-44       Mol_Num
 O2            2  0  0  0  0  0  0  0  0  0  0  1              -6.11650E+07   1.54458E-41       Mol_Num
 NaCl(S)       0  0  0  0  1  1  0  0  0  0  0  2              -4.32620E+08   0.00000E+00       Mol_Num
 OH            1  1  0  0  0  0  0  0  0  0  0  1              -2.26793E+08   3.72100E-14       Mol_Num
   --- Subroutine vcs_dfe called for all species
   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to reevaluate the components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C        N       Na       Cl cn_NaCl_       Fe        E       Si       Ca
   ---         Na+ |         0        0        0        0        1        0        1        0       -1        0        0
   ---         Cl- |         0        0        0        0        0        1       -1        0        1        0        0
   ---          N2 |         0        0        0        2        0        0        0        0        0        0        0
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---         OH- |         1        1        0        0        0        0       -1        0        1        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0
   ---          H+ |         0        1        0        0        0        0        1        0       -1        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        1        1        0        0        0        0        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---     NaCl(S) |         0        0        0        0        1        1        0        0        0        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0

   ---   Na+         (        5) remains                          as component   0
   ---   Cl-         (        5) remains                          as component   1
   ---   N2          (        4) remains                          as component   2
   ---   H2O(L)      (      1.9) remains                          as component   3
   ---   OH-         (     0.12) remains                          as component   4
   ---   OH          (  3.7e-14) replaces CO2         (        0) as component   5
   ---   CO2         (        0) replaces H2          (  1.9e-14) as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:         5         5         4      1.88      0.12  3.72e-14         0
   ---   NonComponent|   Moles  |       Na+       Cl-       N2        H2O(L)    OH-       OH        CO2       
   ---   7 H+        |      0.12|      -0.00      -0.00      -0.00      -1.00       1.00      -0.00      -0.00
   ---   8 H2O       |  2.26e-20|      -0.00      -0.00      -0.00      -1.00      -0.00      -0.00      -0.00
   ---   9 NaCl      |  3.55e-44|      -1.00      -1.00      -0.00       0.00      -0.00      -0.00      -0.00
   ---  10 O2        |  1.54e-41|      -0.00      -0.00      -0.00       2.00      -0.00      -4.00      -0.00
   ---  11 NaCl(S)   |         0|      -1.00      -1.00      -0.00       0.00      -0.00      -0.00      -0.00
   ---  12 H2        |  1.86e-14|      -0.00      -0.00      -0.00      -2.00      -0.00       2.00      -0.00
   -----------------------------------------------------------------------------
   --- Subroutine vcs_deltag called for all noncomponents
  --- MAJOR-MINOR decision is reavaluated: All species are minor except for:
  ---      Major Species          : H+
  ---      Major Species          : NaCl(S)
  ---      Major Species          : H2
  ---
   --- Element Abundance check passed

==============================================================================================================
 Iteration =   0, Iterations since last evaluation of optimal basis =   0 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-01  -1.2035E-01   1.2434E+02 | Delta damped from -21.851 to -0.120348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-20   0.0000E+00   3.1714E+01 | Skipped: IC =   0 and DG >0:   3.171E+01
   --- NaCl            3.5467E-44   3.4897E-42  -9.8394E+01 | Normal Calc
   --- O2              1.5446E-41   3.4057E-40  -2.2050E+01 | Normal Calc
   --- NaCl(S)         0.0000E+00   5.0000E+00  -7.0667E+01 | Delta damped from 176.667 to 5 due to component 0 (       Na+) going neg
   --- H2              1.8605E-14  -1.8605E-14   1.0089E+02 | Delta damped from -6.2569e-13 to -1.8605e-14 due to component 5 (        OH) going neg
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.2035E-01 1.2035E-02-1.0831E-01 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 2.2562E-20 2.2562E-30-2.2562E-20 | minor species alternative calc
   --- NaCl          0 3.5467E-44 4.0000E-04 4.0000E-04 | minor species alternative calc
   --- O2            0 1.5446E-41 5.8190E-32 5.8190E-32 | minor species alternative calc
   --- NaCl(S)       1 0.0000E+00 5.0000E+00 5.0000E+00 | Normal Major Calc
   --- H2            1 1.8605E-14 1.8605E-15-1.6744E-14 | initial nonpos kmoles=   0.000E+00
   --- Na+           c 5.0000E+00-4.0000E-04-5.0004E+00 |
   --- Cl-           c 5.0000E+00-4.0000E-04-5.0004E+00 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.8797E+00 1.9880E+00 1.0831E-01 |
   --- OH-           c 1.2035E-01 1.2035E-02-1.0831E-01 |
   --- OH            c 3.7210E-14 3.7210E-15-3.3489E-14 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Reduction in step size due to component Na+ going negative =   9.899E-01
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -9.9366705866970E+02
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2120969437601E+03
   --- subroutine FORCE: Beginning Slope = -363.144
   --- subroutine FORCE: End Slope       = 15.486
   --- subroutine FORCE produced no adjustments (al = 0.9591)
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   Na+               5.000000E+00  5.000000E-02 -6.616232E+01 -1.040248E+02
   ---   Cl-               5.000000E+00  5.000000E-02 -3.768766E+01 -7.425510E+01
   ---   N2                4.000000E+00  4.000000E+00 -2.304522E+01 -2.304532E+01
   ---   H2O(L)            1.879652E+00  1.986873E+00 -1.986003E+02 -1.237784E+02
   ---   OH-               1.203483E-01  1.312655E-02 -9.196202E+01 -9.312026E+01
   ---   OH                3.720999E-14  4.058542E-15 -1.237955E+02 -1.260114E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082000E+02 -5.082001E+02
   ---   H+            1   1.203483E-01  1.312655E-02  1.770351E+01 -1.122536E+00  1.243418E+02  2.953561E+01
   ---   H2O           0   2.256174E-20  2.274041E-22 -1.668861E+02 -1.714835E+02  3.171423E+01 -4.770505E+01
   ---   NaCl          0   3.546703E-44  3.959683E-04 -2.022440E+02 -1.100305E+02 -9.839401E+01  6.824935E+01
   ---   O2            0   1.544580E-41  5.760392E-32 -1.200312E+02 -9.799178E+01 -2.204965E+01  1.584970E+02
   ---   NaCl(S)       1   0.000000E+00  4.949604E+00 -1.745168E+02 -1.745168E+02 -7.066681E+01  3.763099E+00
   ---   H2            1   1.860500E-14  2.029271E-15 -4.871886E+01 -5.093472E+01  1.008907E+02 -5.540070E+01
   ---                                                        Norms of Delta G():  8.347351E+01  7.747127E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.1131265E+00
   ---                  air =   4.0003960E+00
   ---              NaCl(S) =   4.9496040E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -9.9366705866970E+02
   --- Total New Dimensionless Gibbs Free Energy = -1.2120969437601E+03
   --- Increment counter increased, step is accepted:    1
   --- Normal element abundance check - passed
   --- Get a new basis because NaCl(S) is better than comp Na+ and share nonzero stoic: -1.0     
   -----------------------------------------------------------------------------
   --- Subroutine BASOPT called to reevaluate the components

   ---       Formula Matrix used in BASOPT calculation
   ---      Active |     1    1    1    1    1    1    1    1    0    1    1  
   ---     Species |         O        H        C        N       Na       Cl cn_NaCl_       Fe        E       Si       Ca
   ---         Na+ |         0        0        0        0        1        0        1        0       -1        0        0
   ---         Cl- |         0        0        0        0        0        1       -1        0        1        0        0
   ---          N2 |         0        0        0        2        0        0        0        0        0        0        0
   ---      H2O(L) |         1        2        0        0        0        0        0        0       -0        0        0
   ---         OH- |         1        1        0        0        0        0       -1        0        1        0        0
   ---          OH |         1        1        0        0        0        0        0        0        0        0        0
   ---         CO2 |         2        0        1        0        0        0        0        0        0        0        0
   ---          H+ |         0        1        0        0        0        0        1        0       -1        0        0
   ---         H2O |         1        2        0        0        0        0        0        0        0        0        0
   ---        NaCl |         0        0        0        0        1        1        0        0        0        0        0
   ---          O2 |         2        0        0        0        0        0        0        0        0        0        0
   ---     NaCl(S) |         0        0        0        0        1        1        0        0        0        0        0
   ---          H2 |         0        2        0        0        0        0        0        0        0        0        0

   ---   NaCl(S)     (      4.9) replaces Na+         (     0.05) as component   0
   ---   N2          (        4) replaces Cl-         (     0.05) as component   1
   ---   H2O(L)      (        2) replaces Cl-         (     0.05) as component   2
   ---   Cl-         (     0.05) remains                          as component   3
   ---   OH-         (    0.013) remains                          as component   4
   ---   OH          (  4.1e-15) remains                          as component   5
   ---   CO2         (        0) remains                          as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:      4.95         4      1.99      0.05    0.0131  4.06e-15         0
   ---   NonComponent|   Moles  |       NaCl(S)   N2        H2O(L)    Cl-       OH-       OH        CO2       
   ---   7 H+        |    0.0131|      -0.00      -0.00      -1.00      -0.00       1.00      -0.00       0.00
   ---   8 H2O       |  2.27e-22|      -0.00      -0.00      -1.00      -0.00       0.00      -0.00       0.00
   ---   9 NaCl      |  0.000396|      -1.00      -0.00       0.00      -0.00       0.00      -0.00       0.00
   ---  10 O2        |  5.76e-32|      -0.00      -0.00       2.00      -0.00       0.00      -4.00       0.00
   ---  11 Na+       |      0.05|      -1.00      -0.00       0.00       1.00       0.00      -0.00       0.00
   ---  12 H2        |  2.03e-15|      -0.00      -0.00      -2.00      -0.00       0.00       2.00       0.00
   -----------------------------------------------------------------------------
   --- Subroutine vcs_deltag called for all noncomponents
  --- MAJOR-MINOR decision is reavaluated: All species are minor except for:
  ---      Major Species          : H+
  ---      Major Species          : Na+
  ---      Major Species          : H2
  ---
   --- Element Abundance check passed

==============================================================================================================
 Iteration =   1, Iterations since last evaluation of optimal basis =   0 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-02  -1.3127E-02   2.9536E+01 | Delta damped from -0.206003 to -0.0131265 due to component 4 (       OH-) going neg
   --- H2O             2.2740E-22   1.0848E-20  -4.7705E+01 | Normal Calc
   --- NaCl            3.9597E-04   0.0000E+00   6.4486E+01 | Skipped: IC =   0 and DG >0:   6.449E+01
   --- O2              5.7604E-32   0.0000E+00   1.5850E+02 | Skipped: IC =   0 and DG >0:   1.585E+02
   --- Na+             5.0000E-02   9.0116E-02  -3.7631E+00 | Normal calc: diag adjusted from 38.1071 to 41.7586 due to act coeff
   --- H2              2.0293E-15   3.7474E-14  -5.5401E+01 | Normal Calc
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-02 1.3127E-03-1.1814E-02 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 2.2740E-22 4.0004E-04 4.0004E-04 | minor species alternative calc
   --- NaCl          0 3.9597E-04 3.9597E-14-3.9597E-04 | minor species alternative calc
   --- O2            0 5.7604E-32 5.7604E-42-5.7604E-32 | minor species alternative calc
   --- Na+           1 5.0000E-02 1.4012E-01 9.0116E-02 | Normal Major Calc
   --- H2            1 2.0293E-15 3.9504E-14 3.7474E-14 | Normal Major Calc
   --- NaCl(S)       c 4.9496E+00 4.8599E+00-8.9720E-02 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9869E+00 1.9983E+00 1.1414E-02 |
   --- Cl-           c 5.0000E-02 1.4012E-01 9.0116E-02 |
   --- OH-           c 1.3127E-02 1.3127E-03-1.1814E-02 |
   --- OH            c 4.0585E-15 7.9007E-14 7.4949E-14 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2120969437601E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2126835377536E+03
   --- subroutine FORCE: Beginning Slope = -0.732663
   --- subroutine FORCE: End Slope       = -0.455873
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.949604E+00  4.859884E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   H2O(L)            1.986873E+00  1.998287E+00 -1.237784E+02 -1.238730E+02
   ---   Cl-               5.000000E-02  1.401155E-01 -7.425510E+01 -7.310540E+01
   ---   OH-               1.312655E-02  1.312655E-03 -9.312026E+01 -9.542293E+01
   ---   OH                4.058542E-15  7.900722E-14 -1.260114E+02 -1.230427E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.312655E-02  1.312655E-03 -1.122536E+00 -2.548359E+00  2.953561E+01  2.590174E+01
   ---   H2O           0   2.274041E-22  4.000396E-04 -1.714835E+02 -1.294721E+02 -4.770505E+01 -5.599064E+00
   ---   NaCl          0   3.959683E-04  3.959685E-14 -1.100305E+02 -1.330564E+02  6.448625E+01  4.146040E+01
   ---   O2            0   5.760392E-32  5.760388E-42 -9.799178E+01 -1.210176E+02  1.584970E+02  1.234070E+02
   ---   Na+           1   5.000000E-02  1.401155E-01 -1.040248E+02 -1.028675E+02 -3.763099E+00 -1.456082E+00
   ---   H2            1   2.029271E-15  3.950361E-14 -5.093472E+01 -4.796600E+01 -5.540070E+01 -4.630530E+01
   ---                                                        Norms of Delta G():  7.693210E+01  5.744100E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.2811437E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.8598845E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2120969437601E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2126835377536E+03
   --- Increment counter increased, step is accepted:    2
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   2, Iterations since last evaluation of optimal basis =   1 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-03  -1.3127E-03   2.5902E+01 | Delta damped from -0.0170301 to -0.00131265 due to component 4 (       OH-) going neg
   --- H2O             4.0004E-04   2.2400E-03  -5.5991E+00 | Normal calc: diag adjusted from 2499.56 to 2499.6 due to act coeff
   --- NaCl            3.9597E-14   0.0000E+00   4.1460E+01 | Skipped: IC =   0 and DG >0:   4.146E+01
   --- O2              5.7604E-42   0.0000E+00   1.2341E+02 | Skipped: IC =   0 and DG >0:   1.234E+02
   --- Na+             1.4012E-01   7.3291E-02  -1.4561E+00 | Normal calc: diag adjusted from 12.5204 to 19.867 due to act coeff
   --- H2              3.9504E-14   6.0974E-13  -4.6305E+01 | Normal calc: diag adjusted from 7.59424e+13 to 7.59424e+13 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-03 1.3127E-04-1.1814E-03 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 4.0004E-04 4.0004E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 3.9597E-14 3.9597E-14 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-42 5.7604E-42 0.0000E+00 | minor species not considered
   --- Na+           1 1.4012E-01 2.1341E-01 7.3291E-02 | Normal Major Calc
   --- H2            1 3.9504E-14 6.4925E-13 6.0974E-13 | Normal Major Calc
   --- NaCl(S)       c 4.8599E+00 4.7866E+00-7.3291E-02 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9983E+00 1.9995E+00 1.1814E-03 |
   --- Cl-           c 1.4012E-01 2.1341E-01 7.3291E-02 |
   --- OH-           c 1.3127E-03 1.3127E-04-1.1814E-03 |
   --- OH            c 7.9007E-14 1.2985E-12 1.2195E-12 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2126835377536E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127705059142E+03
   --- subroutine FORCE: Beginning Slope = -0.137318
   --- subroutine FORCE: End Slope       = -0.037668
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.859884E+00  4.786593E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   H2O(L)            1.998287E+00  1.999469E+00 -1.238730E+02 -1.239862E+02
   ---   Cl-               1.401155E-01  2.134070E-01 -7.310540E+01 -7.245606E+01
   ---   OH-               1.312655E-03  1.312655E-04 -9.542293E+01 -9.770873E+01
   ---   OH                7.900722E-14  1.298492E-12 -1.230427E+02 -1.202433E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.312655E-03  1.312655E-04 -2.548359E+00 -4.143415E+00  2.590174E+01  2.213408E+01
   ---   H2O           0   4.000396E-04  4.000396E-04 -1.294721E+02 -1.294721E+02 -5.599064E+00 -5.485865E+00
   ---   NaCl          0   3.959685E-14  3.959685E-14 -1.330564E+02 -1.330564E+02  4.146040E+01  4.146040E+01
   ---   O2            0   5.760388E-42  5.760388E-42 -1.210176E+02 -1.210176E+02  1.234070E+02  1.119829E+02
   ---   Na+           1   1.401155E-01  2.134070E-01 -1.028675E+02 -1.022179E+02 -1.456082E+00 -1.571674E-01
   ---   H2            1   3.950361E-14  6.492458E-13 -4.796600E+01 -4.516658E+01 -4.630530E+01 -3.768065E+01
   ---                                                        Norms of Delta G():  5.744100E+01  5.195994E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4265453E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7865930E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2126835377536E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127705059142E+03
   --- Increment counter increased, step is accepted:    3
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   3, Iterations since last evaluation of optimal basis =   2 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-04  -1.3127E-04   2.2134E+01 | Delta damped from -0.00145263 to -0.000131265 due to component 4 (       OH-) going neg
   --- H2O             4.0004E-04   2.1946E-03  -5.4859E+00 | Normal calc: diag adjusted from 2499.59 to 2499.69 due to act coeff
   --- NaCl            3.9597E-14   0.0000E+00   4.1460E+01 | Skipped: IC =   0 and DG >0:   4.146E+01
   --- O2              5.7604E-42   0.0000E+00   1.1198E+02 | Skipped: IC =   0 and DG >0:   1.120E+02
   --- Na+             2.1341E-01   9.5232E-03  -1.5717E-01 | Normal calc: diag adjusted from 7.72333 to 16.5036 due to act coeff
   --- H2              6.4925E-13   8.1547E-12  -3.7681E+01 | Normal calc: diag adjusted from 4.62075e+12 to 4.62075e+12 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-04 1.3127E-05-1.1814E-04 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 4.0004E-04 4.0004E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 3.9597E-14 3.9597E-14 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-42 5.7604E-42 0.0000E+00 | minor species not considered
   --- Na+           1 2.1341E-01 2.2293E-01 9.5232E-03 | Normal Major Calc
   --- H2            1 6.4925E-13 8.8039E-12 8.1547E-12 | Normal Major Calc
   --- NaCl(S)       c 4.7866E+00 4.7771E+00-9.5232E-03 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9995E+00 1.9996E+00 1.1814E-04 |
   --- Cl-           c 2.1341E-01 2.2293E-01 9.5232E-03 |
   --- OH-           c 1.3127E-04 1.3127E-05-1.1814E-04 |
   --- OH            c 1.2985E-12 1.7608E-11 1.6309E-11 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2127705059142E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127737105506E+03
   --- subroutine FORCE: Beginning Slope = -0.00411164
   --- subroutine FORCE: End Slope       = -0.00210899
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.786593E+00  4.777070E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   H2O(L)            1.999469E+00  1.999587E+00 -1.239862E+02 -1.240030E+02
   ---   Cl-               2.134070E-01  2.229302E-01 -7.245606E+01 -7.237877E+01
   ---   OH-               1.312655E-04  1.312655E-05 -9.770873E+01 -1.000073E+02
   ---   OH                1.298492E-12  1.760783E-11 -1.202433E+02 -1.176361E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.312655E-04  1.312655E-05 -4.143415E+00 -6.353376E+00  2.213408E+01  1.764226E+01
   ---   H2O           0   4.000396E-04  4.000396E-04 -1.294721E+02 -1.294721E+02 -5.485865E+00 -5.469115E+00
   ---   NaCl          0   3.959685E-14  3.959685E-14 -1.330564E+02 -1.330564E+02  4.146040E+01  4.146040E+01
   ---   O2            0   5.760388E-42  5.760388E-42 -1.210176E+02 -1.210176E+02  1.119829E+02  1.015209E+02
   ---   Na+           1   2.134070E-01  2.229302E-01 -1.022179E+02 -1.021406E+02 -1.571674E-01 -2.598987E-03
   ---   H2            1   6.492458E-13  8.803914E-12 -4.516658E+01 -4.255944E+01 -3.768065E+01 -2.982573E+01
   ---                                                        Norms of Delta G():  5.195994E+01  4.700388E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4454736E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7770698E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127705059142E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127737105506E+03
   --- Increment counter increased, step is accepted:    4
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents

==============================================================================================================
 Iteration =   4, Iterations since last evaluation of optimal basis =   3 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-05  -1.3127E-05   1.7642E+01 | Delta damped from -0.00011579 to -1.31265e-05 due to component 4 (       OH-) going neg
   --- H2O             4.0004E-04   2.1879E-03  -5.4691E+00 | Normal calc: diag adjusted from 2499.59 to 2499.7 due to act coeff
   --- NaCl            3.9597E-14   0.0000E+00   4.1460E+01 | Skipped: IC =   0 and DG >0:   4.146E+01
   --- O2              5.7604E-42   0.0000E+00   1.0152E+02 | Skipped: IC =   0 and DG >0:   1.015E+02
   --- Na+             2.2293E-01   1.5986E-04  -2.5990E-03 | Normal calc: diag adjusted from 7.33574 to 16.2581 due to act coeff
   --- H2              8.8039E-12   8.7528E-11  -2.9826E+01 | Normal calc: diag adjusted from 3.40758e+11 to 3.40758e+11 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-05 1.3127E-06-1.1814E-05 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 4.0004E-04 4.0004E-02 3.9604E-02 | minor species alternative calc
   --- NaCl          0 3.9597E-14 3.9597E-24-3.9597E-14 | minor species alternative calc
   --- O2            0 5.7604E-42 5.7604E-52-5.7604E-42 | minor species alternative calc
   --- Na+           1 2.2293E-01 2.2309E-01 1.5986E-04 | Normal Major Calc
   --- H2            1 8.8039E-12 9.6332E-11 8.7528E-11 | Normal Major Calc
   --- NaCl(S)       c 4.7771E+00 4.7769E+00-1.5986E-04 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9996E+00 1.9600E+00-3.9592E-02 |
   --- Cl-           c 2.2293E-01 2.2309E-01 1.5986E-04 |
   --- OH-           c 1.3127E-05 1.3127E-06-1.1814E-05 |
   --- OH            c 1.7608E-11 1.9266E-10 1.7506E-10 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2127737105506E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128459024040E+03
   --- subroutine FORCE: Beginning Slope = -0.216807
   --- subroutine FORCE: End Slope       = -0.0344129
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.777070E+00  4.776910E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.305518E+01
   ---   H2O(L)            1.999587E+00  1.959995E+00 -1.240030E+02 -1.240115E+02
   ---   Cl-               2.229302E-01  2.230901E-01 -7.237877E+01 -7.234102E+01
   ---   OH-               1.312655E-05  1.312655E-06 -1.000073E+02 -1.022879E+02
   ---   OH                1.760783E-11  1.926633E-10 -1.176361E+02 -1.152534E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082100E+02
   ---   H+            1   1.312655E-05  1.312655E-06 -6.353376E+00 -8.590552E+00  1.764226E+01  1.313303E+01
   ---   H2O           0   4.000396E-04  4.000396E-02 -1.294721E+02 -1.248768E+02 -5.469115E+00 -8.653044E-01
   ---   NaCl          0   3.959685E-14  3.959684E-24 -1.330564E+02 -1.560921E+02  4.146040E+01  1.842469E+01
   ---   O2            0   5.760388E-42  5.760391E-52 -1.210176E+02 -1.440533E+02  1.015209E+02  6.893721E+01
   ---   Na+           1   2.229302E-01  2.230901E-01 -1.021406E+02 -1.021029E+02 -2.598987E-03  7.288686E-02
   ---   H2            1   8.803914E-12  9.633163E-11 -4.255944E+01 -4.017669E+01 -2.982573E+01 -2.266050E+01
   ---                                                        Norms of Delta G():  4.700388E+01  3.103368E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4061775E+00
   ---                  air =   4.0400040E+00
   ---              NaCl(S) =   4.7769099E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127737105506E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128459024040E+03
   --- Increment counter increased, step is accepted:    5
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   ---   Noncomponent turned from minor to major: H2O
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   5, Iterations since last evaluation of optimal basis =   4 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-06  -1.3127E-06   1.3133E+01 | Delta damped from -8.61956e-06 to -1.31265e-06 due to component 4 (       OH-) going neg
   --- H2O             4.0004E-02   3.4663E-02  -8.6530E-01 | Normal calc: diag adjusted from 24.8446 to 24.9634 due to act coeff
   --- NaCl            3.9597E-24   0.0000E+00   1.8425E+01 | Skipped: IC =   0 and DG >0:   1.842E+01
   --- O2              5.7604E-52   0.0000E+00   6.8937E+01 | Skipped: IC =   0 and DG >0:   6.894E+01
   --- Na+             2.2309E-01  -4.4243E-03   7.2887E-02 | Normal calc: diag adjusted from 7.3026 to 16.4741 due to act coeff
   --- H2              9.6332E-11   7.2764E-10  -2.2660E+01 | Normal calc: diag adjusted from 3.11424e+10 to 3.11424e+10 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-06 1.3127E-07-1.1814E-06 | initial nonpos kmoles=   0.000E+00
   --- H2O           1 4.0004E-02 7.4667E-02 3.4663E-02 | Normal Major Calc
   --- NaCl          0 3.9597E-24 3.9597E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-52 5.7604E-52 0.0000E+00 | minor species not considered
   --- Na+           1 2.2309E-01 2.1867E-01-4.4243E-03 | Normal Major Calc
   --- H2            1 9.6332E-11 8.2397E-10 7.2764E-10 | Normal Major Calc
   --- NaCl(S)       c 4.7769E+00 4.7813E+00 4.4243E-03 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9600E+00 1.9253E+00-3.4662E-02 |
   --- Cl-           c 2.2309E-01 2.1867E-01-4.4243E-03 |
   --- OH-           c 1.3127E-06 1.3127E-07-1.1814E-06 |
   --- OH            c 1.9266E-10 1.6479E-09 1.4553E-09 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128459024040E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128644454499E+03
   --- subroutine FORCE: Beginning Slope = -0.0303319
   --- subroutine FORCE: End Slope       = -0.00898719
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.776910E+00  4.781334E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.305518E+01 -2.306372E+01
   ---   H2O(L)            1.959995E+00  1.925333E+00 -1.240115E+02 -1.240106E+02
   ---   Cl-               2.230901E-01  2.186658E-01 -7.234102E+01 -7.234505E+01
   ---   OH-               1.312655E-06  1.312655E-07 -1.022879E+02 -1.045728E+02
   ---   OH                1.926633E-10  1.647945E-09 -1.152534E+02 -1.131156E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082100E+02 -5.082185E+02
   ---   H+            1   1.312655E-06  1.312655E-07 -8.590552E+00 -1.088013E+01  1.313303E+01  8.557572E+00
   ---   H2O           1   4.000396E-02  7.466684E-02 -1.248768E+02 -1.242613E+02 -8.653044E-01 -2.507065E-01
   ---   NaCl          0   3.959684E-24  3.959684E-24 -1.560921E+02 -1.561006E+02  1.842469E+01  1.841615E+01
   ---   O2            0   5.760391E-52  5.760391E-52 -1.440533E+02 -1.440619E+02  6.893721E+01  6.037931E+01
   ---   Na+           1   2.230901E-01  2.186658E-01 -1.021029E+02 -1.021069E+02  7.288686E-02  6.483724E-02
   ---   H2            1   9.633163E-11  8.239725E-10 -4.017669E+01 -3.803890E+01 -2.266050E+01 -1.624894E+01
   ---                                                        Norms of Delta G():  3.103368E+01  2.683947E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3626648E+00
   ---                  air =   4.0746668E+00
   ---              NaCl(S) =   4.7813342E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128459024040E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128644454499E+03
   --- Increment counter increased, step is accepted:    6
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   6, Iterations since last evaluation of optimal basis =   5 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-07  -1.3127E-07   8.5576E+00 | Delta damped from -5.61657e-07 to -1.31265e-07 due to component 4 (       OH-) going neg
   --- H2O             7.4667E-02   1.8760E-02  -2.5071E-01 | Normal calc: diag adjusted from 13.2435 to 13.3639 due to act coeff
   --- NaCl            3.9597E-24   0.0000E+00   1.8416E+01 | Skipped: IC =   0 and DG >0:   1.842E+01
   --- O2              5.7604E-52   0.0000E+00   6.0379E+01 | Skipped: IC =   0 and DG >0:   6.038E+01
   --- Na+             2.1867E-01  -3.8634E-03   6.4837E-02 | Normal calc: diag adjusted from 7.45338 to 16.7825 due to act coeff
   --- H2              8.2397E-10   4.4629E-09  -1.6249E+01 | Normal calc: diag adjusted from 3.6409e+09 to 3.6409e+09 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-07 1.3127E-08-1.1814E-07 | initial nonpos kmoles=   0.000E+00
   --- H2O           1 7.4667E-02 9.3427E-02 1.8760E-02 | Normal Major Calc
   --- NaCl          0 3.9597E-24 3.9597E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-52 5.7604E-52 0.0000E+00 | minor species not considered
   --- Na+           1 2.1867E-01 2.1480E-01-3.8634E-03 | Normal Major Calc
   --- H2            1 8.2397E-10 5.2869E-09 4.4629E-09 | Normal Major Calc
   --- NaCl(S)       c 4.7813E+00 4.7852E+00 3.8634E-03 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9253E+00 1.9066E+00-1.8760E-02 |
   --- Cl-           c 2.1867E-01 2.1480E-01-3.8634E-03 |
   --- OH-           c 1.3127E-07 1.3127E-08-1.1814E-07 |
   --- OH            c 1.6479E-09 1.0574E-08 8.9258E-09 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128644454499E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128672368083E+03
   --- subroutine FORCE: Beginning Slope = -0.00495483
   --- subroutine FORCE: End Slope       = -0.000784316
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.781334E+00  4.785198E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306372E+01 -2.306831E+01
   ---   H2O(L)            1.925333E+00  1.906573E+00 -1.240106E+02 -1.240072E+02
   ---   Cl-               2.186658E-01  2.148024E-01 -7.234505E+01 -7.235975E+01
   ---   OH-               1.312655E-07  1.312655E-08 -1.045728E+02 -1.068664E+02
   ---   OH                1.647945E-09  1.057373E-08 -1.131156E+02 -1.112613E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082185E+02 -5.082231E+02
   ---   H+            1   1.312655E-07  1.312655E-08 -1.088013E+01 -1.319060E+01  8.557572E+00  3.950186E+00
   ---   H2O           1   7.466684E-02  9.342684E-02 -1.242613E+02 -1.240417E+02 -2.507065E-01 -3.448382E-02
   ---   NaCl          0   3.959684E-24  3.959684E-24 -1.561006E+02 -1.561052E+02  1.841615E+01  1.841156E+01
   ---   O2            0   5.760391E-52  5.760391E-52 -1.440619E+02 -1.440665E+02  6.037931E+01  5.296437E+01
   ---   Na+           1   2.186658E-01  2.148024E-01 -1.021069E+02 -1.021216E+02  6.483724E-02  3.543137E-02
   ---   H2            1   8.239725E-10  5.286865E-09 -3.803890E+01 -3.618465E+01 -1.624894E+01 -1.069284E+01
   ---                                                        Norms of Delta G():  2.683947E+01  2.336005E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3361779E+00
   ---                  air =   4.0934269E+00
   ---              NaCl(S) =   4.7851976E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128644454499E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128672368083E+03
   --- Increment counter increased, step is accepted:    7
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   7, Iterations since last evaluation of optimal basis =   6 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-08  -1.3127E-08   3.9502E+00 | Delta damped from -2.59261e-08 to -1.31265e-08 due to component 4 (       OH-) going neg
   --- H2O             9.3427E-02   3.2303E-03  -3.4484E-02 | Normal calc: diag adjusted from 10.5557 to 10.675 due to act coeff
   --- NaCl            3.9597E-24   0.0000E+00   1.8412E+01 | Skipped: IC =   0 and DG >0:   1.841E+01
   --- O2              5.7604E-52   0.0000E+00   5.2964E+01 | Skipped: IC =   0 and DG >0:   5.296E+01
   --- Na+             2.1480E-01  -2.0852E-03   3.5431E-02 | Normal calc: diag adjusted from 7.59869 to 16.9922 due to act coeff
   --- H2              5.2869E-09   1.8844E-08  -1.0693E+01 | Normal calc: diag adjusted from 5.67444e+08 to 5.67444e+08 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-08 1.3127E-09-1.1814E-08 | initial nonpos kmoles=   0.000E+00
   --- H2O           1 9.3427E-02 9.6657E-02 3.2303E-03 | Normal Major Calc
   --- NaCl          0 3.9597E-24 3.9597E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-52 5.7604E-52 0.0000E+00 | minor species not considered
   --- Na+           1 2.1480E-01 2.1272E-01-2.0852E-03 | Normal Major Calc
   --- H2            1 5.2869E-09 2.4131E-08 1.8844E-08 | Normal Major Calc
   --- NaCl(S)       c 4.7852E+00 4.7873E+00 2.0852E-03 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9066E+00 1.9033E+00-3.2304E-03 |
   --- Cl-           c 2.1480E-01 2.1272E-01-2.0852E-03 |
   --- OH-           c 1.3127E-08 1.3127E-09-1.1814E-08 |
   --- OH            c 1.0574E-08 4.8261E-08 3.7688E-08 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128672368083E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673430854E+03
   --- subroutine FORCE: Beginning Slope = -0.000185523
   --- subroutine FORCE: End Slope       = -2.75894e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.785198E+00  4.787283E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306831E+01 -2.306910E+01
   ---   H2O(L)            1.906573E+00  1.903343E+00 -1.240072E+02 -1.240039E+02
   ---   Cl-               2.148024E-01  2.127172E-01 -7.235975E+01 -7.237442E+01
   ---   OH-               1.312655E-08  1.312655E-09 -1.068664E+02 -1.091681E+02
   ---   OH                1.057373E-08  4.826147E-08 -1.112613E+02 -1.097438E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082231E+02 -5.082239E+02
   ---   H+            1   1.312655E-08  1.312655E-09 -1.319060E+01 -1.550908E+01  3.950186E+00 -6.732505E-01
   ---   H2O           1   9.342684E-02  9.665719E-02 -1.240417E+02 -1.240085E+02 -3.448382E-02 -4.572594E-03
   ---   NaCl          0   3.959684E-24  3.959684E-24 -1.561052E+02 -1.561060E+02  1.841156E+01  1.841077E+01
   ---   O2            0   5.760391E-52  5.760391E-52 -1.440665E+02 -1.440673E+02  5.296437E+01  4.690027E+01
   ---   Na+           1   2.148024E-01  2.127172E-01 -1.021216E+02 -1.021363E+02  3.543137E-02  6.095652E-03
   ---   H2            1   5.286865E-09  2.413074E-08 -3.618465E+01 -3.466717E+01 -1.069284E+01 -6.147009E+00
   ---                                                        Norms of Delta G():  2.336005E+01  2.072370E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3287772E+00
   ---                  air =   4.0966573E+00
   ---              NaCl(S) =   4.7872828E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128672368083E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673430854E+03
   --- Increment counter increased, step is accepted:    8
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents

==============================================================================================================
 Iteration =   8, Iterations since last evaluation of optimal basis =   7 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.3127E-09   4.4187E-10  -6.7325E-01 | Normal calc: diag adjusted from 1.52363e+09 to 1.52363e+09 due to act coeff
   --- H2O             9.6657E-02   4.4330E-04  -4.5726E-03 | Normal calc: diag adjusted from 10.1977 to 10.3149 due to act coeff
   --- NaCl            3.9597E-24   0.0000E+00   1.8411E+01 | Skipped: IC =   0 and DG >0:   1.841E+01
   --- O2              5.7604E-52   0.0000E+00   4.6900E+01 | Skipped: IC =   0 and DG >0:   4.690E+01
   --- Na+             2.1272E-01  -3.5717E-04   6.0957E-03 | Normal calc: diag adjusted from 7.68451 to 17.0663 due to act coeff
   --- H2              2.4131E-08   4.9444E-08  -6.1470E+00 | Normal calc: diag adjusted from 1.24323e+08 to 1.24323e+08 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.3127E-09 1.7545E-09 4.4187E-10 | Normal Major Calc
   --- H2O           1 9.6657E-02 9.7100E-02 4.4330E-04 | Normal Major Calc
   --- NaCl          0 3.9597E-24 3.9991E-32-3.9597E-24 | minor species alternative calc
   --- O2            0 5.7604E-52 5.7604E-62-5.7604E-52 | minor species alternative calc
   --- Na+           1 2.1272E-01 2.1236E-01-3.5717E-04 | Normal Major Calc
   --- H2            1 2.4131E-08 7.3575E-08 4.9444E-08 | Normal Major Calc
   --- NaCl(S)       c 4.7873E+00 4.7876E+00 3.5717E-04 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9033E+00 1.9029E+00-4.4340E-04 |
   --- Cl-           c 2.1272E-01 2.1236E-01-3.5717E-04 |
   --- OH-           c 1.3127E-09 1.7545E-09 4.4187E-10 |
   --- OH            c 4.8261E-08 1.4715E-07 9.8888E-08 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673430854E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673456972E+03
   --- subroutine FORCE: Beginning Slope = -4.50847e-06
   --- subroutine FORCE: End Slope       = -7.4645e-07
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787283E+00  4.787640E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306910E+01 -2.306921E+01
   ---   H2O(L)            1.903343E+00  1.902899E+00 -1.240039E+02 -1.240033E+02
   ---   Cl-               2.127172E-01  2.123601E-01 -7.237442E+01 -7.237704E+01
   ---   OH-               1.312655E-09  1.754527E-09 -1.091681E+02 -1.088779E+02
   ---   OH                4.826147E-08  1.471494E-07 -1.097438E+02 -1.086291E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082239E+02 -5.082240E+02
   ---   H+            1   1.312655E-09  1.754527E-09 -1.550908E+01 -1.522184E+01 -6.732505E-01 -9.635679E-02
   ---   H2O           1   9.665719E-02  9.710049E-02 -1.240085E+02 -1.240040E+02 -4.572594E-03 -6.915951E-04
   ---   NaCl          0   3.959684E-24  3.999130E-32 -1.561060E+02 -1.745169E+02  1.841077E+01 -1.082493E-04
   ---   O2            0   5.760391E-52  5.760390E-62 -1.440673E+02 -1.670932E+02  4.690027E+01  1.941666E+01
   ---   Na+           1   2.127172E-01  2.123601E-01 -1.021363E+02 -1.021389E+02  6.095652E-03  8.432280E-04
   ---   H2            1   2.413074E-08  7.357469E-08 -3.466717E+01 -3.355247E+01 -6.147009E+00 -2.804063E+00
   ---                                                        Norms of Delta G():  2.072370E+01  8.009148E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3276195E+00
   ---                  air =   4.0971007E+00
   ---              NaCl(S) =   4.7876399E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673430854E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673456972E+03
   --- Increment counter increased, step is accepted:    9
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =   9, Iterations since last evaluation of optimal basis =   8 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.7545E-09   8.4530E-11  -9.6357E-02 | Normal calc: diag adjusted from 1.13991e+09 to 1.13991e+09 due to act coeff
   --- H2O             9.7100E-02   6.7359E-05  -6.9160E-04 | Normal calc: diag adjusted from 10.1504 to 10.2672 due to act coeff
   --- NaCl            3.9991E-32   4.3290E-36  -1.0825E-04 | Normal Calc
   --- O2              5.7604E-62   0.0000E+00   1.9417E+01 | Skipped: IC =   0 and DG >0:   1.942E+01
   --- Na+             2.1236E-01  -4.9374E-05   8.4323E-04 | Normal calc: diag adjusted from 7.69947 to 17.0785 due to act coeff
   --- H2              7.3575E-08   6.8769E-08  -2.8041E+00 | Normal calc: diag adjusted from 4.07749e+07 to 4.07749e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.7545E-09 1.8391E-09 8.4530E-11 | Normal Major Calc
   --- H2O           1 9.7100E-02 9.7168E-02 6.7359E-05 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-62 5.7604E-62 0.0000E+00 | minor species not considered
   --- Na+           1 2.1236E-01 2.1231E-01-4.9374E-05 | Normal Major Calc
   --- H2            1 7.3575E-08 1.4234E-07 6.8769E-08 | Normal Major Calc
   --- NaCl(S)       c 4.7876E+00 4.7877E+00 4.9374E-05 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9029E+00 1.9028E+00-6.7497E-05 |
   --- Cl-           c 2.1236E-01 2.1231E-01-4.9374E-05 |
   --- OH-           c 1.7545E-09 1.8391E-09 8.4530E-11 |
   --- OH            c 1.4715E-07 2.8469E-07 1.3754E-07 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673456972E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458650E+03
   --- subroutine FORCE: Beginning Slope = -2.81061e-07
   --- subroutine FORCE: End Slope       = -6.94015e-08
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787640E+00  4.787689E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306921E+01 -2.306923E+01
   ---   H2O(L)            1.902899E+00  1.902832E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123601E-01  2.123107E-01 -7.237704E+01 -7.237740E+01
   ---   OH-               1.754527E-09  1.839058E-09 -1.088779E+02 -1.088308E+02
   ---   OH                1.471494E-07  2.846881E-07 -1.086291E+02 -1.079692E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.754527E-09  1.839058E-09 -1.522184E+01 -1.517518E+01 -9.635679E-02 -2.704330E-03
   ---   H2O           1   9.710049E-02  9.716785E-02 -1.240040E+02 -1.240034E+02 -6.915951E-04 -9.436152E-05
   ---   NaCl          0   3.999130E-32  3.999130E-32 -1.745169E+02 -1.745169E+02 -1.082493E-04 -1.247403E-04
   ---   O2            0   5.760390E-62  5.760390E-62 -1.670932E+02 -1.670932E+02  1.941666E+01  1.677709E+01
   ---   Na+           1   2.123601E-01  2.123107E-01 -1.021389E+02 -1.021393E+02  8.432280E-04  1.285966E-04
   ---   H2            1   7.357469E-08  1.423440E-07 -3.355247E+01 -3.289254E+01 -2.804063E+00 -8.244342E-01
   ---                                                        Norms of Delta G():  8.009148E+00  6.857481E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274532E+00
   ---                  air =   4.0971683E+00
   ---              NaCl(S) =   4.7876893E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673456972E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458650E+03
   --- Increment counter increased, step is accepted:   10
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  10, Iterations since last evaluation of optimal basis =   9 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8391E-09   2.4867E-12  -2.7043E-03 | Normal calc: diag adjusted from 1.08751e+09 to 1.08751e+09 due to act coeff
   --- H2O             9.7168E-02   9.1970E-06  -9.4362E-05 | Normal calc: diag adjusted from 10.1433 to 10.26 due to act coeff
   --- NaCl            3.9991E-32   4.9885E-36  -1.2474E-04 | Normal Calc
   --- O2              5.7604E-62   0.0000E+00   1.6777E+01 | Skipped: IC =   0 and DG >0:   1.678E+01
   --- Na+             2.1231E-01  -7.5290E-06   1.2860E-04 | Normal calc: diag adjusted from 7.70154 to 17.0802 due to act coeff
   --- H2              1.4234E-07   3.9118E-08  -8.2443E-01 | Normal calc: diag adjusted from 2.10757e+07 to 2.10757e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8391E-09 1.8415E-09 2.4867E-12 | Normal Major Calc
   --- H2O           1 9.7168E-02 9.7177E-02 9.1970E-06 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-62 5.7604E-62 0.0000E+00 | minor species not considered
   --- Na+           1 2.1231E-01 2.1230E-01-7.5290E-06 | Normal Major Calc
   --- H2            1 1.4234E-07 1.8146E-07 3.9118E-08 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 7.5290E-06 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-9.2752E-06 |
   --- Cl-           c 2.1231E-01 2.1230E-01-7.5290E-06 |
   --- OH-           c 1.8391E-09 1.8415E-09 2.4867E-12 |
   --- OH            c 2.8469E-07 3.6292E-07 7.8236E-08 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458650E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458835E+03
   --- subroutine FORCE: Beginning Slope = -3.40861e-08
   --- subroutine FORCE: End Slope       = -4.02327e-09
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787689E+00  4.787697E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902832E+00  1.902823E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123107E-01  2.123031E-01 -7.237740E+01 -7.237746E+01
   ---   OH-               1.839058E-09  1.841544E-09 -1.088308E+02 -1.088294E+02
   ---   OH                2.846881E-07  3.629236E-07 -1.079692E+02 -1.077264E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.839058E-09  1.841544E-09 -1.517518E+01 -1.517389E+01 -2.704330E-03 -7.372741E-05
   ---   H2O           1   9.716785E-02  9.717705E-02 -1.240034E+02 -1.240033E+02 -9.436152E-05 -1.436473E-05
   ---   NaCl          0   3.999130E-32  3.999130E-32 -1.745169E+02 -1.745169E+02 -1.247403E-04 -1.270136E-04
   ---   O2            0   5.760390E-62  5.760390E-62 -1.670932E+02 -1.670932E+02  1.677709E+01  1.580592E+01
   ---   Na+           1   2.123107E-01  2.123031E-01 -1.021393E+02 -1.021393E+02  1.285966E-04  1.767514E-05
   ---   H2            1   1.423440E-07  1.814618E-07 -3.289254E+01 -3.264974E+01 -8.244342E-01 -9.607098E-02
   ---                                                        Norms of Delta G():  6.857481E+00  6.452861E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274289E+00
   ---                  air =   4.0971776E+00
   ---              NaCl(S) =   4.7876969E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458650E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458835E+03
   --- Increment counter increased, step is accepted:   11
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  11, Iterations since last evaluation of optimal basis =  10 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8415E-09   6.7886E-14  -7.3727E-05 | Normal calc: diag adjusted from 1.08604e+09 to 1.08604e+09 due to act coeff
   --- H2O             9.7177E-02   1.4002E-06  -1.4365E-05 | Normal calc: diag adjusted from 10.1423 to 10.2591 due to act coeff
   --- NaCl            3.9991E-32   5.0794E-36  -1.2701E-04 | Normal Calc
   --- O2              5.7604E-62   0.0000E+00   1.5806E+01 | Skipped: IC =   0 and DG >0:   1.581E+01
   --- Na+             2.1230E-01  -1.0348E-06   1.7675E-05 | Normal calc: diag adjusted from 7.70186 to 17.0805 due to act coeff
   --- H2              1.8146E-07   5.8111E-09  -9.6071E-02 | Normal calc: diag adjusted from 1.65324e+07 to 1.65324e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8415E-09 1.8416E-09 6.7886E-14 | Normal Major Calc
   --- H2O           1 9.7177E-02 9.7178E-02 1.4002E-06 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.7604E-62 5.7604E-62 0.0000E+00 | minor species not considered
   --- Na+           1 2.1230E-01 2.1230E-01-1.0348E-06 | Normal Major Calc
   --- H2            1 1.8146E-07 1.8727E-07 5.8111E-09 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 1.0348E-06 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-1.4118E-06 |
   --- Cl-           c 2.1230E-01 2.1230E-01-1.0348E-06 |
   --- OH-           c 1.8415E-09 1.8416E-09 6.7886E-14 |
   --- OH            c 3.6292E-07 3.7455E-07 1.1622E-08 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458835E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -5.96679e-10
   --- subroutine FORCE: End Slope       = -1.43262e-11
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787697E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902823E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123031E-01  2.123021E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841544E-09  1.841612E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.629236E-07  3.745458E-07 -1.077264E+02 -1.076949E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841544E-09  1.841612E-09 -1.517389E+01 -1.517386E+01 -7.372741E-05 -9.547559E-06
   ---   H2O           1   9.717705E-02  9.717845E-02 -1.240033E+02 -1.240033E+02 -1.436473E-05 -1.973950E-06
   ---   NaCl          0   3.999130E-32  3.999130E-32 -1.745169E+02 -1.745169E+02 -1.270136E-04 -1.273596E-04
   ---   O2            0   5.760390E-62  5.760390E-62 -1.670932E+02 -1.670932E+02  1.580592E+01  1.567984E+01
   ---   Na+           1   2.123031E-01  2.123021E-01 -1.021393E+02 -1.021393E+02  1.767514E-05  2.690509E-06
   ---   H2            1   1.814618E-07  1.872729E-07 -3.264974E+01 -3.261822E+01 -9.607098E-02 -1.510575E-03
   ---                                                        Norms of Delta G():  6.452861E+00  6.401269E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274254E+00
   ---                  air =   4.0971790E+00
   ---              NaCl(S) =   4.7876979E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458835E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   12
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents

==============================================================================================================
 Iteration =  12, Iterations since last evaluation of optimal basis =  11 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   8.7915E-15  -9.5476E-06 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7178E-02   1.9241E-07  -1.9740E-06 | Normal calc: diag adjusted from 10.1422 to 10.2589 due to act coeff
   --- NaCl            3.9991E-32   5.0933E-36  -1.2736E-04 | Normal Calc
   --- O2              5.7604E-62   0.0000E+00   1.5680E+01 | Skipped: IC =   0 and DG >0:   1.568E+01
   --- Na+             2.1230E-01  -1.5752E-07   2.6905E-06 | Normal calc: diag adjusted from 7.7019 to 17.0805 due to act coeff
   --- H2              1.8727E-07   9.4297E-11  -1.5106E-03 | Normal calc: diag adjusted from 1.60194e+07 to 1.60194e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-09 1.8416E-09 8.7915E-15 | Normal Major Calc
   --- H2O           1 9.7178E-02 9.7179E-02 1.9241E-07 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9996E-32 5.0936E-36 | minor species alternative calc
   --- O2            0 5.7604E-62 8.9286E-69-5.7604E-62 | minor species alternative calc
   --- Na+           1 2.1230E-01 2.1230E-01-1.5752E-07 | Normal Major Calc
   --- H2            1 1.8727E-07 1.8737E-07 9.4297E-11 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 1.5752E-07 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-1.9260E-07 |
   --- Cl-           c 2.1230E-01 2.1230E-01-1.5752E-07 |
   --- OH-           c 1.8416E-09 1.8416E-09 8.7915E-15 |
   --- OH            c 3.7455E-07 3.7473E-07 1.8859E-10 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -9.46064e-13
   --- subroutine FORCE: End Slope       = -1.1568e-13
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123021E-01  2.123020E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841612E-09  1.841621E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.745458E-07  3.747344E-07 -1.076949E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841612E-09  1.841621E-09 -1.517386E+01 -1.517386E+01 -9.547559E-06 -1.507876E-06
   ---   H2O           1   9.717845E-02  9.717864E-02 -1.240033E+02 -1.240033E+02 -1.973950E-06 -3.002171E-07
   ---   NaCl          0   3.999130E-32  3.999639E-32 -1.745169E+02 -1.745168E+02 -1.273596E-04 -4.703145E-08
   ---   O2            0   5.760390E-62  8.928592E-69 -1.670932E+02 -1.827731E+02  1.567984E+01 -2.012933E-03
   ---   Na+           1   2.123021E-01  2.123020E-01 -1.021393E+02 -1.021393E+02  2.690509E-06  3.670431E-07
   ---   H2            1   1.872729E-07  1.873672E-07 -3.261822E+01 -3.261772E+01 -1.510575E-03 -1.039615E-06
   ---                                                        Norms of Delta G():  6.401269E+00  8.217768E-04
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274249E+00
   ---                  air =   4.0971792E+00
   ---              NaCl(S) =   4.7876980E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   13
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  13, Iterations since last evaluation of optimal basis =  12 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   1.3885E-15  -1.5079E-06 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   2.9264E-08  -3.0022E-07 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   1.8811E-39  -4.7031E-08 | Normal Calc
   --- O2              8.9286E-69   1.7973E-71  -2.0129E-03 | Normal Calc
   --- Na+             2.1230E-01  -2.1489E-08   3.6704E-07 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   6.4930E-14  -1.0396E-06 | Normal calc: diag adjusted from 1.60113e+07 to 1.60113e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-09 1.8416E-09 1.3885E-15 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 2.9264E-08 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9286E-69 8.9286E-69 0.0000E+00 | minor species not considered
   --- Na+           1 2.1230E-01 2.1230E-01-2.1489E-08 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 6.4930E-14 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 2.1489E-08 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-2.9264E-08 |
   --- Cl-           c 2.1230E-01 2.1230E-01-2.1489E-08 |
   --- OH-           c 1.8416E-09 1.8416E-09 1.3885E-15 |
   --- OH            c 3.7473E-07 3.7473E-07 1.2986E-13 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -1.66731e-14
   --- subroutine FORCE: End Slope       = -2.39685e-15
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123020E-01  2.123019E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841621E-09  1.841622E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747344E-07  3.747345E-07 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841621E-09  1.841622E-09 -1.517386E+01 -1.517386E+01 -1.507876E-06 -1.983690E-07
   ---   H2O           1   9.717864E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -3.002171E-07 -4.095189E-08
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -4.703145E-08 -5.417400E-08
   ---   O2            0   8.928592E-69  8.928592E-69 -1.827731E+02 -1.827731E+02 -2.012933E-03 -2.014228E-03
   ---   Na+           1   2.123020E-01  2.123019E-01 -1.021393E+02 -1.021393E+02  3.670431E-07  5.576918E-08
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -1.039615E-06 -9.088683E-08
   ---                                                        Norms of Delta G():  8.217768E-04  8.223052E-04
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+00
   ---                  air =   4.0971792E+00
   ---              NaCl(S) =   4.7876981E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   14
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  14, Iterations since last evaluation of optimal basis =  13 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   1.8266E-16  -1.9837E-07 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   3.9918E-09  -4.0952E-08 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   2.1668E-39  -5.4174E-08 | Normal Calc
   --- O2              8.9286E-69   1.7984E-71  -2.0142E-03 | Normal Calc
   --- Na+             2.1230E-01  -3.2651E-09   5.5769E-08 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   5.6764E-15  -9.0887E-08 | Normal calc: diag adjusted from 1.60113e+07 to 1.60113e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-09 1.8416E-09 1.8266E-16 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 3.9918E-09 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9286E-69 8.9286E-69 0.0000E+00 | minor species not considered
   --- Na+           1 2.1230E-01 2.1230E-01-3.2651E-09 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 5.6764E-15 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 3.2651E-09 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-3.9919E-09 |
   --- Cl-           c 2.1230E-01 2.1230E-01-3.2651E-09 |
   --- OH-           c 1.8416E-09 1.8416E-09 1.8266E-16 |
   --- OH            c 3.7473E-07 3.7473E-07 1.1353E-14 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -3.45565e-16
   --- subroutine FORCE: End Slope       = -4.96771e-17
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841622E-09  1.841623E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-07  3.747345E-07 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841622E-09  1.841623E-09 -1.517386E+01 -1.517386E+01 -1.983690E-07 -3.125595E-08
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -4.095189E-08 -6.222294E-09
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -5.417400E-08 -5.514829E-08
   ---   O2            0   8.928592E-69  8.928592E-69 -1.827731E+02 -1.827731E+02 -2.014228E-03 -2.014336E-03
   ---   Na+           1   2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02  5.576918E-08  7.607326E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -9.088683E-08 -1.366988E-08
   ---                                                        Norms of Delta G():  8.223052E-04  8.223491E-04
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+00
   ---                  air =   4.0971792E+00
   ---              NaCl(S) =   4.7876981E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   15
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  15, Iterations since last evaluation of optimal basis =  14 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   2.8781E-17  -3.1256E-08 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   6.0653E-10  -6.2223E-09 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   2.2057E-39  -5.5148E-08 | Normal Calc
   --- O2              8.9286E-69   1.7985E-71  -2.0143E-03 | Normal Calc
   --- Na+             2.1230E-01  -4.4538E-10   7.6073E-09 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   8.5376E-16  -1.3670E-08 | Normal calc: diag adjusted from 1.60113e+07 to 1.60113e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-09 1.8416E-09 2.8781E-17 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 6.0653E-10 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9286E-69 8.9286E-69 0.0000E+00 | minor species not considered
   --- Na+           1 2.1230E-01 2.1230E-01-4.4538E-10 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 8.5376E-16 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 4.4538E-10 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-6.0653E-10 |
   --- Cl-           c 2.1230E-01 2.1230E-01-4.4538E-10 |
   --- OH-           c 1.8416E-09 1.8416E-09 2.8781E-17 |
   --- OH            c 3.7473E-07 3.7473E-07 1.7075E-15 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -7.16216e-18
   --- subroutine FORCE: End Slope       = -1.02961e-18
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841623E-09  1.841623E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-07  3.747345E-07 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841623E-09  1.841623E-09 -1.517386E+01 -1.517386E+01 -3.125595E-08 -4.111385E-09
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -6.222294E-09 -8.487717E-10
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -5.514829E-08 -5.529634E-08
   ---   O2            0   8.928592E-69  8.928592E-69 -1.827731E+02 -1.827731E+02 -2.014336E-03 -2.014352E-03
   ---   Na+           1   2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02  7.607326E-09  1.155868E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -1.366988E-08 -1.883734E-09
   ---                                                        Norms of Delta G():  8.223491E-04  8.223558E-04
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+00
   ---                  air =   4.0971792E+00
   ---              NaCl(S) =   4.7876981E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   16
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species:  MAJOR SPECIES CONVERGENCE achieved
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents
   --- Equilibrium check for minor species: O2 failed

==============================================================================================================
 Iteration =  16, Iterations since last evaluation of optimal basis =  15 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-09 1.8416E-09 0.0000E+00 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 0.0000E+00 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 2.2117E-39 | minor species alternative calc
   --- O2            0 8.9286E-69 8.9466E-69 1.8003E-71 | minor species alternative calc
   --- Na+           1 2.1230E-01 2.1230E-01 0.0000E+00 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 0.0000E+00 | Normal Major Calc
   --- NaCl(S)       c 4.7877E+00 4.7877E+00-2.2117E-39 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00 3.6007E-71 |
   --- Cl-           c 2.1230E-01 2.1230E-01 0.0000E+00 |
   --- OH-           c 1.8416E-09 1.8416E-09 0.0000E+00 |
   --- OH            c 3.7473E-07 3.7473E-07-7.2014E-71 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -1.22296e-46
   --- subroutine FORCE: End Slope       = 1.02338e-84
   --- subroutine FORCE produced no adjustments (al = 1)
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   Cl-               2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01
   ---   OH-               1.841623E-09  1.841623E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-07  3.747345E-07 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841623E-09  1.841623E-09 -1.517386E+01 -1.517386E+01 -4.111385E-09 -4.111385E-09
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -8.487717E-10 -8.487717E-10
   ---   NaCl          0   3.999639E-32  3.999640E-32 -1.745168E+02 -1.745168E+02 -5.529634E-08  0.000000E+00
   ---   O2            0   8.928592E-69  8.946596E-69 -1.827731E+02 -1.827711E+02 -2.014352E-03  5.684342E-14
   ---   Na+           1   2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02  1.155868E-09  1.155868E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -1.883734E-09 -1.883734E-09
   ---                                                        Norms of Delta G():  8.223558E-04  1.936853E-09
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+00
   ---                  air =   4.0971792E+00
   ---              NaCl(S) =   4.7876981E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- Increment counter increased, step is accepted:   17
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   --- Equilibrium check for major species:  MAJOR SPECIES CONVERGENCE achieved
   --- Equilibrium check for minor species:  CONVERGENCE achieved
   --- Subroutine vcs_deltag called for all noncomponents
   --- Equilibrium check for major species:  MAJOR SPECIES CONVERGENCE achieved
   --- Equilibrium check for minor species:  CONVERGENCE achieved
   --- Check the Full Element Abundances:  passed
   --- Start rechecking deleted species in multispec phases
   --- Subroutine vcs_dfe called for all species




--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
		 VCS_TP REPORT
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
		Temperature  =           3e+02 Kelvin
		Pressure     =      1.0132e+05 Pa 
		total Volume =          100.41 m**3


--------------------------------------------------------------------------------
 Species                 Equilibrium kmoles   Mole Fraction    ChemPot/RT    SpecUnkType
--------------------------------------------------------------------------------
 NaCl(S)                   4.7876981E+00      1.0000000E+00     -1.7452E+02     0
 N2                        4.0000000E+00      9.7628143E-01     -2.3069E+01     0
 H2O(L)                    1.9028209E+00      8.1756495E-01     -1.2400E+02     0
 Cl-                       2.1230193E-01      9.1217525E-02     -7.2377E+01     0
 OH-                       1.8416226E-09      7.9127051E-10     -1.0883E+02     0
 OH                        3.7473450E-07      9.1461582E-08     -1.0769E+02     0
 CO2                       0.0000000E+00      0.0000000E+00     -5.0822E+02     0
 Na+                       2.1230193E-01      9.1217525E-02     -1.0214E+02  KMolNum 
 H2O                       9.7178674E-02      2.3718434E-02     -1.2400E+02  KMolNum 
 H2                        1.8736725E-07      4.5730791E-08     -3.2618E+01  KMolNum 
 H+                        1.8416226E-09      7.9127051E-10     -1.5174E+01  KMolNum 
 NaCl                      3.9996395E-32      9.7619345E-33     -1.7452E+02  KMolNum 
 O2                        8.9465959E-69      2.1835989E-69     -1.8277E+02  KMolNum 
--------------------------------------------------------------------------------


-------------------------------------------------------------------------------------------------------------------
               |ComponentID|          0          1          2          3          4          5          6 |           |
               | Components|    NaCl(S)         N2     H2O(L)        Cl-        OH-         OH        CO2 |           |
 NonComponent  |   Moles   |       4.79          4        1.9      0.212   1.84e-09   3.75e-07          0 | DG/RT Rxn |
-------------------------------------------------------------------------------------------------------------------
   7 H+        |  1.84e-09 |       0.00       0.00      -1.00       0.00       1.00       0.00       0.00 | -4.11e-09 |
   8 H2O       |    0.0972 |       0.00       0.00      -1.00       0.00       0.00       0.00       0.00 | -8.49e-10 |
   9 NaCl      |     4e-32 |      -1.00       0.00       0.00       0.00       0.00       0.00       0.00 |         0 |
  10 O2        |  8.95e-69 |       0.00       0.00       2.00       0.00       0.00      -4.00       0.00 |  5.68e-14 |
  11 Na+       |     0.212 |      -1.00       0.00       0.00       1.00       0.00       0.00       0.00 |  1.16e-09 |
  12 H2        |  1.87e-07 |       0.00       0.00      -2.00       0.00       0.00       2.00       0.00 | -1.88e-09 |
-------------------------------------------------------------------------------------------------------------------




------------------------------------------------------------------------------------------------------------------------------------------------------------------------
                  | ElementID |          0          1          2          3          4          5          6          7          8          9         10 |                     |
                  | Element   |          O          H          C          N         Na         Cl cn_NaCl_el         Fe          E         Si         Ca |                     |
    PhaseName     |KMolTarget |          2          4          0          8          5          5          0          0          0          0          0 |     Gibbs Total     |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
   0 NaCl_electro | 2.327e+00 |        1.9       3.81          0          0      0.212      0.212          0          0          0          0          0 | -2.73006234970E+02 |
   1 air          | 4.097e+00 |     0.0972      0.194          0          8      4e-32      4e-32          0          0          0          0          0 | -1.04327434116E+02 |
   2 NaCl(S)      | 4.788e+00 |          0          0          0          0       4.79       4.79          0          0          0          0          0 | -8.35533676798E+02 |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    TOTAL         | 1.121e+01 |          2          4          0          8          5          5          0          0          0          0          0 | -1.21286734588E+03 |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------


	Total Dimensionless Gibbs Free Energy = G/RT =  -1.2128673E+03

Elemental Abundances (kmol):          Actual                    Target         Type      ElActive
                          O   2.000000000000E+00    2.000000000000E+00     0       1
                          H   4.000000000000E+00    4.000000000000E+00     0       1
                          C   0.000000000000E+00    0.000000000000E+00     0       1
                          N   8.000000000000E+00    8.000000000000E+00     0       1
                          Na  5.000000000000E+00    5.000000000000E+00     0       1
                          Cl  5.000000000000E+00    5.000000000000E+00     0       1
                          cn  0.000000000000E+00    0.000000000000E+00     2       1
                          Fe  0.000000000000E+00    0.000000000000E+00     0       1
                          E   0.000000000000E+00    0.000000000000E+00     1       0
                          Si  0.000000000000E+00    0.000000000000E+00     0       1
                          Ca  0.000000000000E+00    0.000000000000E+00     0       1


---------------------------------------------------------------------------------------------
Chemical Potentials of the Species: (dimensionless)
		(RT = 2.47896e+06 J/kmol)
    Name        TKMoles     StandStateChemPot      ln(AC)       ln(X_i)      |   F z_i phi   |    ChemPot    | (-lnMnaught)|  (MolNum ChemPot)|
---------------------------------------------------------------------------------------------------------------------------------------------------
 NaCl(S)       4.7876981E+00 -1.7451679E+02   0.0000000E+00   0.0000000E+00  | 0.0000000E+00 |   -1.7452E+02 |             | -8.3553367679750E+02 |
 N2            4.0000000E+00 -2.3045225E+01   0.0000000E+00  -2.4004385E-02  | 0.0000000E+00 |   -2.3069E+01 |             | -9.2276916077697E+01 |
 H2O(L)        1.9028209E+00 -1.2371551E+02  -8.6317426E-02  -2.0142493E-01  | 0.0000000E+00 |   -1.2400E+02 |             | -2.3595598432279E+02 |
 Cl-           2.1230193E-01 -7.4214051E+01   2.1455966E-01  -2.3945082E+00  |-0.0000000E+00 |   -7.2377E+01 |(4.01653E+00)| -1.5365875462838E+01 |
 OH-           1.8416226E-09 -9.1483483E+01  -4.0506365E-01  -2.0957381E+01  |-0.0000000E+00 |   -1.0883E+02 |(4.01653E+00)| -2.0042267118620E-07 |
 OH            3.7473450E-07 -9.1487045E+01   0.0000000E+00  -1.6207347E+01  | 0.0000000E+00 |   -1.0769E+02 |             | -4.0356803750687E-05 |
 CO2           0.0000000E+00 -1.8445182E+02   0.0000000E+00  -3.2377221E+02  | 0.0000000E+00 |   -5.0822E+02 |             | -0.0000000000000E+00 |
 Na+           2.1230193E-01 -1.0397591E+02   2.1455966E-01  -2.3945082E+00  | 0.0000000E+00 |   -1.0214E+02 |(4.01653E+00)| -2.1684374956081E+01 |
 H2O           9.7178674E-02 -1.2026175E+02   0.0000000E+00  -3.7415027E+00  | 0.0000000E+00 |   -1.2400E+02 |             | -1.2050471570195E+01 |
 H2            1.8736725E-07 -1.5717224E+01   0.0000000E+00  -1.6900494E+01  | 0.0000000E+00 |   -3.2618E+01 |             | -6.1114920650939E-06 |
 H+            1.8416226E-09  0.0000000E+00   1.7669883E+00  -2.0957381E+01  | 0.0000000E+00 |   -1.5174E+01 |(4.01653E+00)| -2.7944519907221E-08 |
 NaCl          3.9996395E-32 -1.0080997E+02   0.0000000E+00  -7.3706817E+01  | 0.0000000E+00 |   -1.7452E+02 |             | -6.9800423137655E-30 |
 O2            8.9465959E-69 -2.4673669E+01   0.0000000E+00  -1.5809740E+02  | 0.0000000E+00 |   -1.8277E+02 |             | -1.6351788748965E-66 |
                                                                                                                             -1.2128673458838E+03
---------------------------------------------------------------------------------------------------------------------------------------------------


Counters:         Iterations          Time (seconds)
    vcs_basopt:       2                 NA     
    vcs_TP:          17                 NA     
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------

TCounters:   Num_Calls   Total_Its       Total_Time (seconds)
    vcs_basopt:       2          2             NA     
    vcs_TP:           1         17             NA     
    vcs_inest:        1                        NA     
    vcs_TotalTime:                             NA     

 Results from vcs:

Temperature = 298.15 Kelvin
Pressure    = 101325 Pa

-------------------------------------------------------------
 Name             Mole_Number(kmol)  Mole_Fraction     Chem_Potential (J/kmol)
-------------------------------------------------------------
H2O(L)              1.903e+00         8.176e-01       -3.074e+08
Cl-                 2.123e-01         9.122e-02       -1.794e+08
H+                  1.842e-09         7.913e-10       -3.762e+07
Na+                 2.123e-01         9.122e-02       -2.532e+08
OH-                 1.842e-09         7.913e-10       -2.698e+08
O2                  8.947e-69         2.184e-69       -4.531e+08
H2                  1.874e-07         4.573e-08       -8.086e+07
CO2                 0.000e+00         0.000e+00       -1.000e+300
H2O                 9.718e-02         2.372e-02       -3.074e+08
NaCl                4.000e-32         9.762e-33       -4.326e+08
N2                  4.000e+00         9.763e-01       -5.719e+07
OH                  3.747e-07         9.146e-08       -2.670e+08
NaCl(S)             4.788e+00         1.000e+00       -4.326e+08
-------------------------------------------------------------
*************** NaCl_electrolyte *****************
Moles: 2.32742
 
  NaCl_electrolyte:
 
       temperature          298.15  K
          pressure          101325  Pa
           density         1216.41  kg/m^3
  mean mol. weight         20.0596  amu
         potential               0  V
                pH            6.59  
 
                          1 kg            1 kmol
                       -----------      ------------
          enthalpy    -1.35106e+07        -2.71e+08     J
   internal energy    -1.35107e+07        -2.71e+08     J
           entropy          3304.6        6.629e+04     J/K
    Gibbs function    -1.44958e+07       -2.908e+08     J
 heat capacity c_p         3053.16        6.125e+04     J/K
 heat capacity c_v    <not implemented>       
 
                           X           Molalities         Chem.Pot.    ChemPotSS    ActCoeffMolal
                                                          (J/kmol)      (J/kmol)                 
                     -------------    ------------     ------------  ------------    ------------
            H2O(L)      0.817565          55.5084     -3.07399e+08   -3.06686e+08       0.917303
               Cl-     0.0912175           6.1932     -1.79421e+08   -1.83974e+08        1.01322
                H+   7.91271e-10      5.37232e-08     -3.76154e+07              0        4.78537
               Na+     0.0912175           6.1932     -2.53199e+08   -2.57752e+08        1.01322
               OH-   7.91271e-10      5.37232e-08     -2.69784e+08   -2.26784e+08       0.545262

*************** air *****************
Moles: 4.09718
 
  air:
 
       temperature          298.15  K
          pressure          101325  Pa
           density         1.13533  kg/m^3
  mean mol. weight         27.7763  amu
 
                          1 kg            1 kmol
                       -----------      ------------
          enthalpy         -206497       -5.736e+06     J
   internal energy         -295744       -8.215e+06     J
           entropy         6929.48        1.925e+05     J/K
    Gibbs function    -2.27252e+06       -6.312e+07     J
 heat capacity c_p         1052.34        2.923e+04     J/K
 heat capacity c_v         753.001        2.092e+04     J/K
 
                           X                 Y          Chem. Pot. / RT    
                     -------------     ------------     ------------
                O2     2.1836e-69      2.51554e-69         -182.771
                H2    4.57308e-08      3.31893e-09         -32.6177
               CO2              0                0     
               H2O      0.0237184        0.0153834         -124.003
              NaCl    9.76193e-33      2.05395e-32         -174.517
                N2       0.976281         0.984617         -23.0692
                OH    9.14616e-08      5.60016e-08         -107.694

*************** NaCl(S) *****************
Moles: 4.7877
 
  NaCl(S):
 
       temperature          298.15  K
          pressure          101325  Pa
           density            2165  kg/m^3
  mean mol. weight         58.4425  amu
 
                          1 kg            1 kmol
                       -----------      ------------
          enthalpy    -7.03462e+06       -4.111e+08     J
   internal energy    -7.03467e+06       -4.111e+08     J
           entropy         1233.85        7.211e+04     J/K
    Gibbs function     -7.4025e+06       -4.326e+08     J
 heat capacity c_p         864.119         5.05e+04     J/K
 heat capacity c_v         864.119         5.05e+04     J/K
 
                           X                 Y          Chem. Pot. / RT    
                     -------------     ------------     ------------
           NaCl(S)              1                1         -174.517

NUMBER OF SUCCESSES =  1
NUMBER OF FAILURES  =  0
