Hi, I'm looking at the fixed resolution extraction example, and I'm wondering: what ordering (row or column) will the final will have? Thanks, Rick http://yt.enzotools.org/doc/cookbook/recipes.html#extractfixed resolutiondata
Hi, The end of that question should read: "what ordering (row or column) will the final array have?" My apologies for the extra noise. Rick On Dec 2, 2009, at 3:20 PM, Rick Wagner wrote:
Hi,
I'm looking at the fixed resolution extraction example, and I'm wondering: what ordering (row or column) will the final will have?
Thanks, Rick
http://yt.enzotools.org/doc/cookbook/recipes.html#extractfixed resolutiondata _______________________________________________ ytusers mailing list ytusers@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/ytusersspacepope.org
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? Matt
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 _______________________________________________ ytusers mailing list ytusers@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/ytusersspacepope.org
participants (2)

Matthew Turk

Rick Wagner