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
All,
This isn't going to fix it. 'print profile[overdensity_field]' gives this for the bad halo:
+ virial = {}
+ if (len(profile[overdensity_
field]) < 2):
+ for field in fields:
+ virial[field] = 0.0
+ return virial
[ 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