[Yt-dev] HaloProfiler Problems

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