data = thermo.dat verify time end = 20 days T = 10 decouple Fe+++ suppress Hematite, Goethite, Magnetite swap Fe(OH)3(ppd) for Fe+++ 1 free ug Fe(OH)3(ppd) pH = 7; fix pH swap O2(g) for O2(aq) f O2(g) = 0.2; fix f O2(g) Na+ = 1 mmolal Cl- = 1 mmolal Fe++ = 30 mmolal kinetic redox-1 rxn = \ "Fe++ + 1/4 O2(aq) + 5/2 H2O -> Fe(OH)3(ppd) + 2 H+" \ mpower(Fe++) = 1, mpower(O2(aq)) = 1, apower(OH-) = 2 \ rate_con = 7e12 suffix _Redoxkin