From 2afb1f4d04e4a34f8023882f6ba0777528059985 Mon Sep 17 00:00:00 2001 From: yuanxue2870 Date: Mon, 11 Dec 2023 19:49:25 +0000 Subject: [PATCH] add comment to smc pert unit & fix bug in stc pert --- lndp_apply_perts.F90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/lndp_apply_perts.F90 b/lndp_apply_perts.F90 index 705ef015..c88cbba5 100644 --- a/lndp_apply_perts.F90 +++ b/lndp_apply_perts.F90 @@ -232,6 +232,8 @@ subroutine lndp_apply_perts(blksz, lsm, lsm_noah, lsm_ruc, lsm_noahmp, iopt_dveg ! perturb total soil moisture ! factor of sldepth*1000 converts from mm to m3/m3 + ! NOTE: smc in the pertlist specified in input.nml + ! is in the unit of mm/hour pert = sfc_wts(nb,i,v)*smc_vertscale(k)*lndp_prt_list(v)/(zslayer(k)*1000.) pert = pert*tfactor_state @@ -248,7 +250,7 @@ subroutine lndp_apply_perts(blksz, lsm, lsm_noah, lsm_ruc, lsm_noahmp, iopt_dveg if (do_pert_state) then do k=1,lsoil pert = sfc_wts(nb,i,v)*stc_vertscale(k)*lndp_prt_list(v) - pert = tfactor_state + pert = pert*tfactor_state call apply_pert('stc',pert,print_flag, stc(nb,i,k),ierr) enddo endif