Merge pull request #946 from GEOS-ESM/bugfix/rreichle/hcorr_double_co…

…unting bug fix: remove double-counting of heat corr (StieglitzSnow.F90)
GEOS-ESM · Jul 10, 2024 · f5fd2a3 · f5fd2a3
2 parents 8bfa81a + b02a43c
commit f5fd2a3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Shared/StieglitzSnow.F90 b/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