add canopy heat storage and gvf impact on thermal conductivity

wzzheng90 · Mar 7, 2022 · b90d4e2 · b90d4e2
1 parent 10fa17e
commit b90d4e2
Show file tree

Hide file tree

Showing 5 changed files with 88 additions and 35 deletions.
diff --git a/physics/module_sf_noahmp_glacier.f90 b/physics/module_sf_noahmp_glacier.f90
@@ -1152,7 +1152,8 @@ subroutine glacier_flux (nsoil     ,nsnow   ,emg     ,isnow   ,df      ,dzsnso
 ! the following only applies to opt_sfc =3, opt_sfc = 1 still done its old way
 
         snwd      = snowh*1000.0
-        zlvli     = zlvl - zpd
+!       zlvli     = zlvl - zpd
+        zlvli     = zlvl 
 
 !       fv     = ustarx                                   ! the input maybe too high for glacial
         fv     = ur*vkc/log(zlvli/z0m)