# React script, saved Tue May 22 2012 by bfarrell data = thermo.dat verify time start = 0 day, end = 20 day temperature = 10 C decouple Fe+++ H2O = 1 free kg swap Fe(OH)3(ppd) for Fe+++ Fe(OH)3(ppd) = 1 free ug pH = 7 swap O2(g) for O2(aq) O2(g) = .2 fugacity Na+ = 1 mmolal balance on Cl- Cl- = 1 mmolal Fe++ = 30 mmolal fix pH fix fugacity of O2(g) kinetic redox-1 rxn = "Fe++ + .25*O2(aq) + 2.5*H2O -> Fe(OH)3(ppd) + 2*H+" rate_con = 7e12 mpower(Fe++) = 1 mpower(O2(aq)) = 1 apower(OH-) = 2 suppress Hematite Goethite Magnetite