Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16308)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16309)
@@ -3712,5 +3712,5 @@
 
 	/*Return: */
-	return base*(surface - bed + rho_water/rho_ice * bathymetry);
+	return base*(surface - bed + min( rho_water/rho_ice * bathymetry, 0.) );
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16308)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16309)
@@ -2565,7 +2565,7 @@
 	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.));
 }
 /*}}}*/