Index: /issm/trunk-jpl/src/m/classes/slr.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/slr.py	(revision 23032)
+++ /issm/trunk-jpl/src/m/classes/slr.py	(revision 23033)
@@ -161,4 +161,22 @@
 		if (size(self.love_h,0) != size(self.love_k,0) | size(self.love_h,0) != size(self.love_l,0)):
 			error('slr error message: love numbers should be provided at the same level of accuracy')
+
+		#cross check that whereever we have an ice load, the mask is <0 on each vertex: 
+		pos=np.where(self.deltathickness)
+		maskpos=md.mask.ice_levelset[md.mesh.elements[pos,:]] 
+		els=np.where(maskpos>0)
+		if len(els[1])>0:
+			warnings.warn('slr checkconsistency fail: there are elements with ice loads where some vertices are not on the ice!')
+		
+		#check that  if geodetic is requested, we are a mesh3dsurface model (planet), or if we are not, 
+		#a coupler to a planet model is provided. 
+		if self.geodetic:
+			if md.transient.iscoupler: 
+				#we are good; 
+			else
+				if domaintype(md.mesh)=='mesh3dsurface':
+					#we are good
+				else
+					error('model is requesting geodetic computations without being a mesh3dsurface, or being coupled to one!');
 		return md
 	# }}}