_RL SurfHeatFluxConvergToSnowMelt (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
SurfHeatFluxConvergToSnowMelt (I,J) = 0.0
SurfHeatFluxConvergToSnowMelt(I,J) =
SurfHeatFluxConvergToSnowMelt(I,J) =F_ia_net(I,J)
& F_ia_net(I,J)-SurfHeatFluxConvergToSnowMelt(I,J)