Unknown Cantera EOS to VCSnonideal: 45012
vcs_Cantera_convert: Species Type 4 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
    5 PHASE1 SPECIES         7 PHASE2 SPECIES       1 SINGLE SPECIES 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  1              -2.57752E+08   5.00000E+00       Mol_Num
 Cl-           0  0  0  0  0  1 -1  0  1  0  0  1              -1.83974E+08   5.00000E+00       Mol_Num
 N2            0  0  0  2  0  0  0  0  0  0  0  2              -5.71282E+07   4.00000E+00       Mol_Num
 H2O(L)        1  2  0  0  0  0  0  0 -0  0  0  1              -3.06686E+08   1.87965E+00       Mol_Num
 OH-           1  1  0  0  0  0 -1  0  1  0  0  1              -2.26784E+08   1.20348E-01       Mol_Num
 CO2           2  0  1  0  0  0  0  0  0  0  0  2              -4.57249E+08   0.00000E+00       Mol_Num
 H2            0  2  0  0  0  0  0  0  0  0  0  2              -3.89624E+07   1.86050E-14       Mol_Num
 H+            0  1  0  0  0  0  1  0 -1  0  0  1               0.00000E+00   1.20348E-01       Mol_Num
 H2O           1  2  0  0  0  0  0  0  0  0  0  2              -2.98124E+08   2.25617E-20       Mol_Num
 NaCl          0  0  0  0  1  1  0  0  0  0  0  2              -2.49904E+08   3.54670E-44       Mol_Num
 O2            2  0  0  0  0  0  0  0  0  0  0  2              -6.11650E+07   1.54458E-41       Mol_Num
 NaCl(S)       0  0  0  0  1  1  0  0  0  0  0  0              -4.32620E+08   0.00000E+00       Mol_Num
 OH            1  1  0  0  0  0  0  0  0  0  0  2              -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
   -----------------------------------------------------------------------------
  --- MAJOR-MINOR decision is reavaluated: All species are minor except for:
  ---      Major Species          : H+
  ---      Zeroed-SS Phase Species: NaCl(S)
  ---      Major Species          : H2
  ---
   --- Element Abundance check passed
   --- Subroutine vcs_deltag called for all noncomponents

==============================================================================================================
 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
   --- Zeroed species changed to major: NaCl(S)     
   --- NaCl(S)       1 0.0000E+00 5.0000E-02 5.0000E-02 | Born:IC=-1 to IC=1:DG=-7.0667E+01
   --- H2            1 1.8605E-14 1.8605E-15-1.6744E-14 | initial nonpos kmoles=   0.000E+00
   --- Na+           c 5.0000E+00 4.9496E+00-5.0400E-02 |
   --- Cl-           c 5.0000E+00 4.9496E+00-5.0400E-02 |
   --- 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
   --- 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.0101404441423E+03
   --- subroutine FORCE: Beginning Slope = -17.0406
   --- subroutine FORCE: End Slope       = -16.2472
   --- 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
   ---   Na+               5.000000E+00  4.949600E+00 -6.616232E+01 -6.651394E+01
   ---   Cl-               5.000000E+00  4.949600E+00 -3.768766E+01 -3.686990E+01
   ---   N2                4.000000E+00  4.000000E+00 -2.304522E+01 -2.304532E+01
   ---   H2O(L)            1.879652E+00  1.987965E+00 -1.986003E+02 -1.994452E+02
   ---   OH-               1.203483E-01  1.203483E-02 -9.196202E+01 -9.656442E+01
   ---   OH                3.720999E-14  3.720999E-15 -1.237955E+02 -1.260982E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082000E+02 -5.082001E+02
   ---   H+            1   1.203483E-01  1.203483E-02  1.770351E+01  1.425959E+01  1.243418E+02  1.171404E+02
   ---   H2O           0   2.256174E-20  2.256173E-30 -1.668861E+02 -1.899120E+02  3.171423E+01  9.533195E+00
   ---   NaCl          0   3.546703E-44  4.000000E-04 -2.022440E+02 -1.100204E+02 -9.839401E+01 -6.636575E+00
   ---   O2            0   1.544580E-41  5.819043E-32 -1.200312E+02 -9.798165E+01 -2.204965E+01  7.520803E+00
   ---   NaCl(S)       1   0.000000E+00  5.000000E-02 -1.745168E+02 -1.745168E+02 -7.066681E+01 -7.113295E+01
   ---   H2            1   1.860500E-14  1.860500E-15 -4.871886E+01 -5.102155E+01  1.008907E+02  9.567244E+01
   ---                                                        Norms of Delta G():  8.347351E+01  6.846711E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   1.1911235E+01
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   5.0000000E-02
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -9.9366705866970E+02
   --- Total New Dimensionless Gibbs Free Energy = -1.0101404441423E+03
   --- Increment counter increased, step is accepted:    1
   --- 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 =   1, Iterations since last evaluation of optimal basis =   1 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-02  -1.2035E-02   1.1714E+02 | Delta damped from -0.84723 to -0.0120348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-30   0.0000E+00   9.5332E+00 | Skipped: IC =   0 and DG >0:   9.533E+00
   --- NaCl            4.0000E-04   2.6545E-03  -6.6366E+00 | Normal calc: diag adjusted from 2499.82 to 2500.15 due to act coeff
   --- O2              5.8190E-32   0.0000E+00   7.5208E+00 | Skipped: IC =   0 and DG >0:   7.521E+00
   --- NaCl(S)         5.0000E-02   4.9496E+00  -7.1133E+01 | Delta damped from 176.04 to 4.9496 due to component 0 (       Na+) going neg
   --- H2              1.8605E-15  -1.8605E-15   9.5672E+01 | Delta damped from -5.93329e-14 to -1.8605e-15 due to component 5 (        OH) going neg
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 1.2035E-02 1.2035E-03-1.0831E-02 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 2.2562E-30 2.2562E-30 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-04 4.0000E-04 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-32 5.8190E-32 0.0000E+00 | minor species not considered
   --- NaCl(S)       1 5.0000E-02 2.5248E+00 2.4748E+00 | Line Search reduced step size from 4.9496 to 2.4748
   --- H2            1 1.8605E-15 1.8605E-16-1.6744E-15 | initial nonpos kmoles=   0.000E+00
   --- Na+           c 4.9496E+00 2.4748E+00-2.4748E+00 |
   --- Cl-           c 4.9496E+00 2.4748E+00-2.4748E+00 |
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- H2O(L)        c 1.9880E+00 1.9988E+00 1.0831E-02 |
   --- OH-           c 1.2035E-02 1.2035E-03-1.0831E-02 |
   --- OH            c 3.7210E-15 3.7210E-16-3.3489E-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.0101404441423E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.1714270977885E+03
   --- subroutine FORCE: Beginning Slope = -177.309
   --- subroutine FORCE: End Slope       = -100.647
   --- 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
   ---   Na+               4.949600E+00  2.474800E+00 -6.651394E+01 -8.195274E+01
   ---   Cl-               4.949600E+00  2.474800E+00 -3.686990E+01 -5.219880E+01
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   H2O(L)            1.987965E+00  1.998797E+00 -1.994452E+02 -1.533335E+02
   ---   OH-               1.203483E-02  1.203483E-03 -9.656442E+01 -9.615896E+01
   ---   OH                3.720999E-15  3.720999E-16 -1.260982E+02 -1.284008E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.203483E-02  1.203483E-03  1.425959E+01  1.214346E+01  1.171404E+02  6.931801E+01
   ---   H2O           0   2.256173E-30  2.256173E-30 -1.899120E+02 -1.899120E+02  9.533195E+00 -3.657852E+01
   ---   NaCl          0   4.000000E-04  4.000000E-04 -1.100204E+02 -1.100204E+02 -6.636575E+00  2.413113E+01
   ---   O2            0   5.819043E-32  5.819043E-32 -9.798165E+01 -9.798165E+01  7.520803E+00  1.089546E+02
   ---   NaCl(S)       1   5.000000E-02  2.524800E+00 -1.745168E+02 -1.745168E+02 -7.113295E+01 -4.036525E+01
   ---   H2            1   1.860500E-15  1.860500E-16 -5.102155E+01 -5.332413E+01  9.567244E+01 -3.458737E+00
   ---                                                        Norms of Delta G():  6.846711E+01  5.807710E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   6.9508035E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   2.5248000E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.0101404441423E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.1714270977885E+03
   --- Increment counter increased, step is accepted:    2
   --- 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

   ---   N2          (        4) replaces Na+         (      2.5) as component   0
   ---   NaCl(S)     (      2.5) replaces Cl-         (      2.5) as component   1
   ---   Na+         (      2.5) remains                          as component   2
   ---   H2O(L)      (        2) remains                          as component   3
   ---   OH-         (   0.0012) remains                          as component   4
   ---   OH          (  3.7e-16) 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      2.52      2.47         2    0.0012  3.72e-16         0
   ---   NonComponent|   Moles  |       N2        NaCl(S)   Na+       H2O(L)    OH-       OH        CO2       
   ---   7 H+        |    0.0012|      -0.00      -0.00       0.00      -1.00       1.00      -0.00      -0.00
   ---   8 H2O       |  2.26e-30|      -0.00      -0.00       0.00      -1.00      -0.00      -0.00      -0.00
   ---   9 NaCl      |    0.0004|      -0.00      -1.00       0.00       0.00      -0.00      -0.00      -0.00
   ---  10 O2        |  5.82e-32|      -0.00      -0.00       0.00       2.00      -0.00      -4.00      -0.00
   ---  11 Cl-       |      2.47|      -0.00      -1.00       1.00       0.00      -0.00      -0.00      -0.00
   ---  12 H2        |  1.86e-16|      -0.00      -0.00       0.00      -2.00      -0.00       2.00      -0.00
   -----------------------------------------------------------------------------
  --- MAJOR-MINOR decision is reavaluated: All species are minor except for:
  ---      Major Species          : H+
  ---      Major Species          : Cl-
  ---      Major Species          : H2
  ---
   --- Element Abundance check passed
   --- Subroutine vcs_deltag called for all noncomponents
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents

