Changeset 13984 for issm/trunk-jpl/src/m/geometry/FlagElements.py
- Timestamp:
- 11/19/12 14:11:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/m/geometry/FlagElements.py ¶
r13741 r13984 24 24 if isinstance(region,(str,unicode)): 25 25 if not region: 26 flag=numpy.zeros(md.mesh.numberofelements, 'bool')26 flag=numpy.zeros(md.mesh.numberofelements,bool) 27 27 invert=0 28 28 elif strcmpi(region,'all'): 29 flag=numpy.ones(md.mesh.numberofelements, 'bool')29 flag=numpy.ones(md.mesh.numberofelements,bool) 30 30 invert=0 31 31 else: … … 44 44 xlim,ylim=basinzoom('basin',region) 45 45 flag_nodes=numpy.logical_and(numpy.logical_and(md.mesh.x<xlim[1],md.mesh.x>xlim[0]),numpy.logical_and(md.mesh.y<ylim[1],md.mesh.y>ylim[0])).astype(float) 46 flag=numpy.prod(flag_nodes[md.mesh.elements],axis=1) 46 flag=numpy.prod(flag_nodes[md.mesh.elements],axis=1).astype(bool) 47 47 else: 48 48 #ok, flag elements 49 [flag,dum]=ContourToMesh(md.mesh.elements[:,0:3].copy(),md.mesh.x.reshape(-1,1),md.mesh.y.reshape(-1,1),region,'element',1) 49 [flag,dum]=ContourToMesh(md.mesh.elements[:,0:3].astype(float),md.mesh.x.reshape(-1,1),md.mesh.y.reshape(-1,1),region,'element',1) 50 flag=flag.astype(bool) 50 51 51 52 if invert:
Note:
See TracChangeset
for help on using the changeset viewer.