[Yt-dev] HaloProfiler Problems

[Stephen Skory]

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 at physics.ucsd.edu           o__  Stephen Skory
http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
________________________________(_)_\(_)_______________