_RL c1Frz, c2Frz, aFrz, hFrz
hFrz = hThinIce+c2Frz/aFrz
hFrz = hThickIce
xxFrz = aFrz*(hIce(i,j)-hFrz)