On Mon, Dec 12, 2011 at 2:00 PM, Nicholas Earl <
earlnich@msu.edu> wrote:
> Hey yt users,
>
> I am getting some weird exiting happening with the clump finder and Ranger,
> and I was wondering if anyone else has experienced the issue, or know what I
> might be doing wrong? I see that it gets to a point where the min and max
> densities are the same, but I didn't think it would cause the script to
> crash...
>
> Here is what I'm getting:
>
>> Finding clumps: min: 1.000000e-22, max: 1.000000e-21, step: 10.000000
>> Wiping out existing children clumps.
>> Joining...
>> Finished joining in 1.06e+00 seconds
>> Investigating 2 children.
>> Finding clumps: min: 1.000000e-21, max: 1.000000e-21, step: 10.000000
>> Traceback (most recent call last):
>> File
>> "/scratch/01198/earln/run_files/clump_analysis/find_clumps_plot.py", line
>> 69, in <module>
>> amods.level_sets.find_clumps(master_clump, c_min, c_max, step)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py",
>> line 179, in find_clumps
>> find_clumps(clump, min_val*d_clump, max_val, d_clump)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py",
>> line 188, in find_clumps
>> elif (child._isValid()):
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py",
>> line 143, in _isValid
>> self.function_value = eval(self.function)
>> File "<string>", line 1, in <module>
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py",
>> line 93, in __call__
>> return self._call_func_unlazy(args, kwargs)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py",
>> line 116, in _call_func_unlazy
>> retval = self.func(self._data_source, *args, **kwargs)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py",
>> line 320, in _IsBound
>> thermal = (data["ThermalEnergy"] * data["CellMass"]).sum()
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 286, in __getitem__
>> self.get_data(key)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 2277, in get_data
>> if self._generate_field(field):
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 2347, in _generate_field
>> self[field] = self.pf.field_info[field](self)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/field_info_container.py",
>> line 309, in __call__
>> dd = self._function(self, data)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/frontends/enzo/fields.py",
>> line 131, in _ThermalEnergy
>> return data["GasEnergy"]
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 286, in __getitem__
>> self.get_data(key)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 2277, in get_data
>> if self._generate_field(field):
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py",
>> line 2347, in _generate_field
>> self[field] = self.pf.field_info[field](self)
>> File
>> "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/field_info_container.py",
>> line 310, in __call__
>> dd *= self._convert_function(data)
>> TypeError: unsupported operand type(s) for *=: 'NoneType' and 'float'
>> TACC: Cleaning up after job: 2259004
>> TACC: Done.
>
>
> I've also attempted lowering the step size to 5.0, but with the same result.
> It also seems to be the case that I do not get this behavior when I use
> function=True parameter, but I kept keeping jobs that would never end, they
> would just be outputting Getting field Density from 1 the entire time. So
> I've had to rerun the clump finder a few times, putting in some minimum
> conditions, and still it'd happen.
>
> Any ideas?
>
> Thanks!
> Nick
>