Hi Daniel, I am making some small progress, will send update in about a week. 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" I guess it is something to do with mismatched array sizes, but I can't seem to find the correct combination. Thanks again for your help till now. David
