Maybe you could try increasing the value of PROFILE_RADIUS_THRESHOLD in HaloProfiler.py. Right now, it's set to 2, meaning that it will allow profiles for halos whose outer radii are only 2 times the inner radii. That's probably too generous. Even halos of that size should probably be ignored.
On Tue, Mar 24, 2009 at 11:02 AM, Matthew Turk firstname.lastname@example.org:
So, how do I make it skip that halo? I see there is a len(sphere._grids)
== 0 test, but this case there are two grids so it passes that.
I think the problem might be with the bounds. If no points within the sphere satisfy the radius bounds, it will fail. I will commit the changes in that patch I sent you and then we can insert a try/except block for the exception I added; that is the pragmatic solution, I think. But I really, *really* think we need to figure out why the radius bounds are not giving us any points.
-Matt _______________________________________________ Yt-dev mailing list Ytemail@example.com http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org