newIceFrac=MIN(newIce/iceThick,iceMaskmax-iceFraction)
compact = iceFraction + newIceFrac
snowThick = snowThick*iceFraction/compact
_RL iceFraction
iceFraction = compact
IF (iceFraction.EQ.0. _d 0) THEN