Changeset 13467

Timestamp:

09/27/12 11:00:28 (12 years ago)

Author:

jschierm

Message:

FIX: Fixed call to ContourToMesh.

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py (modified) (1 diff)

issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py ¶

@@ -25 +25 @@
                 if not os.path.exists(icefrontfile):
                         raise IOError("SetIceShelfBC error message: ice front file '%s' not found." % icefrontfile)
-                #nodeonicefront= numpy.logical_and( md.mesh.vertexonboundary, nodeinsideicefront.astype(float)) )
-                nodeinsideicefront=ContourToMesh(md.mesh.elements,md.mesh.x.reshape(-1,1),md.mesh.y.reshape(-1,1),icefrontfile,'node',2)
-                nodeonicefront= numpy.bitwise_and( map(int,md.mesh.vertexonboundary), map(int,nodeinsideicefront[0].ravel()) )
+                [nodeinsideicefront,dum]=ContourToMesh(md.mesh.elements,md.mesh.x.reshape(-1,1),md.mesh.y.reshape(-1,1),icefrontfile,'node',2)
+                nodeonicefront=numpy.logical_and(md.mesh.vertexonboundary,nodeinsideicefront).astype(float)
         else:
                 nodeonicefront=numpy.zeros((md.mesh.numberofvertices,1))