From d98defadb2d13d64204188a6ebba2049a99de87c Mon Sep 17 00:00:00 2001 From: Michael Barlage Date: Mon, 17 Apr 2023 17:26:53 -0400 Subject: [PATCH] modify a few terms in vegetation diag3 --- physics/module_sf_noahmplsm.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/physics/module_sf_noahmplsm.F90 b/physics/module_sf_noahmplsm.F90 index 65b3d77b3..e28b5447e 100644 --- a/physics/module_sf_noahmplsm.F90 +++ b/physics/module_sf_noahmplsm.F90 @@ -4339,15 +4339,15 @@ subroutine vege_flux(parameters,nsnow ,nsoil ,isnow ,vegtyp ,veg , & fhi = fh2/fh wrk = 1.0 - fhi if(thsfc_loc) then ! Use local potential temperature - t2mv = tv*wrk + sfctmp*prslkix*fhi - (grav+grav)/cp + t2mv = tah*wrk + sfctmp*prslkix*fhi - (grav+grav)/cp else ! Use potential temperature referenced to 1000 hPa - t2mv = tv*wrk + sfctmp*fhi - (grav+grav)/cp + t2mv = tah*wrk + sfctmp*fhi - (grav+grav)/cp endif - if(evg >= 0.) then ! for evaporation>0, use inferred qsurf to deduce q2v + if((evc+tr)/fveg+evg >= 0.) then ! for evaporation>0, use inferred qsurf to deduce q2v q2v = qsfc*wrk + max(qmin,qair)*fhi else ! for dew formation, use saturated q at tskin - qss = fpvs(tv) + qss = fpvs(tah) qss = ep_2 * qss / (psfc + epsm1 * qss) q2v= qss*wrk + max(qmin,qair)*fhi endif