Index: /issm/trunk-jpl/src/m/mesh/bamg.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/bamg.py	(revision 13909)
+++ /issm/trunk-jpl/src/m/mesh/bamg.py	(revision 13910)
@@ -10,4 +10,5 @@
 from MatlabFuncs import *
 from BamgMesher import *
+from ContourToNodes import *
 
 def bamg(md,*args):
@@ -91,5 +92,5 @@
 			#Checks that all holes are INSIDE the principle domain outline
 			if i:
-				flags=ContourToNodes(domaini['x'],domaini['y'],domain[0],0)
+				flags=ContourToNodes(domaini['x'].reshape(-1,1),domaini['y'].reshape(-1,1),domainfile,0)
 				if numpy.any(numpy.logical_not(flags)):
 					raise RuntimeError("bamg error message: All holes should be strictly inside the principal domain")
@@ -117,5 +118,5 @@
 
 				#detect whether all points of the rift are inside the domain
-				flags=ContourToNodes(rifti['x'],rifti['y'],domain[0],0)
+				flags=ContourToNodes(rifti['x'].reshape(-1,1),rifti['y'].reshape(-1,1),domain[0],0)
 				if numpy.all(numpy.logical_not(flags)):
 					raise RuntimeError("one rift has all its points outside of the domain outline")