From bff4e2919682fc11c7bcc5764d2cafe601249f61 Mon Sep 17 00:00:00 2001 From: drnimbusrain Date: Wed, 15 May 2024 02:18:10 +0000 Subject: [PATCH] Updated rca calculation with rsmin and rsmax limits. --- physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 b/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 index f2f6b6941..db3ae472f 100644 --- a/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 +++ b/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 @@ -1063,14 +1063,14 @@ subroutine noahmpdrv_run & zorl (i) = z0_total * 100.0 ! convert to cm ztmax (i) = z0h_total - ! total stomatal/canopy resistance Based on Bonan et al. (2011) conductance (1/Rs) equation + !LAI-scale canopy resistance based on weighted sunlit shaded fraction if(rs_sunlit .le. 0.0 .or. rs_shaded .le. 0.0 .or. & lai_sunlit .eq. 0.0 .or. lai_shaded .eq. 0.0) then - rca(i) = 0.0 - else + rca(i) = parameters%rsmax + else !calculate LAI-scale canopy conductance (1/Rs) rca(i) = ((1.0/(rs_sunlit+leaf_air_resistance)*lai_sunlit) + & ((1.0/(rs_shaded+leaf_air_resistance))*lai_shaded)) - rca(i) = 1.0/rca(i) !resistance + rca(i) = max((1.0/rca(i)),parameters%rsmin) !resistance end if smc (i,:) = soil_moisture_vol