<div dir="ltr">Hi, <div><br></div><div>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:</div>
<div><br></div><div><div>    sp = pf.h.sphere(center, radius)</div><div>    warmhot = sp.cut_region(['grid["Temperature"] >= 1e5', 'grid["Temperature"] < 1e6'])</div><div>    p0 = np.array([center[0], center[1]-250./pf['kpc'], center[2]])</div>
<div>    p1 = np.array([center[0], center[1]+250./pf['kpc'], center[2]])</div><div>    ray = warmhot.hierarchy.ray(p0, p1)</div><div>    tempt = warmhot['Temperature']</div><div>    ray_tempt = ray['Temperature']</div>
<div><br></div><div> when I examine the temperature range for warmnot, it is </div><div>     (100000.1171875, 999999.5625)</div><div>but for ray_tempt, it is: </div><div>     (6607.2802734375, 86239664.0)</div></div><div>
which is strange since the 'ray' data should be a subset of 'warmhot'... Is there any explanation for it? </div><div><br></div><div>Thanks! </div><div>- Yong</div><div><br></div></div>