Changeset 16962

Timestamp:

11/27/13 00:31:42 (11 years ago)

Author:

Eric.Larour

Message:

CHG: fixed section dealing with required vertices.

File:

• issm/trunk-jpl/src/m/mesh/bamg.py (modified) (1 diff)

issm/trunk-jpl/src/m/mesh/bamg.py

@@ -245 +245 @@
 
                         #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