_RL stuff
stuff = hguess*hguess*hguess
h3po4g = (pt*hguess*hguess*hguess) / stuff
h2po4g = (pt*k1plocal*hguess*hguess) / stuff
hpo4g = (pt*k1plocal*k2plocal*hguess) / stuff
po4g = (pt*k1plocal*k2plocal*k3plocal) / stuff
stuff = (1.0 _d 0-gamm)*(1.0 _d 0-gamm)*k1local*k1local
hnew = 0.5 _d 0*( (gamm-1.0 _d 0)*k1local + sqrt(stuff) )