[yt-users] How to find densest point in subvolume of Gadget data

[Joseph Smidt]

Hey everyone,

    I am trying to find the densest point in Gadget data in the
subvolume with left edge [700]*3 and right edge [800]*3 where 700 and
800 are in code units.

     This is what I do [1],  which was inspired by this page [2],
except I try and take the all_data confined to a region.  The output
[3] says however where you will note the center is not inside the
region. It found the densest point in the whole data set, not the
region.

    So how would I do this instead for a subvolume of Gadget data?  Thanks.


[1] http://pastebin.com/kyUgEi3a

[2] http://yt-project.org/doc/cookbook/gadget_notebook.html

[3]  http://pastebin.com/E4aLb5VC



-- 
------------------------------------------------------------------------
Joseph Smidt <josephsmidt at gmail.com>

Theoretical Division
P.O. Box 1663, Mail Stop B283
Los Alamos, NM 87545
Office: 505-665-9752
Fax:    505-667-1931
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org