newIceFrac=MIN(newIce/iceThick,iceMaskmax-iceFraction)
& iceThick, newIce, newIceFrac
_RL newIce
newIce = esurp*deltaTice/qicAv
IF (newIce.GE.himin0*iceMaskmin) THEN
iceThick = MAX(himin0,newIce/iceMaskmax)
newIceFrac = MIN(newIce/himin0,iceMaskmax)