From a08e189ca35242dc755ef96e49212d65c87ffcdd Mon Sep 17 00:00:00 2001 From: einaraz Date: Tue, 27 May 2025 15:55:41 -0600 Subject: [PATCH] fix units of evaporation for imp_scheme==2 in slucm --- phys/module_sf_urban.F | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/phys/module_sf_urban.F b/phys/module_sf_urban.F index 048d27dcac..a25e7d899c 100644 --- a/phys/module_sf_urban.F +++ b/phys/module_sf_urban.F @@ -1032,7 +1032,8 @@ SUBROUTINE urban(LSOLAR, & ! L IF (IMP_SCHEME==2) then IF (FLXHUMRP <= 0.) FLXHUMRP = 0. ! Compute water retention depth from previous time step - DrelR = DrelRP+(RAIN1-FLXHUMRP)*DELT/porimp(IMPR) + ! Convert kinematic water flux to evaporation in m/s: multiply flux by rho_air/who_water in SI units + DrelR = DrelRP+(RAIN1-FLXHUMRP*RHOO/1000.)*DELT/porimp(IMPR) IF (RAIN > 0. .AND. DrelR < DrelRP) DrelR = DrelRP IF (DrelR <= 0.) then @@ -1284,9 +1285,10 @@ SUBROUTINE urban(LSOLAR, & ! L IF (FLXHUMBP <= 0.) FLXHUMBP = 0. IF (FLXHUMGP <= 0.) FLXHUMGP = 0. ! Compute water retention from previous time step for wall and ground - DrelB = DrelBP+(RAIN1-FLXHUMBP)*DELT/porimp(IMPB) + ! Convert kinematic water flux to evaporation in m/s: multiply flux by rho_air/who_water in SI units + DrelB = DrelBP+(RAIN1-FLXHUMBP*RHOO/1000.)*DELT/porimp(IMPB) IF (RAIN > 0. .AND. DrelB < DrelBP) DrelB = DrelBP - DrelG = DrelGP+(RAIN1-FLXHUMGP)*DELT/porimp(IMPG) + DrelG = DrelGP+(RAIN1-FLXHUMGP*RHOO/1000.)*DELT/porimp(IMPG) IF (RAIN > 0. .AND. DrelG < DrelGP) DrelG = DrelGP IF (DrelB <= 0.) then