Okay, so this could in fact be a problem with the Profile.py module.
Here's the info I'd need -- what are the radii it's trying to
calculate the profile over, and what's the number of bins it's trying
to do? How many bins are you generating?
-Matt
On Sun, Jan 25, 2009 at 3:15 PM, Stephen Skory
Hi,
I'm having some trouble with the HaloProfiler. I'm running with the trunk r1147 sample_halo_profile.par and runHaloProfiler.py. Here's the error I'm getting:
File "runHaloProfiler.py", line 8, in <module> q.makeProfiles() File "/share/home/00649/tg457850/yt/lib/python2.6/site-packages/yt-1.5dev.dev_r1147-py2.6-linux-x86_64.egg/yt/extensions/HaloProfiler.py", line 107, in makeProfiles virial = self._CalculateVirialQuantities(profile) File "/share/home/00649/tg457850/yt/lib/python2.6/site-packages/yt-1.5dev.dev_r1147-py2.6-linux-x86_64.egg/yt/extensions/HaloProfiler.py", line 276, in _CalculateVirialQuantities if (overDensity[1] <= self.haloProfilerParameters['VirialOverdensity']): IndexError: list index out of range
For some of the haloes, profile['CellVolume'] and profile['TotalMassMsun'] is being set to all zeros. If I put this 'if' below in HaloProfiler.py before the call to CalculateVirialQuantities(), say line 104, I don't get the error above and the script runs to completion, making profiles and images.
if profile['CellVolume'][0] == 0.0 : continue
profile['CellVolume'] being zero gives nans in _AddActualOverdensity() in HaloProfiler.py.
If I run this in a unaltered enzo dataset (with no .yt, HopAnalysis.out etc... files), I get this problem for every dataset I've tried. I've tried two runs of 64^3 root grid with 3 levels cosmology, one is 1 mpc/h and the other 64 mpc/h.
I've done some sleuthing, but I don't undertstand enough of how Profiles.py to make sense of it. 'args' in _lazy_add_fields() is returning empty from _get_bins(), and _get_bins() is returning empty because self._get_field() is returning nothing.
Those of you who understand this better, how can I help you?
Thanks!
_______________________________________________________ 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