fix units

GEOS-ESM · Jul 5, 2024 · d95a4a1 · d95a4a1
1 parent afdc6ad
commit d95a4a1
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/...GridComp/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatch_GridComp/GEOS_CatchGridComp.F90 b/...GridComp/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatch_GridComp/GEOS_CatchGridComp.F90
@@ -5673,13 +5673,13 @@ subroutine Driver ( RC )
 
         SUMEV = EVPICE+EVPSOI+EVPVEG+EVPINT
         where(SUMEV/=0.)                          ! apportion residual based on (non-zero) evap/dewfall flux
-           EVPICE = EVPICE - EVACC*EVPICE/SUMEV
-           EVPSOI = EVPSOI - EVACC*EVPSOI/SUMEV
-           EVPINT = EVPINT - EVACC*EVPINT/SUMEV
-           EVPVEG = EVPVEG - EVACC*EVPVEG/SUMEV
+           EVPICE = EVPICE - EVACC*MAPL_ALHS*EVPICE/SUMEV 
+           EVPSOI = EVPSOI - EVACC*MAPL_ALHL*EVPSOI/SUMEV
+           EVPINT = EVPINT - EVACC*MAPL_ALHL*EVPINT/SUMEV
+           EVPVEG = EVPVEG - EVACC*MAPL_ALHL*EVPVEG/SUMEV
         elsewhere                                 ! apportion residual based on ASNOW (for simplicity, add snow-free portion to EVPSOI)
-           EVPICE = EVPICE - EVACC*ASNOW           
-           EVPSOI = EVPSOI - EVACC*(1.-ASNOW)
+           EVPICE = EVPICE - EVACC*MAPL_ALHS*ASNOW
+           EVPSOI = EVPSOI - EVACC*MAPL_ALHL*(1.-ASNOW)
         endwhere
 
 !        endif