Merge pull request NOAA-GSL#39 from barlage/NoahMP_p8c_mynn

add snow density limits for each layer
zhanglikate · May 23, 2022 · d33124c · d33124c
2 parents 0a96734 + 3407bf1
commit d33124c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/physics/module_sf_noahmp_glacier.f90 b/physics/module_sf_noahmp_glacier.f90
@@ -2823,6 +2823,7 @@ subroutine compact_glacier (nsnow  ,nsoil  ,dt     ,stc    ,snice , & !in
            ! the change in dz due to compaction
 
            dzsnso(j) = dzsnso(j)*(1.+pdzdtc)
+           dzsnso(j) = min(max(dzsnso(j),(snliq(j)+snice(j))/500.0),(snliq(j)+snice(j))/50.0)  ! limit adjustment to a reasonable density
         end if
 
         ! pressure of overlying snow

diff --git a/physics/module_sf_noahmplsm.f90 b/physics/module_sf_noahmplsm.f90
@@ -7619,7 +7619,7 @@ subroutine compact (parameters,nsnow  ,nsoil  ,dt     ,stc    ,snice  , & !in
            ! the change in dz due to compaction
 
            dzsnso(j) = dzsnso(j)*(1.+pdzdtc)
-           dzsnso(j) = max(dzsnso(j),snice(j)/denice + snliq(j)/denh2o)
+           dzsnso(j) = min(max(dzsnso(j),(snliq(j)+snice(j))/500.0),(snliq(j)+snice(j))/50.0)  ! limit adjustment to a reasonable density
         end if
 
         ! pressure of overlying snow
@@ -7773,7 +7773,7 @@ subroutine snowh2o (parameters,nsnow  ,nsoil  ,dt     ,qsnfro ,qsnsub , & !in
    end do
 
    do j = isnow+1, 0
-     dzsnso(j) = max(dzsnso(j),snliq(j)/denh2o + snice(j)/denice)
+     dzsnso(j) = min(max(dzsnso(j),(snliq(j)+snice(j))/500.0),(snliq(j)+snice(j))/50.0)  ! limit adjustment to a reasonable density
    end do
 
 ! liquid water from snow bottom to soil