[Yt-dev] HaloProfiler Problems

[Matthew Turk]

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 <stephenskory at yahoo.com> wrote:
> 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 at physics.ucsd.edu           o__  Stephen Skory
> http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
> ________________________________(_)_\(_)_______________
> _______________________________________________
> Yt-dev mailing list
> Yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>