enforced hpond >= 0 in flush_ponds

ESCOMP · Jan 14, 2025 · 7ae6a71 · 7ae6a71
1 parent 250b4a0
commit 7ae6a71
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/columnphysics/icepack_therm_mushy.F90 b/columnphysics/icepack_therm_mushy.F90
@@ -3225,6 +3225,7 @@ subroutine flush_pond(w, hpond, apond, dt, flpnd, expnd, &
           else
                hpond = hpond + dhpond
           endif
+          hpond = max(hpond, c0)
 
           !-------------------------------------------------------------          
           ! exponential decay of pond (macro-flaw drainage)
@@ -3252,6 +3253,7 @@ subroutine flush_pond(w, hpond, apond, dt, flpnd, expnd, &
           else
                hpond = hpond + dhpond
           endif
+          hpond = max(hpond, c0)
        endif
     endif