[yt-users] add_ray problem

Elizabeth Tasker taskere at mcmaster.ca
Thu Oct 28 11:47:12 PDT 2010


Hi,

I've got some (different) problems with add_ray. I have output from the 
Sedov blast test (so unitless sytems) which is run with a boxsize of 10. 
I can analyse the results using:

ray = pf.h.ray([5., 5., 5.], [10.,10.,10.], "Density")
p = pc.add_ray([5.,5.,5.], [10.,10.,10.], "Density")
oray = pf.h.ortho_ray(0, [5., 5.], "Density")

but only ortho_ray gives the correct result, and it's limited to going 
along the axes.

I've attached the results from the 3 techniques. fig0010.png shows the 
analytical solution (red line), the ortho_ray solution (green x) and the 
pf.h.ray solution (blue circles). DD0010_LineQueryPlot shows the results 
from pc.add_ray.

Does anyone know why this isn't behaving? A slice of the shock confirms 
it looks nice and spherical.

Elizabeth
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fig0010.png
Type: image/png
Size: 27594 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20101028/ad3edca4/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DD0010_LineQueryPlot_0_t_Density.png
Type: image/png
Size: 30967 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20101028/ad3edca4/attachment-0001.png>


More information about the yt-users mailing list