newIceFrac=MIN(newIce/iceThick,iceMaskmax-iceFraction)
qleft= iceThick*newIceFrac*qicAv/deltaTice
fresh=-(rhoi*iceThick)*newIceFrac/deltaTice
fsalt=-(rhoi*iceThick*saltice)*newIceFrac/deltaTice
& iceThick, newIce, newIceFrac
U sst, compact, iceThick, snowThick, qicen,
_RL iceThick
iceThick = MAX(himin0,newIce/iceMaskmax)