fwflxFac = 0. _d 0
IF ( tLoc(I,J) .GT. thetaFreeze ) fwflxFac = dFac
& *shiTransCoeffS(i,j,bi,bj) * fwflxFac
& - freshWaterFlux*SHELFICEheatCapacity_Cp*fwflxFac )
_RL cFac, rFac, dFac, fwflxFac