adjust denominator for flux limiting by -puny (fbot+fside<0) (#474)

CICE-Consortium · Oct 27, 2023 · 71093e6 · 71093e6
1 parent 6ad0f44
commit 71093e6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
@@ -638,7 +638,7 @@ subroutine frzmlt_bottom_lateral (dt,       ncat,     &
       ! FYI: fside is not yet correct for fsd, may need to adjust fbot further
       !-----------------------------------------------------------------
 
-         xtmp = frzmlt/(fbot + fside + puny)
+         xtmp = frzmlt/(fbot + fside - puny)
          xtmp = min(xtmp, c1)
          fbot  = fbot  * xtmp
          rside = rside * xtmp