Merge pull request NCAR#70 from barlage/hr2_accumulate

hr2 accumulate
mkavulich · May 1, 2023 · 4fcd86a · 4fcd86a
2 parents 2ca652b + 2fc95fe
commit 4fcd86a
Show file tree

Hide file tree

Showing 5 changed files with 472 additions and 146 deletions.
diff --git a/physics/module_sf_noahmp_glacier.F90 b/physics/module_sf_noahmp_glacier.F90
@@ -1583,7 +1583,7 @@ subroutine sfcdif1_glacier(iter   ,zlvl ,zpd    ,z0h    ,z0m , & !in
 #endif
 
 ! outputs
-    real (kind=kind_phys),                intent(out) :: fv     !< friction velocity (m/s)
+    real (kind=kind_phys),              intent(inout) :: fv     !< friction velocity (m/s)
     real (kind=kind_phys),                intent(out) :: cm     !< drag coefficient for momentum
     real (kind=kind_phys),                intent(out) :: ch     !< drag coefficient for heat
     real (kind=kind_phys),                intent(out) :: ch2    !< drag coefficient for heat
@@ -2646,9 +2646,9 @@ subroutine snowwater_glacier (nsnow   ,nsoil   ,imelt    ,dt     ,sfctmp , & !in
 
 !to obtain equilibrium state of snow in glacier region
 
-   if(sneqv > 2000.) then   ! 2000 mm -> maximum water depth
+   if(sneqv > 100.) then   ! 100 mm -> maximum water depth
       bdsnow      = snice(0) / dzsnso(0)
-      snoflow     = (sneqv - 2000.)
+      snoflow     = (sneqv - 100.)
       snice(0)    = snice(0)  - snoflow 
       dzsnso(0)   = dzsnso(0) - snoflow/bdsnow
       snoflow     = snoflow / dt