[yt-users] Covering Grid Array Order
Rick Wagner
rwagner at physics.ucsd.edu
Wed Dec 2 16:12:57 PST 2009
Hi Matt,
On Dec 2, 2009, at 3:27 PM, Matthew Turk wrote:
> Hi Rick,
>
>> I'm looking at the fixed resolution extraction example, and I'm
>> wondering:
>> what ordering (row or column) will the final will have?
>
> The "cube" object in that example can also return the x, y, z fields.
> If you add this code in to the example:
>
> for xi in [0, -1]:
> for yi in [0, -1]:
> for zi in [0, -1]:
> print xi, yi, zi,
> print cube['x'][xi,yi,zi],
> print cube['y'][xi,yi,zi],
> print cube['z'][xi,yi,zi]
>
> the results give back:
>
> 0 0 0 0.00390625 0.00390625 0.00390625
> 0 0 -1 0.00390625 0.00390625 0.99609375
> 0 -1 0 0.00390625 0.99609375 0.00390625
> 0 -1 -1 0.00390625 0.99609375 0.99609375
> -1 0 0 0.99609375 0.00390625 0.00390625
> -1 0 -1 0.99609375 0.00390625 0.99609375
> -1 -1 0 0.99609375 0.99609375 0.00390625
> -1 -1 -1 0.99609375 0.99609375 0.99609375
>
> So the set of array indices [0,0,0] is the lowest value in all three
> coordinates and the set of array indices [-1, -1, -1] is the maximum
> in all three coordinates. These correspond to the coordinate system
> Enzo uses. Does that answer your question?
I think it does, since the field referenced by cube['x'] is the same
type of object (i.e., Numpy ndarray) as cube["Density"]. After
writing out the x, y and z fields to a test file everything seems to
be in order (some pun intended).
Thanks,
Rick
>
> -Matt
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
More information about the yt-users
mailing list