Changeset 27238

Timestamp:

08/25/22 19:33:28 (3 years ago)

Author:

schlegel

Message:

CHG: tighten GEMB constraints on ice depth

File:

• issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp

@@ -1853 +1853 @@
                         }
                         // if reaches impermeable ice layer all liquid water runs off (R)
-                        else if (d[i] >= dIce-Dtol || (d[i] >= dPHC-Dtol && depthice>0.1)){  // dPHC = pore hole close off [kg m-3]
+                        else if (d[i] >= dIce-Dtol || (d[i] >= dPHC-Dtol && depthice>0.1+Dtol)){  // dPHC = pore hole close off [kg m-3]
                                 // _printf_("ICE LAYER");
                                 // no water freezes in this cell