source: issm/oecreview/Archive/16133-16554/ISSM-16308-16309.diff

../trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -2564 +2564 @@
         surface_input->GetInputAverage(&surface);
         bed_input->GetInputAverage(&bed);
         bathymetry_input->GetInputAverage(&bathymetry);
-
+        
         /*Return: */
-        return base*(surface-bed+rho_water/rho_ice*bathymetry);
+        return base*(surface-bed+min(rho_water/rho_ice*bathymetry,0.));
 }
 /*}}}*/
 /*FUNCTION Tria::MassFlux {{{*/

../trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -3711 +3711 @@
         bathymetry_input->GetInputAverage(&bathymetry);
 
         /*Return: */
-        return base*(surface - bed + rho_water/rho_ice * bathymetry);
+        return base*(surface - bed + min( rho_water/rho_ice * bathymetry, 0.) );
 }
 /*}}}*/
 /*FUNCTION Penta::MinVel{{{*/