I t,s,diclocal,pt,sit,ta,
diclocal=diclocal*permil
fn = k1local*x*diclocal/b +
& 2.0*diclocal*k12/b +
df = ((k1local*diclocal*b) - k1local*x*diclocal*db)/b2 -
& 2.0*diclocal*k12*db/b2 -
co2star=diclocal*htotal2/(htotal2 + k1local*htotal
diclocal=diclocal/permil
_RL pCO2surfloc, diclocal, pHlocal
gamm = diclocal/cag
co2s = diclocal/