zdum = (pressc*10.0d0 - 10.0d0)/10.0d0
_RL pressc
pressc = 1.0d0
pressc = 1.0d0 + 0.1d0*cdepth
& exp( (24.2-0.085*t)*(pressc-1.0)/(83.143*tk) )
& exp( (16.4-0.040*t)*(pressc-1.0)/(83.143*tk) )
pfactor = - (dv/(bigR*tk))*pressc
& + (0.5*dk/(bigR*tk))*pressc*pressc