set nsoil to max and remove factorsm

HelinWei-NOAA · Aug 25, 2023 · 1336b80 · 1336b80
1 parent 6bf9907
commit 1336b80
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 9 deletions.
diff --git 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
@@ -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)