It seems that ortho_ray produces a 'dx' field (also 'CellLength') that is not the same order as the other data fields. This is again for a 1D simulation, this time with a few levels of AMR on. ray = pf.h.ortho_ray(0, [0.5, 0.5]) print ray['x'] yt INFO 2010-04-23 16:54:42,181 Getting field x from 6 [ 0.005 0.015 0.025 0.035 0.045 0.055 0.065 0.075 0.085 0.095 0.105 0.115 0.125 0.135 0.145 0.155 0.165 0.175 0.185 0.195 0.205 0.215 0.225 0.235 0.245 0.255 0.265 0.275 0.285 0.295 0.305 0.315 0.325 0.335 0.345 0.355 0.365 0.375 0.385 0.395 0.405 0.415 0.425 0.435 0.445 0.455 0.465 0.475 0.485 0.495 0.505 0.515 0.525 0.535 0.545 0.555 0.5625 0.5675 0.5725 0.5775 0.5825 0.5875 0.595 0.605 0.615 0.625 0.635 0.645 0.655 0.665 0.675 0.685 0.695 0.705 0.7125 0.7175 0.72125 0.72375 0.725625 0.726875 0.728125 0.7290625 0.7296875 0.7303125 0.7309375 0.7315625 0.7321875 0.7328125 0.7334375 0.7340625 0.7346875 0.7353125 0.7359375 0.736875 0.738125 0.739375 0.74125 0.74375 0.7475 0.7525 0.7575 0.765 0.775 0.785 0.795 0.805 0.815 0.825 0.835 0.845 0.855 0.865 0.875 0.885 0.895 0.905 0.915 0.925 0.935 0.945 0.955 0.965 0.975 0.985 0.995 ] print ray['dx'] yt INFO 2010-04-23 16:54:42,232 Getting field dx from 6 [ 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.0025 0.0025 0.0025 0.0025 0.00125 0.00125 0.00125 0.00125 0.00125 0.00125 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625] Any ideas how to fix this bug? Or some other way to access the ray's 'dx' values in the correct order? I suppose I could simply calculate the dx's myself from the 'x' field... Thanks, Mike -- ********************************************************************* * * * Dr. Michael Kuhlen Theoretical Astrophysics Center * * email: mqk@astro.berkeley.edu UC Berkeley * * cell phone: (831) 588-1468 601 Campbell Hall * * skype username: mikekuhlen Berkeley, CA 94720 * * * *********************************************************************