==============================================================================================================
 Iteration =   2, Iterations since last evaluation of optimal basis =   0 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-03  -1.2035E-03   6.9318E+01 | Delta damped from -0.0403055 to -0.00120348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-30   8.2527E-29  -3.6579E+01 | Normal Calc
   --- NaCl            4.0000E-04   0.0000E+00   6.4496E+01 | Skipped: IC =   0 and DG >0:   6.450E+01
   --- O2              5.8190E-32   0.0000E+00   1.0895E+02 | Skipped: IC =   0 and DG >0:   1.090E+02
   --- Cl-             2.4748E+00  -1.7109E+00   4.0365E+01 | Normal calc: diag adjusted from 0.232673 to 23.5929 due to act coeff
   --- H2              1.8605E-16   2.1450E-16  -3.4587E+00 | Normal calc: diag adjusted from 1.61247e+16 to 1.61247e+16 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.2035E-03 1.2035E-04-1.0831E-03 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 2.2562E-30 1.7347E-14 1.7347E-14 | minor species alternative calc
   --- NaCl          0 4.0000E-04 4.0000E-14-4.0000E-04 | minor species alternative calc
   --- O2            0 5.8190E-32 5.8190E-42-5.8190E-32 | minor species alternative calc
   --- Cl-           1 2.4748E+00 7.6389E-01-1.7109E+00 | Normal Major Calc
   --- H2            1 1.8605E-16 4.0055E-16 2.1450E-16 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 2.5248E+00 4.2361E+00 1.7113E+00 |
   --- Na+           c 2.4748E+00 7.6389E-01-1.7109E+00 |
   --- H2O(L)        c 1.9988E+00 1.9999E+00 1.0831E-03 |
   --- OH-           c 1.2035E-03 1.2035E-04-1.0831E-03 |
   --- OH            c 3.7210E-16 8.0110E-16 4.2900E-16 |
   --- 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.1714270977885E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2106126373154E+03
   --- subroutine FORCE: Beginning Slope = -69.1622
   --- subroutine FORCE: End Slope       = -13.5717
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304522E+01
   ---   NaCl(S)           2.524800E+00  4.236110E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.474800E+00  7.638904E-01 -8.195274E+01 -9.818712E+01
   ---   H2O(L)            1.998797E+00  1.999880E+00 -1.533335E+02 -1.259446E+02
   ---   OH-               1.203483E-03  1.203483E-04 -9.615896E+01 -9.745474E+01
   ---   OH                3.720999E-16  8.010985E-16 -1.284008E+02 -1.276339E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082000E+02
   ---   H+            1   1.203483E-03  1.203483E-04  1.214346E+01  8.605827E-01  6.931801E+01  2.935049E+01
   ---   H2O           0   2.256173E-30  1.734686E-14 -1.899120E+02 -1.533334E+02 -3.657852E+01 -2.738877E+01
   ---   NaCl          0   4.000000E-04  4.000001E-14 -1.100204E+02 -1.330462E+02  6.449638E+01  4.147063E+01
   ---   O2            0   5.819043E-32  5.819045E-42 -9.798165E+01 -1.210074E+02  1.089546E+02  1.376389E+02
   ---   Cl-           1   2.474800E+00  7.638904E-01 -5.219880E+01 -6.842548E+01  4.036525E+01  7.904190E+00
   ---   H2            1   1.860500E-16  4.005493E-16 -5.332413E+01 -5.255721E+01 -3.458737E+00 -5.593571E+01
   ---                                                        Norms of Delta G():  6.300155E+01  6.515002E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   3.5279012E+00
   ---                  air =   4.0000000E+00
   ---              NaCl(S) =   4.2361096E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.1714270977885E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2106126373154E+03
   --- Increment counter increased, step is accepted:    3
   --- 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 =   3, Iterations since last evaluation of optimal basis =   1 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-04  -1.2035E-04   2.9350E+01 | Delta damped from -0.00176449 to -0.000120348 due to component 4 (       OH-) going neg
   --- H2O             1.7347E-14   4.7511E-13  -2.7389E+01 | Normal calc: diag adjusted from 5.76473e+13 to 5.76473e+13 due to act coeff
   --- NaCl            4.0000E-14   0.0000E+00   4.1471E+01 | Skipped: IC =   0 and DG >0:   4.147E+01
   --- O2              5.8190E-42   0.0000E+00   1.3764E+02 | Skipped: IC =   0 and DG >0:   1.376E+02
   --- Cl-             7.6389E-01  -5.3361E-01   7.9042E+00 | Normal calc: diag adjusted from 1.48436 to 14.8128 due to act coeff
   --- H2              4.0055E-16   7.4683E-15  -5.5936E+01 | Normal calc: diag adjusted from 7.48972e+15 to 7.48972e+15 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.2035E-04 1.2035E-05-1.0831E-04 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 1.7347E-14 1.7347E-14 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-14 4.0000E-14 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-42 5.8190E-42 0.0000E+00 | minor species not considered
   --- Cl-           1 7.6389E-01 2.3029E-01-5.3361E-01 | Normal Major Calc
   --- H2            1 4.0055E-16 7.8689E-15 7.4683E-15 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.2361E+00 4.7697E+00 5.3361E-01 |
   --- Na+           c 7.6389E-01 2.3029E-01-5.3361E-01 |
   --- H2O(L)        c 1.9999E+00 2.0000E+00 1.0831E-04 |
   --- OH-           c 1.2035E-04 1.2035E-05-1.0831E-04 |
   --- OH            c 8.0110E-16 1.5738E-14 1.4937E-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.2106126373154E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127707282422E+03
   --- subroutine FORCE: Beginning Slope = -4.2209
   --- subroutine FORCE: End Slope       = -0.0635798
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.236110E+00  4.769715E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               7.638904E-01  2.302853E-01 -9.818712E+01 -1.020815E+02
   ---   H2O(L)            1.999880E+00  1.999988E+00 -1.259446E+02 -1.240164E+02
   ---   OH-               1.203483E-04  1.203483E-05 -9.745474E+01 -1.000912E+02
   ---   OH                8.010985E-16  1.573777E-14 -1.276339E+02 -1.246561E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082000E+02 -5.082000E+02
   ---   H+            1   1.203483E-04  1.203483E-05  8.605827E-01 -6.369306E+00  2.935049E+01  1.755590E+01
   ---   H2O           0   1.734686E-14  1.734686E-14 -1.533334E+02 -1.533334E+02 -2.738877E+01 -2.931704E+01
   ---   NaCl          0   4.000001E-14  4.000001E-14 -1.330462E+02 -1.330462E+02  4.147063E+01  4.147063E+01
   ---   O2            0   5.819045E-42  5.819045E-42 -1.210074E+02 -1.210074E+02  1.376389E+02  1.295841E+02
   ---   Cl-           1   7.638904E-01  2.302853E-01 -6.842548E+01 -7.231967E+01  7.904190E+00  1.155879E-01
   ---   H2            1   4.005493E-16  7.868885E-15 -5.255721E+01 -4.957938E+01 -5.593571E+01 -5.085875E+01
   ---                                                        Norms of Delta G():  6.515002E+01  6.091821E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4605826E+00
   ---                  air =   4.0000000E+00
   ---              NaCl(S) =   4.7697147E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2106126373154E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127707282422E+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

