thetaFreeze = a0*saltFreeze + eps4
& eps1 * ( thetaFreeze - tLoc ) /
& (thetaFreeze - thetaIce))
& FreshWaterFlux(I,J) ) * ( thetaFreeze - tLoc )
_RL thetaFreeze, saltFreeze
thetaFreeze =
& SHELFICEheatTransCoeff * ( tLoc(I,J) - thetaFreeze )
& * ( thetaFreeze - tLoc(I,J) )
& ( eps3*( thetaFreeze - SHELFICEthetaSurface )
& - HeatCapacity_Cp*( thetaFreeze - rFac*tLoc(I,J) ) )