Ok, this last commit should fix the problem.
+ virial = {}
+ if (len(profile[overdensity_field]) < 2):
+ for field in fields:
+ virial[field] = 0.0
+ return virial
I don't have data to test it, so if you find that it doesn't work, feel free to respond to this thread.
On Sun, Jan 25, 2009 at 4:47 PM, Stephen Skory
<stephenskory@yahoo.com> wrote:
> Thanks for finding that!
You're welcome. I guess you're just used to running your code on actual high-resolution science-quality datasets, while I'm just playing around with crummy rough test cases.