From 81a326afa210c402144e9dddc56a45f85c745a70 Mon Sep 17 00:00:00 2001 From: helin wei Date: Thu, 24 Mar 2022 20:18:12 +0000 Subject: [PATCH] put a upper/lower limit on cwpc --- physics/module_sf_noahmplsm.f90 | 1 + physics/sfcsub.F | 1 - 2 files changed, 1 insertion(+), 1 deletion(-) diff --git a/physics/module_sf_noahmplsm.f90 b/physics/module_sf_noahmplsm.f90 index cdc43635b..98364b19c 100644 --- 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 index cdc91cca9..78e5201be 100644 --- 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