Hey yt users,<div><br></div><div><div>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...</div>
<div><br></div><div>Here is what I'm getting:</div><div><br></div><blockquote class="gmail_quote" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<font face="'courier new', monospace">Finding clumps: min: 1.000000e-22, max: 1.000000e-21, step: 10.000000<br></font><font face="'courier new', monospace">Wiping out existing children clumps.<br>
</font><font face="'courier new', monospace">Joining...<br></font><font face="'courier new', monospace">Finished joining in 1.06e+00 seconds<br></font><font face="'courier new', monospace">Investigating 2 children.<br>
</font><font face="'courier new', monospace">Finding clumps: min: 1.000000e-21, max: 1.000000e-21, step: 10.000000<br></font><font face="'courier new', monospace">Traceback (most recent call last):<br>
</font><font face="'courier new', monospace"> File "/scratch/01198/earln/run_files/clump_analysis/find_clumps_plot.py", line 69, in <module><br></font><font face="'courier new', monospace"> amods.level_sets.find_clumps(master_clump, c_min, c_max, step)<br>
</font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py", line 179, in find_clumps<br></font><font face="'courier new', monospace"> find_clumps(clump, min_val*d_clump, max_val, d_clump)<br>
</font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py", line 188, in find_clumps<br></font><font face="'courier new', monospace"> elif (child._isValid()):<br>
</font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/analysis_modules/level_sets/clump_handling.py", line 143, in _isValid<br></font><font face="'courier new', monospace"> self.function_value = eval(self.function)<br>
</font><font face="'courier new', monospace"> File "<string>", line 1, in <module><br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py", line 93, in __call__<br>
</font><font face="'courier new', monospace"> return self._call_func_unlazy(args, kwargs)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py", line 116, in _call_func_unlazy<br>
</font><font face="'courier new', monospace"> retval = self.func(self._data_source, *args, **kwargs)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/derived_quantities.py", line 320, in _IsBound<br>
</font><font face="'courier new', monospace"> thermal = (data["ThermalEnergy"] * data["CellMass"]).sum()<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 286, in __getitem__<br>
</font><font face="'courier new', monospace"> self.get_data(key)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 2277, in get_data<br>
</font><font face="'courier new', monospace"> if self._generate_field(field):<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 2347, in _generate_field<br>
</font><font face="'courier new', monospace"> self[field] = self.pf.field_info[field](self)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/field_info_container.py", line 309, in __call__<br>
</font><font face="'courier new', monospace"> dd = self._function(self, data)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/frontends/enzo/fields.py", line 131, in _ThermalEnergy<br>
</font><font face="'courier new', monospace"> return data["GasEnergy"]<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 286, in __getitem__<br>
</font><font face="'courier new', monospace"> self.get_data(key)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 2277, in get_data<br>
</font><font face="'courier new', monospace"> if self._generate_field(field):<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/data_containers.py", line 2347, in _generate_field<br>
</font><font face="'courier new', monospace"> self[field] = self.pf.field_info[field](self)<br></font><font face="'courier new', monospace"> File "/share/home/00688/tg458377/local/src/yt-hg/yt/data_objects/field_info_container.py", line 310, in __call__<br>
</font><font face="'courier new', monospace"> dd *= self._convert_function(data)<br></font><font face="'courier new', monospace">TypeError: unsupported operand type(s) for *=: 'NoneType' and 'float'<br>
</font><font face="'courier new', monospace">TACC: Cleaning up after job: 2259004<br></font><font face="'courier new', monospace">TACC: Done.</font></blockquote>
<div><br></div><div>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<font class="Apple-style-span" face="'courier new', monospace"> function=True</font> parameter, but I kept keeping jobs that would never end, they would just be outputting <font class="Apple-style-span" face="'courier new', monospace">Getting field Density from 1</font> 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.</div>
<div><br></div><div>Any ideas?</div><div><br></div><div>Thanks!</div><div>Nick</div>
</div>