==============================================================================================================
 Iteration =   4, Iterations since last evaluation of optimal basis =   2 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-05  -1.2035E-05   1.7556E+01 | Delta damped from -0.00010564 to -1.20348e-05 due to component 4 (       OH-) going neg
   --- H2O             1.7347E-14   5.0856E-13  -2.9317E+01 | Normal calc: diag adjusted from 5.76473e+13 to 5.76473e+13 due to act coeff
   --- NaCl            4.0000E-14   0.0000E+00   4.1471E+01 | Skipped: IC =   0 and DG >0:   4.147E+01
   --- O2              5.8190E-42   0.0000E+00   1.2958E+02 | Skipped: IC =   0 and DG >0:   1.296E+02
   --- Cl-             2.3029E-01  -7.1861E-03   1.1559E-01 | Normal calc: diag adjusted from 7.05925 to 16.0849 due to act coeff
   --- H2              7.8689E-15   1.3340E-13  -5.0859E+01 | Normal calc: diag adjusted from 3.81248e+14 to 3.81248e+14 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.2035E-05 1.2035E-06-1.0831E-05 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 1.7347E-14 1.7347E-14 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-14 4.0000E-14 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-42 5.8190E-42 0.0000E+00 | minor species not considered
   --- Cl-           1 2.3029E-01 2.2314E-01-7.1495E-03 | Line Search reduced step size from -0.00718611 to -0.00714954
   --- H2            1 7.8689E-15 1.4127E-13 1.3340E-13 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7697E+00 4.7769E+00 7.1495E-03 |
   --- Na+           c 2.3029E-01 2.2314E-01-7.1495E-03 |
   --- H2O(L)        c 2.0000E+00 2.0000E+00 1.0831E-05 |
   --- OH-           c 1.2035E-05 1.2035E-06-1.0831E-05 |
   --- OH            c 1.5738E-14 2.8254E-13 2.6680E-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.2127707282422E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713152646E+03
   --- subroutine FORCE: Beginning Slope = -0.00101655
   --- subroutine FORCE: End Slope       = -0.00013844
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.769715E+00  4.776864E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.302853E-01  2.231357E-01 -1.020815E+02 -1.021394E+02
   ---   H2O(L)            1.999988E+00  1.999999E+00 -1.240164E+02 -1.240032E+02
   ---   OH-               1.203483E-05  1.203483E-06 -1.000912E+02 -1.023969E+02
   ---   OH                1.573777E-14  2.825389E-13 -1.246561E+02 -1.217683E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082000E+02 -5.082000E+02
   ---   H+            1   1.203483E-05  1.203483E-06 -6.369306E+00 -8.741424E+00  1.755590E+01  1.286494E+01
   ---   H2O           0   1.734686E-14  1.734686E-14 -1.533334E+02 -1.533334E+02 -2.931704E+01 -2.933017E+01
   ---   NaCl          0   4.000001E-14  4.000001E-14 -1.330462E+02 -1.330462E+02  4.147063E+01  4.147063E+01
   ---   O2            0   5.819045E-42  5.819045E-42 -1.210074E+02 -1.210074E+02  1.295841E+02  1.180593E+02
   ---   Cl-           1   2.302853E-01  2.231357E-01 -7.231967E+01 -7.237753E+01  1.155879E-01 -1.264886E-04
   ---   H2            1   7.868885E-15  1.412694E-13 -4.957938E+01 -4.669163E+01 -5.085875E+01 -4.222175E+01
   ---                                                        Norms of Delta G():  6.091821E+01  5.547711E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462727E+00
   ---                  air =   4.0000000E+00
   ---              NaCl(S) =   4.7768643E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127707282422E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713152646E+03
   --- Increment counter increased, step is accepted:    5
   --- 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 =   5, Iterations since last evaluation of optimal basis =   3 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-06  -1.2035E-06   1.2865E+01 | Delta damped from -7.74136e-06 to -1.20348e-06 due to component 4 (       OH-) going neg
   --- H2O             1.7347E-14   5.0879E-13  -2.9330E+01 | Normal calc: diag adjusted from 5.76473e+13 to 5.76473e+13 due to act coeff
   --- NaCl            4.0000E-14   0.0000E+00   4.1471E+01 | Skipped: IC =   0 and DG >0:   4.147E+01
   --- O2              5.8190E-42   0.0000E+00   1.1806E+02 | Skipped: IC =   0 and DG >0:   1.181E+02
   --- Cl-             2.2314E-01   7.7835E-06  -1.2649E-04 | Normal calc: diag adjusted from 7.32801 to 16.2508 due to act coeff
   --- H2              1.4127E-13   1.9882E-12  -4.2222E+01 | Normal calc: diag adjusted from 2.1236e+13 to 2.1236e+13 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.2035E-06 1.2035E-07-1.0831E-06 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 1.7347E-14 4.0000E-04 4.0000E-04 | minor species alternative calc
   --- NaCl          0 4.0000E-14 4.0000E-24-4.0000E-14 | minor species alternative calc
   --- O2            0 5.8190E-42 5.8190E-52-5.8190E-42 | minor species alternative calc
   --- Cl-           1 2.2314E-01 2.2314E-01 7.7835E-06 | Normal Major Calc
   --- H2            1 1.4127E-13 2.1295E-12 1.9882E-12 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7769E+00-7.7835E-06 |
   --- Na+           c 2.2314E-01 2.2314E-01 7.7835E-06 |
   --- H2O(L)        c 2.0000E+00 1.9996E+00-3.9892E-04 |
   --- OH-           c 1.2035E-06 1.2035E-07-1.0831E-06 |
   --- OH            c 2.8254E-13 4.2590E-12 3.9764E-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.2127713152646E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127739151271E+03
   --- subroutine FORCE: Beginning Slope = -0.011746
   --- subroutine FORCE: End Slope       = -0.00219649
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304522E+01 -2.304532E+01
   ---   NaCl(S)           4.776864E+00  4.776856E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231357E-01  2.231435E-01 -1.021394E+02 -1.021390E+02
   ---   H2O(L)            1.999999E+00  1.999600E+00 -1.240032E+02 -1.240033E+02
   ---   OH-               1.203483E-06  1.203483E-07 -1.023969E+02 -1.046992E+02
   ---   OH                2.825389E-13  4.258968E-12 -1.217683E+02 -1.190554E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082000E+02 -5.082001E+02
   ---   H+            1   1.203483E-06  1.203483E-07 -8.741424E+00 -1.104330E+01  1.286494E+01  8.260791E+00
   ---   H2O           0   1.734686E-14  4.000000E-04 -1.533334E+02 -1.294722E+02 -2.933017E+01 -5.468858E+00
   ---   NaCl          0   4.000001E-14  4.000004E-24 -1.330462E+02 -1.560721E+02  4.147063E+01  1.844468E+01
   ---   O2            0   5.819045E-42  5.819042E-52 -1.210074E+02 -1.440334E+02  1.180593E+02  8.418172E+01
   ---   Cl-           1   2.231357E-01  2.231435E-01 -7.237753E+01 -7.237711E+01 -1.264886E-04  7.124068E-04
   ---   H2            1   1.412694E-13  2.129484E-12 -4.669163E+01 -4.397876E+01 -4.222175E+01 -3.408297E+01
   ---                                                        Norms of Delta G():  5.547711E+01  3.804945E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4458872E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7768565E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713152646E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127739151271E+03
   --- Increment counter increased, step is accepted:    6
   --- 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 =   6, Iterations since last evaluation of optimal basis =   4 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-07  -1.2035E-07   8.2608E+00 | Delta damped from -4.97086e-07 to -1.20348e-07 due to component 4 (       OH-) going neg
   --- H2O             4.0000E-04   2.1876E-03  -5.4689E+00 | Normal calc: diag adjusted from 2499.84 to 2499.95 due to act coeff
   --- NaCl            4.0000E-24   0.0000E+00   1.8445E+01 | Skipped: IC =   0 and DG >0:   1.844E+01
   --- O2              5.8190E-52   0.0000E+00   8.4182E+01 | Skipped: IC =   0 and DG >0:   8.418E+01
   --- Cl-             2.2314E-01  -4.3833E-05   7.1241E-04 | Normal calc: diag adjusted from 7.32744 to 16.2528 due to act coeff
   --- H2              2.1295E-12   2.4193E-11  -3.4083E+01 | Normal calc: diag adjusted from 1.40879e+12 to 1.40879e+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.2035E-07 1.2035E-08-1.0831E-07 | initial nonpos kmoles=   0.000E+00
   --- H2O           0 4.0000E-04 4.0000E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-24 4.0000E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-52 5.8190E-52 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E-01 2.2310E-01-4.3832E-05 | Line Search reduced step size from -4.3833e-05 to -4.38315e-05
   --- H2            1 2.1295E-12 2.6323E-11 2.4193E-11 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7769E+00 4.3832E-05 |
   --- Na+           c 2.2314E-01 2.2310E-01-4.3832E-05 |
   --- H2O(L)        c 1.9996E+00 1.9996E+00 1.0827E-07 |
   --- OH-           c 1.2035E-07 1.2035E-08-1.0831E-07 |
   --- OH            c 4.2590E-12 5.2645E-11 4.8386E-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.2127739151271E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127739158769E+03
   --- subroutine FORCE: Beginning Slope = -9.26806e-07
   --- subroutine FORCE: End Slope       = -3.96482e-07
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   NaCl(S)           4.776856E+00  4.776900E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231435E-01  2.230997E-01 -1.021390E+02 -1.021393E+02
   ---   H2O(L)            1.999600E+00  1.999600E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.203483E-07  1.203483E-08 -1.046992E+02 -1.070018E+02
   ---   OH                4.258968E-12  5.264505E-11 -1.190554E+02 -1.165409E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.203483E-07  1.203483E-08 -1.104330E+01 -1.334631E+01  8.260791E+00  3.655096E+00
   ---   H2O           0   4.000000E-04  4.000000E-04 -1.294722E+02 -1.294722E+02 -5.468858E+00 -5.468938E+00
   ---   NaCl          0   4.000004E-24  4.000004E-24 -1.560721E+02 -1.560721E+02  1.844468E+01  1.844468E+01
   ---   O2            0   5.819042E-52  5.819042E-52 -1.440334E+02 -1.440334E+02  8.418172E+01  7.412370E+01
   ---   Cl-           1   2.231435E-01  2.230997E-01 -7.237711E+01 -7.237747E+01  7.124068E-04 -1.295812E-06
   ---   H2            1   2.129484E-12  2.632253E-11 -4.397876E+01 -4.146421E+01 -3.408297E+01 -2.653949E+01
   ---                                                        Norms of Delta G():  3.804945E+01  3.312135E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4457994E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7769003E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127739151271E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127739158769E+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 =   5 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-08  -1.2035E-08   3.6551E+00 | Delta damped from -2.19942e-08 to -1.20348e-08 due to component 4 (       OH-) going neg
   --- H2O             4.0000E-04   2.1876E-03  -5.4689E+00 | Normal calc: diag adjusted from 2499.84 to 2499.95 due to act coeff
   --- NaCl            4.0000E-24   0.0000E+00   1.8445E+01 | Skipped: IC =   0 and DG >0:   1.844E+01
   --- O2              5.8190E-52   0.0000E+00   7.4124E+01 | Skipped: IC =   0 and DG >0:   7.412E+01
   --- Cl-             2.2310E-01   7.9724E-08  -1.2958E-06 | Normal calc: diag adjusted from 7.32915 to 16.2538 due to act coeff
   --- H2              2.6323E-11   2.3286E-10  -2.6539E+01 | Normal calc: diag adjusted from 1.13971e+11 to 1.13971e+11 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.2035E-08 3.4381E-09-8.5968E-09 | Line Search reduced step size from -1.08314e-08 to -8.59678e-09
   --- H2O           0 4.0000E-04 4.0000E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-24 4.0000E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-52 5.8190E-52 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2310E-01 2.2310E-01 7.9724E-08 | Normal Major Calc
   --- H2            1 2.6323E-11 2.5918E-10 2.3286E-10 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7769E+00-7.9724E-08 |
   --- Na+           c 2.2310E-01 2.2310E-01 7.9724E-08 |
   --- H2O(L)        c 1.9996E+00 1.9996E+00 8.1311E-09 |
   --- OH-           c 1.2035E-08 3.4381E-09-8.5968E-09 |
   --- OH            c 5.2645E-11 5.1837E-10 4.6572E-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.2127739158769E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127739159049E+03
   --- subroutine FORCE: Beginning Slope = -3.76022e-08
   --- subroutine FORCE: End Slope       = -1.44626e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   NaCl(S)           4.776900E+00  4.776900E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.230997E-01  2.230998E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.999600E+00  1.999600E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.203483E-08  3.438056E-09 -1.070018E+02 -1.082547E+02
   ---   OH                5.264505E-11  5.183693E-10 -1.165409E+02 -1.142538E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   1.203483E-08  3.438056E-09 -1.334631E+01 -1.459921E+01  3.655096E+00  1.149299E+00
   ---   H2O           0   4.000000E-04  4.000000E-04 -1.294722E+02 -1.294722E+02 -5.468938E+00 -5.468938E+00
   ---   NaCl          0   4.000004E-24  4.000004E-24 -1.560721E+02 -1.560721E+02  1.844468E+01  1.844468E+01
   ---   O2            0   5.819042E-52  5.819042E-52 -1.440334E+02 -1.440334E+02  7.412370E+01  6.497524E+01
   ---   Cl-           1   2.230997E-01  2.230998E-01 -7.237747E+01 -7.237746E+01 -1.295812E-06 -1.020104E-07
   ---   H2            1   2.632253E-11  2.591847E-10 -4.146421E+01 -3.917710E+01 -2.653949E+01 -1.967814E+01
   ---                                                        Norms of Delta G():  3.312135E+01  2.881101E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4457996E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7769002E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127739158769E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127739159049E+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

