data = thermo.dat verify surface_data = FeOH+.dat decouple Fe+++ swap Fe(OH)3(ppd) for Fe+++ 1 free gram Fe(OH)3(ppd) 100 mmolal Na+ 100 mmolal Cl- 10 mmolal Ca++ pH = 7 fix pH kinetic >(w)FeOCa+ 1500 umolal, rate_con = 2.e-4, mpower(>(w)FeOCa+) = 1 precip = off time end 12 hours suffix _Sorbkin