On Wed, Jan 9, 2013 at 3:33 PM, Fausto Arinos de A. Barbuto < [email protected]> wrote:
> > > Leftface = ((mesh.facesLeft & (Y < H/2))) > Rightface = mesh.facesRight > Bottomface = mesh.facesBottom > Step1 = ((mesh.facesTop & (X > 0.0)) | (mesh.facesTop & (X < 2.0))) > Step2 = ((mesh.facesTop & (X > 2.0))) > Midface = ((mesh.facesTop & (Y > H/2)) | (mesh.facesTop & (Y < H))) > > Hi Fausto, Try this: Leftface = mesh.facesLeft Rightface = mesh.facesRight BottomFace = mesh.facesBottom Step1 = (X < 2) & (Y == 1) Step2 = mesh.facesTop Midface = (X == 2) & (Y > 1) from fipy import numerix print numerix.sum(Midface) print numerix.sum(Step1) print numerix.sum(Step2) Seems like the correct number of faces. Cheers, -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
