Merge pull request NCAR#219 from HelinWei-NOAA/land_upgrade_hr4

grantfirl · web-flow · commit 2a50cccd916c · 2024-07-23T15:33:49.000-04:00
Land surface upgrades for HR4
diff --git a/physics/Interstitials/UFS_SCM_NEPTUNE/GFS_phys_time_vary.fv3.F90 b/physics/Interstitials/UFS_SCM_NEPTUNE/GFS_phys_time_vary.fv3.F90
@@ -593,8 +593,10 @@ subroutine GFS_phys_time_vary_init (
 
                  isnow = nint(snowxy(ix))+1 ! snowxy <=0.0, dzsno >= 0.0
 
+! using stc and tgxy to linearly interpolate the snow temp for each layer
+
                  do is = isnow,0
-                   tsnoxy(ix,is)  = tgxy(ix)
+                   tsnoxy(ix,is) =  tgxy(ix) + (( sum(dzsno(isnow:is)) -0.5*dzsno(is) )/snd)*(stc(ix,1)-tgxy(ix))
                    snliqxy(ix,is) = zero
                    snicexy(ix,is) = one * dzsno(is) * weasd(ix)/snd
                  enddo
diff --git a/physics/SFC_Layer/UFS/sfc_diag_post.F90 b/physics/SFC_Layer/UFS/sfc_diag_post.F90
@@ -14,15 +14,14 @@ module sfc_diag_post
 !!
 #endif
       subroutine sfc_diag_post_run (im, lsm, lsm_noahmp, opt_diag, dry, lssav, dtf, con_eps, con_epsm1, pgr,&
-                 vegtype,t2mmp,q2mp, t2m, q2m, u10m, v10m, tmpmin, tmpmax, spfhmin, spfhmax,                  &
+                 t2mmp,q2mp, t2m, q2m, u10m, v10m, tmpmin, tmpmax, spfhmin, spfhmax,                  &
                          wind10mmax, u10mmax, v10mmax, dpt2m, errmsg, errflg)
 
         use machine,               only: kind_phys, kind_dbl_prec
 
         implicit none
 
         integer,                             intent(in) :: im, lsm, lsm_noahmp,opt_diag
-        integer,              dimension(:),  intent(in) :: vegtype    !  vegetation type (integer index)
         logical,                             intent(in) :: lssav
         real(kind=kind_phys),                intent(in) :: dtf, con_eps, con_epsm1
         logical             , dimension(:),  intent(in) :: dry
@@ -42,17 +41,6 @@ subroutine sfc_diag_post_run (im, lsm, lsm_noahmp, opt_diag, dry, lssav, dtf, co
         errflg = 0
 
         if (lsm == lsm_noahmp) then
-!     over shrublands use opt_diag=2
-          do i=1, im
-           if(dry(i)) then
-             if (vegtype(i) == 6 .or. vegtype(i) == 7  & 
-                .or. vegtype(i) == 16) then
-              t2m(i) = t2mmp(i)
-              q2m(i) = q2mp(i)
-             endif
-           endif
-          enddo
-             
          if (opt_diag == 2 .or. opt_diag == 3) then
           do i=1,im
             if(dry(i)) then
diff --git a/physics/SFC_Layer/UFS/sfc_diag_post.meta b/physics/SFC_Layer/UFS/sfc_diag_post.meta
@@ -82,13 +82,6 @@
   type = real
   kind = kind_phys
   intent = in
-[vegtype]
-  standard_name = vegetation_type_classification
-  long_name = vegetation type at each grid cell
-  units = index
-  dimensions = (horizontal_loop_extent)
-  type = integer
-  intent= in
 [t2mmp]
   standard_name = temperature_at_2m_from_noahmp
   long_name = 2 meter temperature from noahmp
diff --git a/physics/SFC_Models/Land/Noahmp/module_sf_noahmplsm.F90 b/physics/SFC_Models/Land/Noahmp/module_sf_noahmplsm.F90
@@ -1989,7 +1989,7 @@ subroutine energy (parameters,ice    ,vegtyp ,ist    ,nsnow  ,nsoil  , & !in
 
   real (kind=kind_phys), parameter                   :: mpe    = 1.e-6
   real (kind=kind_phys), parameter                   :: psiwlt = -150.  !metric potential for wilting point (m)
-  real (kind=kind_phys), parameter                   :: z0     = 0.002  ! bare-soil roughness length (m) (i.e., under the canopy)
+  real (kind=kind_phys), parameter                   :: z0     = 0.015  ! bare-soil roughness length (m) (i.e., under the canopy)
 
 ! ---------------------------------------------------------------------------------------------------
 ! initialize fluxes from veg. fraction
@@ -2629,10 +2629,10 @@ subroutine csnow (parameters,isnow   ,nsnow   ,nsoil   ,snice   ,snliq   ,dzsnso
 ! thermal conductivity of snow
 
   do iz = isnow+1, 0
-!     tksno(iz) = 3.2217e-6*bdsnoi(iz)**2.           ! stieglitz(yen,1965)
+!    tksno(iz) = 3.2217e-6*bdsnoi(iz)**2.           ! stieglitz(yen,1965)
 !    tksno(iz) = 2e-2+2.5e-6*bdsnoi(iz)*bdsnoi(iz)   ! anderson, 1976
-!    tksno(iz) = 0.35                                ! constant
-    tksno(iz) = 2.576e-6*bdsnoi(iz)**2. + 0.074    ! verseghy (1991)
+     tksno(iz) = 0.35                                ! constant
+!   tksno(iz) = 2.576e-6*bdsnoi(iz)**2. + 0.074    ! verseghy (1991)
 !    tksno(iz) = 2.22*(bdsnoi(iz)/1000.)**1.88      ! douvill(yen, 1981)
   enddo
 
@@ -4056,11 +4056,6 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
           
         end if
 
-! prepare for longwave rad.
-
-        air = -emv*(1.+(1.-emv)*(1.-emg))*lwdn - emv*emg*sb*tg**4  
-        cir = (2.-emv*(1.-emg))*emv*sb
-!
        if(opt_sfc == 4) then
 
         gdx  = sqrt(garea1)
@@ -4207,6 +4202,11 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
         end if
      end if
 
+! prepare for longwave rad.
+
+        air = -emv*(1.+(1.-emv)*(1.-emg))*lwdn - emv*emg*sb*tg**4  
+        cir = (2.-emv*(1.-emg))*emv*sb
+
 ! prepare for sensible heat flux above veg.
 
         cah  = 1./rahc
@@ -4269,7 +4269,7 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
 
 ! update vegetation surface temperature
         tv  = tv + dtv
-!        tah = ata + bta*tv               ! canopy air t; update here for consistency
+        tah = ata + bta*tv               ! canopy air t; update here for consistency
 
 ! for computing m-o length in the next iteration
         h  = rhoair*cpair*(tah - sfctmp) /rahc        
@@ -4282,15 +4282,7 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
            qfx = (qsfc-qair)*rhoair*caw
         endif
 
-
-        if (liter == 1) then
-           exit loop1 
-        endif
-        if (iter >= 5 .and. abs(dtv) <= 0.01 .and. liter == 0) then
-           liter = 1
-        endif
-
-     end do loop1 ! end stability iteration
+! after canopy balance, do the under-canopy ground balance
 
 ! under-canopy fluxes and tg
 
@@ -4300,8 +4292,6 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
         cev = rhoair*cpair / (gammag*(rawg+rsurf))  ! barlage: change to ground v3.6
         cgh = 2.*df(isnow+1)/dzsnso(isnow+1)
 
-     loop2: do iter = 1, niterg
-
         t = tdc(tg)
         call esat(t, esatw, esati, dsatw, dsati)
         if (t .gt. 0.) then
@@ -4327,7 +4317,14 @@ subroutine vege_flux(parameters,nsnow   ,nsoil   ,isnow   ,vegtyp  ,veg     , &
         gh  = gh  + cgh*dtg
         tg  = tg  + dtg
 
-     end do loop2
+        if (liter == 1) then
+           exit loop1 
+        endif
+        if (iter >= 5 .and. abs(dtv) <= 0.01  .and. abs(dtg) <= 0.01 .and. liter == 0) then
+           liter = 1   ! if conditions are met, then do one final loop
+        endif
+
+     end do loop1
      
 !     tah = (cah*sfctmp + cvh*tv + cgh*tg)/(cah + cvh + cgh)
 
@@ -5824,7 +5821,8 @@ subroutine thermalz0(parameters,    fveg,          z0m, z0mg,       zlvl,
 
       if (opt_trs == z0heqz0m) then
 
-        z0m_out = exp(fveg * log(z0m)      + (1.0 - fveg) * log(z0mg))
+!       z0m_out = exp(fveg * log(z0m)      + (1.0 - fveg) * log(z0mg))
+        z0m_out = fveg * z0m      + (1.0 - fveg) * z0mg
         z0h_out = z0m_out
 
       elseif (opt_trs == chen09) then
@@ -5841,7 +5839,7 @@ subroutine thermalz0(parameters,    fveg,          z0m, z0mg,       zlvl,
         endif
 
         z0h_out = exp( fveg        * log(z0m * exp(-czil*0.4*258.2*sqrt(ustarx*z0m))) + &
-                      (1.0 - fveg) * log(max(z0m/exp(kb_sigma_f0),1.0e-6)) )
+                      (1.0 - fveg) * log(max(z0mg/exp(kb_sigma_f0),1.0e-6)) )
 
       elseif (opt_trs == tessel) then
 
@@ -5880,15 +5878,15 @@ subroutine thermalz0(parameters,    fveg,          z0m, z0mg,       zlvl,
 
         z0h_out = z0m_out
 
-      elseif (opt_trs == chen09 .or. opt_trs == tessel) then
+      elseif (opt_trs == tessel) then
 
         if (vegtyp <= 5) then
           z0h_out = z0m_out
         else
           z0h_out = z0m_out * 0.01
         endif
 
-      elseif (opt_trs == blumel99) then
+      elseif (opt_trs == chen09 .or. opt_trs == blumel99) then
 
         reyn = ustarx*z0m_out/viscosity                      ! Blumel99 eqn 36c
         if (reyn > 2.0) then
diff --git a/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90 b/physics/SFC_Models/Land/Noahmp/noahmpdrv.F90
@@ -663,6 +663,7 @@ subroutine noahmpdrv_run                                       &
   real (kind=kind_phys) :: precip_freeze_frac_in  ! used for penman calculation
  
   real (kind=kind_phys) :: virtfac1               ! virtual factor
+  real (kind=kind_phys) :: tflux                  ! surface flux temp
   real (kind=kind_phys) :: tvs1                   ! surface virtual temp
   real (kind=kind_phys) :: vptemp                 ! virtual potential temp
 
@@ -944,7 +945,8 @@ subroutine noahmpdrv_run                                       &
         t2mmp(i)               = temperature_bare_2m
         q2mp(i)                = spec_humidity_bare_2m
 
-        tskin(i)               = temperature_ground
+        tskin(i)               = temperature_radiative
+        tflux                  = temperature_ground
         surface_temperature    = temperature_ground
         vegetation_fraction    = vegetation_frac
         ch_vegetated           = 0.0
@@ -1038,7 +1040,8 @@ subroutine noahmpdrv_run                                       &
          q2mp(i)  = spec_humidity_veg_2m * vegetation_fraction + &
                    spec_humidity_bare_2m * (1-vegetation_fraction)
 
-         tskin(i) = surface_temperature
+         tskin(i)               = temperature_radiative
+         tflux                  = surface_temperature
 
       endif          ! glacial split ends
 
@@ -1194,9 +1197,9 @@ subroutine noahmpdrv_run                                       &
        endif
 
        if(thsfc_loc) then ! Use local potential temperature
-              tvs1   = tskin(i) * virtfac1
+              tvs1   = tflux * virtfac1
          else ! Use potential temperature referenced to 1000 hPa
-              tvs1   = tskin(i)/prsik1(i) * virtfac1
+              tvs1   = tflux/prsik1(i) * virtfac1
        endif
 
       z0_total  = max(min(z0_total,forcing_height),1.0e-6)
diff --git a/physics/SFC_Models/Land/Noahmp/noahmptable.tbl b/physics/SFC_Models/Land/Noahmp/noahmptable.tbl
@@ -217,7 +217,7 @@
  !---------------------------------------------------------------------------------------------------------------------------------------------------------------------
  ch2op =   0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,    0.1,
  dleaf =  0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,   0.04,
- z0mvt =  1.09,   1.10,   0.85,   0.80,   0.80,   0.20,   0.06,   0.60,   0.50,   0.12,   0.30,   0.15,   1.00,   0.14,   0.00,   0.00,   0.00,   0.30,   0.20,   0.03,
+ z0mvt =  1.00,   1.50,   0.75,   0.90,   0.85,   0.20,   0.10,   0.90,   0.60,   0.20,   0.30,   0.25,   1.00,   0.25,   0.00,  0.015,   0.00,   0.30,   0.10,   0.05,
  hvt   =  20.0,   20.0,   18.0,   16.0,   16.0,   1.10,   1.10,   13.0,   10.0,   1.00,   5.00,   2.00,   15.0,   1.50,   0.00,   0.00,   0.00,   4.00,   2.00,   0.50,
  hvb   =  8.50,   8.00,   7.00,   11.5,   10.0,   0.10,   0.10,   0.10,   0.10,   0.05,   0.10,   0.10,   1.00,   0.10,   0.00,   0.00,   0.00,   0.30,   0.20,   0.10,
  z0mhvt=  0.0545, 0.055,  0.047,  0.050,  0.050,  0.182,  0.0545, 0.046,  0.050,  0.120,  0.060,  0.075,  0.067,  0.093,  0.000,  0.000,  0.000,  0.075,  0.100,  0.060,
@@ -226,32 +226,34 @@
 !mfsno =  2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,   2.50,
 ! c. he 12/17/2020: optimized mfsno values dependent on land type based on evaluation with snotel swe and modis scf, surface albedo
  mfsno =  1.00,   1.00,   1.00,   1.00,   1.00,   2.00,   2.00,   2.00,   2.00,   2.00,   3.00,   3.00,   4.00,   4.00,   2.50,   3.00,   3.00,   3.50,   3.50,   3.50,
+!mfsno =  1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,   1.97,
 ! c. he 12/17/2020: optimized snow cover factor (m) in scf formulation to replace original constant 2.5*z0,z0=0.002m, based on evaluation with snotel swe and modis scf, surface albedo
 ! scffac = 0.008, 0.008,  0.008,  0.008,  0.008,  0.016,  0.016,  0.020,  0.020,  0.020,  0.020,  0.014,  0.042,  0.026,  0.030,  0.016,  0.030,  0.030,  0.030,  0.030,
- scffac = 0.005, 0.005,  0.005,  0.005,  0.005,  0.008,  0.008,  0.010,  0.010,  0.010,  0.010,  0.007,  0.021,  0.013,  0.015,  0.008,  0.015,  0.015,  0.015,  0.015,
+  scffac = 0.005, 0.005,  0.005,  0.005,  0.005,  0.008,  0.008,  0.010,  0.010,  0.010,  0.010,  0.007,  0.021,  0.013,  0.015,  0.008,  0.015,  0.015,  0.015,  0.015,
+! scffac = 0.059, 0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,  0.059,
  cbiom  =  0.02,  0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,   0.02,
 
  ! row 1:  vis
  ! row 2:  near ir
  rhol_vis=0.07,   0.10,   0.07,   0.10,   0.10,   0.07,   0.07,   0.07,   0.10,   0.11,  0.105,   0.11,   0.00,   0.11,   0.00,   0.00,   0.00,   0.10,   0.10,   0.10,
- rhol_nir=0.35,   0.45,   0.35,   0.45,   0.45,   0.35,   0.35,   0.35,   0.45,   0.58,  0.515,   0.58,   0.00,   0.58,   0.00,   0.00,   0.00,   0.45,   0.45,   0.45,
+ rhol_nir=0.35,   0.45,   0.35,   0.45,   0.45,   0.35,   0.35,   0.35,   0.45,   0.35,  0.515,   0.35,   0.00,   0.35,   0.00,   0.00,   0.00,   0.45,   0.45,   0.45,
 
  ! row 1:  vis
  ! row 2:  near ir
- rhos_vis=0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.36,   0.26,   0.36,   0.00,   0.36,   0.00,   0.00,   0.00,   0.16,   0.16,   0.16,
- rhos_nir=0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.58,  0.485,   0.58,   0.00,   0.58,   0.00,   0.00,   0.00,   0.39,   0.39,   0.39,
+ rhos_vis=0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.16,   0.31,   0.26,   0.31,   0.00,   0.31,   0.00,   0.00,   0.00,   0.16,   0.16,   0.16,
+ rhos_nir=0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.39,   0.53,  0.485,   0.53,   0.00,   0.53,   0.00,   0.00,   0.00,   0.39,   0.39,   0.39,
 
  ! row 1:  vis
  ! row 2:  near ir
- taul_vis=0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.07,   0.06,   0.07,   0.00,   0.07,   0.00,   0.00,   0.00,   0.05,   0.05,   0.05,
- taul_nir=0.10,   0.25,   0.10,   0.25,   0.25,   0.10,   0.10,   0.10,   0.25,   0.25,   0.25,   0.25,   0.00,   0.25,   0.00,   0.00,   0.00,   0.25,   0.25,   0.25,
+ taul_vis=0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.05,   0.06,   0.05,   0.00,   0.05,   0.00,   0.00,   0.00,   0.05,   0.05,   0.05,
+ taul_nir=0.10,   0.25,   0.10,   0.25,   0.25,   0.10,   0.10,   0.10,   0.25,   0.34,   0.25,   0.34,   0.00,   0.34,   0.00,   0.00,   0.00,   0.25,   0.25,   0.25,
 
  ! row 1:  vis
  ! row 2:  near ir
- taus_vis=0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.220, 0.1105,  0.220,  0.000,  0.220,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,
- taus_nir=0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.380, 0.1905,  0.380,  0.000,  0.380,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,
+ taus_vis=0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.120, 0.1105,  0.120,  0.000,  0.120,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,
+ taus_nir=0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.001,  0.250, 0.1905,  0.250,  0.000,  0.250,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,
 
- xl    = 0.010,  0.010,  0.010,  0.250,  0.250,  0.010,  0.010,  0.010,  0.010,  -0.30, -0.025,  -0.30,  0.000,  -0.30,  0.000,  0.000,  0.000,  0.250,  0.250,  0.250,
+ xl    = 0.010,  0.10,  0.010,  0.250,  0.250,  0.010,  0.010,  0.010,  0.010,  -0.30, -0.025,  -0.30,  0.000,  -0.30,  0.000,  0.000,  0.000,  0.250,  0.250,  0.250,
 ! make cwpvt vegetation dependent according to j. goudriaan, crop micrometeorology: a simulation study (simulation monographs), 1977). c. he, 12/17/2020
 ! cwpvt =  0.18,   0.67,   0.18,   0.67,   0.29,    1.0,    2.0,    1.3,    1.0,    5.0,   1.17,   1.67,   1.67,   1.67,   0.18,   0.18,   0.18,   0.67,    1.0,   0.18,
  cwpvt =  0.09,   0.335,  0.09,   0.335,  0.145,   0.5,    1.0,    0.65,   0.5,    2.5,   0.585,  0.835,  0.835,  0.835,  0.09,   0.09,   0.09,   0.335,   0.5,   0.09,
@@ -335,10 +337,10 @@
  !--------------------------------------------------------------------------------------------------------------------------------------------------------------------------
  !                1       2       3       4       5       6       7       8       9      10      11      12      13      14      15      16      17      18      19      20  soil color index for soil albedo
  !--------------------------------------------------------------------------------------------------------------------------------------------------------------------------
- albsat_vis =   0.25,   0.23,   0.21,   0.20,   0.19,   0.18,   0.17,   0.16,   0.15,   0.14,   0.13,   0.12,   0.11,   0.10,   0.09,   0.08,   0.07,   0.06,   0.05,   0.04  ! saturated soil albedos
- albsat_nir =   0.50,   0.46,   0.42,   0.40,   0.38,   0.36,   0.34,   0.32,   0.30,   0.28,   0.26,   0.24,   0.22,   0.20,   0.18,   0.16,   0.14,   0.12,   0.10,   0.08  ! saturated soil albedos
- albdry_vis =   0.36,   0.34,   0.32,   0.31,   0.30,   0.29,   0.28,   0.27,   0.26,   0.25,   0.24,   0.23,   0.22,   0.20,   0.18,   0.16,   0.14,   0.12,   0.10,   0.08  ! dry soil albedos
- albdry_nir =   0.61,   0.57,   0.53,   0.51,   0.49,   0.48,   0.45,   0.43,   0.41,   0.39,   0.37,   0.35,   0.33,   0.31,   0.29,   0.27,   0.25,   0.23,   0.21,   0.16  ! dry soil albedos
+ albsat_vis =   0.21,   0.20,   0.18,   0.17,   0.16,   0.15,   0.14,   0.13,   0.13,   0.12,   0.11,   0.10,   0.10,   0.09,   0.08,   0.08,   0.08,   0.07,   0.07,   0.06  ! saturated soil albedos
+ albsat_nir =   0.42,   0.40,   0.36,   0.34,   0.32,   0.30,   0.28,   0.26,   0.26,   0.24,   0.22,   0.20,   0.20,   0.18,   0.16,   0.16,   0.16,   0.14,   0.14,   0.13  ! saturated soil albedos
+ albdry_vis =   0.31,   0.30,   0.28,   0.27,   0.26,   0.24,   0.23,   0.22,   0.22,   0.22,   0.20,   0.19,   0.20,   0.18,   0.16,   0.16,   0.16,   0.14,   0.14,   0.13  ! dry soil albedos
+ albdry_nir =   0.52,   0.50,   0.46,   0.44,   0.42,   0.40,   0.38,   0.37,   0.36,   0.34,   0.32,   0.30,   0.30,   0.28,   0.27,   0.27,   0.27,   0.26,   0.25,   0.25  ! dry soil albedos
  albice     =   0.80,   0.55                                                   ! albedo land ice: 1=vis, 2=nir
  alblak     =   0.60,   0.40                                                   ! albedo frozen lakes: 1=vis, 2=nir
  omegas     =   0.8 ,   0.4                                                    ! two-stream parameter omega for snow
@@ -397,7 +399,7 @@
   class_sno_age = 3600.0  ! snow aging e-folding time (s) in class albedo scheme
   class_alb_new = 0.84  ! fresh snow albedo in class scheme
   psiwlt = -150.0  !metric potential for wilting point (m)
-  z0soil = 0.002  ! bare-soil roughness length (m) (i.e., under the canopy)
+  z0soil = 0.015  ! bare-soil roughness length (m) (i.e., under the canopy)
   z0lake = 0.01   ! lake surface roughness length (m)
 /
 
