cFac = 0. _d 0
IF ( SHELFICEconserve ) cFac = 1. _d 0
& - cFac * shelfIceFreshWaterFlux(I,J,bi,bj)*mass2rUnit
& * ( cFac*sLoc(I,J) + (1. _d 0-cFac)*convertFW2SaltLoc )
print *, 'ml-shelfice: cFac = ', cFac
& - cFac*freshWaterFlux*( SHELFICElatentHeat
& - cFac*shelfIceFreshWaterFlux(I,J,bi,bj)*mass2rUnit )
_RL cFac, rFac, dFac, fwflxFac