I've just committed a fix for this that should take care of cases like this.  The problem turned out to be that r_max for the halo was actually smaller than r_min, so the profiling itself failed.  I've cut the minimum radius by a factor of 2 to be a little more generous, and also now require that r_max / r_min be at least 2 as well.  There is also an addition check to catch halos that make it through to the profiling but create no valid data.  I've tested both of these and they seem to work.  Let me know if this comes up again.

Britton

On Mon, Jan 26, 2009 at 10:50 AM, Stephen Skory <stephenskory@yahoo.com> wrote:
All,

+        virial = {}
+        if (len(profile[overdensity_
field]) < 2):
+            for field in fields:
+                virial[field] = 0.0
+            return virial


This isn't going to fix it. 'print profile[overdensity_field]' gives this for the bad halo:


[ NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
 NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
 NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
 NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
 NaN  NaN  NaN  NaN]

Again, sorry for not following up more quickly on this.


 _______________________________________________________
sskory@physics.ucsd.edu           o__  Stephen Skory
http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
________________________________(_)_\(_)_______________
_______________________________________________
Yt-dev mailing list
Yt-dev@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org