tLoc = theta(I,J,K,bi,bj)
& ( eps1 * ( tLoc - thetaFreeze ) )*recip_horiVertRatio
_RL tLoc, sLoc, pLoc
& SHELFICEheatTransCoeff * ( tLoc - thetaFreeze )
I tLoc, sLoc, pLoc,
_RL sLoc, tLoc, pLoc
t1 = tLoc
& -tAlpha*(tLoc-tRef(1))
_RL FUNCTION STATE1(sLoc,tLoc,I,J,KREF,bi,bj,mythid)
_RL tLoc,sLoc
& - tAlpha*(tLoc-tRef(KREF)) ) + dRho