==============================================================================================================
 Iteration =   8, Iterations since last evaluation of optimal basis =   6 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              3.4381E-09  -1.9757E-09   1.1493E+00 | Normal calc: diag adjusted from 5.81724e+08 to 5.81724e+08 due to act coeff
   --- H2O             4.0000E-04   2.1876E-03  -5.4689E+00 | Normal calc: diag adjusted from 2499.84 to 2499.95 due to act coeff
   --- NaCl            4.0000E-24   0.0000E+00   1.8445E+01 | Skipped: IC =   0 and DG >0:   1.844E+01
   --- O2              5.8190E-52   0.0000E+00   6.4975E+01 | Skipped: IC =   0 and DG >0:   6.498E+01
   --- Cl-             2.2310E-01   6.2761E-09  -1.0201E-07 | Normal calc: diag adjusted from 7.32914 to 16.2538 due to act coeff
   --- H2              2.5918E-10   1.7001E-09  -1.9678E+01 | Normal calc: diag adjusted from 1.15748e+10 to 1.15748e+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 3.4381E-09 2.1099E-09-1.3281E-09 | Line Search reduced step size from -1.97568e-09 to -1.32811e-09
   --- H2O           0 4.0000E-04 4.0000E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 4.0000E-24 4.0000E-24 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-52 5.8190E-52 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2310E-01 2.2310E-01 6.2761E-09 | Normal Major Calc
   --- H2            1 2.5918E-10 1.9593E-09 1.7001E-09 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7769E+00-6.2761E-09 |
   --- Na+           c 2.2310E-01 2.2310E-01 6.2761E-09 |
   --- H2O(L)        c 1.9996E+00 1.9996E+00-2.0721E-09 |
   --- OH-           c 3.4381E-09 2.1099E-09-1.3281E-09 |
   --- OH            c 5.1837E-10 3.9186E-09 3.4002E-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.2127739159049E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127739159325E+03
   --- subroutine FORCE: Beginning Slope = -3.4981e-08
   --- subroutine FORCE: End Slope       = -2.33674e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.304532E+01
   ---   NaCl(S)           4.776900E+00  4.776900E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.230998E-01  2.230998E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.999600E+00  1.999600E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               3.438056E-09  2.109948E-09 -1.082547E+02 -1.087430E+02
   ---   OH                5.183693E-10  3.918551E-09 -1.142538E+02 -1.122310E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082001E+02
   ---   H+            1   3.438056E-09  2.109948E-09 -1.459921E+01 -1.508745E+01  1.149299E+00  1.728134E-01
   ---   H2O           0   4.000000E-04  4.000000E-04 -1.294722E+02 -1.294722E+02 -5.468938E+00 -5.468938E+00
   ---   NaCl          0   4.000004E-24  4.000004E-24 -1.560721E+02 -1.560721E+02  1.844468E+01  1.844468E+01
   ---   O2            0   5.819042E-52  5.819042E-52 -1.440334E+02 -1.440334E+02  6.497524E+01  5.688408E+01
   ---   Cl-           1   2.230998E-01  2.230998E-01 -7.237746E+01 -7.237746E+01 -1.020104E-07 -9.723792E-09
   ---   H2            1   2.591847E-10  1.959275E-09 -3.917710E+01 -3.715431E+01 -1.967814E+01 -1.360977E+01
   ---                                                        Norms of Delta G():  2.881101E+01  2.513686E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4457996E+00
   ---                  air =   4.0004000E+00
   ---              NaCl(S) =   4.7769002E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127739159049E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2127739159325E+03
   --- Increment counter increased, step is accepted:    9
   --- 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 =   9, Iterations since last evaluation of optimal basis =   7 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              2.1099E-09  -1.8231E-10   1.7281E-01 | Normal calc: diag adjusted from 9.47891e+08 to 9.47891e+08 due to act coeff
   --- H2O             4.0000E-04   2.1876E-03  -5.4689E+00 | Normal calc: diag adjusted from 2499.84 to 2499.95 due to act coeff
   --- NaCl            4.0000E-24   0.0000E+00   1.8445E+01 | Skipped: IC =   0 and DG >0:   1.844E+01
   --- O2              5.8190E-52   0.0000E+00   5.6884E+01 | Skipped: IC =   0 and DG >0:   5.688E+01
   --- Cl-             2.2310E-01   5.9825E-10  -9.7238E-09 | Normal calc: diag adjusted from 7.32914 to 16.2538 due to act coeff
   --- H2              1.9593E-09   8.8884E-09  -1.3610E+01 | Normal calc: diag adjusted from 1.53118e+09 to 1.53118e+09 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     KMoles  Tent_KMoles Rxn_Adj   |    Comment 
   --- H+            1 2.1099E-09 1.9356E-09-1.7432E-10 | Line Search reduced step size from -1.82314e-10 to -1.74319e-10
   --- H2O           0 4.0000E-04 4.0000E-02 3.9600E-02 | minor species alternative calc
   --- NaCl          0 4.0000E-24 3.9052E-32-4.0000E-24 | minor species alternative calc
   --- O2            0 5.8190E-52 5.8190E-62-5.8190E-52 | minor species alternative calc
   --- Cl-           1 2.2310E-01 2.2310E-01 5.9825E-10 | Normal Major Calc
   --- H2            1 1.9593E-09 1.0848E-08 8.8884E-09 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7769E+00-5.9825E-10 |
   --- Na+           c 2.2310E-01 2.2310E-01 5.9825E-10 |
   --- H2O(L)        c 1.9996E+00 1.9600E+00-3.9600E-02 |
   --- OH-           c 2.1099E-09 1.9356E-09-1.7432E-10 |
   --- OH            c 3.9186E-09 2.1695E-08 1.7777E-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.2127739159325E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128459130107E+03
   --- subroutine FORCE: Beginning Slope = -0.21657
   --- subroutine FORCE: End Slope       = -0.0342694
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.304532E+01 -2.305517E+01
   ---   NaCl(S)           4.776900E+00  4.776900E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.230998E-01  2.230998E-01 -1.021393E+02 -1.021028E+02
   ---   H2O(L)            1.999600E+00  1.960000E+00 -1.240033E+02 -1.240115E+02
   ---   OH-               2.109948E-09  1.935630E-09 -1.087430E+02 -1.088073E+02
   ---   OH                3.918551E-09  2.169542E-08 -1.122310E+02 -1.105295E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082001E+02 -5.082100E+02
   ---   H+            1   2.109948E-09  1.935630E-09 -1.508745E+01 -1.510984E+01  1.728134E-01  9.438799E-02
   ---   H2O           0   4.000000E-04  4.000000E-02 -1.294722E+02 -1.248769E+02 -5.468938E+00 -8.653882E-01
   ---   NaCl          0   4.000004E-24  3.905164E-32 -1.560721E+02 -1.745266E+02  1.844468E+01 -9.850335E-03
   ---   O2            0   5.819042E-52  5.819040E-62 -1.440334E+02 -1.670691E+02  5.688408E+01  2.702580E+01
   ---   Cl-           1   2.230998E-01  2.230998E-01 -7.237746E+01 -7.234095E+01 -9.723792E-09  7.302295E-02
   ---   H2            1   1.959275E-09  1.084771E-08 -3.715431E+01 -3.545278E+01 -1.360977E+01 -8.488729E+00
   ---                                                        Norms of Delta G():  2.513686E+01  1.157019E+01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4061995E+00
   ---                  air =   4.0400000E+00
   ---              NaCl(S) =   4.7769002E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127739159325E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128459130107E+03
   --- Increment counter increased, step is accepted:   10
   --- 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 =  10, Iterations since last evaluation of optimal basis =   8 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9356E-09  -9.1350E-11   9.4388E-02 | Normal calc: diag adjusted from 1.03326e+09 to 1.03326e+09 due to act coeff
   --- H2O             4.0000E-02   3.4663E-02  -8.6539E-01 | Normal calc: diag adjusted from 24.8471 to 24.9659 due to act coeff
   --- NaCl            3.9052E-32   3.8467E-34  -9.8503E-03 | Normal Calc
   --- O2              5.8190E-62   0.0000E+00   2.7026E+01 | Skipped: IC =   0 and DG >0:   2.703E+01
   --- Cl-             2.2310E-01  -4.4327E-03   7.3023E-02 | Normal calc: diag adjusted from 7.30223 to 16.4738 due to act coeff
   --- H2              1.0848E-08   3.0694E-08  -8.4887E+00 | Normal calc: diag adjusted from 2.76556e+08 to 2.76556e+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.9356E-09 1.8465E-09-8.9177E-11 | Line Search reduced step size from -9.13501e-11 to -8.91771e-11
   --- H2O           1 4.0000E-02 7.4663E-02 3.4663E-02 | Normal Major Calc
   --- NaCl          0 3.9052E-32 3.9052E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-62 5.8190E-62 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2310E-01 2.1868E-01-4.4182E-03 | Line Search reduced step size from -0.00443267 to -0.00441821
   --- H2            1 1.0848E-08 4.1542E-08 3.0694E-08 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+00 4.7813E+00 4.4182E-03 |
   --- Na+           c 2.2310E-01 2.1868E-01-4.4182E-03 |
   --- H2O(L)        c 1.9600E+00 1.9253E+00-3.4663E-02 |
   --- OH-           c 1.9356E-09 1.8465E-09-8.9177E-11 |
   --- OH            c 2.1695E-08 8.3084E-08 6.1389E-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.2128459130107E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128644445665E+03
   --- subroutine FORCE: Beginning Slope = -0.0303197
   --- subroutine FORCE: End Slope       = -0.00897876
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.305517E+01 -2.306372E+01
   ---   NaCl(S)           4.776900E+00  4.781318E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.230998E-01  2.186816E-01 -1.021028E+02 -1.021068E+02
   ---   H2O(L)            1.960000E+00  1.925337E+00 -1.240115E+02 -1.240106E+02
   ---   OH-               1.935630E-09  1.846453E-09 -1.088073E+02 -1.088368E+02
   ---   OH                2.169542E-08  8.308425E-08 -1.105295E+02 -1.091952E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082100E+02 -5.082185E+02
   ---   H+            1   1.935630E-09  1.846453E-09 -1.510984E+01 -1.514393E+01  9.438799E-02  2.984791E-02
   ---   H2O           1   4.000000E-02  7.466279E-02 -1.248769E+02 -1.242613E+02 -8.653882E-01 -2.507309E-01
   ---   NaCl          0   3.905164E-32  3.905164E-32 -1.745266E+02 -1.745352E+02 -9.850335E-03 -1.839366E-02
   ---   O2            0   5.819040E-62  5.819040E-62 -1.670691E+02 -1.670776E+02  2.702580E+01  2.168222E+01
   ---   Cl-           1   2.230998E-01  2.186816E-01 -7.234095E+01 -7.234492E+01  7.302295E-02  6.509339E-02
   ---   H2            1   1.084771E-08  4.154213E-08 -3.545278E+01 -3.411857E+01 -8.488729E+00 -4.487900E+00
   ---                                                        Norms of Delta G():  1.157019E+01  9.039988E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3627003E+00
   ---                  air =   4.0746629E+00
   ---              NaCl(S) =   4.7813184E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128459130107E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128644445665E+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 =   9 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8465E-09  -2.7556E-11   2.9848E-02 | Normal calc: diag adjusted from 1.08316e+09 to 1.08316e+09 due to act coeff
   --- H2O             7.4663E-02   1.8761E-02  -2.5073E-01 | Normal calc: diag adjusted from 13.2443 to 13.3646 due to act coeff
   --- NaCl            3.9052E-32   7.1830E-34  -1.8394E-02 | Normal Calc
   --- O2              5.8190E-62   0.0000E+00   2.1682E+01 | Skipped: IC =   0 and DG >0:   2.168E+01
   --- Cl-             2.1868E-01  -3.8787E-03   6.5093E-02 | Normal calc: diag adjusted from 7.45274 to 16.7821 due to act coeff
   --- H2              4.1542E-08   6.2146E-08  -4.4879E+00 | Normal calc: diag adjusted from 7.22158e+07 to 7.22159e+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.8465E-09 1.8191E-09-2.7350E-11 | Line Search reduced step size from -2.75564e-11 to -2.73502e-11
   --- H2O           1 7.4663E-02 9.3424E-02 1.8761E-02 | Normal Major Calc
   --- NaCl          0 3.9052E-32 3.9052E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-62 5.8190E-62 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1868E-01 2.1481E-01-3.8674E-03 | Line Search reduced step size from -0.00387873 to -0.00386742
   --- H2            1 4.1542E-08 1.0369E-07 6.2146E-08 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7813E+00 4.7852E+00 3.8674E-03 |
   --- Na+           c 2.1868E-01 2.1481E-01-3.8674E-03 |
   --- H2O(L)        c 1.9253E+00 1.9066E+00-1.8761E-02 |
   --- OH-           c 1.8465E-09 1.8191E-09-2.7350E-11 |
   --- OH            c 8.3084E-08 2.0738E-07 1.2429E-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.2128644445665E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128672367281E+03
   --- subroutine FORCE: Beginning Slope = -0.00495593
   --- subroutine FORCE: End Slope       = -0.000785067
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.306372E+01 -2.306831E+01
   ---   NaCl(S)           4.781318E+00  4.785186E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.186816E-01  2.148141E-01 -1.021068E+02 -1.021215E+02
   ---   H2O(L)            1.925337E+00  1.906576E+00 -1.240106E+02 -1.240072E+02
   ---   OH-               1.846453E-09  1.819102E-09 -1.088368E+02 -1.088427E+02
   ---   OH                8.308425E-08  2.073755E-07 -1.091952E+02 -1.082852E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082185E+02 -5.082231E+02
   ---   H+            1   1.846453E-09  1.819102E-09 -1.514393E+01 -1.516678E+01  2.984791E-02 -2.260815E-03
   ---   H2O           1   7.466279E-02  9.342356E-02 -1.242613E+02 -1.240417E+02 -2.507309E-01 -3.449644E-02
   ---   NaCl          0   3.905164E-32  3.905164E-32 -1.745352E+02 -1.745398E+02 -1.839366E-02 -2.298739E-02
   ---   O2            0   5.819040E-62  5.819040E-62 -1.670776E+02 -1.670822E+02  2.168222E+01  1.804396E+01
   ---   Cl-           1   2.186816E-01  2.148141E-01 -7.234492E+01 -7.235965E+01  6.509339E-02  3.562514E-02
   ---   H2            1   4.154213E-08  1.036878E-07 -3.411857E+01 -3.320849E+01 -4.487900E+00 -1.764320E+00
   ---                                                        Norms of Delta G():  9.039988E+00  7.401582E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3362045E+00
   ---                  air =   4.0934239E+00
   ---              NaCl(S) =   4.7851859E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128644445665E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128672367281E+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

