Index: /issm/trunk-jpl/src/m/classes/geometry.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/geometry.m	(revision 23827)
+++ /issm/trunk-jpl/src/m/classes/geometry.m	(revision 23828)
@@ -66,4 +66,12 @@
 				if strcmp(solution,'TransientSolution') & md.transient.isgroundingline,
 					md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+					if any(self.bed>self.base),
+						md = checkmessage(md,['base<bed on one or more vertex']);
+					end 
+					pos = find(md.mask.groundedice_levelset>0);
+					if any(abs(self.bed(pos)-self.base(pos))>10^-9),
+						md = checkmessage(md,['equality base=bed on grounded ice violated']);
+					end 
+					md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
 				end
 			end