Index: /issm/trunk-jpl/src/m/mesh/bamg.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/bamg.py	(revision 16961)
+++ /issm/trunk-jpl/src/m/mesh/bamg.py	(revision 16962)
@@ -245,11 +245,12 @@
 
 			#recover RequiredVertices
-			requiredvertices=float(options.getfieldvalue('RequiredVertices'))    #for some reason, it is of class "single"
+			requiredvertices=options.getfieldvalue('RequiredVertices')    #for some reason, it is of class "single"
 			if numpy.size(requiredvertices,axis=1)==2:
 				requiredvertices=numpy.hstack((requiredvertices,4.*numpy.ones((numpy.size(requiredvertices,axis=0),1))))
-	
+			
+
 			#only keep those inside
-			flags=ContourToNodes(requiredvertices[:,0],domain[0],0)
-			requiredvertices=requiredvertices[numpy.nonzero(flags),:]
+			flags=ContourToNodes(requiredvertices[:,0],requiredvertices[:,1],domainfile,0)[0]
+			requiredvertices=requiredvertices[numpy.nonzero(flags)[0],:]
 
 			#Add all points to bamg_geometry