==============================================================================================================
 Iteration =  12, Iterations since last evaluation of optimal basis =  10 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8191E-09   2.0563E-12  -2.2608E-03 | Normal calc: diag adjusted from 1.09944e+09 to 1.09944e+09 due to act coeff
   --- H2O             9.3424E-02   3.2314E-03  -3.4496E-02 | Normal calc: diag adjusted from 10.5561 to 10.6753 due to act coeff
   --- NaCl            3.9052E-32   8.9770E-34  -2.2987E-02 | Normal Calc
   --- O2              5.8190E-62   0.0000E+00   1.8044E+01 | Skipped: IC =   0 and DG >0:   1.804E+01
   --- Cl-             2.1481E-01  -2.0966E-03   3.5625E-02 | Normal calc: diag adjusted from 7.5982 to 16.9918 due to act coeff
   --- H2              1.0369E-07   6.0979E-08  -1.7643E+00 | Normal calc: diag adjusted from 2.8933e+07 to 2.8933e+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.8191E-09 1.8212E-09 2.0563E-12 | Normal Major Calc
   --- H2O           1 9.3424E-02 9.6655E-02 3.2314E-03 | Normal Major Calc
   --- NaCl          0 3.9052E-32 3.9052E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.8190E-62 5.8190E-62 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1481E-01 2.1272E-01-2.0932E-03 | Line Search reduced step size from -0.0020966 to -0.00209322
   --- H2            1 1.0369E-07 1.6467E-07 6.0979E-08 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7852E+00 4.7873E+00 2.0932E-03 |
   --- Na+           c 2.1481E-01 2.1272E-01-2.0932E-03 |
   --- H2O(L)        c 1.9066E+00 1.9033E+00-3.2315E-03 |
   --- OH-           c 1.8191E-09 1.8212E-09 2.0563E-12 |
   --- OH            c 2.0738E-07 3.2933E-07 1.2196E-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.2128672367281E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673433895E+03
   --- subroutine FORCE: Beginning Slope = -0.000186151
   --- subroutine FORCE: End Slope       = -2.77347e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306831E+01 -2.306910E+01
   ---   NaCl(S)           4.785186E+00  4.787279E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.148141E-01  2.127209E-01 -1.021215E+02 -1.021362E+02
   ---   H2O(L)            1.906576E+00  1.903345E+00 -1.240072E+02 -1.240039E+02
   ---   OH-               1.819102E-09  1.821159E-09 -1.088427E+02 -1.088407E+02
   ---   OH                2.073755E-07  3.293344E-07 -1.082852E+02 -1.078234E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082231E+02 -5.082239E+02
   ---   H+            1   1.819102E-09  1.821159E-09 -1.516678E+01 -1.518162E+01 -2.260815E-03 -1.836595E-02
   ---   H2O           1   9.342356E-02  9.665497E-02 -1.240417E+02 -1.240085E+02 -3.449644E-02 -4.588517E-03
   ---   NaCl          0   3.905164E-32  3.905164E-32 -1.745398E+02 -1.745406E+02 -2.298739E-02 -2.377653E-02
   ---   O2            0   5.819040E-62  5.819040E-62 -1.670822E+02 -1.670830E+02  1.804396E+01  1.620278E+01
   ---   Cl-           1   2.148141E-01  2.127209E-01 -7.235965E+01 -7.237439E+01  3.562514E-02  6.155026E-03
   ---   H2            1   1.036878E-07  1.646672E-07 -3.320849E+01 -3.274673E+01 -1.764320E+00 -3.856738E-01
   ---                                                        Norms of Delta G():  7.401582E+00  6.616641E+00
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3287865E+00
   ---                  air =   4.0966555E+00
   ---              NaCl(S) =   4.7872791E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128672367281E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673433895E+03
   --- Increment counter increased, step is accepted:   13
   --- 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 =  13, Iterations since last evaluation of optimal basis =  11 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8212E-09   1.6724E-11  -1.8366E-02 | Normal calc: diag adjusted from 1.0982e+09 to 1.0982e+09 due to act coeff
   --- H2O             9.6655E-02   4.4483E-04  -4.5885E-03 | Normal calc: diag adjusted from 10.198 to 10.3151 due to act coeff
   --- NaCl            3.9052E-32   9.2851E-34  -2.3777E-02 | Normal Calc
   --- O2              5.8190E-62   0.0000E+00   1.6203E+01 | Skipped: IC =   0 and DG >0:   1.620E+01
   --- Cl-             2.1272E-01  -3.6066E-04   6.1550E-03 | Normal calc: diag adjusted from 7.68436 to 17.0662 due to act coeff
   --- H2              1.6467E-07   2.1169E-08  -3.8567E-01 | Normal calc: diag adjusted from 1.82186e+07 to 1.82186e+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.8212E-09 1.8379E-09 1.6724E-11 | Normal Major Calc
   --- H2O           1 9.6655E-02 9.7100E-02 4.4483E-04 | Normal Major Calc
   --- NaCl          0 3.9052E-32 3.9991E-32 9.3964E-34 | minor species alternative calc
   --- O2            0 5.8190E-62 5.3466E-69-5.8190E-62 | minor species alternative calc
   --- Cl-           1 2.1272E-01 2.1236E-01-3.6055E-04 | Line Search reduced step size from -0.000360656 to -0.000360555
   --- H2            1 1.6467E-07 1.8584E-07 2.1169E-08 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7873E+00 4.7876E+00 3.6055E-04 |
   --- Na+           c 2.1272E-01 2.1236E-01-3.6055E-04 |
   --- H2O(L)        c 1.9033E+00 1.9029E+00-4.4488E-04 |
   --- OH-           c 1.8212E-09 1.8379E-09 1.6724E-11 |
   --- OH            c 3.2933E-07 3.7167E-07 4.2339E-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.2128673433895E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458315E+03
   --- subroutine FORCE: Beginning Slope = -4.26851e-06
   --- subroutine FORCE: End Slope       = -6.16713e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306910E+01 -2.306921E+01
   ---   NaCl(S)           4.787279E+00  4.787640E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.127209E-01  2.123604E-01 -1.021362E+02 -1.021389E+02
   ---   H2O(L)            1.903345E+00  1.902900E+00 -1.240039E+02 -1.240033E+02
   ---   OH-               1.821159E-09  1.837882E-09 -1.088407E+02 -1.088314E+02
   ---   OH                3.293344E-07  3.716730E-07 -1.078234E+02 -1.077026E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082239E+02 -5.082240E+02
   ---   H+            1   1.821159E-09  1.837882E-09 -1.518162E+01 -1.517543E+01 -1.836595E-02 -3.524920E-03
   ---   H2O           1   9.665497E-02  9.709980E-02 -1.240085E+02 -1.240040E+02 -4.588517E-03 -6.981099E-04
   ---   NaCl          0   3.905164E-32  3.999128E-32 -1.745406E+02 -1.745169E+02 -2.377653E-02 -1.085940E-04
   ---   O2            0   5.819040E-62  5.346571E-69 -1.670830E+02 -1.832859E+02  1.620278E+01 -4.822514E-01
   ---   Cl-           1   2.127209E-01  2.123604E-01 -7.237439E+01 -7.237704E+01  6.155026E-03  8.477354E-04
   ---   H2            1   1.646672E-07  1.858365E-07 -3.274673E+01 -3.262590E+01 -3.856738E-01 -2.436303E-02
   ---                                                        Norms of Delta G():  6.616641E+00  1.971352E-01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3276206E+00
   ---                  air =   4.0971004E+00
   ---              NaCl(S) =   4.7876396E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673433895E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458315E+03
   --- Increment counter increased, step is accepted:   14
   --- 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 =  14, Iterations since last evaluation of optimal basis =  12 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8379E-09   3.2392E-12  -3.5249E-03 | Normal calc: diag adjusted from 1.08821e+09 to 1.08821e+09 due to act coeff
   --- H2O             9.7100E-02   6.7993E-05  -6.9811E-04 | Normal calc: diag adjusted from 10.1505 to 10.2673 due to act coeff
   --- NaCl            3.9991E-32   4.3428E-36  -1.0859E-04 | Normal Calc
   --- O2              5.3466E-69   2.5784E-69  -4.8225E-01 | Normal Calc
   --- Cl-             2.1236E-01  -4.9638E-05   8.4774E-04 | Normal calc: diag adjusted from 7.69946 to 17.0785 due to act coeff
   --- H2              1.8584E-07   1.5092E-09  -2.4363E-02 | Normal calc: diag adjusted from 1.61432e+07 to 1.61432e+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.8379E-09 1.8411E-09 3.2392E-12 | Normal Major Calc
   --- H2O           1 9.7100E-02 9.7168E-02 6.7993E-05 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.3466E-69 5.3466E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1236E-01 2.1231E-01-4.9636E-05 | Line Search reduced step size from -4.96377e-05 to -4.96358e-05
   --- H2            1 1.8584E-07 1.8735E-07 1.5092E-09 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7876E+00 4.7877E+00 4.9636E-05 |
   --- Na+           c 2.1236E-01 2.1231E-01-4.9636E-05 |
   --- H2O(L)        c 1.9029E+00 1.9028E+00-6.7997E-05 |
   --- OH-           c 1.8379E-09 1.8411E-09 3.2392E-12 |
   --- OH            c 3.7167E-07 3.7469E-07 3.0184E-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.2128673458315E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458827E+03
   --- subroutine FORCE: Beginning Slope = -8.95817e-08
   --- subroutine FORCE: End Slope       = -1.28811e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306921E+01 -2.306923E+01
   ---   NaCl(S)           4.787640E+00  4.787689E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123604E-01  2.123107E-01 -1.021389E+02 -1.021393E+02
   ---   H2O(L)            1.902900E+00  1.902832E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.837882E-09  1.841122E-09 -1.088314E+02 -1.088297E+02
   ---   OH                3.716730E-07  3.746914E-07 -1.077026E+02 -1.076945E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.837882E-09  1.841122E-09 -1.517543E+01 -1.517406E+01 -3.524920E-03 -4.603026E-04
   ---   H2O           1   9.709980E-02  9.716780E-02 -1.240040E+02 -1.240034E+02 -6.981099E-04 -9.484509E-05
   ---   NaCl          0   3.999128E-32  3.999128E-32 -1.745169E+02 -1.745169E+02 -1.085940E-04 -1.251905E-04
   ---   O2            0   5.346571E-69  5.346571E-69 -1.832859E+02 -1.832859E+02 -4.822514E-01 -5.143942E-01
   ---   Cl-           1   2.123604E-01  2.123107E-01 -7.237704E+01 -7.237740E+01  8.477354E-04  1.295801E-04
   ---   H2            1   1.858365E-07  1.873457E-07 -3.262590E+01 -3.261783E+01 -2.436303E-02 -3.084553E-04
   ---                                                        Norms of Delta G():  1.971352E-01  2.100007E-01
   ---           Phase_Name    KMoles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274533E+00
   ---                  air =   4.0971684E+00
   ---              NaCl(S) =   4.7876893E+00
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458315E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458827E+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 =  13 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8411E-09   4.2374E-13  -4.6030E-04 | Normal calc: diag adjusted from 1.08629e+09 to 1.08629e+09 due to act coeff
   --- H2O             9.7168E-02   9.2441E-06  -9.4845E-05 | Normal calc: diag adjusted from 10.1433 to 10.26 due to act coeff
   --- NaCl            3.9991E-32   5.0065E-36  -1.2519E-04 | Normal Calc
   --- O2              5.3466E-69   2.7502E-69  -5.1439E-01 | Normal Calc
   --- Cl-             2.1231E-01  -7.5866E-06   1.2958E-04 | Normal calc: diag adjusted from 7.70154 to 17.0802 due to act coeff
   --- H2              1.8735E-07   1.9263E-11  -3.0846E-04 | Normal calc: diag adjusted from 1.60132e+07 to 1.60132e+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.8411E-09 1.8415E-09 4.2374E-13 | Normal Major Calc
   --- H2O           1 9.7168E-02 9.7177E-02 9.2441E-06 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.3466E-69 5.3466E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1231E-01 2.1230E-01-7.5865E-06 | Line Search reduced step size from -7.58657e-06 to -7.58652e-06
   --- H2            1 1.8735E-07 1.8736E-07 1.9263E-11 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 7.5865E-06 |
   --- Na+           c 2.1231E-01 2.1230E-01-7.5865E-06 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-9.2442E-06 |
   --- OH-           c 1.8411E-09 1.8415E-09 4.2374E-13 |
   --- OH            c 3.7469E-07 3.7473E-07 3.8525E-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.2128673458827E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458837E+03
   --- subroutine FORCE: Beginning Slope = -1.85983e-09
   --- subroutine FORCE: End Slope       = -2.67342e-10
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787689E+00  4.787697E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123107E-01  2.123031E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902832E+00  1.902823E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841122E-09  1.841545E-09 -1.088297E+02 -1.088294E+02
   ---   OH                3.746914E-07  3.747299E-07 -1.076945E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841122E-09  1.841545E-09 -1.517406E+01 -1.517389E+01 -4.603026E-04 -7.275328E-05
   ---   H2O           1   9.716780E-02  9.717704E-02 -1.240034E+02 -1.240033E+02 -9.484509E-05 -1.446242E-05
   ---   NaCl          0   3.999128E-32  3.999128E-32 -1.745169E+02 -1.745169E+02 -1.251905E-04 -1.274467E-04
   ---   O2            0   5.346571E-69  5.346571E-69 -1.832859E+02 -1.832859E+02 -5.143942E-01 -5.147737E-01
   ---   Cl-           1   2.123107E-01  2.123031E-01 -7.237740E+01 -7.237746E+01  1.295801E-04  1.761666E-05
   ---   H2            1   1.873457E-07  1.873649E-07 -3.261783E+01 -3.261773E+01 -3.084553E-04 -3.176905E-05
   ---                                                        Norms of Delta G():  2.100007E-01  2.101555E-01
   ---           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.2128673458827E+03
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458837E+03
   --- Increment counter increased, step is accepted:   16
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  16, Iterations since last evaluation of optimal basis =  14 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8415E-09   6.6989E-14  -7.2753E-05 | Normal calc: diag adjusted from 1.08604e+09 to 1.08604e+09 due to act coeff
   --- H2O             9.7177E-02   1.4097E-06  -1.4462E-05 | Normal calc: diag adjusted from 10.1423 to 10.2591 due to act coeff
   --- NaCl            3.9991E-32   5.0968E-36  -1.2745E-04 | Normal Calc
   --- O2              5.3466E-69   2.7523E-69  -5.1477E-01 | Normal Calc
   --- Cl-             2.1230E-01  -1.0314E-06   1.7617E-05 | Normal calc: diag adjusted from 7.70186 to 17.0805 due to act coeff
   --- H2              1.8736E-07   1.9841E-12  -3.1769E-05 | Normal calc: diag adjusted from 1.60115e+07 to 1.60115e+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.6989E-14 | Normal Major Calc
   --- H2O           1 9.7177E-02 9.7178E-02 1.4097E-06 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9991E-32 0.0000E+00 | minor species not considered
   --- O2            0 5.3466E-69 5.3466E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E-01 2.1230E-01-1.0314E-06 | Line Search reduced step size from -1.03139e-06 to -1.03139e-06
   --- H2            1 1.8736E-07 1.8737E-07 1.9841E-12 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 1.0314E-06 |
   --- Na+           c 2.1230E-01 2.1230E-01-1.0314E-06 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-1.4097E-06 |
   --- OH-           c 1.8415E-09 1.8416E-09 6.6989E-14 |
   --- OH            c 3.7473E-07 3.7473E-07 3.9683E-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.2128673458837E+03
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+03
   --- subroutine FORCE: Beginning Slope = -3.85578e-11
   --- subroutine FORCE: End Slope       = -5.54189e-12
   --- 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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787697E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123031E-01  2.123021E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902823E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841545E-09  1.841612E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747299E-07  3.747339E-07 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.082240E+02 -5.082240E+02
   ---   H+            1   1.841545E-09  1.841612E-09 -1.517389E+01 -1.517386E+01 -7.275328E-05 -9.509607E-06
   ---   H2O           1   9.717704E-02  9.717845E-02 -1.240033E+02 -1.240033E+02 -1.446242E-05 -1.965642E-06
   ---   NaCl          0   3.999128E-32  3.999128E-32 -1.745169E+02 -1.745169E+02 -1.274467E-04 -1.277908E-04
   ---   O2            0   5.346571E-69  5.346571E-69 -1.832859E+02 -1.832859E+02 -5.147737E-01 -5.148117E-01
   ---   Cl-           1   2.123031E-01  2.123021E-01 -7.237746E+01 -7.237746E+01  1.761666E-05  2.686530E-06
   ---   H2            1   1.873649E-07  1.873669E-07 -3.261773E+01 -3.261772E+01 -3.176905E-05 -4.363969E-06
   ---                                                        Norms of Delta G():  2.101555E-01  2.101710E-01
   ---           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.2128673458837E+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
   --- Equilibrium check for major species: H+ failed
   --- Subroutine vcs_dfe called for components and minors
   --- Subroutine vcs_deltag called for minor noncomponents

