bug fix: remove double-counting of heat corr (StieglitzSnow.F90)

In the rare circumstance that more snow would need to be sublimated than is available in the top layer, the residual energy accounting term was added twice, once here and again ~30 lines further down.

Author: gmao-rreichle

GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Shared/StieglitzSnow.F90

@@ -671,7 +671,7 @@ subroutine StieglitzSnow_snowrt(tile_lon, tile_lat,                            &
           if(-dw > wesn(1) ) then
              dw = -wesn(1)
              evap = -dw/(dts*areasc)
-             hcorr=hcorr+(lhflux-evap*alhv)*areasc
+             !hcorr=hcorr+(lhflux-evap*alhv)*areasc  ! removed, was double-counting corr, see "Store excess heat in hcorr" below; koster+reichle, 31 May 2024
              lhflux=evap*alhv
           endif
           wesn(1)  = wesn(1) + dw