[yt-users] 'dx' field from ortho_ray

[Matthew Turk]

Hi Mike,

You're right, this is a bug.  I've committed a fix in r1695.
Essentially, fields that can be generated just from existing ray data
(i.e., things like pressure, where you don't need to touch the actual
grids if you already have the constituent parts in the ray's .data
object) are automatically sorted, by definition; so, they don't need
to be resorted.  However, the logic that was in place, meant to ensure
that only *new* fields read from disk were sorted, was not allowing
for something like dx -- which needs the grid object but is not found
on disk -- to be sorted.  So this would apply to any field using the
validator ValidateGridType or ValidateSpatial -- which dx is!

Let me know if this doesn't work, and thanks very much for the bug report!

-Matt

On Fri, Apr 23, 2010 at 5:04 PM, Michael Kuhlen <mqk at astro.berkeley.edu> wrote:
> 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 at astro.berkeley.edu   UC Berkeley                      *
> *  cell phone: (831) 588-1468      601 Campbell Hall                *
> *  skype username: mikekuhlen      Berkeley, CA 94720               *
> *                                                                   *
> *********************************************************************
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>