==============================================================================================================
 Iteration =  17, Iterations since last evaluation of optimal basis =  15 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   8.7565E-15  -9.5096E-06 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7178E-02   1.9160E-07  -1.9656E-06 | Normal calc: diag adjusted from 10.1422 to 10.2589 due to act coeff
   --- NaCl            3.9991E-32   5.1105E-36  -1.2779E-04 | Normal Calc
   --- O2              5.3466E-69   2.7525E-69  -5.1481E-01 | Normal Calc
   --- Cl-             2.1230E-01  -1.5729E-07   2.6865E-06 | Normal calc: diag adjusted from 7.7019 to 17.0805 due to act coeff
   --- H2              1.8737E-07   2.7255E-13  -4.3640E-06 | Normal calc: diag adjusted from 1.60114e+07 to 1.60114e+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.7565E-15 | Normal Major Calc
   --- H2O           1 9.7178E-02 9.7179E-02 1.9160E-07 | Normal Major Calc
   --- NaCl          0 3.9991E-32 3.9996E-32 5.1108E-36 | minor species alternative calc
   --- O2            0 5.3466E-69 8.9465E-69 3.6000E-69 | minor species alternative calc
   --- Cl-           1 2.1230E-01 2.1230E-01-1.5729E-07 | Line Search reduced step size from -1.57286e-07 to -1.57286e-07
   --- H2            1 1.8737E-07 1.8737E-07 2.7255E-13 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 1.5729E-07 |
   --- Na+           c 2.1230E-01 2.1230E-01-1.5729E-07 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-1.9160E-07 |
   --- OH-           c 1.8416E-09 1.8416E-09 8.7565E-15 |
   --- OH            c 3.7473E-07 3.7473E-07 5.4511E-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 = -7.9918e-13
   --- subroutine FORCE: End Slope       = -1.14864e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123021E-01  2.123020E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841612E-09  1.841621E-09 -1.088294E+02 -1.088294E+02
   ---   OH                3.747339E-07  3.747344E-07 -1.076944E+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.509607E-06 -1.507395E-06
   ---   H2O           1   9.717845E-02  9.717864E-02 -1.240033E+02 -1.240033E+02 -1.965642E-06 -2.997442E-07
   ---   NaCl          0   3.999128E-32  3.999639E-32 -1.745169E+02 -1.745168E+02 -1.277908E-04 -4.676497E-08
   ---   O2            0   5.346571E-69  8.946542E-69 -1.832859E+02 -1.827711E+02 -5.148117E-01 -5.160322E-06
   ---   Cl-           1   2.123021E-01  2.123020E-01 -7.237746E+01 -7.237746E+01  2.686530E-06  3.651417E-07
   ---   H2            1   1.873669E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -4.363969E-06 -6.583027E-07
   ---                                                        Norms of Delta G():  2.101710E-01  2.219605E-06
   ---           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:   18
   --- 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 =  18, Iterations since last evaluation of optimal basis =  16 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   1.3880E-15  -1.5074E-06 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   2.9218E-08  -2.9974E-07 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   1.8704E-39  -4.6765E-08 | Normal Calc
   --- O2              8.9465E-69   4.6167E-74  -5.1603E-06 | Normal Calc
   --- Cl-             2.1230E-01  -2.1378E-08   3.6514E-07 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   4.1115E-14  -6.5830E-07 | 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.3880E-15 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 2.9218E-08 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9465E-69 8.9465E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E-01 2.1230E-01-2.1378E-08 | Line Search reduced step size from -2.13777e-08 to -2.13777e-08
   --- H2            1 1.8737E-07 1.8737E-07 4.1115E-14 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 2.1378E-08 |
   --- Na+           c 2.1230E-01 2.1230E-01-2.1378E-08 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-2.9218E-08 |
   --- OH-           c 1.8416E-09 1.8416E-09 1.3880E-15 |
   --- OH            c 3.7473E-07 3.7473E-07 8.2230E-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 = -1.65639e-14
   --- subroutine FORCE: End Slope       = -2.38067e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123020E-01  2.123019E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   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.507395E-06 -1.970826E-07
   ---   H2O           1   9.717864E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -2.997442E-07 -4.073974E-08
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -4.676497E-08 -5.389626E-08
   ---   O2            0   8.946542E-69  8.946542E-69 -1.827711E+02 -1.827711E+02 -5.160322E-06 -5.947611E-06
   ---   Cl-           1   2.123020E-01  2.123019E-01 -7.237746E+01 -7.237746E+01  3.651417E-07  5.568121E-08
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -6.583027E-07 -9.044825E-08
   ---                                                        Norms of Delta G():  2.219605E-06  2.429978E-06
   ---           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:   19
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  19, Iterations since last evaluation of optimal basis =  17 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   1.8148E-16  -1.9708E-07 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   3.9712E-09  -4.0740E-08 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   2.1557E-39  -5.3896E-08 | Normal Calc
   --- O2              8.9465E-69   5.3211E-74  -5.9476E-06 | Normal Calc
   --- Cl-             2.1230E-01  -3.2599E-09   5.5681E-08 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   5.6490E-15  -9.0448E-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.8148E-16 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 3.9712E-09 | Normal Major Calc
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9465E-69 8.9465E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E-01 2.1230E-01-3.2599E-09 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 5.6490E-15 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 3.2599E-09 |
   --- Na+           c 2.1230E-01 2.1230E-01-3.2599E-09 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-3.9712E-09 |
   --- OH-           c 1.8416E-09 1.8416E-09 1.8148E-16 |
   --- OH            c 3.7473E-07 3.7473E-07 1.1298E-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.43302e-16
   --- subroutine FORCE: End Slope       = -4.93417e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   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.970826E-07 -3.124197E-08
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -4.073974E-08 -6.212474E-09
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -5.389626E-08 -5.486550E-08
   ---   O2            0   8.946542E-69  8.946542E-69 -1.827711E+02 -1.827711E+02 -5.947611E-06 -6.054564E-06
   ---   Cl-           1   2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01  5.568121E-08  7.567920E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -9.044825E-08 -1.364396E-08
   ---                                                        Norms of Delta G():  2.429978E-06  2.471909E-06
   ---           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:   20
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  20, Iterations since last evaluation of optimal basis =  18 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species        KMoles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-09   2.8768E-17  -3.1242E-08 | Normal calc: diag adjusted from 1.086e+09 to 1.086e+09 due to act coeff
   --- H2O             9.7179E-02   6.0557E-10  -6.2125E-09 | Normal calc: diag adjusted from 10.1421 to 10.2589 due to act coeff
   --- NaCl            3.9996E-32   2.1944E-39  -5.4865E-08 | Normal Calc
   --- O2              8.9465E-69   5.4167E-74  -6.0546E-06 | Normal Calc
   --- Cl-             2.1230E-01  -4.4307E-10   7.5679E-09 | Normal calc: diag adjusted from 7.70191 to 17.0805 due to act coeff
   --- H2              1.8737E-07   8.5214E-16  -1.3644E-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.8768E-17 | Normal Major Calc
   --- H2O           1 9.7179E-02 9.7179E-02 6.0557E-10 | Line Search reduced step size from 6.05571e-10 to 6.05569e-10
   --- NaCl          0 3.9996E-32 3.9996E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.9465E-69 8.9465E-69 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E-01 2.1230E-01-4.4307E-10 | Normal Major Calc
   --- H2            1 1.8737E-07 1.8737E-07 8.5214E-16 | Normal Major Calc
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00 4.4307E-10 |
   --- Na+           c 2.1230E-01 2.1230E-01-4.4307E-10 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00-6.0557E-10 |
   --- OH-           c 1.8416E-09 1.8416E-09 2.8768E-17 |
   --- OH            c 3.7473E-07 3.7473E-07 1.7043E-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.11524e-18
   --- subroutine FORCE: End Slope       = -1.02266e-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
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   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.124197E-08 -4.084725E-09
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -6.212474E-09 -8.443806E-10
   ---   NaCl          0   3.999639E-32  3.999639E-32 -1.745168E+02 -1.745168E+02 -5.486550E-08 -5.501332E-08
   ---   O2            0   8.946542E-69  8.946542E-69 -1.827711E+02 -1.827711E+02 -6.054564E-06 -6.070882E-06
   ---   Cl-           1   2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01  7.567920E-09  1.154049E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -1.364396E-08 -1.874611E-09
   ---                                                        Norms of Delta G():  2.471909E-06  2.478530E-06
   ---           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:   21
   --- 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 =  21, Iterations since last evaluation of optimal basis =  19 (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.2003E-39 | minor species alternative calc
   --- O2            0 8.9465E-69 8.9466E-69 5.4314E-74 | minor species alternative calc
   --- Cl-           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
   --- N2            c 4.0000E+00 4.0000E+00 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+00 4.7877E+00-2.2003E-39 |
   --- Na+           c 2.1230E-01 2.1230E-01 0.0000E+00 |
   --- H2O(L)        c 1.9028E+00 1.9028E+00 1.0863E-73 |
   --- OH-           c 1.8416E-09 1.8416E-09 0.0000E+00 |
   --- OH            c 3.7473E-07 3.7473E-07-2.1725E-73 |
   --- 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.21048e-46
   --- subroutine FORCE: End Slope       = 6.25373e-53
   --- subroutine FORCE produced no adjustments (al = 0.999999)
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_KMoles Final_KMoles Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+00  4.000000E+00 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+00  4.787698E+00 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E-01  2.123019E-01 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+00  1.902821E+00 -1.240033E+02 -1.240033E+02
   ---   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.084725E-09 -4.084725E-09
   ---   H2O           1   9.717867E-02  9.717867E-02 -1.240033E+02 -1.240033E+02 -8.443806E-10 -8.443806E-10
   ---   NaCl          0   3.999639E-32  3.999640E-32 -1.745168E+02 -1.745168E+02 -5.501332E-08  2.842171E-14
   ---   O2            0   8.946542E-69  8.946596E-69 -1.827711E+02 -1.827711E+02 -6.070882E-06  0.000000E+00
   ---   Cl-           1   2.123019E-01  2.123019E-01 -7.237746E+01 -7.237746E+01  1.154049E-09  1.154049E-09
   ---   H2            1   1.873672E-07  1.873672E-07 -3.261772E+01 -3.261772E+01 -1.874611E-09 -1.874611E-09
   ---                                                        Norms of Delta G():  2.478530E-06  1.925442E-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:   22
   --- 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 
		Volume      =          100.41 m**3


--------------------------------------------------------------------------------
 Species                 Equilibrium kmoles   Mole Fraction    ChemPot/RT    SpecUnkType
--------------------------------------------------------------------------------
 N2                        4.0000000E+00      9.7628143E-01     -2.3069E+01     0
 NaCl(S)                   4.7876981E+00      1.0000000E+00     -1.7452E+02     0
 Na+                       2.1230193E-01      9.1217525E-02     -1.0214E+02     0
 H2O(L)                    1.9028209E+00      8.1756495E-01     -1.2400E+02     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
 Cl-                       2.1230193E-01      9.1217525E-02     -7.2377E+01  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|         N2    NaCl(S)        Na+     H2O(L)        OH-         OH        CO2 |           |
 NonComponent  |   Moles   |          4       4.79      0.212        1.9   1.84e-09   3.75e-07          0 | DG/RT Rxn |
-------------------------------------------------------------------------------------------------------------------
   7 H+        |  1.84e-09 |       0.00       0.00       0.00      -1.00       1.00       0.00       0.00 | -4.08e-09 |
   8 H2O       |    0.0972 |       0.00       0.00       0.00      -1.00       0.00       0.00       0.00 | -8.44e-10 |
   9 NaCl      |     4e-32 |       0.00      -1.00       0.00       0.00       0.00       0.00       0.00 |  2.84e-14 |
  10 O2        |  8.95e-69 |       0.00       0.00       0.00       2.00       0.00      -4.00       0.00 |         0 |
  11 Cl-       |     0.212 |       0.00      -1.00       1.00       0.00       0.00       0.00       0.00 |  1.15e-09 |
  12 H2        |  1.87e-07 |       0.00       0.00       0.00      -2.00       0.00       2.00       0.00 | -1.87e-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)
