# React script, saved Tue May 22 2012 by bfarrell data = thermo.dat verify surface_data = FeOH+.dat time start = 0 day, end = 12 hr temperature = 25 C decouple Fe+++ H2O = 1 free kg swap Fe(OH)3(ppd) for Fe+++ Fe(OH)3(ppd) = 1 free g Na+ = 100 mmolal balance on Cl- Cl- = 100 mmolal Ca++ = 10 mmolal pH = 7 fix pH kinetic >(w)FeOCa+ 1500 umolal rate_con = 2e-4 mpower(>(w)FeOCa+) = 1 precip = off