_RL hguess
hguess = 10.0**(-phguess)
bohg = btlocal*kblocal/(hguess+kblocal)
stuff = hguess*hguess*hguess
& + (k1plocal*hguess*hguess)
& + (k1plocal*k2plocal*hguess)
h3po4g = (pt*hguess*hguess*hguess) / stuff
h2po4g = (pt*k1plocal*hguess*hguess) / stuff
hpo4g = (pt*k1plocal*k2plocal*hguess) / stuff
siooh3g = sit*ksilocal / (ksilocal + hguess)
cag = ta - bohg - (kwlocal/hguess) + hguess