_RL htotal
htotal = x
htotal2=htotal*htotal
co2star=diclocal*htotal2/(htotal2 + k1local*htotal
phlocal=-log10(htotal)