add comment to smc pert unit & fix bug in stc pert

NOAA-PSL · Dec 11, 2023 · 2afb1f4 · 2afb1f4
1 parent 62b8914
commit 2afb1f4
Showing 1 changed file with 3 additions and 1 deletion.
diff --git 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