Index: ../trunk-jpl/src/m/mesh/bamg.py =================================================================== --- ../trunk-jpl/src/m/mesh/bamg.py (revision 20906) +++ ../trunk-jpl/src/m/mesh/bamg.py (revision 20907) @@ -248,12 +248,10 @@ 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],requiredvertices[:,1],domainfile,0)[0] + flags=ContourToNodes(requiredvertices[:,0],requiredvertices[:,1],domainfile,0) requiredvertices=requiredvertices[numpy.nonzero(flags)[0],:] - #Add all points to bamg_geometry nods=numpy.size(requiredvertices,axis=0) bamg_geometry.Vertices=numpy.vstack((bamg_geometry.Vertices,requiredvertices))