hp.add_halo_filter(HP.VirialFilter,must_be_virialized=True,
overdensity_field='ActualOverdensity',
virial_overdensity=200,
virial_filters=[['TotalMassMsun','>=','1e14']],
virial_quantities=['TotalMassMsun','RadiusMpc'])
If so, please post it. The number of halos > 10^14 Msun is fairly close to 500 in a 200 Mpc box. If you want the halo profiler to use all halos, set virial_filters=[['TotalMassMsun','>=','1e14']] accordingly.
Sam
On Sun, Feb 28, 2010 at 11:08 AM, Agarwal, Shankar
<sagarwal@ku.edu> wrote:
Hi Stephen,
I read http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2009-March/000110.html
Do you remember the cosmology parameters for your run: particles, grid, box, refinement, AMR ON or OFF ?
I have run a 200 Mpc/h box, 512^3 particle, 512^3 grid, AMR-OFF simulation and I get about 30,000 halos using HOP. When I try to finding virial halos, I get only 500. I may have to increase the grid to 1024^3.
Let me know how you got 40,000 virial halos.
Thanks
Shankar
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
--
Samuel W. Skillman
DOE Computational Science Graduate Fellow
Center for Astrophysics and Space Astronomy
University of Colorado at Boulder
samuel.skillman[at]
colorado.edu