freshWaterFlux = rUnit2mass
& - freshWaterFlux*SHELFICEheatCapacity_Cp*fwflxFac )
& - cFac*freshWaterFlux*( SHELFICElatentHeat
shelfIceFreshWaterFlux(I,J,bi,bj) = freshWaterFlux
_RL freshWaterFlux, convertFW2SaltLoc