_RL co2star
co2star=diclocal*htotal2/(htotal2 + k1local*htotal
fco2 = co2star / k0local
pCO2surfloc = co2star / fflocal