put a upper/lower limit on cwpc

HelinWei-NOAA · Mar 24, 2022 · 81a326a · 81a326a
1 parent c8d6545
commit 81a326a
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 1 deletion.
diff --git a/physics/module_sf_noahmplsm.f90 b/physics/module_sf_noahmplsm.f90
@@ -4688,6 +4688,7 @@ subroutine ragrb(parameters,iter   ,vai    ,rhoair ,hg     ,tah    , & !in
 
        cwpc = (cwp * vai * hcan * fhg)**0.5
 !       cwpc = (cwp*fhg)**0.5
+       cwpc = max(min(cwpc,5.0),1.0)
 
        tmp1 = exp( -cwpc*z0hg/hcan )
        tmp2 = exp( -cwpc*(z0h+zpd)/hcan )

diff --git a/physics/sfcsub.F b/physics/sfcsub.F
@@ -34,7 +34,6 @@ module sfccyc_module
       integer, parameter :: kpdalf(2)=(/214,217/)
 !
       real (kind=kind_io8), parameter :: ten=10.0, one=1.0, zero=0.0
-!     integer, parameter :: xdata=5000, ydata=2500, mdata=xdata*ydata
       integer, parameter :: xdata=7200, ydata=3600, mdata=xdata*ydata
       integer            :: veg_type_landice
       integer            :: soil_type_landice