& *ICEFRONTheatTransCoeff(I,J)
eps1 = rUnit2mass*HeatCapacity_Cp*ICEFRONTheatTransCoeff(I,J)
& ( -rUnit2mass*ICEFRONTheatTransCoeff(I,J) +
_RS ICEFRONTheatTransCoeff ( 1:sNx, 1:sNy )
ICEFRONTheatTransCoeff(I,J) = 1.0 _d -02
ICEFRONTheatTransCoeff(I,J) = max
& (ICEFRONTheatTransCoeff(I,J),1. _d -04)