# React script, saved Tue May 22 2012 by bfarrell data = thermo.dat verify surface_data = FeOH.dat temperature = 25 C decouple Fe+++ H2O = 1 free kg swap Fe(OH)3(ppd) for Fe+++ Fe(OH)3(ppd) = 1 free g pH = 4 Na+ = 100 mmolal balance on Cl- Cl- = 100 mmolal Ca++ = 10 mmolal SO4-- = 10 mmolal Hg++ = .1 mmolal Zn++ = .1 mmolal Cr+++ = .1 mmolal As(OH)4- = .1 mmolal slide pH to 12 suppress Hematite Goethite precip = off