-------------------------------------------------------------------------------------------------------------------
 N2            4.0000000E+00 -2.3045225E+01   0.0000000E+00  -2.4004385E-02  | 0.0000000E+00 |   -2.3069E+01 |
 NaCl(S)       4.7876981E+00 -1.7451679E+02   0.0000000E+00   0.0000000E+00  | 0.0000000E+00 |   -1.7452E+02 |
 Na+           2.1230193E-01 -1.0397591E+02   2.1455966E-01  -2.3945082E+00  | 0.0000000E+00 |   -1.0214E+02 | ( 4.0165350E+00)
 H2O(L)        1.9028209E+00 -1.2371551E+02  -8.6317426E-02  -2.0142493E-01  | 0.0000000E+00 |   -1.2400E+02 |
 OH-           1.8416226E-09 -9.1483483E+01  -4.0506365E-01  -2.0957381E+01  |-0.0000000E+00 |   -1.0883E+02 | ( 4.0165350E+00)
 OH            3.7473450E-07 -9.1487045E+01   0.0000000E+00  -1.6207347E+01  | 0.0000000E+00 |   -1.0769E+02 |
 CO2           0.0000000E+00 -1.8445182E+02   0.0000000E+00  -3.2377221E+02  | 0.0000000E+00 |   -5.0822E+02 |
 Cl-           2.1230193E-01 -7.4214051E+01   2.1455966E-01  -2.3945082E+00  |-0.0000000E+00 |   -7.2377E+01 | ( 4.0165350E+00)
 H2O           9.7178674E-02 -1.2026175E+02   0.0000000E+00  -3.7415027E+00  | 0.0000000E+00 |   -1.2400E+02 |
 H2            1.8736725E-07 -1.5717224E+01   0.0000000E+00  -1.6900494E+01  | 0.0000000E+00 |   -3.2618E+01 |
 H+            1.8416226E-09  0.0000000E+00   1.7669883E+00  -2.0957381E+01  | 0.0000000E+00 |   -1.5174E+01 | ( 4.0165350E+00)
 NaCl          3.9996395E-32 -1.0080997E+02   0.0000000E+00  -7.3706817E+01  | 0.0000000E+00 |   -1.7452E+02 |
 O2            8.9465959E-69 -2.4673669E+01   0.0000000E+00  -1.5809740E+02  | 0.0000000E+00 |   -1.8277E+02 |
-------------------------------------------------------------------------------------------------------------------


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

TCounters:   Num_Calls   Total_Its       Total_Time (seconds)
    vcs_basopt:       2          2             NA     
    vcs_TP:           1         22             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.61        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
