Daniel,
It worked beautifully, thanks!
Fausto
>________________________________
> From: Daniel Wheeler <[email protected]>
>To: Fausto Arinos de A. Barbuto <[email protected]>; Multiple recipients
>of list <[email protected]>
>Sent: Thursday, January 10, 2013 1:31:42 PM
>Subject: Re: Geometries & FiPy
>
>
>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 ]
