[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