_RL xxFrz, tmpFrz
xxFrz = aFrz*(hIce(i,j)-hFrz)
xxFrz = xxFrz**powerLaw
tmpFrz = xxFrz
xxFrz = xxFrz*tmpFrz
xxFrz = fracEnFreez-xxFrz
enFrcFrz(i,j) = MAX( 0. _d 0, MIN( xxFrz, 1. _d 0 ) )