Shankar,

r_min is defined in your email:
r_min = 2 * self.pf.h.get_smallest_dx() * self.pf['mpc']

it's twice the smallest dx in units of mpc, dx being the smallest cell size.  it is being rejected because your halo has something like 3 cells (radial) in it, which is almost certainly not resolved.  

sam

On Tue, Dec 8, 2009 at 2:11 PM, Agarwal, Shankar <sagarwal@ku.edu> wrote:
Hi,

While running HaloProfiler, I am seeing this...

P001 yt.lagos   ERROR      2009-12-08 15:53:22,403 Skipping halo with r_max / r_min = 1.885682.



I looked at yt/extensions/HaloProfiler.py...

          r_min = 2 * self.pf.h.get_smallest_dx() * self.pf['mpc']
           if (halo['r_max'] / r_min < PROFILE_RADIUS_THRESHOLD):
               mylog.error("Skipping halo with r_max / r_min = %f." % (halo['r_max']/r_min))


r_max would be the distance to the furthest particle in the halo as found by HaloFinder.
But what is r_min ? And what is the basis for skipping this halo ?


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