add extra condition in line 2655 of module_sf_noahmp_glacier.F90 to a…

…void FPE errors in some tests
joeolson42 · Nov 9, 2023 · 793ec64 · 793ec64
1 parent 94bc140
commit 793ec64
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/physics/module_sf_noahmp_glacier.F90 b/physics/module_sf_noahmp_glacier.F90
@@ -2652,7 +2652,7 @@ subroutine snowwater_glacier (nsnow   ,nsoil   ,imelt    ,dt     ,sfctmp , & !in
 
 !to obtain equilibrium state of snow in glacier region
 
-   if(sneqv > mwd) then   ! 100 mm -> maximum water depth
+   if(sneqv > mwd .and. isnow /= 0) then   ! 100 mm -> maximum water depth
       bdsnow      = snice(0) / dzsnso(0)
       snoflow     = (sneqv - mwd)
       snice(0)    = snice(0)  - snoflow