Added roughness length over land to the perturbed variables.

NOAA-PSL · Jan 12, 2021 · ffc88da · ffc88da
1 parent aeb543d
commit ffc88da
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 2 deletions.
diff --git a/compns_stochy.F90 b/compns_stochy.F90
@@ -294,7 +294,7 @@ subroutine compns_stochy (me,sz_nml,input_nml_file,fn_nml,nlunit,deltim,iret)
             'land perturbations will be applied to selected paramaters, using newer scheme designed for DA ens spread'
                do k =1,n_var_lndp
                    select case (lndp_var_list(k))
-                   case('vgf','smc','stc','alb', 'sal','emi') 
+                   case('vgf','smc','stc','alb', 'sal','emi','zol','zoi','zow') 
                        if (me==0) print*, 'land perturbation will be applied to ', lndp_var_list(k)
                    case default
                       print*, 'ERROR: land perturbation requested for new parameter - will need to be coded in lndp_apply_pert', lndp_var_list(k)

diff --git a/lndp_apply_perts.F90 b/lndp_apply_perts.F90
@@ -21,7 +21,7 @@ subroutine lndp_apply_perts(blksz, lsm, lsm_noah, lsm_ruc, lsoil,           &
                 n_var_lndp, lndp_var_list, lndp_prt_list,                       &
                 sfc_wts, xlon, xlat, stype, smcmax, smcmin, param_update_flag,  &
                 smc, slc, stc, vfrac, alvsf, alnsf, alvwf, alnwf, facsf, facwf, &
-                snoalb, semis, ierr)
+                snoalb, semis, zorll, ierr)
 
         implicit none
 
@@ -55,6 +55,7 @@ subroutine lndp_apply_perts(blksz, lsm, lsm_noah, lsm_ruc, lsoil,           &
         real(kind=kind_dbl_prec),     intent(inout) :: facsf(:,:)
         real(kind=kind_dbl_prec),     intent(inout) :: facwf(:,:)
         real(kind=kind_dbl_prec),     intent(inout) :: semis(:,:)
+        real(kind=kind_dbl_prec),     intent(inout) :: zorll(:,:)
 
         ! intent(out)
         integer,                        intent(out) :: ierr
@@ -220,6 +221,16 @@ subroutine lndp_apply_perts(blksz, lsm, lsm_noah, lsm_ruc, lsoil,           &
                          pert = pert*factor
                          call apply_pert ('semis',pert,print_flag, semis(nb,i), ierr,p,min_bound, max_bound)
                      endif
+                case('zol')  ! land roughness length
+                     if (param_update_flag .or. lndp_each_step) then
+                         p =5.
+                         min_bound=0.
+                         max_bound=300.
+
+                         pert = sfc_wts(nb,i,v)*lndp_prt_list(v)
+                         pert = pert*factor
+                         call apply_pert ('zol',pert,print_flag, zorll(nb,i), ierr,p,min_bound, max_bound)
+                     endif
                 case default
                     print*, &
                      'ERROR: unrecognised lndp_prt_list option in lndp_apply_pert, exiting', trim(lndp_var_list(v))