Thanks Jonothan - this works. On Fri, May 1, 2009 at 4:42 PM, Jonathan Guyer <[email protected]> wrote:
> > > On May 1, 2009, at 9:01 AM, david wende wrote: > > In meantime, how do I select a subset of faces from 3D mesh? >> >> [1] m3D = Grid3D(nx=nx, ny=ny, nz=nz) >> [2] x, y, z = m3D.getCellCenters() >> [3] facesRight = m3D.getFacesRight() >> [4] facesLeft = (m3D.getFacesLeft() & (x < 2) ) >> >> Line 3 is OK. >> Line 4 I get an error message: >> "ValueError: shape mismatch: objects cannot be broadcast to a single >> shape" >> > > Cells get selected with Cell centers. Faces get selected with Face centers. > > Try: > > X, Y, Z = m3D.getFaceCenters() > facesLeft = (m3D.getFacesLeft() & (X < 2) ) > > > x, y, and z (or X, Y, and Z) are not mathematical abstractions of the > abscissae, they are arrays of coordinates representing the discrete points > at the centers of the Cells or the centers of the Faces. > > -- David Wende home +972-8-9353488 work +972-2-5886116 mobile +972-54-234-6479
