Solution 1 -units mg/L pH 5.7 Na 1.6 K 3.2 Mg 0.2 Ca 3.3 Cl 5.9 S(6) 7.1 Al 0.001 Fe(2) 0.001 Mn 0.001 C(4) 0.8EQUILIBRIUM_PHASES# K-mica 0 0# Quartz 0 0.001 Chalcedony 0 0 Chlorite(14A) 0 0 Gibbsite 0 0 Kaolinite 0 0Kinetics 1Quartz -parms 2.740e-15 -m0 1810.70 -m 1810.70 -tol 1e-3K-mica# Should be modified as Muscovite -parms 1.196e-14 -m0 77.98 -m 77.98 -tol 1e-3K-feldspar -parms 1.269e-13 -m0 293.44 -m 293.44 -tol 1e-3Albite -parms 8.298e-14 # m^2/mol Albite, fraction adjusts lab rate to field rate -m0 327.15 -m 327.15 -tol 1e-3 -steps 1 5 10 15 20 25 35 45 55 100 200 300 days -bad_step_max 5000INCREMENTAL_REACTIONS trueRATESQuartz -start 10 REM M = Current number of moles of Quartz 20 SI_Q = SI("Quartz") 30 IF ( ABS(SI_Q) < 1.e-1 ) THEN GOTO 100 40 Rate = m0*PARM(1)*(1-SR("Quartz")) 50 Moles = Rate * Time 100 SAVE moles -endK-mica -start 10 REM M = current number of moles of K-mica 20 si_km = SI("K-mica") 30 IF ( abs(si_km) < 1.e-1 ) Then goto 100 40 Rate = m0*parm(1)*(1-SR("K-mica")) 50 moles = Rate * Time 100 SAVE moles -endK-feldspar -start 10 REM M = Current number of moles of K-feldspar 20 SI_Kf = SI("K-feldspar") 30 IF ( M <= 0 or ABS(SI_Kf) < 1.e-1 ) THEN GOTO 100 40 Rate = m0*parm(1)*(1-SR("K-feldspar")) 50 Moles = Rate * TIME 100 SAVE moles -endAlbite -start 10 REM M = current number of moles of Albite 20 si_Al = SI("Albite") 30 IF ( M <= 0 or ABS(si_Al) < 1.e-1 ) THEN GOTO 100 40 Rate = m0*PARM(1)*(1-SR("Albite")) 50 Moles = Rate * TIME 100 SAVE Moles -endEND