[Yt-dev] HaloProfiler Problems

[Matthew Turk]

OH!  And what's the radius of the halo, and how many cells does it get
if you make a sphere from the HOPGroup?  THAT is the key thing,
actually...

On Sun, Jan 25, 2009 at 3:20 PM, Matthew Turk <matthewturk at gmail.com> wrote:
> 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
>>
>