TITLE Example from GWB User's Guide
SURFACE 1
        -equilibrate with solution 1
        Hfo_sOH        6.262e-5    600.    1.0
        Hfo_wOH        2.505e-3
EQUILIBRIUM_PHASES 1
    Hematite    0.  0.006262
SOLUTION 1
     -units mg/kgw
     ph     5.3
     Ca     2.4
     Mg     3.5
     Ba     .04
     Na     0.4
     Zn     1.2
     S(6)     15.
     Cl      8     charge
SELECTED_OUTPUT
     -file  sel_HFO.out
     -reset false
USER_PUNCH
     -headings  Ba  Ca SO4  Zn
     -start
     10 sorbed_Ba = mol("Hfo_sOHBa+2")  + mol("Hfo_wOBa+")
     20 sorbed_Ca = mol("Hfo_sOHCa+2")  + mol("Hfo_wOCa+")
     30 sorbed_SO4= mol("Hfo_wOHSO4-2") + mol("Hfo_wSO4-")
     40 sorbed_Zn = mol("Hfo_sOZn+")    + mol("Hfo_wOZn+")
     50 frac_Ba   = sorbed_Ba/ (sorbed_Ba + TOT("Ba"))
     60 frac_Ca   = sorbed_Ca/ (sorbed_Ca + TOT("Ca"))
     70 frac_SO4  = sorbed_SO4/(sorbed_SO4 + TOT("S(6)"))
     80 frac_Zn   = sorbed_Zn/ (sorbed_Zn + TOT("Zn"))
     90 PUNCH frac_Ba  frac_Ca  frac_SO4  frac_Zn
     -end
END