Index: ../trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16308)
+++ ../trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16309)
@@ -2564,9 +2564,9 @@
 	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 {{{*/
Index: ../trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16308)
+++ ../trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16309)
@@ -3711,7 +3711,7 @@
 	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{{{*/