CALL DIAGNOSTICS_SCALE_FILL(surfaceForcingTice,tmpFac,1,
_RL tmpFac
tmpFac = HeatCapacity_Cp*recip_horiVertRatio*rhoConst
tmpFac = 0.
& tmpFac = (0.5+abEps) / implicDiv2Dflow
& + tmpFac*(etaN(i,j,bi,bj)-etaH(i,j,bi,bj))