& PotSnowMeltRateFromSurf(I,J),
_RL PotSnowMeltRateFromSurf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
PotSnowMeltRateFromSurf (I,J) = 0.0
PotSnowMeltRateFromSurf(I,J)= - F_ia_net(I,J)*QS
& PotSnowMeltRateFromSurf(I,J)* SEAICE_deltaTtherm
& PotSnowMeltRateFromSurf(I,J)