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
