newIceFrac=MIN(newIce/iceThick,iceMaskmax-iceFraction)
compact = iceFraction + newIceFrac
sst=(1. _d 0-newIceFrac)*sst+newIceFrac*Tf
qleft= iceThick*newIceFrac*qicAv/deltaTice
fresh=-(rhoi*iceThick)*newIceFrac/deltaTice
fsalt=-(rhoi*iceThick*saltice)*newIceFrac/deltaTice
& iceThick, newIce, newIceFrac
_RL newIceFrac
newIceFrac = 0. _d 0
newIceFrac = MIN(newIce/himin0,iceMaskmax)
compact = newIceFrac