Index: /issm/trunk-jpl/src/m/classes/groundingline.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/groundingline.m	(revision 17036)
+++ /issm/trunk-jpl/src/m/classes/groundingline.m	(revision 17037)
@@ -39,6 +39,5 @@
 					md = checkmessage(md,['bathymetry not equal to bed on grounded ice!']);
 				end
-				pos=find(md.mask.groundedice_levelset<0. & md.mask.ice_levelset>=0);
-				if any(md.geometry.bathymetry(pos)-md.geometry.bed(pos)>10^-9),
+				if any(md.geometry.bathymetry - md.geometry.bed > 10^-9),
 					md = checkmessage(md,['bathymetry superior to bed on floating ice!']);
 				end
Index: /issm/trunk-jpl/src/m/classes/groundingline.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/groundingline.py	(revision 17036)
+++ /issm/trunk-jpl/src/m/classes/groundingline.py	(revision 17037)
@@ -50,6 +50,5 @@
 			if any(numpy.abs(md.geometry.bed[pos]-md.geometry.bathymetry[pos])>10**-10):
 				md.checkmessage("bathymetry not equal to bed on grounded ice!")
-			pos=numpy.nonzero(numpy.logical_and(md.mask.groundedice_levelset<0,md.mask.ice_levelset>=0))[0]
-			if any(md.geometry.bathymetry[pos]-md.geometry.bed[pos]>10**-9):
+			if any(md.geometry.bathymetry - md.geometry.bed > 10**-9):
 				md.checkmessage("bathymetry superior to bed on floating ice!")