DATABASE C:\PROGRA~1\USGS\PHREEQ~1.7\phreeqc.dat
SOLUTION_SPREAD
    -units    mg/l
       Description	 Number	 temp	   pH	  Ca	   Mg	   Na	    K	 Alkalinity	     Cl	  S(6)	    Br	    Si	 Fe(2)	   Mn	   Sr	   Ba	   N(5)	  N(-3)
                 	       	     	     	    	     	     	     	    as HCO3	       	      	      	      	      	     	     	     	 as NO3	 as NH4
        Background	      1	 16.6	 7.01	 166	 52.6	 91.8	  2.5	        626	 180.98	 113.6	 0.793	 18.60	  0.13	 1.94	 1.06	 0.14	   0.05	    2.1
 Most contaminated	      2	 16.6	 6.78	 514	  232	  606	 14.4	       2642	 1026.6	   0.1	 7.536	 35.90	 19.30	 0.90	 9.87	 7.12	   3.48	   15.0
END
RATES
    Sulfate_reduction
-start
10 if (m <= 0) then goto 100
20 mSO4 = tot("S(6)")
30 rate = 9.606e-12*mSO4/(9.259e-10 + mSO4)
40 moles = rate * time
50 if (moles > m) then moles = m
100 save moles
-end
END
USE solution 2
SELECTED_OUTPUT
    -file                 OK-3-norman-sulf-kinetics-sel
    -totals               S(6)  S(-2)  C(4)
    -kinetic_reactants    Sulfate_reduction
KINETICS 1
Sulfate_reduction
    -formula  CH2O  2
    -m        0.01325
    -m0       1
    -tol      1e-008
-steps       1 3150000 31500000 315000000 630700000 1260000000 1892000000 2523000000 3150000000
-step_divide 1
-runge_kutta 3
END