From 1336b80140557089269af4afab02a540e3e76c27 Mon Sep 17 00:00:00 2001 From: Helin Wei Date: Fri, 25 Aug 2023 12:51:37 -0400 Subject: [PATCH] set nsoil to max and remove factorsm --- physics/module_sf_noahmplsm.F90 | 3 ++- physics/noahmpdrv.F90 | 16 ++++++++-------- 2 files changed, 10 insertions(+), 9 deletions(-) diff --git a/physics/module_sf_noahmplsm.F90 b/physics/module_sf_noahmplsm.F90 index dc7ac54d1..fe9cfb530 100644 --- a/physics/module_sf_noahmplsm.F90 +++ b/physics/module_sf_noahmplsm.F90 @@ -198,7 +198,7 @@ module module_sf_noahmplsm real (kind=kind_phys), parameter :: denice = 917. !< density of ice (kg/m3) integer, private, parameter :: mband = 2 - integer, private, parameter :: nsoil = 9 + integer, private, parameter :: nsoil = 20 integer, private, parameter :: nstage = 8 type noahmp_parameters ! define a noahmp parameters type @@ -747,6 +747,7 @@ subroutine noahmp_sflx (parameters, & canhs = 0. ! -------------------------------------------------------------------------------------------------- + ! re-process atmospheric forcing call atm (parameters,ep_2, epsm1, sfcprs ,sfctmp ,q2 , & diff --git a/physics/noahmpdrv.F90 b/physics/noahmpdrv.F90 index 4bf1eb5fa..b545fb29f 100644 --- a/physics/noahmpdrv.F90 +++ b/physics/noahmpdrv.F90 @@ -2015,14 +2015,14 @@ subroutine noahmpsoilinit (lsm_cold_start, im, lsoil_lsm, lsoil, & ! in soilm(i,k,1) = max(0.02,soilm(i,k,1)*smtotn(i,1)/(0.9*smtotr(i,1))) enddo - if( soilm(i,2,1) > soilm(i,1,1) .and. soilm(i,3,1) > soilm(i,2,1)) then - ! typical for daytime, no recent precip - factorsm(1) = 0.75 - factorsm(2) = 0.8 - factorsm(3) = 0.85 - factorsm(4) = 0.9 - factorsm(5) = 0.95 - endif +! if( soilm(i,2,1) > soilm(i,1,1) .and. soilm(i,3,1) > soilm(i,2,1)) then +! ! typical for daytime, no recent precip +! factorsm(1) = 0.75 +! factorsm(2) = 0.8 +! factorsm(3) = 0.85 +! factorsm(4) = 0.9 +! factorsm(5) = 0.95 +! endif do k=1,lsoil_lsm soilm(i,k,1) = factorsm(k) * soilm(i,k,1)