Yes, we know. Britton's just committed a fix... It's been tested on
"RD0005-mine."
On Mon, Jan 26, 2009 at 9:50 AM, 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@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