[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