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_gMols
          H2O(L)          0   NaCl_electrolyte             2000      
             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             4000      
              OH          1                air             0         
         NaCl(S)          2            NaCl(S)             5000      

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

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


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

             Phase IDs of species
            species     phaseID        phaseName    Initial_Estimated_gMols
          H2O(L)          0   NaCl_electrolyte             2000      
             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             4000      
              OH          1                air             0         
         NaCl(S)          2            NaCl(S)             5000      

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

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


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

	Solve a constant T, P problem:
		T    = 298.15 K
		Pres = 0.999951 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       Tmoles
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_gmol    ElemType ElActive
                          O   2.000000000000E+03    0         1
                          H   4.000000000000E+03    0         1
                          C   0.000000000000E+00    0         1
                          Fe  0.000000000000E+00    0         1
                          Si  0.000000000000E+00    0         1
                          N   8.000000000000E+03    0         1
                          Na  5.000000000000E+03    0         1
                          Cl  5.000000000000E+03    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)         5000  -       -174.5
 ---                        N2         4000  -       -23.05
 ---                    H2O(L)         2000  -       -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 = -1.2122e+06
 --- 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+03   0.000000000000E+00   2.000000000000E+03
   ---       H    4.000000000000E+03   0.000000000000E+00   4.000000000000E+03
   ---       C    0.000000000000E+00   0.000000000000E+00   0.000000000000E+00
   ---       N    8.000000000000E+03   0.000000000000E+00   8.000000000000E+03
   ---       Na   5.000000000000E+03   0.000000000000E+00   5.000000000000E+03
   ---       Cl   5.000000000000E+03   0.000000000000E+00   5.000000000000E+03
   ---       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+03   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)     (    5e+03) remains                          as component   0
   ---   N2          (    4e+03) remains                          as component   1
   ---   H2O(L)      (    2e+03) 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:     5e+03     4e+03     2e+03         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         4000  -       -23.05
 ---                    H2O(L)         2000  -       -123.7
 ---                        O2            0  -       -24.67
 ---                        H+            0  -            0
 ---                       Na+         5000  -         -104
 ---                       CO2            0  -       -184.5
 ---                        H2            0  -       -15.72
 ---                       H2O            0  -       -120.3
 ---                      NaCl            0  -       -100.8
 ---                       Cl-         5000  -       -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 = -1.23056e+06
   -----------------------------------------------------------------------------
   --- 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+         (    5e+03) replaces NaCl(S)     (        0) as component   0
   ---   Cl-         (    5e+03) replaces N2          (    4e+03) as component   1
   ---   N2          (    4e+03) replaces H2O(L)      (    2e+03) as component   2
   ---   H2O(L)      (    2e+03) 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:     5e+03     5e+03     4e+03     2e+03         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+         5000  -         -104
 ---                       Cl-         5000  -       -74.21
 ---                        N2         4000  -       -23.05
 ---                    H2O(L)         2000  -       -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 = -1.23056e+06
   --- vcs_inest:  Mole Numbers returned from linear programming (vcs_inest initial guess):
   --- vcs_inest:      SPECIES          MOLE_NUMBER      -SS_ChemPotential
   --- vcs_inest:      Na+                     5000           104
   --- vcs_inest:      Cl-                     5000          74.2
   --- vcs_inest:      N2                      4000            23
   --- vcs_inest:      H2O(L)                  2000           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                2e+03        2e+03
   --- vcs_inest:         H                4e+03        4e+03
   --- vcs_inest:         C                    0            0
   --- vcs_inest:         N                8e+03        8e+03
   --- vcs_inest:         Na               5e+03        5e+03
   --- vcs_inest:         Cl               5e+03        5e+03
   --- 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+         (    5e+03) remains                          as component   0
   ---   Cl-         (    5e+03) remains                          as component   1
   ---   N2          (    4e+03) remains                          as component   2
   ---   H2O(L)      (    2e+03) 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:     5e+03     5e+03     4e+03     2e+03     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* =         -169.95 ff =         -91.483
   --- vcs_inest: CO2         fe* =         -261.82 ff =         -184.45
   --- vcs_inest: H2          fe* =         -93.089 ff =         -15.717
   --- vcs_inest: H+          fe* =               0 ff =               0 dg* =         -44.446
   --- 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* =          40.163
   --- vcs_inest: NaCl(S)     fe* =         -174.52 ff =         -174.52 dg* =          5.4241
   --- vcs_inest: OH          fe* =         -91.487 ff =         -91.487 dg* =         -12.524
   --- vcs_inest: direction (Na+         ) = -4.30402e-31
   --- vcs_inest: direction (Cl-         ) = -4.30402e-31
   --- vcs_inest: direction (N2          ) = 0
   --- vcs_inest: direction (H2O(L)      ) = -2.81169e+23
   --- vcs_inest: direction (OH-         ) = 2.81169e+23
   --- vcs_inest: direction (CO2         ) = 0
   --- vcs_inest: direction (H2          ) = 1.37454e+09
   --- vcs_inest: direction (H+          ) = 2.81169e+23
   --- vcs_inest: direction (H2O         ) = 52.7108
   --- vcs_inest: direction (NaCl        ) = 4.30402e-31
   --- vcs_inest: direction (O2          ) = 3.60859e-14
   --- vcs_inest: direction (NaCl(S)     ) = 0 (ssPhase doesn't exist -> stability not checked)
   --- vcs_inest: direction (OH          ) = 2.74908e+09
   --- 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+          5000
   --- vcs_inest:      Cl-          5000
   --- vcs_inest:      N2           4000
   --- vcs_inest:      H2O(L)       1879.65
   --- vcs_inest:      OH-          120.348
   --- vcs_inest:      CO2          1e-30
   --- vcs_inest:      H2           5.88342e-13
   --- vcs_inest:      H+           120.348
   --- vcs_inest:      H2O          2.25617e-20
   --- vcs_inest:      NaCl         1.84225e-52
   --- vcs_inest:      O2           1.54458e-35
   --- vcs_inest:      NaCl(S)      0
   --- vcs_inest:      OH           1.17668e-12
   --- 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+03   2.000000000000E+03   2.000000000000E+03
   ---       H    4.000000000000E+03   4.000000000000E+03   4.000000000000E+03
   ---       C    0.000000000000E+00   1.000000000000E-30   0.000000000000E+00
   ---       N    8.000000000000E+03   8.000000000000E+03   8.000000000000E+03
   ---       Na   5.000000000000E+03   5.000000000000E+03   5.000000000000E+03
   ---       Cl   5.000000000000E+03   5.000000000000E+03   5.000000000000E+03
   ---       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:           1.371114852090E-13   0.000000000000E+00
   --- vcs_inest: Initial guess now satisfies element abundances
   --- vcs_inest: Total Dimensionless Gibbs Free Energy =  -1.0103528E+06
VCS CALCULATION METHOD

 MultiPhase Object


   13 SPECIES      11 ELEMENTS               7 COMPONENTS
    5 PHASE1 SPECIES         7 PHASE2 SPECIES       1 SINGLE SPECIES PHASES

 PRESSURE            101325.000 ATM
 TEMPERATURE            298.150 K
 PHASE1 INERTS            0.000
 PHASE2 INERTS            0.000

 ELEMENTAL ABUNDANCES             CORRECT          FROM ESTIMATE           Type

                          O   2.000000000000E+03  2.000000000000E+03       0
                          H   4.000000000000E+03  4.000000000000E+03       0
                          C   0.000000000000E+00  0.000000000000E+00       0
                          N   8.000000000000E+03  8.000000000000E+03       0
                          Na  5.000000000000E+03  5.000000000000E+03       0
                          Cl  5.000000000000E+03  5.000000000000E+03       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
 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+03       Mol_Num
 Cl-           0  0  0  0  0  1 -1  0  1  0  0  1              -1.83974E+08   5.00000E+03       Mol_Num
 N2            0  0  0  2  0  0  0  0  0  0  0  2              -5.71282E+07   4.00000E+03       Mol_Num
 H2O(L)        1  2  0  0  0  0  0  0 -0  0  0  1              -3.06686E+08   1.87965E+03       Mol_Num
 OH-           1  1  0  0  0  0 -1  0  1  0  0  1              -2.26784E+08   1.20348E+02       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   3.60968E-13       Mol_Num
 H+            0  1  0  0  0  0  1  0 -1  0  0  1               0.00000E+00   1.20348E+02       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   1.84225E-52       Mol_Num
 O2            2  0  0  0  0  0  0  0  0  0  0  2              -6.11650E+07   1.54458E-35       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   1.17668E-12       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+         (    5e+03) remains                          as component   0
   ---   Cl-         (    5e+03) remains                          as component   1
   ---   N2          (    4e+03) remains                          as component   2
   ---   H2O(L)      (  1.9e+03) remains                          as component   3
   ---   OH-         (  1.2e+02) remains                          as component   4
   ---   OH          (  1.2e-12) replaces CO2         (        0) as component   5
   ---   CO2         (        0) replaces H2          (  3.6e-13) as component   6
   ---   Total number of components found =   7 (ne = 11)
    ---                Components:       0         1         2         3         4         5         6  
   ---          Components Moles:     5e+03     5e+03     4e+03  1.88e+03       120  1.18e-12         0
   ---   NonComponent|   Moles  |       Na+       Cl-       N2        H2O(L)    OH-       OH        CO2       
   ---   7 H+        |       120|      -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      |  1.84e-52|      -1.00      -1.00      -0.00       0.00      -0.00      -0.00      -0.00
   ---  10 O2        |  1.54e-35|      -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        |  3.61e-13|      -0.00      -0.00      -0.00      -2.00      -0.00       2.00      -0.00
   -----------------------------------------------------------------------------
   --- Minor species changed to major: H2O         
   --- Minor species changed to major: NaCl        
   --- Minor species changed to major: O2          
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E+02  -1.2035E+02   1.2434E+02 | Delta damped from -21851 to -120.348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-20  -2.2562E-20   2.4806E+01 | Delta damped from -5.59677e-19 to -2.25617e-20 due to H2O going negative
   --- NaCl            1.8422E-52   2.2913E-50  -1.2438E+02 | Normal Calc
   --- O2              1.5446E-35   2.0487E-35  -1.3264E+00 | Normal Calc
   --- NaCl(S)         0.0000E+00   5.0000E+03  -7.0667E+01 | Delta damped from 176667 to 5000 due to component 0 (       Na+) going neg
   --- H2              3.6097E-13  -3.6097E-13   9.0041E+01 | Delta damped from -5.88342e-13 to -3.60968e-13 due to H2 going negative
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- Major species changed to minor: H+          
   --- H+            0 1.2035E+02 1.2035E+01-1.0831E+02 | initial nonpos moles=   0.000E+00
   --- Major species changed to minor: H2O         
   --- H2O           0 2.2562E-20 2.2562E-21-2.0306E-20 | initial nonpos moles=   0.000E+00
   --- NaCl          1 1.8422E-52 2.3098E-50 2.2913E-50 | Normal Major Calc
   --- O2            1 1.5446E-35 3.5933E-35 2.0487E-35 | Normal Major Calc
   --- Zeroed species changed to major: NaCl(S)     
   --- NaCl(S)       1 0.0000E+00 5.0000E+01 5.0000E+01 | Born:IC=-1 to IC=1:DG=-7.0667E+01
   --- Major species changed to minor: H2          
   --- H2            0 3.6097E-13 3.6097E-14-3.2487E-13 | initial nonpos moles=   0.000E+00
   --- Na+           c 5.0000E+03 4.9500E+03-5.0000E+01 |
   --- Cl-           c 5.0000E+03 4.9500E+03-5.0000E+01 |
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- H2O(L)        c 1.8797E+03 1.9880E+03 1.0831E+02 |
   --- OH-           c 1.2035E+02 1.2035E+01-1.0831E+02 |
   --- OH            c 1.1767E-12 5.2694E-13-6.4974E-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 = -9.9366705866972E+05
   --- Total tentative Dimensionless Gibbs Free Energy = -1.0101365895000E+06
   --- subroutine FORCE: Beginning Slope = -17001.2
   --- subroutine FORCE: End Slope       = -16244.5
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   Na+               5.000000E+03  4.950000E+03 -6.616232E+01 -6.651386E+01
   ---   Cl-               5.000000E+03  4.950000E+03 -3.768766E+01 -3.686981E+01
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   H2O(L)            1.879652E+03  1.987965E+03 -1.986003E+02 -1.994452E+02
   ---   OH-               1.203483E+02  1.203483E+01 -9.196202E+01 -9.656442E+01
   ---   OH                1.176683E-12  5.269410E-13 -1.272494E+02 -1.280528E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E+02  1.203483E+01  1.770351E+01  1.425959E+01  1.243418E+02  1.171404E+02
   ---   H2O           0   2.256174E-20  2.256174E-21 -1.737938E+02 -1.760964E+02  2.480647E+01  2.334881E+01
   ---   NaCl          1   1.842246E-52  2.309761E-50 -2.282275E+02 -2.233961E+02 -1.243775E+02 -1.200125E+02
   ---   O2            1   1.544580E-35  3.593286E-35 -1.131234E+02 -1.122791E+02 -1.326384E+00  1.041560E+00
   ---   NaCl(S)       1   0.000000E+00  5.000000E+01 -1.745168E+02 -1.745168E+02 -7.066681E+01 -7.113311E+01
   ---   H2            0   3.609680E-13  3.609680E-14 -5.266126E+01 -5.496385E+01  9.004051E+01  8.782103E+01
   ---                                                        Norms of Delta G():  8.626399E+01  8.310993E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   1.1912035E+04
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   5.0000000E+01
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -9.9366705866972E+05
   --- Total New Dimensionless Gibbs Free Energy = -1.0101365895000E+06
   --- Increment counter increased, step is accepted:    1
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   ---   Noncomponent turned from minor to major: H+
   ---   Noncomponent turned from major to minor: NaCl
   ---   Noncomponent turned from major to minor: O2
   ---   Noncomponent turned from minor to major: H2
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E+01  -1.2035E+01   1.1714E+02 | Delta damped from -847.23 to -12.0348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-21   0.0000E+00   2.3349E+01 | Skipped: IC =   0 and DG >0:   2.335E+01
   --- NaCl            2.3098E-50   2.7720E-48  -1.2001E+02 | Normal Calc
   --- O2              3.5933E-35   0.0000E+00   1.0416E+00 | Skipped: IC =   0 and DG >0:   1.042E+00
   --- NaCl(S)         5.0000E+01   4.9500E+03  -7.1133E+01 | Delta damped from 176054 to 4950 due to component 0 (       Na+) going neg
   --- H2              3.6097E-14  -3.6097E-14   8.7821E+01 | Delta damped from -2.6347e-13 to -3.60968e-14 due to H2 going negative
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- Major species changed to minor: H+          
   --- H+            0 1.2035E+01 1.2035E+00-1.0831E+01 | initial nonpos moles=   0.000E+00
   --- H2O           0 2.2562E-21 2.2562E-21 0.0000E+00 | minor species not considered
   --- NaCl          0 2.3098E-50 2.3098E-50 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-35 3.5933E-35 0.0000E+00 | minor species not considered
   --- NaCl(S)       1 5.0000E+01 2.5250E+03 2.4750E+03 | Line Search reduced step size from 4950 to 2475
   --- Major species changed to minor: H2          
   --- H2            0 3.6097E-14 3.6097E-15-3.2487E-14 | initial nonpos moles=   0.000E+00
   --- Na+           c 4.9500E+03 2.4750E+03-2.4750E+03 |
   --- Cl-           c 4.9500E+03 2.4750E+03-2.4750E+03 |
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- H2O(L)        c 1.9880E+03 1.9988E+03 1.0831E+01 |
   --- OH-           c 1.2035E+01 1.2035E+00-1.0831E+01 |
   --- OH            c 5.2694E-13 4.6197E-13-6.4974E-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.0101365895000E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.1714444228373E+06
   --- subroutine FORCE: Beginning Slope = -177323
   --- subroutine FORCE: End Slope       = -100667
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   Na+               4.950000E+03  2.475000E+03 -6.651386E+01 -8.195038E+01
   ---   Cl-               4.950000E+03  2.475000E+03 -3.686981E+01 -5.219644E+01
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   H2O(L)            1.987965E+03  1.998797E+03 -1.994452E+02 -1.533394E+02
   ---   OH-               1.203483E+01  1.203483E+00 -9.656442E+01 -9.615932E+01
   ---   OH                5.269410E-13  4.619667E-13 -1.280528E+02 -1.281844E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E+01  1.203483E+00  1.425959E+01  1.214385E+01  1.171404E+02  6.932389E+01
   ---   H2O           0   2.256174E-21  2.256174E-21 -1.760964E+02 -1.760964E+02  2.334881E+01 -2.275706E+01
   ---   NaCl          0   2.309761E-50  2.309761E-50 -2.233961E+02 -2.233961E+02 -1.200125E+02 -8.924931E+01
   ---   O2            0   3.593286E-35  3.593286E-35 -1.122791E+02 -1.122791E+02  1.041560E+00  9.377968E+01
   ---   NaCl(S)       1   5.000000E+01  2.525000E+03 -1.745168E+02 -1.745168E+02 -7.113311E+01 -4.036997E+01
   ---   H2            0   3.609680E-14  3.609680E-15 -5.496385E+01 -5.726643E+01  8.782103E+01 -6.956487E+00
   ---                                                        Norms of Delta G():  8.310993E+01  6.293102E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   6.9512035E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   2.5250000E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.0101365895000E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.1714444228373E+06
   --- 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          (    4e+03) replaces Na+         (  2.5e+03) as component   0
   ---   NaCl(S)     (  2.5e+03) replaces Cl-         (  2.5e+03) as component   1
   ---   Na+         (  2.5e+03) remains                          as component   2
   ---   H2O(L)      (    2e+03) remains                          as component   3
   ---   OH-         (      1.2) remains                          as component   4
   ---   OH          (  4.6e-13) 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:     4e+03  2.52e+03  2.48e+03     2e+03       1.2  4.62e-13         0
   ---   NonComponent|   Moles  |       N2        NaCl(S)   Na+       H2O(L)    OH-       OH        CO2       
   ---   7 H+        |       1.2|      -0.00      -0.00       0.00      -1.00       1.00      -0.00      -0.00
   ---   8 H2O       |  2.26e-21|      -0.00      -0.00       0.00      -1.00      -0.00      -0.00      -0.00
   ---   9 NaCl      |  2.31e-50|      -0.00      -1.00       0.00       0.00      -0.00      -0.00      -0.00
   ---  10 O2        |  3.59e-35|      -0.00      -0.00       0.00       2.00      -0.00      -4.00      -0.00
   ---  11 Cl-       |  2.48e+03|      -0.00      -1.00       1.00       0.00      -0.00      -0.00      -0.00
   ---  12 H2        |  3.61e-15|      -0.00      -0.00       0.00      -2.00      -0.00       2.00      -0.00
   -----------------------------------------------------------------------------
   --- Minor species changed to major: H2O         
   --- Minor species changed to major: NaCl        
   --- Minor species changed to major: O2          
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E+00  -1.2035E+00   6.9324E+01 | Delta damped from -40.3088 to -1.20348 due to component 4 (       OH-) going neg
   --- H2O             2.2562E-21   5.1344E-20  -2.2757E+01 | Normal Calc
   --- NaCl            2.3098E-50   1.1290E-48  -4.8879E+01 | Normal Calc
   --- O2              3.5933E-35  -3.5933E-35   9.3780E+01 | Delta damped from -3.36977e-33 to -3.59329e-35 due to O2 going negative
   --- Cl-             2.4750E+03  -1.7110E+03   4.0370E+01 | Normal calc: diag adjusted from 0.000232641 to 0.023594 due to act coeff
   --- H2              3.6097E-15   2.4350E-14  -6.9565E+00 | Normal calc: diag adjusted from 2.85692e+14 to 2.85692e+14 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- Major species changed to minor: H+          
   --- H+            0 1.2035E+00 1.2035E-01-1.0831E+00 | initial nonpos moles=   0.000E+00
   --- H2O           1 2.2562E-21 5.3600E-20 5.1344E-20 | Normal Major Calc
   --- NaCl          1 2.3098E-50 1.1521E-48 1.1290E-48 | Normal Major Calc
   --- Major species changed to minor: O2          
   --- O2            0 3.5933E-35 3.5933E-36-3.2340E-35 | initial nonpos moles=   0.000E+00
   --- Cl-           1 2.4750E+03 7.6397E+02-1.7110E+03 | Normal Major Calc
   --- H2            1 3.6097E-15 2.7959E-14 2.4350E-14 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 2.5250E+03 4.2360E+03 1.7110E+03 |
   --- Na+           c 2.4750E+03 7.6397E+02-1.7110E+03 |
   --- H2O(L)        c 1.9988E+03 1.9999E+03 1.0831E+00 |
   --- OH-           c 1.2035E+00 1.2035E-01-1.0831E+00 |
   --- OH            c 4.6197E-13 5.1067E-13 4.8699E-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.1714444228373E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2106119897651E+06
   --- subroutine FORCE: Beginning Slope = -69149.2
   --- subroutine FORCE: End Slope       = -13558.2
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           2.525000E+03  4.236028E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.475000E+03  7.639723E+02 -8.195038E+01 -9.818651E+01
   ---   H2O(L)            1.998797E+03  1.999880E+03 -1.533394E+02 -1.259451E+02
   ---   OH-               1.203483E+00  1.203483E-01 -9.615932E+01 -9.745470E+01
   ---   OH                4.619667E-13  5.106660E-13 -1.281844E+02 -1.280842E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E+00  1.203483E-01  1.214385E+01  8.612743E-01  6.932389E+01  2.935168E+01
   ---   H2O           1   2.256174E-21  5.360007E-20 -1.760964E+02 -1.729285E+02 -2.275706E+01 -4.698343E+01
   ---   NaCl          1   2.309761E-50  1.152094E-48 -2.233961E+02 -2.194865E+02 -4.887934E+01 -4.496974E+01
   ---   O2            0   3.593286E-35  3.593286E-36 -1.122791E+02 -1.145817E+02  9.377968E+01  1.458647E+02
   ---   Cl-           1   2.475000E+03  7.639723E+02 -5.219644E+01 -6.842487E+01  4.036997E+01  7.905404E+00
   ---   H2            1   3.609680E-15  2.795933E-14 -5.726643E+01 -5.521930E+01 -6.956487E+00 -5.949740E+01
   ---                                                        Norms of Delta G():  5.505396E+01  7.067554E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   3.5280650E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.2360277E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.1714444228373E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2106119897651E+06
   --- Increment counter increased, step is accepted:    3
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   ---   Noncomponent turned from minor to major: H+
   ---   Noncomponent turned from major to minor: H2O
   ---   Noncomponent turned from major to minor: NaCl
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-01  -1.2035E-01   2.9352E+01 | Delta damped from -1.76456 to -0.120348 due to component 4 (       OH-) going neg
   --- H2O             5.3600E-20   2.5183E-18  -4.6983E+01 | Normal Calc
   --- NaCl            1.1521E-48   5.1809E-47  -4.4970E+01 | Normal Calc
   --- O2              3.5933E-36   0.0000E+00   1.4586E+02 | Skipped: IC =   0 and DG >0:   1.459E+02
   --- Cl-             7.6397E+02  -5.3368E+02   7.9054E+00 | Normal calc: diag adjusted from 0.00148413 to 0.0148131 due to act coeff
   --- H2              2.7959E-14   1.3646E-12  -5.9497E+01 | Normal calc: diag adjusted from 4.35991e+13 to 4.35991e+13 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- Major species changed to minor: H+          
   --- H+            0 1.2035E-01 1.2035E-02-1.0831E-01 | initial nonpos moles=   0.000E+00
   --- H2O           0 5.3600E-20 5.3600E-20 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1521E-48 1.1521E-48 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-36 3.5933E-36 0.0000E+00 | minor species not considered
   --- Cl-           1 7.6397E+02 2.3030E+02-5.3368E+02 | Normal Major Calc
   --- H2            1 2.7959E-14 1.3926E-12 1.3646E-12 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.2360E+03 4.7697E+03 5.3368E+02 |
   --- Na+           c 7.6397E+02 2.3030E+02-5.3368E+02 |
   --- H2O(L)        c 1.9999E+03 2.0000E+03 1.0831E-01 |
   --- OH-           c 1.2035E-01 1.2035E-02-1.0831E-01 |
   --- OH            c 5.1067E-13 3.2400E-12 2.7293E-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.2106119897651E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127707270704E+06
   --- subroutine FORCE: Beginning Slope = -4222.11
   --- subroutine FORCE: End Slope       = -63.6751
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.236028E+03  4.769705E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               7.639723E+02  2.302954E+02 -9.818651E+01 -1.020814E+02
   ---   H2O(L)            1.999880E+03  1.999988E+03 -1.259451E+02 -1.240164E+02
   ---   OH-               1.203483E-01  1.203483E-02 -9.745470E+01 -1.000912E+02
   ---   OH                5.106660E-13  3.239958E-12 -1.280842E+02 -1.262366E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-01  1.203483E-02  8.612743E-01 -6.369208E+00  2.935168E+01  1.755602E+01
   ---   H2O           0   5.360007E-20  5.360007E-20 -1.729285E+02 -1.729285E+02 -4.698343E+01 -4.891215E+01
   ---   NaCl          0   1.152094E-48  1.152094E-48 -2.194865E+02 -2.194865E+02 -4.496974E+01 -4.496974E+01
   ---   O2            0   3.593286E-36  3.593286E-36 -1.145817E+02 -1.145817E+02  1.458647E+02  1.423317E+02
   ---   Cl-           1   7.639723E+02  2.302954E+02 -6.842487E+01 -7.231959E+01  7.905404E+00  1.157509E-01
   ---   H2            1   2.795933E-14  1.392605E-12 -5.521930E+01 -5.131112E+01 -5.949740E+01 -5.575146E+01
   ---                                                        Norms of Delta G():  7.067554E+01  6.842204E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4606029E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7697046E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2106119897651E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127707270704E+06
   --- Increment counter increased, step is accepted:    4
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: Cl- failed

==============================================================================================================
 Iteration =   4, Iterations since last evaluation of optimal basis =   2 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-02   0.0000E+00   1.7556E+01 | Skipped: IC =   0 and DG >0:   1.756E+01
   --- H2O             5.3600E-20   2.6217E-18  -4.8912E+01 | Normal Calc
   --- NaCl            1.1521E-48   5.1809E-47  -4.4970E+01 | Normal Calc
   --- O2              3.5933E-36   0.0000E+00   1.4233E+02 | Skipped: IC =   0 and DG >0:   1.423E+02
   --- Cl-             2.3030E+02  -7.1963E+00   1.1575E-01 | Normal calc: diag adjusted from 0.00705888 to 0.0160847 due to act coeff
   --- H2              1.3926E-12   2.8552E-11  -5.5751E+01 | Normal calc: diag adjusted from 1.95266e+12 to 1.95266e+12 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            0 1.2035E-02 1.2035E-02 0.0000E+00 | minor species not considered
   --- H2O           0 5.3600E-20 5.3600E-20 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1521E-48 1.1521E-48 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-36 3.5933E-36 0.0000E+00 | minor species not considered
   --- Cl-           1 2.3030E+02 2.2314E+02-7.1597E+00 | Line Search reduced step size from -7.19634 to -7.15967
   --- H2            1 1.3926E-12 2.9944E-11 2.8552E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7697E+03 4.7769E+03 7.1597E+00 |
   --- Na+           c 2.3030E+02 2.2314E+02-7.1597E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.7103E-11 |
   --- OH-           c 1.2035E-02 1.2035E-02 0.0000E+00 |
   --- OH            c 3.2400E-12 6.0343E-11 5.7103E-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.2127707270704E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127711421602E+06
   --- subroutine FORCE: Beginning Slope = -0.828738
   --- subroutine FORCE: End Slope       = -2.2127e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.769705E+03  4.776864E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.302954E+02  2.231357E+02 -1.020814E+02 -1.021393E+02
   ---   H2O(L)            1.999988E+03  1.999988E+03 -1.240164E+02 -1.240033E+02
   ---   OH-               1.203483E-02  1.203483E-02 -1.000912E+02 -1.000943E+02
   ---   OH                3.239958E-12  6.034298E-11 -1.262366E+02 -1.233121E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-02  1.203483E-02 -6.369208E+00 -6.438843E+00  1.755602E+01  1.747012E+01
   ---   H2O           0   5.360007E-20  5.360007E-20 -1.729285E+02 -1.729285E+02 -4.891215E+01 -4.892527E+01
   ---   NaCl          0   1.152094E-48  1.152094E-48 -2.194865E+02 -2.194865E+02 -4.496974E+01 -4.496974E+01
   ---   O2            0   3.593286E-36  3.593286E-36 -1.145817E+02 -1.145817E+02  1.423317E+02  1.306600E+02
   ---   Cl-           1   2.302954E+02  2.231357E+02 -7.231959E+01 -7.237746E+01  1.157509E-01  3.090322E-06
   ---   H2            1   1.392605E-12  2.994412E-11 -5.131112E+01 -4.824296E+01 -5.575146E+01 -4.686058E+01
   ---                                                        Norms of Delta G():  6.842204E+01  6.323122E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462835E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768643E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127707270704E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127711421602E+06
   --- Increment counter increased, step is accepted:    5
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: Cl- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-02   0.0000E+00   1.7470E+01 | Skipped: IC =   0 and DG >0:   1.747E+01
   --- H2O             5.3600E-20   2.6224E-18  -4.8925E+01 | Normal Calc
   --- NaCl            1.1521E-48   5.1809E-47  -4.4970E+01 | Normal Calc
   --- O2              3.5933E-36   0.0000E+00   1.3066E+02 | Skipped: IC =   0 and DG >0:   1.307E+02
   --- Cl-             2.2314E+02  -1.9016E-04   3.0903E-06 | Normal calc: diag adjusted from 0.00732802 to 0.0162509 due to act coeff
   --- H2              2.9944E-11   4.7009E-10  -4.6861E+01 | Normal calc: diag adjusted from 9.96833e+10 to 9.96833e+10 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            0 1.2035E-02 1.2035E-07-1.2035E-02 | minor species alternative calc
   --- H2O           0 5.3600E-20 5.3600E-18 5.3064E-18 | minor species alternative calc
   --- NaCl          0 1.1521E-48 1.1521E-46 1.1406E-46 | minor species alternative calc
   --- O2            0 3.5933E-36 3.5933E-41-3.5933E-36 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02-1.9016E-04 | Line Search reduced step size from -0.000190163 to -0.000190163
   --- H2            1 2.9944E-11 5.0004E-10 4.7009E-10 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.9016E-04 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.9016E-04 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03 1.2035E-02 |
   --- OH-           c 1.2035E-02 1.2035E-07-1.2035E-02 |
   --- OH            c 6.0343E-11 1.0005E-09 9.4019E-10 |
   --- CO2           c 0.0000E+00 0.0000E+00 0.0000E+00 |
   --------------------------------------------------------------------------------
   --- Finished Main Loop
   --- Reduction in step size due to component OH- going negative =   9.900E-01
   --- Subroutine vcs_dfe called for all species using tentative solution
   --- Subroutine vcs_deltag called for all noncomponents
   --- Total Old       Dimensionless Gibbs Free Energy = -1.2127711421602E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713275871E+06
   --- subroutine FORCE: Beginning Slope = -0.208148
   --- subroutine FORCE: End Slope       = -0.0984108
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776864E+03  4.776864E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231357E+02  2.231356E+02 -1.021393E+02 -1.021394E+02
   ---   H2O(L)            1.999988E+03  2.000000E+03 -1.240033E+02 -1.240032E+02
   ---   OH-               1.203483E-02  1.203483E-04 -1.000943E+02 -1.046995E+02
   ---   OH                6.034298E-11  9.911396E-10 -1.233121E+02 -1.205133E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-02  1.203483E-04 -6.438843E+00 -1.104401E+01  1.747012E+01  8.259761E+00
   ---   H2O           0   5.360007E-20  5.306995E-18 -1.729285E+02 -1.683333E+02 -4.892527E+01 -4.433007E+01
   ---   NaCl          0   1.152094E-48  1.140699E-46 -2.194865E+02 -2.148913E+02 -4.496974E+01 -4.037450E+01
   ---   O2            0   3.593286E-36  3.593286E-38 -1.145817E+02 -1.191869E+02  1.306600E+02  1.148597E+02
   ---   Cl-           1   2.231357E+02  2.231356E+02 -7.237746E+01 -7.237754E+01  3.090322E-06 -1.424915E-04
   ---   H2            1   2.994412E-11  4.953424E-10 -4.824296E+01 -4.543705E+01 -4.686058E+01 -3.845710E+01
   ---                                                        Norms of Delta G():  6.323122E+01  5.527986E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462712E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768644E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127711421602E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713275871E+06
   --- Increment counter increased, step is accepted:    6
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Reevaluate major-minor status of noncomponents:
   ---   Noncomponent turned from minor to major: H+
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-04  -1.2035E-04   8.2598E+00 | Delta damped from -0.000497024 to -0.000120348 due to component 4 (       OH-) going neg
   --- H2O             5.3070E-18   2.3526E-16  -4.4330E+01 | Normal Calc
   --- NaCl            1.1407E-46   4.6055E-45  -4.0375E+01 | Normal Calc
   --- O2              3.5933E-38   0.0000E+00   1.1486E+02 | Skipped: IC =   0 and DG >0:   1.149E+02
   --- Cl-             2.2314E+02   8.7683E-03  -1.4249E-04 | Normal calc: diag adjusted from 0.00732802 to 0.0162508 due to act coeff
   --- H2              4.9534E-10   6.3518E-09  -3.8457E+01 | Normal calc: diag adjusted from 6.05456e+09 to 6.05456e+09 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- Major species changed to minor: H+          
   --- H+            0 1.2035E-04 1.2035E-05-1.0831E-04 | initial nonpos moles=   0.000E+00
   --- H2O           0 5.3070E-18 5.3070E-18 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-46 1.1407E-46 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-38 3.5933E-38 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02 8.7683E-03 | Normal Major Calc
   --- H2            1 4.9534E-10 6.8471E-09 6.3518E-09 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-8.7683E-03 |
   --- Na+           c 2.2314E+02 2.2314E+02 8.7683E-03 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03 1.0830E-04 |
   --- OH-           c 1.2035E-04 1.2035E-05-1.0831E-04 |
   --- OH            c 9.9114E-10 1.3695E-08 1.2704E-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.2127713275871E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713283213E+06
   --- subroutine FORCE: Beginning Slope = -0.000896137
   --- subroutine FORCE: End Slope       = -0.000396059
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776864E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231356E+02  2.231443E+02 -1.021394E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240032E+02 -1.240033E+02
   ---   OH-               1.203483E-04  1.203483E-05 -1.046995E+02 -1.070020E+02
   ---   OH                9.911396E-10  1.369465E-08 -1.205133E+02 -1.178874E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-04  1.203483E-05 -1.104401E+01 -1.334651E+01  8.259761E+00  3.654696E+00
   ---   H2O           0   5.306995E-18  5.306995E-18 -1.683333E+02 -1.683333E+02 -4.433007E+01 -4.433005E+01
   ---   NaCl          0   1.140699E-46  1.140699E-46 -2.148913E+02 -2.148913E+02 -4.037450E+01 -4.037450E+01
   ---   O2            0   3.593286E-38  3.593286E-38 -1.191869E+02 -1.191869E+02  1.148597E+02  1.043560E+02
   ---   Cl-           1   2.231356E+02  2.231443E+02 -7.237754E+01 -7.237747E+01 -1.424915E-04 -1.296550E-06
   ---   H2            1   4.953424E-10  6.847096E-09 -4.543705E+01 -4.281071E+01 -3.845710E+01 -3.057892E+01
   ---                                                        Norms of Delta G():  5.527986E+01  5.071795E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462887E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768557E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713275871E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713283213E+06
   --- Increment counter increased, step is accepted:    7
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: Cl- failed

==============================================================================================================
 Iteration =   7, Iterations since last evaluation of optimal basis =   5 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-05   0.0000E+00   3.6547E+00 | Skipped: IC =   0 and DG >0:   3.655E+00
   --- H2O             5.3070E-18   2.3526E-16  -4.4330E+01 | Normal Calc
   --- NaCl            1.1407E-46   4.6055E-45  -4.0375E+01 | Normal Calc
   --- O2              3.5933E-38   0.0000E+00   1.0436E+02 | Skipped: IC =   0 and DG >0:   1.044E+02
   --- Cl-             2.2314E+02   7.9785E-05  -1.2965E-06 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              6.8471E-09   6.9794E-08  -3.0579E+01 | Normal calc: diag adjusted from 4.38132e+08 to 4.38132e+08 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            0 1.2035E-05 1.2035E-05 0.0000E+00 | minor species not considered
   --- H2O           0 5.3070E-18 5.3070E-18 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-46 1.1407E-46 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-38 3.5933E-38 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02 7.9785E-05 | Normal Major Calc
   --- H2            1 6.8471E-09 7.6641E-08 6.9794E-08 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-7.9785E-05 |
   --- Na+           c 2.2314E+02 2.2314E+02 7.9785E-05 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-1.3959E-07 |
   --- OH-           c 1.2035E-05 1.2035E-05 0.0000E+00 |
   --- OH            c 1.3695E-08 1.5328E-07 1.3959E-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.2127713283213E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713283231E+06
   --- subroutine FORCE: Beginning Slope = -2.13432e-06
   --- subroutine FORCE: End Slope       = -1.6285e-06
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231443E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.203483E-05  1.203483E-05 -1.070020E+02 -1.070020E+02
   ---   OH                1.369465E-08  1.532823E-07 -1.178874E+02 -1.154721E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-05  1.203483E-05 -1.334651E+01 -1.334651E+01  3.654696E+00  3.654697E+00
   ---   H2O           0   5.306995E-18  5.306995E-18 -1.683333E+02 -1.683333E+02 -4.433005E+01 -4.433005E+01
   ---   NaCl          0   1.140699E-46  1.140699E-46 -2.148913E+02 -2.148913E+02 -4.037450E+01 -4.037450E+01
   ---   O2            0   3.593286E-38  3.593286E-38 -1.191869E+02 -1.191869E+02  1.043560E+02  9.469493E+01
   ---   Cl-           1   2.231443E+02  2.231444E+02 -7.237747E+01 -7.237746E+01 -1.296550E-06  2.530243E-10
   ---   H2            1   6.847096E-09  7.664092E-08 -4.281071E+01 -4.039541E+01 -3.057892E+01 -2.333307E+01
   ---                                                        Norms of Delta G():  5.071795E+01  4.676212E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462888E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713283213E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713283231E+06
   --- Increment counter increased, step is accepted:    8
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H2 failed

==============================================================================================================
 Iteration =   8, Iterations since last evaluation of optimal basis =   6 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-05   0.0000E+00   3.6547E+00 | Skipped: IC =   0 and DG >0:   3.655E+00
   --- H2O             5.3070E-18   2.3526E-16  -4.4330E+01 | Normal Calc
   --- NaCl            1.1407E-46   4.6055E-45  -4.0375E+01 | Normal Calc
   --- O2              3.5933E-38   0.0000E+00   9.4695E+01 | Skipped: IC =   0 and DG >0:   9.469E+01
   --- Cl-             2.2314E+02  -1.5570E-08   2.5302E-10 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              7.6641E-08   5.9609E-07  -2.3333E+01 | Normal calc: diag adjusted from 3.91435e+07 to 3.91435e+07 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            0 1.2035E-05 1.2035E-05 0.0000E+00 | minor species not considered
   --- H2O           0 5.3070E-18 5.3070E-18 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-46 1.1407E-46 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-38 3.5933E-38 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-1.5569E-08 | Line Search reduced step size from -1.55702e-08 to -1.55693e-08
   --- H2            1 7.6641E-08 6.7273E-07 5.9609E-07 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.5569E-08 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.5569E-08 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-1.1922E-06 |
   --- OH-           c 1.2035E-05 1.2035E-05 0.0000E+00 |
   --- OH            c 1.5328E-07 1.3455E-06 1.1922E-06 |
   --- 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.2127713283231E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713283344E+06
   --- subroutine FORCE: Beginning Slope = -1.39086e-05
   --- subroutine FORCE: End Slope       = -1.00241e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.203483E-05  1.203483E-05 -1.070020E+02 -1.070020E+02
   ---   OH                1.532823E-07  1.345463E-06 -1.154721E+02 -1.132999E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-05  1.203483E-05 -1.334651E+01 -1.334651E+01  3.654697E+00  3.654697E+00
   ---   H2O           0   5.306995E-18  5.306995E-18 -1.683333E+02 -1.683333E+02 -4.433005E+01 -4.433005E+01
   ---   NaCl          0   1.140699E-46  1.140699E-46 -2.148913E+02 -2.148913E+02 -4.037450E+01 -4.037450E+01
   ---   O2            0   3.593286E-38  3.593286E-38 -1.191869E+02 -1.191869E+02  9.469493E+01  8.600608E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.530243E-10  2.161556E-09
   ---   H2            1   7.664092E-08  6.727313E-07 -4.039541E+01 -3.822319E+01 -2.333307E+01 -1.681643E+01
   ---                                                        Norms of Delta G():  4.676212E+01  4.337520E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462888E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713283231E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713283344E+06
   --- Increment counter increased, step is accepted:    9
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H2 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.2035E-05   0.0000E+00   3.6547E+00 | Skipped: IC =   0 and DG >0:   3.655E+00
   --- H2O             5.3070E-18   2.3526E-16  -4.4330E+01 | Normal Calc
   --- NaCl            1.1407E-46   4.6055E-45  -4.0375E+01 | Normal Calc
   --- O2              3.5933E-38   0.0000E+00   8.6006E+01 | Skipped: IC =   0 and DG >0:   8.601E+01
   --- Cl-             2.2314E+02  -1.3301E-07   2.1616E-09 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              6.7273E-07   3.7710E-06  -1.6816E+01 | Normal calc: diag adjusted from 4.45943e+06 to 4.45943e+06 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            0 1.2035E-05 3.1133E-07-1.1724E-05 | minor species alternative calc
   --- H2O           0 5.3070E-18 5.3070E-16 5.2539E-16 | minor species alternative calc
   --- NaCl          0 1.1407E-46 1.1407E-44 1.1293E-44 | minor species alternative calc
   --- O2            0 3.5933E-38 3.5933E-43-3.5933E-38 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02-1.3301E-07 | Normal Major Calc
   --- H2            1 6.7273E-07 4.4437E-06 3.7710E-06 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.3301E-07 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.3301E-07 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03 4.1815E-06 |
   --- OH-           c 1.2035E-05 3.1133E-07-1.1724E-05 |
   --- OH            c 1.3455E-06 8.8874E-06 7.5420E-06 |
   --- 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.2127713283344E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713284057E+06
   --- subroutine FORCE: Beginning Slope = -0.00010626
   --- subroutine FORCE: End Slope       = 7.8913e-07
   --- subroutine FORCE produced no adjustments (al = 0.992628)
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.203483E-05  3.113331E-07 -1.070020E+02 -1.106567E+02
   ---   OH                1.345463E-06  8.887422E-06 -1.132999E+02 -1.114120E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            0   1.203483E-05  3.113331E-07 -1.334651E+01 -1.700121E+01  3.654697E+00 -3.654697E+00
   ---   H2O           0   5.306995E-18  5.306995E-16 -1.683333E+02 -1.637281E+02 -4.433005E+01 -3.972488E+01
   ---   NaCl          0   1.140699E-46  1.140699E-44 -2.148913E+02 -2.102861E+02 -4.037450E+01 -3.576933E+01
   ---   O2            0   3.593286E-38  3.593286E-43 -1.191869E+02 -1.306998E+02  8.600608E+01  6.694156E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.161556E-09 -1.265616E-07
   ---   H2            1   6.727313E-07  4.443711E-06 -3.822319E+01 -3.633529E+01 -1.681643E+01 -1.115273E+01
   ---                                                        Norms of Delta G():  4.337520E+01  3.529973E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462888E+03
   ---                  air =   4.0000000E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713283344E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713284057E+06
   --- 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: H+
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              3.1133E-07   5.6891E-07  -3.6547E+00 | Normal calc: diag adjusted from 6.42399e+06 to 6.42399e+06 due to act coeff
   --- H2O             5.3070E-16   2.1082E-14  -3.9725E+01 | Normal Calc
   --- NaCl            1.1407E-44   4.0802E-43  -3.5769E+01 | Normal Calc
   --- O2              3.5933E-43   0.0000E+00   6.6942E+01 | Skipped: IC =   0 and DG >0:   6.694E+01
   --- Cl-             2.2314E+02   7.7881E-06  -1.2656E-07 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              4.4437E-06   1.6520E-05  -1.1153E+01 | Normal calc: diag adjusted from 675111 to 675111 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 3.1133E-07 8.8025E-07 5.6891E-07 | Normal Major Calc
   --- H2O           0 5.3070E-16 5.3070E-16 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-44 1.1407E-44 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-43 3.5933E-43 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02 7.7881E-06 | Line Search reduced step size from 7.78814e-06 to 7.78813e-06
   --- H2            1 4.4437E-06 2.0964E-05 1.6520E-05 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-7.7881E-06 |
   --- Na+           c 2.2314E+02 2.2314E+02 7.7881E-06 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-3.3609E-05 |
   --- OH-           c 3.1133E-07 8.8025E-07 5.6891E-07 |
   --- OH            c 8.8874E-06 4.1927E-05 3.3040E-05 |
   --- 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.2127713284057E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713285433E+06
   --- subroutine FORCE: Beginning Slope = -0.00018632
   --- subroutine FORCE: End Slope       = -0.000108256
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               3.113331E-07  8.802473E-07 -1.106567E+02 -1.096174E+02
   ---   OH                8.887422E-06  4.192709E-05 -1.114120E+02 -1.098607E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   3.113331E-07  8.802473E-07 -1.700121E+01 -1.596187E+01 -3.654697E+00 -1.576018E+00
   ---   H2O           0   5.306995E-16  5.306995E-16 -1.637281E+02 -1.637281E+02 -3.972488E+01 -3.972488E+01
   ---   NaCl          0   1.140699E-44  1.140699E-44 -2.102861E+02 -2.102861E+02 -3.576933E+01 -3.576933E+01
   ---   O2            0   3.593286E-43  3.593286E-43 -1.306998E+02 -1.306998E+02  6.694156E+01  6.073638E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01 -1.265616E-07  6.670993E-08
   ---   H2            1   4.443711E-06  2.096354E-05 -3.633529E+01 -3.478400E+01 -1.115273E+01 -6.498844E+00
   ---                                                        Norms of Delta G():  3.529973E+01  3.314397E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462888E+03
   ---                  air =   4.0000001E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713284057E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713285433E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              8.8025E-07   6.9364E-07  -1.5760E+00 | Normal calc: diag adjusted from 2.27209e+06 to 2.27209e+06 due to act coeff
   --- H2O             5.3070E-16   2.1082E-14  -3.9725E+01 | Normal Calc
   --- NaCl            1.1407E-44   4.0802E-43  -3.5769E+01 | Normal Calc
   --- O2              3.5933E-43   0.0000E+00   6.0736E+01 | Skipped: IC =   0 and DG >0:   6.074E+01
   --- Cl-             2.2314E+02  -4.1051E-06   6.6710E-08 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              2.0964E-05   4.5413E-05  -6.4988E+00 | Normal calc: diag adjusted from 143106 to 143106 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 8.8025E-07 1.5739E-06 6.9364E-07 | Normal Major Calc
   --- H2O           0 5.3070E-16 5.3070E-16 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-44 1.1407E-44 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-43 3.5933E-43 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-4.1051E-06 | Normal Major Calc
   --- H2            1 2.0964E-05 6.6376E-05 4.5413E-05 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 4.1051E-06 |
   --- Na+           c 2.2314E+02 2.2314E+02-4.1051E-06 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-9.1520E-05 |
   --- OH-           c 8.8025E-07 1.5739E-06 6.9364E-07 |
   --- OH            c 4.1927E-05 1.3275E-04 9.0826E-05 |
   --- 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.2127713285433E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713287458E+06
   --- subroutine FORCE: Beginning Slope = -0.000296225
   --- subroutine FORCE: End Slope       = -0.000138396
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               8.802473E-07  1.573890E-06 -1.096174E+02 -1.090363E+02
   ---   OH                4.192709E-05  1.327530E-04 -1.098607E+02 -1.087081E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   8.802473E-07  1.573890E-06 -1.596187E+01 -1.538076E+01 -1.576018E+00 -4.138125E-01
   ---   H2O           0   5.306995E-16  5.306995E-16 -1.637281E+02 -1.637281E+02 -3.972488E+01 -3.972488E+01
   ---   NaCl          0   1.140699E-44  1.140699E-44 -2.102861E+02 -2.102861E+02 -3.576933E+01 -3.576933E+01
   ---   O2            0   3.593286E-43  3.593286E-43 -1.306998E+02 -1.306998E+02  6.073638E+01  5.612615E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  6.670993E-08  1.729748E-07
   ---   H2            1   2.096354E-05  6.637648E-05 -3.478400E+01 -3.363144E+01 -6.498844E+00 -3.041170E+00
   ---                                                        Norms of Delta G():  3.314397E+01  3.166774E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462887E+03
   ---                  air =   4.0000002E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713285433E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713287458E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.5739E-06   3.2565E-07  -4.1381E-01 | Normal calc: diag adjusted from 1.27074e+06 to 1.27074e+06 due to act coeff
   --- H2O             5.3070E-16   2.1082E-14  -3.9725E+01 | Normal Calc
   --- NaCl            1.1407E-44   4.0802E-43  -3.5769E+01 | Normal Calc
   --- O2              3.5933E-43   0.0000E+00   5.6126E+01 | Skipped: IC =   0 and DG >0:   5.613E+01
   --- Cl-             2.2314E+02  -1.0644E-05   1.7297E-07 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              6.6376E-05   6.7287E-05  -3.0412E+00 | Normal calc: diag adjusted from 45196.7 to 45196.7 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.5739E-06 1.8995E-06 3.2565E-07 | Normal Major Calc
   --- H2O           0 5.3070E-16 5.3070E-16 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-44 1.1407E-44 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-43 3.5933E-43 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-1.0644E-05 | Normal Major Calc
   --- H2            1 6.6376E-05 1.3366E-04 6.7287E-05 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.0644E-05 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.0644E-05 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-1.3490E-04 |
   --- OH-           c 1.5739E-06 1.8995E-06 3.2565E-07 |
   --- OH            c 1.3275E-04 2.6733E-04 1.3457E-04 |
   --- 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.2127713287458E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288717E+06
   --- subroutine FORCE: Beginning Slope = -0.000204767
   --- subroutine FORCE: End Slope       = -6.33441e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.573890E-06  1.899538E-06 -1.090363E+02 -1.088482E+02
   ---   OH                1.327530E-04  2.673278E-04 -1.087081E+02 -1.080081E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.573890E-06  1.899538E-06 -1.538076E+01 -1.519270E+01 -4.138125E-01 -3.769177E-02
   ---   H2O           0   5.306995E-16  5.306995E-16 -1.637281E+02 -1.637281E+02 -3.972488E+01 -3.972488E+01
   ---   NaCl          0   1.140699E-44  1.140699E-44 -2.102861E+02 -2.102861E+02 -3.576933E+01 -3.576933E+01
   ---   O2            0   3.593286E-43  3.593286E-43 -1.306998E+02 -1.306998E+02  5.612615E+01  5.332621E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  1.729748E-07  2.478945E-07
   ---   H2            1   6.637648E-05  1.336639E-04 -3.363144E+01 -3.293146E+01 -3.041170E+00 -9.412140E-01
   ---                                                        Norms of Delta G():  3.166774E+01  3.082770E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462885E+03
   ---                  air =   4.0000004E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713287458E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288717E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8995E-06   3.5798E-08  -3.7692E-02 | Normal calc: diag adjusted from 1.05289e+06 to 1.05289e+06 due to act coeff
   --- H2O             5.3070E-16   2.1082E-14  -3.9725E+01 | Normal Calc
   --- NaCl            1.1407E-44   4.0802E-43  -3.5769E+01 | Normal Calc
   --- O2              3.5933E-43   0.0000E+00   5.3326E+01 | Skipped: IC =   0 and DG >0:   5.333E+01
   --- Cl-             2.2314E+02  -1.5255E-05   2.4789E-07 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.3366E-04   4.1935E-05  -9.4121E-01 | Normal calc: diag adjusted from 22444.4 to 22444.4 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8995E-06 1.9353E-06 3.5798E-08 | Normal Major Calc
   --- H2O           0 5.3070E-16 5.3070E-14 5.2539E-14 | minor species alternative calc
   --- NaCl          0 1.1407E-44 1.1407E-42 1.1293E-42 | minor species alternative calc
   --- O2            0 3.5933E-43 3.5933E-48-3.5933E-43 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02-1.5255E-05 | Normal Major Calc
   --- H2            1 1.3366E-04 1.7560E-04 4.1935E-05 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.5255E-05 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.5255E-05 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-8.3907E-05 |
   --- OH-           c 1.8995E-06 1.9353E-06 3.5798E-08 |
   --- OH            c 2.6733E-04 3.5120E-04 8.3871E-05 |
   --- 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.2127713288717E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288932E+06
   --- subroutine FORCE: Beginning Slope = -3.94716e-05
   --- subroutine FORCE: End Slope       = -5.14065e-06
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.899538E-06  1.935336E-06 -1.088482E+02 -1.088296E+02
   ---   OH                2.673278E-04  3.511986E-04 -1.080081E+02 -1.077353E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.899538E-06  1.935336E-06 -1.519270E+01 -1.517403E+01 -3.769177E-02 -3.507540E-04
   ---   H2O           0   5.306995E-16  5.306995E-14 -1.637281E+02 -1.591230E+02 -3.972488E+01 -3.511971E+01
   ---   NaCl          0   1.140699E-44  1.140699E-42 -2.102861E+02 -2.056810E+02 -3.576933E+01 -3.116416E+01
   ---   O2            0   3.593286E-43  3.593286E-48 -1.306998E+02 -1.422127E+02  5.332621E+01  4.072178E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.478945E-07  1.524955E-07
   ---   H2            1   1.336639E-04  1.755993E-04 -3.293146E+01 -3.265858E+01 -9.412140E-01 -1.225846E-01
   ---                                                        Norms of Delta G():  3.082770E+01  2.537347E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288717E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288932E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9353E-06   3.3941E-10  -3.5075E-04 | Normal calc: diag adjusted from 1.03341e+06 to 1.03341e+06 due to act coeff
   --- H2O             5.3070E-14   1.8638E-12  -3.5120E+01 | Normal Calc
   --- NaCl            1.1407E-42   3.5549E-41  -3.1164E+01 | Normal Calc
   --- O2              3.5933E-48   0.0000E+00   4.0722E+01 | Skipped: IC =   0 and DG >0:   4.072E+01
   --- Cl-             2.2314E+02  -9.3840E-06   1.5250E-07 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.7560E-04   7.1753E-06  -1.2258E-01 | Normal calc: diag adjusted from 17084.3 to 17084.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9353E-06 1.9357E-06 3.3941E-10 | Normal Major Calc
   --- H2O           0 5.3070E-14 5.3070E-14 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-42 1.1407E-42 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-48 3.5933E-48 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-9.3840E-06 | Line Search reduced step size from -9.38401e-06 to -9.38401e-06
   --- H2            1 1.7560E-04 1.8277E-04 7.1753E-06 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 9.3840E-06 |
   --- Na+           c 2.2314E+02 2.2314E+02-9.3840E-06 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-1.4351E-05 |
   --- OH-           c 1.9353E-06 1.9357E-06 3.3941E-10 |
   --- OH            c 3.5120E-04 3.6555E-04 1.4351E-05 |
   --- 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.2127713288932E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- subroutine FORCE: Beginning Slope = -8.79577e-07
   --- subroutine FORCE: End Slope       = -1.74959e-08
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935336E-06  1.935676E-06 -1.088296E+02 -1.088294E+02
   ---   OH                3.511986E-04  3.655491E-04 -1.077353E+02 -1.076952E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935336E-06  1.935676E-06 -1.517403E+01 -1.517386E+01 -3.507540E-04 -1.095025E-07
   ---   H2O           0   5.306995E-14  5.306995E-14 -1.591230E+02 -1.591230E+02 -3.511971E+01 -3.511971E+01
   ---   NaCl          0   1.140699E-42  1.140699E-42 -2.056810E+02 -2.056810E+02 -3.116416E+01 -3.116416E+01
   ---   O2            0   3.593286E-48  3.593286E-48 -1.422127E+02 -1.422127E+02  4.072178E+01  4.056158E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  1.524955E-07  2.602316E-08
   ---   H2            1   1.755993E-04  1.827746E-04 -3.265858E+01 -3.261853E+01 -1.225846E-01 -2.438330E-03
   ---                                                        Norms of Delta G():  2.537347E+01  2.533062E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288932E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06   1.0598E-13  -1.0950E-07 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-14   1.8638E-12  -3.5120E+01 | Normal Calc
   --- NaCl            1.1407E-42   3.5549E-41  -3.1164E+01 | Normal Calc
   --- O2              3.5933E-48   0.0000E+00   4.0562E+01 | Skipped: IC =   0 and DG >0:   4.056E+01
   --- Cl-             2.2314E+02  -1.6014E-06   2.6023E-08 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.8277E-04   1.4855E-07  -2.4383E-03 | Normal calc: diag adjusted from 16413.7 to 16413.7 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 1.0598E-13 | Normal Major Calc
   --- H2O           0 5.3070E-14 5.3070E-14 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-42 1.1407E-42 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-48 3.5933E-48 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-1.6014E-06 | Line Search reduced step size from -1.60137e-06 to -1.60137e-06
   --- H2            1 1.8277E-04 1.8292E-04 1.4855E-07 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 1.6014E-06 |
   --- Na+           c 2.2314E+02 2.2314E+02-1.6014E-06 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-2.9711E-07 |
   --- OH-           c 1.9357E-06 1.9357E-06 1.0598E-13 |
   --- OH            c 3.6555E-04 3.6585E-04 2.9711E-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.2127713288937E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- subroutine FORCE: Beginning Slope = -3.62267e-10
   --- subroutine FORCE: End Slope       = -1.4885e-13
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.655491E-04  3.658463E-04 -1.076952E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -1.095025E-07 -1.846186E-08
   ---   H2O           0   5.306995E-14  5.306995E-14 -1.591230E+02 -1.591230E+02 -3.511971E+01 -3.511971E+01
   ---   NaCl          0   1.140699E-42  1.140699E-42 -2.056810E+02 -2.056810E+02 -3.116416E+01 -3.116416E+01
   ---   O2            0   3.593286E-48  3.593286E-48 -1.422127E+02 -1.422127E+02  4.056158E+01  4.055833E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.602316E-08  5.387051E-10
   ---   H2            1   1.827746E-04  1.829231E-04 -3.261853E+01 -3.261772E+01 -2.438330E-03 -9.961790E-07
   ---                                                        Norms of Delta G():  2.533062E+01  2.532975E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- 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         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06   1.7868E-14  -1.8462E-08 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-14   1.8638E-12  -3.5120E+01 | Normal Calc
   --- NaCl            1.1407E-42   3.5549E-41  -3.1164E+01 | Normal Calc
   --- O2              3.5933E-48   0.0000E+00   4.0558E+01 | Skipped: IC =   0 and DG >0:   4.056E+01
   --- Cl-             2.2314E+02  -3.3150E-08   5.3871E-10 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.8292E-04   6.0741E-11  -9.9618E-07 | Normal calc: diag adjusted from 16400.3 to 16400.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 1.7868E-14 | Normal Major Calc
   --- H2O           0 5.3070E-14 5.3070E-14 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-42 1.1407E-42 0.0000E+00 | minor species not considered
   --- O2            0 3.5933E-48 3.5933E-48 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-3.3149E-08 | Line Search reduced step size from -3.31499e-08 to -3.3149e-08
   --- H2            1 1.8292E-04 1.8292E-04 6.0741E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 3.3149E-08 |
   --- Na+           c 2.2314E+02 2.2314E+02-3.3149E-08 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-1.2150E-10 |
   --- OH-           c 1.9357E-06 1.9357E-06 1.7868E-14 |
   --- OH            c 3.6585E-04 3.6585E-04 1.2148E-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.2127713288937E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- subroutine FORCE: Beginning Slope = -7.83672e-17
   --- subroutine FORCE: End Slope       = -1.43861e-20
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658463E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -1.846186E-08 -3.962413E-10
   ---   H2O           0   5.306995E-14  5.306995E-14 -1.591230E+02 -1.591230E+02 -3.511971E+01 -3.511971E+01
   ---   NaCl          0   1.140699E-42  1.140699E-42 -2.056810E+02 -2.056810E+02 -3.116416E+01 -3.116416E+01
   ---   O2            0   3.593286E-48  3.593286E-48 -1.422127E+02 -1.422127E+02  4.055833E+01  4.055833E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  5.387051E-10  2.131628E-13
   ---   H2            1   1.829231E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -9.961790E-07 -1.203944E-10
   ---                                                        Norms of Delta G():  2.532975E+01  2.532975E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Increment counter increased, step is accepted:   17
   --- 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: H2O failed

==============================================================================================================
 Iteration =  17, Iterations since last evaluation of optimal basis =  15 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 0.0000E+00 | Normal Major Calc
   --- H2O           0 5.3070E-14 5.3070E-12 5.2539E-12 | minor species alternative calc
   --- NaCl          0 1.1407E-42 1.1407E-40 1.1293E-40 | minor species alternative calc
   --- O2            0 3.5933E-48 3.5933E-53-3.5933E-48 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 0.0000E+00 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-1.1293E-40 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.2539E-12 |
   --- OH-           c 1.9357E-06 1.9357E-06 0.0000E+00 |
   --- OH            c 3.6585E-04 3.6585E-04 1.4373E-47 |
   --- 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.2127713288937E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- subroutine FORCE: Beginning Slope = -1.84516e-10
   --- subroutine FORCE: End Slope       = -1.60321e-10
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -3.962413E-10 -3.962271E-10
   ---   H2O           0   5.306995E-14  5.306995E-12 -1.591230E+02 -1.545178E+02 -3.511971E+01 -3.051454E+01
   ---   NaCl          0   1.140699E-42  1.140699E-40 -2.056810E+02 -2.010758E+02 -3.116416E+01 -2.655899E+01
   ---   O2            0   3.593286E-48  3.593286E-53 -1.422127E+02 -1.537257E+02  4.055833E+01  2.904540E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.131628E-13  2.273737E-13
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.203944E-10 -1.203944E-10
   ---                                                        Norms of Delta G():  2.532975E+01  2.033123E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- 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:  MAJOR SPECIES CONVERGENCE achieved
   --- Equilibrium check for minor species: H2O failed

==============================================================================================================
 Iteration =  18, Iterations since last evaluation of optimal basis =  16 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 0.0000E+00 | Normal Major Calc
   --- H2O           0 5.3070E-12 5.3070E-10 5.2539E-10 | minor species alternative calc
   --- NaCl          0 1.1407E-40 1.1407E-38 1.1293E-38 | minor species alternative calc
   --- O2            0 3.5933E-53 3.5933E-58-3.5933E-53 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 0.0000E+00 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-1.1293E-38 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.2539E-10 |
   --- OH-           c 1.9357E-06 1.9357E-06 0.0000E+00 |
   --- OH            c 3.6585E-04 3.6585E-04 1.4373E-52 |
   --- 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.2127713288937E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- subroutine FORCE: Beginning Slope = -1.60321e-08
   --- subroutine FORCE: End Slope       = -1.36126e-08
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -3.962271E-10 -3.950049E-10
   ---   H2O           0   5.306995E-12  5.306995E-10 -1.545178E+02 -1.499126E+02 -3.051454E+01 -2.590937E+01
   ---   NaCl          0   1.140699E-40  1.140699E-38 -2.010758E+02 -1.964706E+02 -2.655899E+01 -2.195382E+01
   ---   O2            0   3.593286E-53  3.593286E-58 -1.537257E+02 -1.652386E+02  2.904540E+01  1.753248E+01
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  2.273737E-13  1.207923E-12
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.203944E-10 -1.205933E-10
   ---                                                        Norms of Delta G():  2.033123E+01  1.560264E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Increment counter increased, step is accepted:   19
   --- 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: H2O failed

==============================================================================================================
 Iteration =  19, Iterations since last evaluation of optimal basis =  17 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 0.0000E+00 | Normal Major Calc
   --- H2O           0 5.3070E-10 5.3070E-08 5.2539E-08 | minor species alternative calc
   --- NaCl          0 1.1407E-38 1.1407E-36 1.1293E-36 | minor species alternative calc
   --- O2            0 3.5933E-58 3.5933E-63-3.5933E-58 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 0.0000E+00 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-1.1293E-36 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.2539E-08 |
   --- OH-           c 1.9357E-06 1.9357E-06 0.0000E+00 |
   --- OH            c 3.6585E-04 3.6585E-04 1.4373E-57 |
   --- 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.2127713288937E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713288948E+06
   --- subroutine FORCE: Beginning Slope = -1.36126e-06
   --- subroutine FORCE: End Slope       = -1.11931e-06
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -3.950049E-10 -2.723510E-10
   ---   H2O           0   5.306995E-10  5.306995E-08 -1.499126E+02 -1.453075E+02 -2.590937E+01 -2.130420E+01
   ---   NaCl          0   1.140699E-38  1.140699E-36 -1.964706E+02 -1.918654E+02 -2.195382E+01 -1.734865E+01
   ---   O2            0   3.593286E-58  3.593286E-63 -1.652386E+02 -1.767515E+02  1.753248E+01  6.019554E+00
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  1.207923E-12  9.644907E-11
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.205933E-10 -1.387264E-10
   ---                                                        Norms of Delta G():  1.560264E+01  1.148245E+01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000005E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288937E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713288948E+06
   --- Increment counter increased, step is accepted:   20
   --- 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: H2O failed

==============================================================================================================
 Iteration =  20, Iterations since last evaluation of optimal basis =  18 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 0.0000E+00 | Normal Major Calc
   --- H2O           0 5.3070E-08 5.3070E-06 5.2539E-06 | minor species alternative calc
   --- NaCl          0 1.1407E-36 1.1407E-34 1.1293E-34 | minor species alternative calc
   --- O2            0 3.5933E-63 3.5933E-65-3.5574E-63 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 0.0000E+00 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-1.1293E-34 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.2539E-06 |
   --- OH-           c 1.9357E-06 1.9357E-06 0.0000E+00 |
   --- OH            c 3.6585E-04 3.6585E-04 1.4229E-62 |
   --- 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.2127713288948E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713289876E+06
   --- subroutine FORCE: Beginning Slope = -0.000111931
   --- subroutine FORCE: End Slope       = -8.77355e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  2.000000E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01 -2.723510E-10  1.199358E-08
   ---   H2O           0   5.306995E-08  5.306995E-06 -1.453075E+02 -1.407023E+02 -2.130420E+01 -1.669903E+01
   ---   NaCl          0   1.140699E-36  1.140699E-34 -1.918654E+02 -1.872603E+02 -1.734865E+01 -1.274348E+01
   ---   O2            0   3.593286E-63  3.593286E-65 -1.767515E+02 -1.813567E+02  6.019554E+00  1.414383E+00
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  9.644907E-11  9.622411E-09
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.387264E-10 -1.953509E-09
   ---                                                        Norms of Delta G():  1.148245E+01  8.595101E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462884E+03
   ---                  air =   4.0000006E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713288948E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713289876E+06
   --- Increment counter increased, step is accepted:   21
   --- 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 =  21, Iterations since last evaluation of optimal basis =  19 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06  -1.1608E-14   1.1994E-08 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-06   8.8622E-05  -1.6699E+01 | Normal calc: diag adjusted from 188431 to 188431 due to act coeff
   --- NaCl            1.1407E-34   1.4536E-33  -1.2743E+01 | Normal Calc
   --- O2              3.5933E-65   0.0000E+00   1.4144E+00 | Skipped: IC =   0 and DG >0:   1.414E+00
   --- Cl-             2.2314E+02  -5.9213E-07   9.6224E-09 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.8292E-04   1.1911E-13  -1.9535E-09 | Normal calc: diag adjusted from 16400.3 to 16400.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06-1.1608E-14 | Normal Major Calc
   --- H2O           0 5.3070E-06 5.3070E-04 5.2539E-04 | minor species alternative calc
   --- NaCl          0 1.1407E-34 1.1407E-32 1.1293E-32 | minor species alternative calc
   --- O2            0 3.5933E-65 8.7344E-66-2.7198E-65 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02-5.9213E-07 | Line Search reduced step size from -5.92128e-07 to -5.92127e-07
   --- H2            1 1.8292E-04 1.8292E-04 1.1911E-13 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 5.9213E-07 |
   --- Na+           c 2.2314E+02 2.2314E+02-5.9213E-07 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-5.2539E-04 |
   --- OH-           c 1.9357E-06 1.9357E-06-1.1608E-14 |
   --- OH            c 3.6585E-04 3.6585E-04 2.3823E-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.2127713289876E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- subroutine FORCE: Beginning Slope = -0.00877355
   --- subroutine FORCE: End Slope       = -0.00635402
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231444E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            2.000000E+03  1.999999E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935676E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935676E-06 -1.517386E+01 -1.517386E+01  1.199358E-08  1.219510E-06
   ---   H2O           0   5.306995E-06  5.306995E-04 -1.407023E+02 -1.360971E+02 -1.669903E+01 -1.209386E+01
   ---   NaCl          0   1.140699E-34  1.140699E-32 -1.872603E+02 -1.826551E+02 -1.274348E+01 -8.138314E+00
   ---   O2            0   3.593286E-65  8.734397E-66 -1.813567E+02 -1.827711E+02  1.414383E+00  1.810204E-07
   ---   Cl-           1   2.231444E+02  2.231444E+02 -7.237746E+01 -7.237746E+01  9.622411E-09  9.525955E-07
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.953509E-09 -1.836251E-07
   ---                                                        Norms of Delta G():  8.595101E+00  5.951101E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462878E+03
   ---                  air =   4.0000011E+03
   ---              NaCl(S) =   4.7768556E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713289876E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- 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: H+ failed

==============================================================================================================
 Iteration =  22, Iterations since last evaluation of optimal basis =  20 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06  -1.1803E-12   1.2195E-06 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-04   6.4182E-03  -1.2094E+01 | Normal calc: diag adjusted from 1884.31 to 1884.31 due to act coeff
   --- NaCl            1.1407E-32   9.2834E-32  -8.1383E+00 | Normal Calc
   --- O2              8.7344E-66   0.0000E+00   1.8102E-07 | Skipped: IC =   0 and DG >0:   1.810E-07
   --- Cl-             2.2314E+02  -5.8619E-05   9.5260E-07 | Normal calc: diag adjusted from 0.00732768 to 0.0162506 due to act coeff
   --- H2              1.8292E-04   1.1196E-11  -1.8363E-07 | Normal calc: diag adjusted from 16400.3 to 16400.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06-1.1803E-12 | Line Search reduced step size from -1.18029e-12 to -1.18029e-12
   --- H2O           0 5.3070E-04 5.3070E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-32 1.1407E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.7344E-66 8.7344E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02-5.8619E-05 | Line Search reduced step size from -5.86192e-05 to -5.86192e-05
   --- H2            1 1.8292E-04 1.8292E-04 1.1196E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 5.8619E-05 |
   --- Na+           c 2.2314E+02 2.2314E+02-5.8619E-05 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-2.1213E-11 |
   --- OH-           c 1.9357E-06 1.9357E-06-1.1803E-12 |
   --- OH            c 3.6585E-04 3.6585E-04 2.2393E-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.2127713358426E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- subroutine FORCE: Beginning Slope = -5.58404e-11
   --- subroutine FORCE: End Slope       = -2.38539e-18
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231444E+02  2.231443E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.999999E+03  1.999999E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935676E-06  1.935675E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935676E-06  1.935675E-06 -1.517386E+01 -1.517386E+01  1.219510E-06 -7.012006E-07
   ---   H2O           0   5.306995E-04  5.306995E-04 -1.360971E+02 -1.360971E+02 -1.209386E+01 -1.209386E+01
   ---   NaCl          0   1.140699E-32  1.140699E-32 -1.826551E+02 -1.826551E+02 -8.138314E+00 -8.138314E+00
   ---   O2            0   8.734397E-66  8.734397E-66 -1.827711E+02 -1.827711E+02  1.810204E-07  1.487533E-07
   ---   Cl-           1   2.231444E+02  2.231443E+02 -7.237746E+01 -7.237746E+01  9.525955E-07  1.421085E-14
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -1.836251E-07 -2.125664E-07
   ---                                                        Norms of Delta G():  5.951101E+00  5.951101E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462877E+03
   ---                  air =   4.0000011E+03
   ---              NaCl(S) =   4.7768557E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- Increment counter increased, step is accepted:   23
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  23, Iterations since last evaluation of optimal basis =  21 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06   6.7865E-13  -7.0120E-07 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-04   6.4182E-03  -1.2094E+01 | Normal calc: diag adjusted from 1884.31 to 1884.31 due to act coeff
   --- NaCl            1.1407E-32   9.2834E-32  -8.1383E+00 | Normal Calc
   --- O2              8.7344E-66   0.0000E+00   1.4875E-07 | Skipped: IC =   0 and DG >0:   1.488E-07
   --- Cl-             2.2314E+02   0.0000E+00   1.4211E-14 | Skipped: superconverged DG =   1.421E-14
   --- H2              1.8292E-04   1.2961E-11  -2.1257E-07 | Normal calc: diag adjusted from 16400.3 to 16400.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 6.7865E-13 | Normal Major Calc
   --- H2O           0 5.3070E-04 5.3070E-04 0.0000E+00 | minor species not considered
   --- NaCl          0 1.1407E-32 1.1407E-32 0.0000E+00 | minor species not considered
   --- O2            0 8.7344E-66 8.7344E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 1.2961E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 0.0000E+00 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 2.0000E+03-2.6601E-11 |
   --- OH-           c 1.9357E-06 1.9357E-06 6.7865E-13 |
   --- OH            c 3.6585E-04 3.6585E-04 2.5922E-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.2127713358426E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- subroutine FORCE: Beginning Slope = -3.23096e-18
   --- subroutine FORCE: End Slope       = -3.85767e-26
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304522E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231443E+02  2.231443E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.999999E+03  1.999999E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935675E-06  1.935675E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935675E-06  1.935675E-06 -1.517386E+01 -1.517386E+01 -7.012006E-07 -5.684342E-14
   ---   H2O           0   5.306995E-04  5.306995E-04 -1.360971E+02 -1.360971E+02 -1.209386E+01 -1.209386E+01
   ---   NaCl          0   1.140699E-32  1.140699E-32 -1.826551E+02 -1.826551E+02 -8.138314E+00 -8.138314E+00
   ---   O2            0   8.734397E-66  8.734397E-66 -1.827711E+02 -1.827711E+02  1.487533E-07 -1.346685E-07
   ---   Cl-           1   2.231443E+02  2.231443E+02 -7.237746E+01 -7.237746E+01  1.421085E-14  8.526513E-14
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261772E+01 -2.125664E-07  0.000000E+00
   ---                                                        Norms of Delta G():  5.951101E+00  5.951101E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462877E+03
   ---                  air =   4.0000011E+03
   ---              NaCl(S) =   4.7768557E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- Increment counter increased, step is accepted:   24
   --- 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: H2O failed

==============================================================================================================
 Iteration =  24, Iterations since last evaluation of optimal basis =  22 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9357E-06 0.0000E+00 | major species is converged
   --- H2O           0 5.3070E-04 5.3070E-02 5.2539E-02 | minor species alternative calc
   --- NaCl          0 1.1407E-32 1.1407E-30 1.1293E-30 | minor species alternative calc
   --- O2            0 8.7344E-66 8.7344E-66 1.1762E-72 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8292E-04 1.8292E-04 0.0000E+00 | major species is converged
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03-1.1293E-30 |
   --- Na+           c 2.2314E+02 2.2314E+02 0.0000E+00 |
   --- H2O(L)        c 2.0000E+03 1.9999E+03-5.2539E-02 |
   --- OH-           c 1.9357E-06 1.9357E-06 0.0000E+00 |
   --- OH            c 3.6585E-04 3.6585E-04-4.7050E-72 |
   --- 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.2127713358426E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127717793882E+06
   --- subroutine FORCE: Beginning Slope = -0.635402
   --- subroutine FORCE: End Slope       = -0.39345
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304522E+01 -2.304524E+01
   ---   NaCl(S)           4.776856E+03  4.776856E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231443E+02  2.231443E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.999999E+03  1.999947E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.935675E-06  1.935675E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.658464E-04  3.658464E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151078E+02
   ---   H+            1   1.935675E-06  1.935675E-06 -1.517386E+01 -1.517377E+01 -5.684342E-14  1.226620E-04
   ---   H2O           0   5.306995E-04  5.306995E-02 -1.360971E+02 -1.314920E+02 -1.209386E+01 -7.488694E+00
   ---   NaCl          0   1.140699E-32  1.140699E-30 -1.826551E+02 -1.780499E+02 -8.138314E+00 -3.533157E+00
   ---   O2            0   8.734397E-66  8.734398E-66 -1.827711E+02 -1.827711E+02 -1.346685E-07  1.814679E-05
   ---   Cl-           1   2.231443E+02  2.231443E+02 -7.237746E+01 -7.237742E+01  8.526513E-14  9.526164E-05
   ---   H2            1   1.829232E-04  1.829232E-04 -3.261772E+01 -3.261773E+01  0.000000E+00 -1.814679E-05
   ---                                                        Norms of Delta G():  5.951101E+00  3.380427E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4462352E+03
   ---                  air =   4.0000536E+03
   ---              NaCl(S) =   4.7768557E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127713358426E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127717793882E+06
   --- Increment counter increased, step is accepted:   25
   --- 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 =  25, Iterations since last evaluation of optimal basis =  23 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9357E-06  -1.1872E-10   1.2266E-04 | Normal calc: diag adjusted from 1.03323e+06 to 1.03323e+06 due to act coeff
   --- H2O             5.3070E-02   3.9743E-01  -7.4887E+00 | Normal calc: diag adjusted from 18.8429 to 18.843 due to act coeff
   --- NaCl            1.1407E-30   4.0303E-30  -3.5332E+00 | Normal Calc
   --- O2              8.7344E-66   0.0000E+00   1.8147E-05 | Skipped: IC =   0 and DG >0:   1.815E-05
   --- Cl-             2.2314E+02  -5.8619E-03   9.5262E-05 | Normal calc: diag adjusted from 0.00732764 to 0.0162509 due to act coeff
   --- H2              1.8292E-04   1.1065E-09  -1.8147E-05 | Normal calc: diag adjusted from 16400.3 to 16400.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9357E-06 1.9356E-06-1.1871E-10 | Line Search reduced step size from -1.18717e-10 to -1.18713e-10
   --- H2O           0 5.3070E-02 5.3070E+00 5.2539E+00 | minor species alternative calc
   --- NaCl          0 1.1407E-30 3.9048E-29 3.7908E-29 | minor species alternative calc
   --- O2            0 8.7344E-66 8.7342E-66-1.5850E-70 | minor species alternative calc
   --- Cl-           1 2.2314E+02 2.2314E+02-5.8619E-03 | Line Search reduced step size from -0.00586195 to -0.00586192
   --- H2            1 1.8292E-04 1.8292E-04 1.1065E-09 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7769E+03 5.8619E-03 |
   --- Na+           c 2.2314E+02 2.2314E+02-5.8619E-03 |
   --- H2O(L)        c 1.9999E+03 1.9947E+03-5.2539E+00 |
   --- OH-           c 1.9357E-06 1.9356E-06-1.1871E-10 |
   --- OH            c 3.6585E-04 3.6585E-04 2.2130E-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.2127717793882E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127919394421E+06
   --- subroutine FORCE: Beginning Slope = -39.345
   --- subroutine FORCE: End Slope       = -15.1512
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304524E+01 -2.304655E+01
   ---   NaCl(S)           4.776856E+03  4.776862E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231443E+02  2.231384E+02 -1.021393E+02 -1.021345E+02
   ---   H2O(L)            1.999947E+03  1.994693E+03 -1.240033E+02 -1.240043E+02
   ---   OH-               1.935675E-06  1.935557E-06 -1.088294E+02 -1.088265E+02
   ---   OH                3.658464E-04  3.658486E-04 -1.076944E+02 -1.076957E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151078E+02 -5.151091E+02
   ---   H+            1   1.935675E-06  1.935557E-06 -1.517377E+01 -1.516555E+01  1.226620E-04  1.222325E-02
   ---   H2O           0   5.306995E-02  5.306995E+00 -1.314920E+02 -1.268881E+02 -7.488694E+00 -2.883781E+00
   ---   NaCl          0   1.140699E-30  3.904826E-29 -1.780499E+02 -1.745181E+02 -3.533157E+00 -1.312602E-03
   ---   O2            0   8.734398E-66  8.734239E-66 -1.827711E+02 -1.827724E+02  1.814679E-05  1.801673E-03
   ---   Cl-           1   2.231443E+02  2.231384E+02 -7.237742E+01 -7.237269E+01  9.526164E-05  9.547129E-03
   ---   H2            1   1.829232E-04  1.829243E-04 -3.261773E+01 -3.261904E+01 -1.814679E-05 -1.825869E-03
   ---                                                        Norms of Delta G():  3.380427E+00  1.177316E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4409695E+03
   ---                  air =   4.0053075E+03
   ---              NaCl(S) =   4.7768616E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127717793882E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127919394421E+06
   --- Increment counter increased, step is accepted:   26
   --- 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 =  26, Iterations since last evaluation of optimal basis =  24 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9356E-06  -1.1829E-08   1.2223E-02 | Normal calc: diag adjusted from 1.03329e+06 to 1.03329e+06 due to act coeff
   --- H2O             5.3070E+00   1.5308E+01  -2.8838E+00 | Normal calc: diag adjusted from 0.188273 to 0.188385 due to act coeff
   --- NaCl            3.9048E-29   5.1255E-32  -1.3126E-03 | Normal Calc
   --- O2              8.7342E-66   0.0000E+00   1.8017E-03 | Skipped: IC =   0 and DG >0:   1.802E-03
   --- Cl-             2.2314E+02  -5.8645E-01   9.5471E-03 | Normal calc: diag adjusted from 0.00732435 to 0.0162796 due to act coeff
   --- H2              1.8292E-04   1.1133E-07  -1.8259E-03 | Normal calc: diag adjusted from 16400.2 to 16400.2 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9356E-06 1.9238E-06-1.1793E-08 | Line Search reduced step size from -1.18294e-08 to -1.17932e-08
   --- H2O           0 5.3070E+00 5.3070E+00 0.0000E+00 | minor species not considered
   --- NaCl          0 3.9048E-29 3.9048E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7342E-66 8.7342E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2314E+02 2.2255E+02-5.8619E-01 | Line Search reduced step size from -0.586447 to -0.586192
   --- H2            1 1.8292E-04 1.8304E-04 1.1133E-07 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7769E+03 4.7774E+03 5.8619E-01 |
   --- Na+           c 2.2314E+02 2.2255E+02-5.8619E-01 |
   --- H2O(L)        c 1.9947E+03 1.9947E+03-2.1087E-07 |
   --- OH-           c 1.9356E-06 1.9238E-06-1.1793E-08 |
   --- OH            c 3.6585E-04 3.6607E-04 2.2266E-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.2127919394421E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- subroutine FORCE: Beginning Slope = -0.00559645
   --- subroutine FORCE: End Slope       = -1.37424e-09
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304655E+01 -2.304655E+01
   ---   NaCl(S)           4.776862E+03  4.777448E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.231384E+02  2.225523E+02 -1.021345E+02 -1.021393E+02
   ---   H2O(L)            1.994693E+03  1.994693E+03 -1.240043E+02 -1.240033E+02
   ---   OH-               1.935557E-06  1.923763E-06 -1.088265E+02 -1.088329E+02
   ---   OH                3.658486E-04  3.660713E-04 -1.076957E+02 -1.076951E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151091E+02 -5.151091E+02
   ---   H+            1   1.935557E-06  1.923763E-06 -1.516555E+01 -1.517737E+01  1.222325E-02 -7.032274E-03
   ---   H2O           0   5.306995E+00  5.306995E+00 -1.268881E+02 -1.268881E+02 -2.883781E+00 -2.884847E+00
   ---   NaCl          0   3.904826E-29  3.904826E-29 -1.745181E+02 -1.745181E+02 -1.312602E-03 -1.312602E-03
   ---   O2            0   8.734239E-66  8.734239E-66 -1.827724E+02 -1.827724E+02  1.801673E-03  1.501115E-03
   ---   Cl-           1   2.231384E+02  2.225523E+02 -7.237269E+01 -7.237746E+01  9.547129E-03  2.080583E-09
   ---   H2            1   1.829243E-04  1.830357E-04 -3.261904E+01 -3.261843E+01 -1.825869E-03 -2.133749E-03
   ---                                                        Norms of Delta G():  1.177316E+00  1.177738E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4397972E+03
   ---                  air =   4.0053075E+03
   ---              NaCl(S) =   4.7774477E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127919394421E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Increment counter increased, step is accepted:   27
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  27, Iterations since last evaluation of optimal basis =  25 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9238E-06   6.7642E-09  -7.0323E-03 | Normal calc: diag adjusted from 1.03963e+06 to 1.03963e+06 due to act coeff
   --- H2O             5.3070E+00   1.5314E+01  -2.8848E+00 | Normal calc: diag adjusted from 0.188272 to 0.188384 due to act coeff
   --- NaCl            3.9048E-29   5.1255E-32  -1.3126E-03 | Normal Calc
   --- O2              8.7342E-66   0.0000E+00   1.5011E-03 | Skipped: IC =   0 and DG >0:   1.501E-03
   --- Cl-             2.2255E+02  -1.2769E-07   2.0806E-09 | Normal calc: diag adjusted from 0.00734717 to 0.0162938 due to act coeff
   --- H2              1.8304E-04   1.3018E-07  -2.1337E-03 | Normal calc: diag adjusted from 16390.2 to 16390.2 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9238E-06 1.9305E-06 6.7642E-09 | Normal Major Calc
   --- H2O           0 5.3070E+00 5.3070E+00 0.0000E+00 | minor species not considered
   --- NaCl          0 3.9048E-29 3.9048E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7342E-66 8.7342E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2255E+02 2.2255E+02-1.2769E-07 | Line Search reduced step size from -1.27692e-07 to -1.27691e-07
   --- H2            1 1.8304E-04 1.8317E-04 1.3018E-07 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7774E+03 4.7774E+03 1.2769E-07 |
   --- Na+           c 2.2255E+02 2.2255E+02-1.2769E-07 |
   --- H2O(L)        c 1.9947E+03 1.9947E+03-2.6713E-07 |
   --- OH-           c 1.9238E-06 1.9305E-06 6.7642E-09 |
   --- OH            c 3.6607E-04 3.6633E-04 2.6037E-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.2127919422407E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- subroutine FORCE: Beginning Slope = -3.25348e-10
   --- subroutine FORCE: End Slope       = -1.82304e-13
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304655E+01 -2.304655E+01
   ---   NaCl(S)           4.777448E+03  4.777448E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.225523E+02  2.225523E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.994693E+03  1.994693E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.923763E-06  1.930528E-06 -1.088329E+02 -1.088294E+02
   ---   OH                3.660713E-04  3.663317E-04 -1.076951E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151091E+02 -5.151091E+02
   ---   H+            1   1.923763E-06  1.930528E-06 -1.517737E+01 -1.517386E+01 -7.032274E-03 -1.233520E-05
   ---   H2O           0   5.306995E+00  5.306995E+00 -1.268881E+02 -1.268881E+02 -2.884847E+00 -2.884847E+00
   ---   NaCl          0   3.904826E-29  3.904826E-29 -1.745181E+02 -1.745181E+02 -1.312602E-03 -1.312602E-03
   ---   O2            0   8.734239E-66  8.734239E-66 -1.827724E+02 -1.827724E+02  1.501115E-03 -1.342872E-03
   ---   Cl-           1   2.225523E+02  2.225523E+02 -7.237746E+01 -7.237746E+01  2.080583E-09  5.544649E-10
   ---   H2            1   1.830357E-04  1.831658E-04 -3.261843E+01 -3.261772E+01 -2.133749E-03 -7.588872E-07
   ---                                                        Norms of Delta G():  1.177738E+00  1.177734E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4397972E+03
   ---                  air =   4.0053075E+03
   ---              NaCl(S) =   4.7774477E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Increment counter increased, step is accepted:   28
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  28, Iterations since last evaluation of optimal basis =  26 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9305E-06   1.1907E-11  -1.2335E-05 | Normal calc: diag adjusted from 1.03599e+06 to 1.03599e+06 due to act coeff
   --- H2O             5.3070E+00   1.5314E+01  -2.8848E+00 | Normal calc: diag adjusted from 0.188272 to 0.188384 due to act coeff
   --- NaCl            3.9048E-29   5.1255E-32  -1.3126E-03 | Normal Calc
   --- O2              8.7342E-66   1.1729E-68  -1.3429E-03 | Normal Calc
   --- Cl-             2.2255E+02  -3.4029E-08   5.5446E-10 | Normal calc: diag adjusted from 0.00734717 to 0.0162938 due to act coeff
   --- H2              1.8317E-04   4.6334E-11  -7.5889E-07 | Normal calc: diag adjusted from 16378.6 to 16378.6 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9305E-06 1.9305E-06 1.1907E-11 | Normal Major Calc
   --- H2O           0 5.3070E+00 5.3070E+00 0.0000E+00 | minor species not considered
   --- NaCl          0 3.9048E-29 3.9048E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7342E-66 8.7342E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2255E+02 2.2255E+02-3.4028E-08 | Line Search reduced step size from -3.40292e-08 to -3.40283e-08
   --- H2            1 1.8317E-04 1.8317E-04 4.6334E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7774E+03 4.7774E+03 3.4028E-08 |
   --- Na+           c 2.2255E+02 2.2255E+02-3.4028E-08 |
   --- H2O(L)        c 1.9947E+03 1.9947E+03-1.0457E-10 |
   --- OH-           c 1.9305E-06 1.9305E-06 1.1907E-11 |
   --- OH            c 3.6633E-04 3.6633E-04 9.2668E-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.2127919422407E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- subroutine FORCE: Beginning Slope = -2.00902e-16
   --- subroutine FORCE: End Slope       = -2.21647e-20
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304655E+01 -2.304655E+01
   ---   NaCl(S)           4.777448E+03  4.777448E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.225523E+02  2.225523E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.994693E+03  1.994693E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.930528E-06  1.930540E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.663317E-04  3.663318E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151091E+02 -5.151091E+02
   ---   H+            1   1.930528E-06  1.930540E-06 -1.517386E+01 -1.517386E+01 -1.233520E-05 -4.458656E-10
   ---   H2O           0   5.306995E+00  5.306995E+00 -1.268881E+02 -1.268881E+02 -2.884847E+00 -2.884847E+00
   ---   NaCl          0   3.904826E-29  3.904826E-29 -1.745181E+02 -1.745181E+02 -1.312602E-03 -1.312602E-03
   ---   O2            0   8.734239E-66  8.734239E-66 -1.827724E+02 -1.827724E+02 -1.342872E-03 -1.343883E-03
   ---   Cl-           1   2.225523E+02  2.225523E+02 -7.237746E+01 -7.237746E+01  5.544649E-10  3.268497E-13
   ---   H2            1   1.831658E-04  1.831659E-04 -3.261772E+01 -3.261772E+01 -7.588872E-07 -1.237481E-10
   ---                                                        Norms of Delta G():  1.177734E+00  1.177734E+00
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.4397972E+03
   ---                  air =   4.0053075E+03
   ---              NaCl(S) =   4.7774477E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Increment counter increased, step is accepted:   29
   --- 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: H2O failed

==============================================================================================================
 Iteration =  29, Iterations since last evaluation of optimal basis =  27 (all species)
   --- vcs_RxnStepSizes not called because all majors have converged
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9305E-06 1.9305E-06 0.0000E+00 | Normal Major Calc
   --- H2O           0 5.3070E+00 9.5000E+01 8.9693E+01 | minor species alternative calc
   --- NaCl          0 3.9048E-29 3.9100E-29 5.1288E-32 | minor species alternative calc
   --- O2            0 8.7342E-66 8.7460E-66 1.1746E-68 | minor species alternative calc
   --- Cl-           1 2.2255E+02 2.2255E+02 0.0000E+00 | major species is converged
   --- H2            1 1.8317E-04 1.8317E-04 0.0000E+00 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7774E+03 4.7774E+03-5.1288E-32 |
   --- Na+           c 2.2255E+02 2.2255E+02 0.0000E+00 |
   --- H2O(L)        c 1.9947E+03 1.9050E+03-8.9693E+01 |
   --- OH-           c 1.9305E-06 1.9305E-06 0.0000E+00 |
   --- OH            c 3.6633E-04 3.6633E-04-4.6983E-68 |
   --- 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.2127919422407E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128664718909E+06
   --- subroutine FORCE: Beginning Slope = -258.75
   --- subroutine FORCE: End Slope       = -0.250825
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.304655E+01 -2.306870E+01
   ---   NaCl(S)           4.777448E+03  4.777448E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.225523E+02  2.225523E+02 -1.021393E+02 -1.020546E+02
   ---   H2O(L)            1.994693E+03  1.905000E+03 -1.240033E+02 -1.240226E+02
   ---   OH-               1.930540E-06  1.930540E-06 -1.088294E+02 -1.087787E+02
   ---   OH                3.663318E-04  3.663318E-04 -1.076944E+02 -1.077165E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151091E+02 -5.151313E+02
   ---   H+            1   1.930540E-06  1.930540E-06 -1.517386E+01 -1.502565E+01 -4.458656E-10  2.182319E-01
   ---   H2O           0   5.306995E+00  9.499962E+01 -1.268881E+02 -1.240254E+02 -2.884847E+00 -2.796499E-03
   ---   NaCl          0   3.904826E-29  3.909955E-29 -1.745181E+02 -1.745389E+02 -1.312602E-03 -2.214639E-02
   ---   O2            0   8.734239E-66  8.745985E-66 -1.827724E+02 -1.827932E+02 -1.343883E-03  2.773939E-02
   ---   Cl-           1   2.225523E+02  2.225523E+02 -7.237746E+01 -7.229274E+01  3.268497E-13  1.694501E-01
   ---   H2            1   1.831659E-04  1.831659E-04 -3.261772E+01 -3.263986E+01 -1.237481E-10 -2.773939E-02
   ---                                                        Norms of Delta G():  1.177734E+00  1.142918E-01
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3501045E+03
   ---                  air =   4.0950002E+03
   ---              NaCl(S) =   4.7774477E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2127919422407E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128664718909E+06
   --- Increment counter increased, step is accepted:   30
   --- 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 =  30, Iterations since last evaluation of optimal basis =  28 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.9305E-06  -2.1065E-07   2.1823E-01 | Normal calc: diag adjusted from 1.03598e+06 to 1.03598e+06 due to act coeff
   --- H2O             9.5000E+01   2.6604E-01  -2.7965E-03 | Normal calc: diag adjusted from 0.0103816 to 0.0105116 due to act coeff
   --- NaCl            3.9100E-29   8.6591E-31  -2.2146E-02 | Normal Calc
   --- O2              8.7460E-66   0.0000E+00   2.7739E-02 | Skipped: IC =   0 and DG >0:   2.774E-02
   --- Cl-             2.2255E+02  -1.0080E+01   1.6945E-01 | Normal calc: diag adjusted from 0.0072846 to 0.0168097 due to act coeff
   --- H2              1.8317E-04   1.6936E-06  -2.7739E-02 | Normal calc: diag adjusted from 16378.6 to 16378.6 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.9305E-06 1.7316E-06-1.9894E-07 | Line Search reduced step size from -2.10653e-07 to -1.98939e-07
   --- H2O           1 9.5000E+01 9.5266E+01 2.6604E-01 | Normal Major Calc
   --- NaCl          0 3.9100E-29 3.9100E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7460E-66 8.7460E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.2255E+02 2.1255E+02-1.0007E+01 | Line Search reduced step size from -10.0805 to -10.0066
   --- H2            1 1.8317E-04 1.8486E-04 1.6936E-06 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7774E+03 4.7875E+03 1.0007E+01 |
   --- Na+           c 2.2255E+02 2.1255E+02-1.0007E+01 |
   --- H2O(L)        c 1.9050E+03 1.9047E+03-2.6604E-01 |
   --- OH-           c 1.9305E-06 1.7316E-06-1.9894E-07 |
   --- OH            c 3.6633E-04 3.6972E-04 3.3873E-06 |
   --- 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.2128664718909E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673273685E+06
   --- subroutine FORCE: Beginning Slope = -1.69637
   --- subroutine FORCE: End Slope       = -0.0103127
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306870E+01 -2.306876E+01
   ---   NaCl(S)           4.777448E+03  4.787454E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.225523E+02  2.125456E+02 -1.020546E+02 -1.021391E+02
   ---   H2O(L)            1.905000E+03  1.904734E+03 -1.240226E+02 -1.240033E+02
   ---   OH-               1.930540E-06  1.731601E-06 -1.087787E+02 -1.088920E+02
   ---   OH                3.663318E-04  3.697190E-04 -1.077165E+02 -1.077074E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151313E+02 -5.151313E+02
   ---   H+            1   1.930540E-06  1.731601E-06 -1.502565E+01 -1.523615E+01  2.182319E-01 -1.248312E-01
   ---   H2O           1   9.499962E+01  9.526566E+01 -1.240254E+02 -1.240227E+02 -2.796499E-03 -1.935727E-02
   ---   NaCl          0   3.909955E-29  3.909955E-29 -1.745389E+02 -1.745390E+02 -2.214639E-02 -2.221136E-02
   ---   O2            0   8.745985E-66  8.745985E-66 -1.827932E+02 -1.827933E+02  2.773939E-02  2.970303E-02
   ---   Cl-           1   2.225523E+02  2.125456E+02 -7.229274E+01 -7.237721E+01  1.694501E-01  5.159423E-04
   ---   H2            1   1.831659E-04  1.848595E-04 -3.263986E+01 -3.263073E+01 -2.773939E-02 -3.890700E-02
   ---                                                        Norms of Delta G():  1.142918E-01  5.604637E-02
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3298252E+03
   ---                  air =   4.0952662E+03
   ---              NaCl(S) =   4.7874544E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128664718909E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673273685E+06
   --- Increment counter increased, step is accepted:   31
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  31, Iterations since last evaluation of optimal basis =  29 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.7316E-06   1.0808E-07  -1.2483E-01 | Normal calc: diag adjusted from 1.155e+06 to 1.155e+06 due to act coeff
   --- H2O             9.5266E+01   1.8497E+00  -1.9357E-02 | Normal calc: diag adjusted from 0.0103486 to 0.0104652 due to act coeff
   --- NaCl            3.9100E-29   8.6845E-31  -2.2211E-02 | Normal Calc
   --- O2              8.7460E-66   0.0000E+00   2.9703E-02 | Skipped: IC =   0 and DG >0:   2.970E-02
   --- Cl-             2.1255E+02  -3.0238E-02   5.1594E-04 | Normal calc: diag adjusted from 0.00769288 to 0.0170625 due to act coeff
   --- H2              1.8486E-04   2.3974E-06  -3.8907E-02 | Normal calc: diag adjusted from 16228.5 to 16228.5 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.7316E-06 1.8397E-06 1.0808E-07 | Normal Major Calc
   --- H2O           1 9.5266E+01 9.7115E+01 1.8497E+00 | Normal Major Calc
   --- NaCl          0 3.9100E-29 3.9100E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7460E-66 8.7460E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1255E+02 2.1252E+02-3.0238E-02 | Line Search reduced step size from -0.0302383 to -0.0302376
   --- H2            1 1.8486E-04 1.8726E-04 2.3974E-06 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7875E+03 4.7875E+03 3.0238E-02 |
   --- Na+           c 2.1255E+02 2.1252E+02-3.0238E-02 |
   --- H2O(L)        c 1.9047E+03 1.9029E+03-1.8497E+00 |
   --- OH-           c 1.7316E-06 1.8397E-06 1.0808E-07 |
   --- OH            c 3.6972E-04 3.7451E-04 4.7949E-06 |
   --- 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.2128673273685E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673454999E+06
   --- subroutine FORCE: Beginning Slope = -0.0358203
   --- subroutine FORCE: End Slope       = -0.000556289
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306876E+01 -2.306921E+01
   ---   NaCl(S)           4.787454E+03  4.787485E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.125456E+02  2.125154E+02 -1.021391E+02 -1.021376E+02
   ---   H2O(L)            1.904734E+03  1.902884E+03 -1.240033E+02 -1.240036E+02
   ---   OH-               1.731601E-06  1.839680E-06 -1.088920E+02 -1.088304E+02
   ---   OH                3.697190E-04  3.745139E-04 -1.077074E+02 -1.076950E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151313E+02 -5.151318E+02
   ---   H+            1   1.731601E-06  1.839680E-06 -1.523615E+01 -1.517284E+01 -1.248312E-01  4.173489E-04
   ---   H2O           1   9.526566E+01  9.711533E+01 -1.240227E+02 -1.240039E+02 -1.935727E-02 -2.431358E-04
   ---   NaCl          0   3.909955E-29  3.909955E-29 -1.745390E+02 -1.745394E+02 -2.221136E-02 -2.266292E-02
   ---   O2            0   8.745985E-66  8.745985E-66 -1.827933E+02 -1.827937E+02  2.970303E-02 -2.115645E-02
   ---   Cl-           1   2.125456E+02  2.125154E+02 -7.237721E+01 -7.237570E+01  5.159423E-04  3.524263E-03
   ---   H2            1   1.848595E-04  1.872570E-04 -3.263073E+01 -3.261829E+01 -3.890700E-02 -9.331480E-04
   ---                                                        Norms of Delta G():  5.604637E-02  1.274577E-02
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3279150E+03
   ---                  air =   4.0971159E+03
   ---              NaCl(S) =   4.7874846E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673273685E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673454999E+06
   --- Increment counter increased, step is accepted:   32
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  32, Iterations since last evaluation of optimal basis =  30 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8397E-06  -3.8389E-10   4.1735E-04 | Normal calc: diag adjusted from 1.08715e+06 to 1.08715e+06 due to act coeff
   --- H2O             9.7115E+01   2.3684E-02  -2.4314E-04 | Normal calc: diag adjusted from 0.0101489 to 0.0102659 due to act coeff
   --- NaCl            3.9100E-29   8.8611E-31  -2.2663E-02 | Normal Calc
   --- O2              8.7460E-66   1.8503E-67  -2.1156E-02 | Normal Calc
   --- Cl-             2.1252E+02  -2.0641E-01   3.5243E-03 | Normal calc: diag adjusted from 0.00769281 to 0.0170744 due to act coeff
   --- H2              1.8726E-04   5.8246E-08  -9.3315E-04 | Normal calc: diag adjusted from 16020.8 to 16020.8 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8397E-06 1.8393E-06-3.8385E-10 | Line Search reduced step size from -3.83894e-10 to -3.83854e-10
   --- H2O           1 9.7115E+01 9.7139E+01 2.3684E-02 | Normal Major Calc
   --- NaCl          0 3.9100E-29 3.9100E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.7460E-66 8.7460E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1252E+02 2.1231E+02-2.0637E-01 | Line Search reduced step size from -0.206406 to -0.206373
   --- H2            1 1.8726E-04 1.8732E-04 5.8246E-08 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7875E+03 4.7877E+03 2.0637E-01 |
   --- Na+           c 2.1252E+02 2.1231E+02-2.0637E-01 |
   --- H2O(L)        c 1.9029E+03 1.9029E+03-2.3684E-02 |
   --- OH-           c 1.8397E-06 1.8393E-06-3.8385E-10 |
   --- OH            c 3.7451E-04 3.7463E-04 1.1649E-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.2128673454999E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458758E+06
   --- subroutine FORCE: Beginning Slope = -0.00073307
   --- subroutine FORCE: End Slope       = -1.86325e-05
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306921E+01 -2.306922E+01
   ---   NaCl(S)           4.787485E+03  4.787691E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.125154E+02  2.123090E+02 -1.021376E+02 -1.021393E+02
   ---   H2O(L)            1.902884E+03  1.902861E+03 -1.240036E+02 -1.240033E+02
   ---   OH-               1.839680E-06  1.839296E-06 -1.088304E+02 -1.088307E+02
   ---   OH                3.745139E-04  3.746304E-04 -1.076950E+02 -1.076947E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.839680E-06  1.839296E-06 -1.517284E+01 -1.517512E+01  4.173489E-04 -2.536731E-03
   ---   H2O           1   9.711533E+01  9.713902E+01 -1.240039E+02 -1.240036E+02 -2.431358E-04 -3.934388E-04
   ---   NaCl          0   3.909955E-29  3.909955E-29 -1.745394E+02 -1.745395E+02 -2.266292E-02 -2.266870E-02
   ---   O2            0   8.745985E-66  8.745985E-66 -1.827937E+02 -1.827937E+02 -2.115645E-02 -2.160639E-02
   ---   Cl-           1   2.125154E+02  2.123090E+02 -7.237570E+01 -7.237744E+01  3.524263E-03  4.513365E-05
   ---   H2            1   1.872570E-04  1.873152E-04 -3.261829E+01 -3.261799E+01 -9.331480E-04 -7.942174E-04
   ---                                                        Norms of Delta G():  1.274577E-02  1.283179E-02
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274786E+03
   ---                  air =   4.0971396E+03
   ---              NaCl(S) =   4.7876910E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673454999E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458758E+06
   --- Increment counter increased, step is accepted:   33
   --- 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 =  33, Iterations since last evaluation of optimal basis =  31 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8393E-06   2.3329E-09  -2.5367E-03 | Normal calc: diag adjusted from 1.08737e+06 to 1.08737e+06 due to act coeff
   --- H2O             9.7139E+01   3.8335E-02  -3.9344E-04 | Normal calc: diag adjusted from 0.0101463 to 0.0102631 due to act coeff
   --- NaCl            3.9100E-29   8.8634E-31  -2.2669E-02 | Normal Calc
   --- O2              8.7460E-66   1.8897E-67  -2.1606E-02 | Normal Calc
   --- Cl-             2.1231E+02  -2.6425E-03   4.5134E-05 | Normal calc: diag adjusted from 0.00770163 to 0.0170801 due to act coeff
   --- H2              1.8732E-04   4.9590E-08  -7.9422E-04 | Normal calc: diag adjusted from 16015.8 to 16015.8 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8393E-06 1.8416E-06 2.3329E-09 | Normal Major Calc
   --- H2O           1 9.7139E+01 9.7177E+01 3.8335E-02 | Normal Major Calc
   --- NaCl          0 3.9100E-29 3.9996E-29 8.9646E-31 | minor species alternative calc
   --- O2            0 8.7460E-66 8.9370E-66 1.9103E-67 | minor species alternative calc
   --- Cl-           1 2.1231E+02 2.1231E+02-2.6425E-03 | Line Search reduced step size from -0.00264247 to -0.00264247
   --- H2            1 1.8732E-04 1.8736E-04 4.9590E-08 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 2.6425E-03 |
   --- Na+           c 2.1231E+02 2.1231E+02-2.6425E-03 |
   --- H2O(L)        c 1.9029E+03 1.9028E+03-3.8335E-02 |
   --- OH-           c 1.8393E-06 1.8416E-06 2.3329E-09 |
   --- OH            c 3.7463E-04 3.7473E-04 9.9179E-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.2128673458758E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458836E+06
   --- subroutine FORCE: Beginning Slope = -1.52019e-05
   --- subroutine FORCE: End Slope       = -3.89075e-07
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306922E+01 -2.306923E+01
   ---   NaCl(S)           4.787691E+03  4.787694E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123090E+02  2.123064E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902861E+03  1.902822E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.839296E-06  1.841629E-06 -1.088307E+02 -1.088294E+02
   ---   OH                3.746304E-04  3.747296E-04 -1.076947E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.839296E-06  1.841629E-06 -1.517512E+01 -1.517381E+01 -2.536731E-03  5.923910E-05
   ---   H2O           1   9.713902E+01  9.717735E+01 -1.240036E+02 -1.240033E+02 -3.934388E-04 -5.113459E-06
   ---   NaCl          0   3.909955E-29  3.999601E-29 -1.745395E+02 -1.745168E+02 -2.266870E-02 -9.356609E-06
   ---   O2            0   8.745985E-66  8.937010E-66 -1.827937E+02 -1.827721E+02 -2.160639E-02 -1.036977E-03
   ---   Cl-           1   2.123090E+02  2.123064E+02 -7.237744E+01 -7.237743E+01  4.513365E-05  7.305597E-05
   ---   H2            1   1.873152E-04  1.873648E-04 -3.261799E+01 -3.261773E+01 -7.942174E-04 -2.194406E-05
   ---                                                        Norms of Delta G():  1.283179E-02  4.251988E-04
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274350E+03
   ---                  air =   4.0971779E+03
   ---              NaCl(S) =   4.7876936E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458758E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458836E+06
   --- Increment counter increased, step is accepted:   34
   --- 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 =  34, Iterations since last evaluation of optimal basis =  32 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06  -5.4548E-11   5.9239E-05 | Normal calc: diag adjusted from 1.086e+06 to 1.086e+06 due to act coeff
   --- H2O             9.7177E+01   4.9844E-04  -5.1135E-06 | Normal calc: diag adjusted from 0.0101423 to 0.010259 due to act coeff
   --- NaCl            3.9996E-29   3.7423E-34  -9.3566E-06 | Normal Calc
   --- O2              8.9370E-66   9.2675E-69  -1.0370E-03 | Normal Calc
   --- Cl-             2.1231E+02  -4.2772E-03   7.3056E-05 | Normal calc: diag adjusted from 0.00770172 to 0.0170804 due to act coeff
   --- H2              1.8736E-04   1.3705E-09  -2.1944E-05 | Normal calc: diag adjusted from 16011.5 to 16011.5 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06-5.4547E-11 | Line Search reduced step size from -5.45482e-11 to -5.45474e-11
   --- H2O           1 9.7177E+01 9.7178E+01 4.9844E-04 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.9370E-66 8.9370E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1231E+02 2.1230E+02-4.2772E-03 | Line Search reduced step size from -0.00427719 to -0.00427717
   --- H2            1 1.8736E-04 1.8737E-04 1.3705E-09 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 4.2772E-03 |
   --- Na+           c 2.1231E+02 2.1230E+02-4.2772E-03 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-4.9844E-04 |
   --- OH-           c 1.8416E-06 1.8416E-06-5.4547E-11 |
   --- OH            c 3.7473E-04 3.7473E-04 2.7410E-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.2128673458836E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -3.15022e-07
   --- subroutine FORCE: End Slope       = -8.12561e-09
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787694E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123064E+02  2.123021E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902822E+03  1.902822E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841629E-06  1.841574E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747296E-04  3.747323E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841629E-06  1.841574E-06 -1.517381E+01 -1.517388E+01  5.923910E-05 -5.255330E-05
   ---   H2O           1   9.717735E+01  9.717785E+01 -1.240033E+02 -1.240033E+02 -5.113459E-06 -8.151115E-06
   ---   NaCl          0   3.999601E-29  3.999601E-29 -1.745168E+02 -1.745168E+02 -9.356609E-06 -9.478264E-06
   ---   O2            0   8.937010E-66  8.937010E-66 -1.827721E+02 -1.827721E+02 -1.036977E-03 -1.049781E-03
   ---   Cl-           1   2.123064E+02  2.123021E+02 -7.237743E+01 -7.237746E+01  7.305597E-05  9.498766E-07
   ---   H2            1   1.873648E-04  1.873662E-04 -3.261773E+01 -3.261772E+01 -2.194406E-05 -1.645528E-05
   ---                                                        Norms of Delta G():  4.251988E-04  4.291910E-04
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274259E+03
   ---                  air =   4.0971784E+03
   ---              NaCl(S) =   4.7876979E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458836E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   35
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  35, Iterations since last evaluation of optimal basis =  33 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06   4.8390E-11  -5.2553E-05 | Normal calc: diag adjusted from 1.08603e+06 to 1.08603e+06 due to act coeff
   --- H2O             9.7178E+01   7.9454E-04  -8.1511E-06 | Normal calc: diag adjusted from 0.0101422 to 0.010259 due to act coeff
   --- NaCl            3.9996E-29   3.7909E-34  -9.4783E-06 | Normal Calc
   --- O2              8.9370E-66   9.3819E-69  -1.0498E-03 | Normal Calc
   --- Cl-             2.1230E+02  -5.5612E-05   9.4988E-07 | Normal calc: diag adjusted from 0.0077019 to 0.0170805 due to act coeff
   --- H2              1.8737E-04   1.0277E-09  -1.6455E-05 | Normal calc: diag adjusted from 16011.4 to 16011.4 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06 4.8390E-11 | Normal Major Calc
   --- H2O           1 9.7178E+01 9.7179E+01 7.9454E-04 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.9370E-66 8.9370E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E+02 2.1230E+02-5.5612E-05 | Line Search reduced step size from -5.56118e-05 to -5.56118e-05
   --- H2            1 1.8737E-04 1.8737E-04 1.0277E-09 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 5.5612E-05 |
   --- Na+           c 2.1230E+02 2.1230E+02-5.5612E-05 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-7.9454E-04 |
   --- OH-           c 1.8416E-06 1.8416E-06 4.8390E-11 |
   --- OH            c 3.7473E-04 3.7473E-04 2.0554E-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.2128673458838E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -6.5292e-09
   --- subroutine FORCE: End Slope       = -1.68437e-10
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123021E+02  2.123020E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902822E+03  1.902821E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841574E-06  1.841623E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747323E-04  3.747344E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841574E-06  1.841623E-06 -1.517388E+01 -1.517386E+01 -5.255330E-05  1.249790E-06
   ---   H2O           1   9.717785E+01  9.717865E+01 -1.240033E+02 -1.240033E+02 -8.151115E-06 -1.060134E-07
   ---   NaCl          0   3.999601E-29  3.999601E-29 -1.745168E+02 -1.745168E+02 -9.478264E-06 -9.672187E-06
   ---   O2            0   8.937010E-66  8.937010E-66 -1.827721E+02 -1.827721E+02 -1.049781E-03 -1.071265E-03
   ---   Cl-           1   2.123021E+02  2.123020E+02 -7.237746E+01 -7.237746E+01  9.498766E-07  1.514158E-06
   ---   H2            1   1.873662E-04  1.873672E-04 -3.261772E+01 -3.261772E+01 -1.645528E-05 -4.558975E-07
   ---                                                        Norms of Delta G():  4.291910E-04  4.373609E-04
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274250E+03
   ---                  air =   4.0971792E+03
   ---              NaCl(S) =   4.7876980E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   36
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  36, Iterations since last evaluation of optimal basis =  34 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06  -1.1508E-12   1.2498E-06 | Normal calc: diag adjusted from 1.086e+06 to 1.086e+06 due to act coeff
   --- H2O             9.7179E+01   1.0334E-05  -1.0601E-07 | Normal calc: diag adjusted from 0.0101421 to 0.0102589 due to act coeff
   --- NaCl            3.9996E-29   3.8685E-34  -9.6722E-06 | Normal Calc
   --- O2              8.9370E-66   9.5739E-69  -1.0713E-03 | Normal Calc
   --- Cl-             2.1230E+02  -8.8648E-05   1.5142E-06 | Normal calc: diag adjusted from 0.0077019 to 0.0170805 due to act coeff
   --- H2              1.8737E-04   2.8473E-11  -4.5590E-07 | Normal calc: diag adjusted from 16011.3 to 16011.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06-1.1508E-12 | Line Search reduced step size from -1.15082e-12 to -1.15082e-12
   --- H2O           1 9.7179E+01 9.7179E+01 1.0334E-05 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.9370E-66 8.9370E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E+02 2.1230E+02-8.8648E-05 | Line Search reduced step size from -8.86484e-05 to -8.86484e-05
   --- H2            1 1.8737E-04 1.8737E-04 2.8473E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 8.8648E-05 |
   --- Na+           c 2.1230E+02 2.1230E+02-8.8648E-05 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-1.0334E-05 |
   --- OH-           c 1.8416E-06 1.8416E-06-1.1508E-12 |
   --- OH            c 3.7473E-04 3.7473E-04 5.6947E-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.2128673458838E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -1.35323e-10
   --- subroutine FORCE: End Slope       = -3.49157e-12
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123020E+02  2.123019E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+03  1.902821E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841623E-06  1.841622E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747344E-04  3.747345E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841623E-06  1.841622E-06 -1.517386E+01 -1.517386E+01  1.249790E-06 -1.089206E-06
   ---   H2O           1   9.717865E+01  9.717866E+01 -1.240033E+02 -1.240033E+02 -1.060134E-07 -1.689380E-07
   ---   NaCl          0   3.999601E-29  3.999601E-29 -1.745168E+02 -1.745168E+02 -9.672187E-06 -9.674709E-06
   ---   O2            0   8.937010E-66  8.937010E-66 -1.827721E+02 -1.827721E+02 -1.071265E-03 -1.071532E-03
   ---   Cl-           1   2.123020E+02  2.123019E+02 -7.237746E+01 -7.237746E+01  1.514158E-06  1.969336E-08
   ---   H2            1   1.873672E-04  1.873672E-04 -3.261772E+01 -3.261772E+01 -4.558975E-07 -3.410482E-07
   ---                                                        Norms of Delta G():  4.373609E-04  4.374693E-04
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+03
   ---                  air =   4.0971792E+03
   ---              NaCl(S) =   4.7876981E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   37
   --- 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 =  37, Iterations since last evaluation of optimal basis =  35 (all species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06   1.0030E-12  -1.0892E-06 | Normal calc: diag adjusted from 1.086e+06 to 1.086e+06 due to act coeff
   --- H2O             9.7179E+01   1.6467E-05  -1.6894E-07 | Normal calc: diag adjusted from 0.0101421 to 0.0102589 due to act coeff
   --- NaCl            3.9996E-29   3.8695E-34  -9.6747E-06 | Normal Calc
   --- O2              8.9370E-66   9.5763E-69  -1.0715E-03 | Normal Calc
   --- Cl-             2.1230E+02  -1.1530E-06   1.9693E-08 | Normal calc: diag adjusted from 0.00770191 to 0.0170805 due to act coeff
   --- H2              1.8737E-04   2.1300E-11  -3.4105E-07 | Normal calc: diag adjusted from 16011.3 to 16011.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Full Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06 1.0030E-12 | Normal Major Calc
   --- H2O           1 9.7179E+01 9.7179E+01 1.6467E-05 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 3.8695E-34 | minor species alternative calc
   --- O2            0 8.9370E-66 8.9466E-66 9.5814E-69 | minor species alternative calc
   --- Cl-           1 2.1230E+02 2.1230E+02-1.1530E-06 | Line Search reduced step size from -1.15297e-06 to -1.15297e-06
   --- H2            1 1.8737E-04 1.8737E-04 2.1300E-11 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 1.1530E-06 |
   --- Na+           c 2.1230E+02 2.1230E+02-1.1530E-06 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-1.6468E-05 |
   --- OH-           c 1.8416E-06 1.8416E-06 1.0030E-12 |
   --- OH            c 3.7473E-04 3.7473E-04 4.2601E-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.2128673458838E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -2.8047e-12
   --- subroutine FORCE: End Slope       = -7.23662e-14
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (all species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E+02  2.123019E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+03  1.902821E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841622E-06  1.841623E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-04  3.747345E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841622E-06  1.841623E-06 -1.517386E+01 -1.517386E+01 -1.089206E-06  2.591246E-08
   ---   H2O           1   9.717866E+01  9.717867E+01 -1.240033E+02 -1.240033E+02 -1.689380E-07 -2.197240E-09
   ---   NaCl          0   3.999601E-29  3.999640E-29 -1.745168E+02 -1.745168E+02 -9.674709E-06 -4.019228E-09
   ---   O2            0   8.937010E-66  8.946592E-66 -1.827721E+02 -1.827711E+02 -1.071532E-03 -4.452816E-07
   ---   Cl-           1   2.123019E+02  2.123019E+02 -7.237746E+01 -7.237746E+01  1.969336E-08  3.138233E-08
   ---   H2            1   1.873672E-04  1.873672E-04 -3.261772E+01 -3.261772E+01 -3.410482E-07 -9.449366E-09
   ---                                                        Norms of Delta G():  4.374693E-04  1.825935E-07
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+03
   ---                  air =   4.0971792E+03
   ---              NaCl(S) =   4.7876981E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   38
   --- 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 =  38, Iterations since last evaluation of optimal basis =  36 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06  -2.3860E-14   2.5912E-08 | Normal calc: diag adjusted from 1.086e+06 to 1.086e+06 due to act coeff
   --- H2O             9.7179E+01   2.1418E-07  -2.1972E-09 | Normal calc: diag adjusted from 0.0101421 to 0.0102589 due to act coeff
   --- NaCl            3.9996E-29   1.6075E-37  -4.0192E-09 | Normal Calc
   --- O2              8.9466E-66   3.9838E-72  -4.4528E-07 | Normal Calc
   --- Cl-             2.1230E+02  -1.8373E-06   3.1382E-08 | Normal calc: diag adjusted from 0.00770191 to 0.0170805 due to act coeff
   --- H2              1.8737E-04   5.9017E-13  -9.4494E-09 | Normal calc: diag adjusted from 16011.3 to 16011.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06-2.3860E-14 | Normal Major Calc
   --- H2O           1 9.7179E+01 9.7179E+01 2.1418E-07 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.9466E-66 8.9466E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E+02 2.1230E+02-1.8373E-06 | Normal Major Calc
   --- H2            1 1.8737E-04 1.8737E-04 5.9017E-13 | Line Search reduced step size from 5.90167e-13 to 5.90165e-13
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 1.8373E-06 |
   --- Na+           c 2.1230E+02 2.1230E+02-1.8373E-06 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-2.1418E-07 |
   --- OH-           c 1.8416E-06 1.8416E-06-2.3860E-14 |
   --- OH            c 3.7473E-04 3.7473E-04 1.1803E-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.2128673458838E+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -5.813e-14
   --- subroutine FORCE: End Slope       = -1.49988e-15
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E+02  2.123019E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+03  1.902821E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841623E-06  1.841623E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-04  3.747345E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841623E-06  1.841623E-06 -1.517386E+01 -1.517386E+01  2.591246E-08 -2.257481E-08
   ---   H2O           1   9.717867E+01  9.717867E+01 -1.240033E+02 -1.240033E+02 -2.197240E-09 -3.501398E-09
   ---   NaCl          0   3.999640E-29  3.999640E-29 -1.745168E+02 -1.745168E+02 -4.019228E-09 -4.071524E-09
   ---   O2            0   8.946592E-66  8.946592E-66 -1.827711E+02 -1.827711E+02 -4.452816E-07 -4.508122E-07
   ---   Cl-           1   2.123019E+02  2.123019E+02 -7.237746E+01 -7.237746E+01  3.138233E-08  4.081784E-10
   ---   H2            1   1.873672E-04  1.873672E-04 -3.261772E+01 -3.261772E+01 -9.449366E-09 -7.068564E-09
   ---                                                        Norms of Delta G():  1.825935E-07  1.843096E-07
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+03
   ---                  air =   4.0971792E+03
   ---              NaCl(S) =   4.7876981E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   39
   --- Normal element abundance check - passed
   --- Check for an optimum basis passed
   --- Equilibrium check for major species: H+ failed

==============================================================================================================
 Iteration =  39, Iterations since last evaluation of optimal basis =  37 (only major species)
   ----------------------------------------------------------------------------------
   --- Subroutine vcs_RxnStepSizes called - Details:
   ----------------------------------------------------------------------------------
   --- Species         Moles     Rxn_Adjustment    DeltaG   | Comment
   --- H+              1.8416E-06   2.0787E-14  -2.2575E-08 | Normal calc: diag adjusted from 1.086e+06 to 1.086e+06 due to act coeff
   --- H2O             9.7179E+01   3.4130E-07  -3.5014E-09 | Normal calc: diag adjusted from 0.0101421 to 0.0102589 due to act coeff
   --- NaCl            3.9996E-29   1.6285E-37  -4.0715E-09 | Normal Calc
   --- O2              8.9466E-66   4.0332E-72  -4.5081E-07 | Normal Calc
   --- Cl-             2.1230E+02  -2.3897E-08   4.0818E-10 | Normal calc: diag adjusted from 0.00770191 to 0.0170805 due to act coeff
   --- H2              1.8737E-04   4.4147E-13  -7.0686E-09 | Normal calc: diag adjusted from 16011.3 to 16011.3 due to act coeff
   ----------------------------------------------------------------------------------
   --- Main Loop Treatment of each non-component species - Major Components Calculation:
   --- Species     IC     Moles   Tent_Moles  Rxn_Adj   |    Comment 
   --- H+            1 1.8416E-06 1.8416E-06 2.0787E-14 | Normal Major Calc
   --- H2O           1 9.7179E+01 9.7179E+01 3.4130E-07 | Normal Major Calc
   --- NaCl          0 3.9996E-29 3.9996E-29 0.0000E+00 | minor species not considered
   --- O2            0 8.9466E-66 8.9466E-66 0.0000E+00 | minor species not considered
   --- Cl-           1 2.1230E+02 2.1230E+02-2.3897E-08 | Line Search reduced step size from -2.38973e-08 to -2.38965e-08
   --- H2            1 1.8737E-04 1.8737E-04 4.4147E-13 | Normal Major Calc
   --- N2            c 4.0000E+03 4.0000E+03 0.0000E+00 |
   --- NaCl(S)       c 4.7877E+03 4.7877E+03 2.3897E-08 |
   --- Na+           c 2.1230E+02 2.1230E+02-2.3897E-08 |
   --- H2O(L)        c 1.9028E+03 1.9028E+03-3.4131E-07 |
   --- OH-           c 1.8416E-06 1.8416E-06 2.0787E-14 |
   --- OH            c 3.7473E-04 3.7473E-04 8.8295E-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+06
   --- Total tentative Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- subroutine FORCE: Beginning Slope = -1.2048e-15
   --- subroutine FORCE: End Slope       = -3.10832e-17
   --- subroutine FORCE produced no adjustments, s2 < 0
   -------------------------------------------------------------------------------------------------------
   --- Summary of the Update  (only major species):
   ---      Species Status Initial_Moles  Final_Moles  Initial_Mu/RT     Mu/RT     Init_Del_G/RT   Delta_G/RT
   ---   N2                4.000000E+03  4.000000E+03 -2.306923E+01 -2.306923E+01
   ---   NaCl(S)           4.787698E+03  4.787698E+03 -1.745168E+02 -1.745168E+02
   ---   Na+               2.123019E+02  2.123019E+02 -1.021393E+02 -1.021393E+02
   ---   H2O(L)            1.902821E+03  1.902821E+03 -1.240033E+02 -1.240033E+02
   ---   OH-               1.841623E-06  1.841623E-06 -1.088294E+02 -1.088294E+02
   ---   OH                3.747345E-04  3.747345E-04 -1.076944E+02 -1.076944E+02
   ---   CO2               0.000000E+00  0.000000E+00 -5.151318E+02 -5.151318E+02
   ---   H+            1   1.841623E-06  1.841623E-06 -1.517386E+01 -1.517386E+01 -2.257481E-08  5.370708E-10
   ---   H2O           1   9.717867E+01  9.717867E+01 -1.240033E+02 -1.240033E+02 -3.501398E-09 -4.553158E-11
   ---   NaCl          0   3.999640E-29  3.999640E-29 -1.745168E+02 -1.745168E+02 -4.071524E-09 -4.154828E-09
   ---   O2            0   8.946592E-66  8.946592E-66 -1.827711E+02 -1.827711E+02 -4.508122E-07 -4.600411E-07
   ---   Cl-           1   2.123019E+02  2.123019E+02 -7.237746E+01 -7.237746E+01  4.081784E-10  6.504308E-10
   ---   H2            1   1.873672E-04  1.873672E-04 -3.261772E+01 -3.261772E+01 -7.068564E-09 -1.958256E-10
   ---                                                        Norms of Delta G():  1.843096E-07  1.878190E-07
   ---           Phase_Name     Moles(after update)
   ---   --------------------------------------------------
   ---     NaCl_electrolyte =   2.3274248E+03
   ---                  air =   4.0971792E+03
   ---              NaCl(S) =   4.7876981E+03
   -------------------------------------------------------------------------------------------------------
   --- Total Old Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Total New Dimensionless Gibbs Free Energy = -1.2128673458838E+06
   --- Increment counter increased, step is accepted:   40
   --- 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:  CONVERGENCE achieved
   --- Subroutine vcs_deltag called for all noncomponents
   --- 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:  CONVERGENCE achieved
   --- Check the Full Element Abundances:  passed
   --- Start rechecking deleted species in multispec phases




--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
		 VCS_TP REPORT
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
		Temperature =           3e+02 Kelvin
		Pressure    =      1.0132e+05 Atmos
		Volume      =      1.0041e+05 cm**3


--------------------------------------------------------------------------------
 Species                 Equilibrium moles    Mole Fraction    ChemPot/RT    SpecUnkType
--------------------------------------------------------------------------------
 N2                        4.0000000E+03      9.7628143E-01     -2.3069E+01     0
 NaCl(S)                   4.7876981E+03      1.0000000E+00     -1.7452E+02     0
 Na+                       2.1230193E+02      9.1217525E-02     -1.0214E+02     0
 H2O(L)                    1.9028209E+03      8.1756495E-01     -1.2400E+02     0
 OH-                       1.8416226E-06      7.9127051E-10     -1.0883E+02     0
 OH                        3.7473450E-04      9.1461582E-08     -1.0769E+02     0
 CO2                       0.0000000E+00      0.0000000E+00     -5.1513E+02     0
 Cl-                       2.1230193E+02      9.1217525E-02     -7.2377E+01   MolNum 
 H2O                       9.7178674E+01      2.3718434E-02     -1.2400E+02   MolNum 
 H2                        1.8736725E-04      4.5730791E-08     -3.2618E+01   MolNum 
 H+                        1.8416226E-06      7.9127051E-10     -1.5174E+01   MolNum 
 NaCl                      3.9996395E-29      9.7619344E-33     -1.7452E+02   MolNum 
 O2                        8.9465918E-66      2.1835979E-69     -1.8277E+02   MolNum 
--------------------------------------------------------------------------------


-------------------------------------------------------------------------------------------------------------------
               |ComponentID|          0          1          2          3          4          5          6 |           |
               | Components|         N2    NaCl(S)        Na+     H2O(L)        OH-         OH        CO2 |           |
 NonComponent  |   Moles   |      4e+03   4.79e+03        212    1.9e+03   1.84e-06   0.000375          0 | DG/RT Rxn |
-------------------------------------------------------------------------------------------------------------------
   7 H+        |  1.84e-06 |       0.00       0.00       0.00      -1.00       1.00       0.00       0.00 |  5.37e-10 |
   8 H2O       |      97.2 |       0.00       0.00       0.00      -1.00       0.00       0.00       0.00 | -4.55e-11 |
   9 NaCl      |     4e-29 |       0.00      -1.00       0.00       0.00       0.00       0.00       0.00 | -4.15e-09 |
  10 O2        |  8.95e-66 |       0.00       0.00       0.00       2.00       0.00      -4.00       0.00 |  -4.6e-07 |
  11 Cl-       |       212 |       0.00      -1.00       1.00       0.00       0.00       0.00       0.00 |   6.5e-10 |
  12 H2        |  0.000187 |       0.00       0.00       0.00      -2.00       0.00       2.00       0.00 | -1.96e-10 |
-------------------------------------------------------------------------------------------------------------------




------------------------------------------------------------------------------------------------------------------------------------------------------------------------
                  | 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     | MolTarget |      2e+03      4e+03          0      8e+03      5e+03      5e+03          0          0          0          0          0 |     Gibbs Total     |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
   0 NaCl_electro | 2.327e+03 |    1.9e+03   3.81e+03          0          0        212        212          0          0          0          0          0 | -2.73006234953E+05 |
   1 air          | 4.097e+03 |       97.2        194          0      8e+03      4e-29      4e-29          0          0          0          0          0 | -1.04327434127E+05 |
   2 NaCl(S)      | 4.788e+03 |          0          0          0          0   4.79e+03   4.79e+03          0          0          0          0          0 | -8.35533676804E+05 |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    TOTAL         | 1.121e+04 |      2e+03      4e+03          0      8e+03      5e+03      5e+03          0          0          0          0          0 | -1.21286734588E+06 |
------------------------------------------------------------------------------------------------------------------------------------------------------------------------


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

Elemental Abundances:              Actual                    Target         Type      ElActive
                          O   2.000000000000E+03    2.000000000000E+03     0       1
                          H   4.000000000000E+03    4.000000000000E+03     0       1
                          C   0.000000000000E+00    0.000000000000E+00     0       1
                          N   8.000000000000E+03    8.000000000000E+03     0       1
                          Na  5.000000000000E+03    5.000000000000E+03     0       1
                          Cl  5.000000000000E+03    5.000000000000E+03     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         TMoles     StandStateChemPot      ln(AC)       ln(X_i)      |   F z_i phi   |    ChemPot    | (-lnMnaught)
-------------------------------------------------------------------------------------------------------------------
 N2            4.0000000E+03 -2.3045225E+01   0.0000000E+00  -2.4004385E-02  | 0.0000000E+00 |   -2.3069E+01 |
 NaCl(S)       4.7876981E+03 -1.7451679E+02   0.0000000E+00   0.0000000E+00  | 0.0000000E+00 |   -1.7452E+02 |
 Na+           2.1230193E+02 -1.0397591E+02   2.1455966E-01  -2.3945082E+00  | 0.0000000E+00 |   -1.0214E+02 | ( 4.0165350E+00)
 H2O(L)        1.9028209E+03 -1.2371551E+02  -8.6317426E-02  -2.0142493E-01  | 0.0000000E+00 |   -1.2400E+02 |
 OH-           1.8416226E-06 -9.1483483E+01  -4.0506365E-01  -2.0957381E+01  |-0.0000000E+00 |   -1.0883E+02 | ( 4.0165350E+00)
 OH            3.7473450E-04 -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.3067997E+02  | 0.0000000E+00 |   -5.1513E+02 |
 Cl-           2.1230193E+02 -7.4214051E+01   2.1455966E-01  -2.3945082E+00  |-0.0000000E+00 |   -7.2377E+01 | ( 4.0165350E+00)
 H2O           9.7178674E+01 -1.2026175E+02   0.0000000E+00  -3.7415027E+00  | 0.0000000E+00 |   -1.2400E+02 |
 H2            1.8736725E-04 -1.5717224E+01   0.0000000E+00  -1.6900494E+01  | 0.0000000E+00 |   -3.2618E+01 |
 H+            1.8416226E-06  0.0000000E+00   1.7669883E+00  -2.0957381E+01  | 0.0000000E+00 |   -1.5174E+01 | ( 4.0165350E+00)
 NaCl          3.9996395E-29 -1.0080997E+02   0.0000000E+00  -7.3706817E+01  | 0.0000000E+00 |   -1.7452E+02 |
 O2            8.9465918E-66 -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:          40                 NA     
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------

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