[yt-users] object.hierarchy.ray( )

[Yong Zheng]

Hi,

I am trying to use the .hierarchy.ray( ) to extract data along a particular
light ray, but when I check the temperature field of that light ray data, I
find that the data points it extracted were not really what I want, the
temperature range goes beyond my limit, here is part of my code:

    sp = pf.h.sphere(center, radius)
    warmhot = sp.cut_region(['grid["Temperature"] >= 1e5',
'grid["Temperature"] < 1e6'])
    p0 = np.array([center[0], center[1]-250./pf['kpc'], center[2]])
    p1 = np.array([center[0], center[1]+250./pf['kpc'], center[2]])
    ray = warmhot.hierarchy.ray(p0, p1)
    tempt = warmhot['Temperature']
    ray_tempt = ray['Temperature']

 when I examine the temperature range for warmnot, it is
     (100000.1171875, 999999.5625)
but for ray_tempt, it is:
     (6607.2802734375, 86239664.0)
which is strange since the 'ray' data should be a subset of 'warmhot'... Is
there any explanation for it?

Thanks!
- Yong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131106/36495532/attachment.htm>