[yt-dev] Issue #949: Halo analysis cookbook notebook fails (yt_analysis/yt)
Nathan Goldbaum
issues-reply at bitbucket.org
Sat Dec 6 12:50:51 PST 2014
New issue 949: Halo analysis cookbook notebook fails
https://bitbucket.org/yt_analysis/yt/issue/949/halo-analysis-cookbook-notebook-fails
Nathan Goldbaum:
>From the current docs build:
```
#!python
INFO:root:Cell returned
INFO:root:Running cell:
# Define a virial radius for the halo.
hc.add_callback("virial_quantities", ["radius"],
profile_storage = "virial_quantities_profiles")
INFO:root:Cell returned
INFO:root:Running cell:
hc.add_callback('delete_attribute','virial_quantities_profiles')
INFO:root:Cell returned
INFO:root:Running cell:
hc.add_callback('sphere', radius_field='radius_200', factor=5,
field_parameters=dict(virial_radius=('quantity', 'radius_200')))
INFO:root:Cell returned
INFO:root:Running cell:
hc.add_callback('profile', 'virial_radius', [('gas','temperature')],
storage='virial_profiles',
weight_field='cell_mass',
accumulation=False, output_dir='profiles')
INFO:root:Cell returned
INFO:root:Running cell:
# Save the profiles
hc.add_callback("save_profiles", storage="virial_profiles", output_dir="profiles")
INFO:root:Cell returned
INFO:root:Running cell:
hc.create()
/home/xarth/yt-x86_64/lib/python2.7/site-packages/numpy/ma/core.py:790: RuntimeWarning: invalid value encountered in less_equal
return umath.less_equal(x, self.critical_value)
INFO:root:Cell raised uncaught exception:
---------------------------------------------------------------------------
YTFieldNotFound Traceback (most recent call last)
<ipython-input-12-76041198ce66> in <module>()
----> 1 hc.create()
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/analysis_modules/halo_analysis/halo_catalog.pyc in create(self, save_halos, save_catalog, njobs, dynamic)
295
296 """
--> 297 self._run(save_halos, save_catalog, njobs=njobs, dynamic=dynamic)
298
299 def load(self, save_halos=True, save_catalog=False, njobs=-1, dynamic=False):
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/utilities/parallel_tools/parallel_analysis_interface.pyc in barrierize(*args, **kwargs)
273 def barrierize(*args, **kwargs):
274 if not parallel_capable:
--> 275 return func(*args, **kwargs)
276 mylog.debug("Entering barrier before %s", func.func_name)
277 comm = _get_comm(args)
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/analysis_modules/halo_analysis/halo_catalog.pyc in _run(self, save_halos, save_catalog, njobs, dynamic)
387 for action_type, action in self.actions:
388 if action_type == "callback":
--> 389 action(new_halo)
390 elif action_type == "filter":
391 halo_filter = action(new_halo)
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/analysis_modules/halo_analysis/halo_callbacks.pyc in __call__(self, halo)
54
55 def __call__(self, halo):
---> 56 self.function(halo, *self.args, **self.kwargs)
57 return True
58
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/analysis_modules/halo_analysis/halo_callbacks.pyc in profile(halo, bin_fields, profile_fields, n_bins, extrema, logs, weight_field, accumulation, fractional, storage, output_dir)
224 my_profile = create_profile(halo.data_object, bin_fields, profile_fields, n_bins=n_bins,
225 extrema=extrema, logs=logs, weight_field=weight_field,
--> 226 accumulation=accumulation, fractional=fractional)
227
228 prof_store = dict([(field, my_profile[field]) \
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/data_objects/profiles.pyc in create_profile(data_source, bin_fields, fields, n_bins, extrema, logs, units, weight_field, accumulation, fractional)
1307
1308 """
-> 1309 bin_fields = data_source._determine_fields(bin_fields)
1310 fields = ensure_list(fields)
1311 if len(bin_fields) == 1:
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/data_objects/data_containers.pyc in _determine_fields(self, fields)
496 else:
497 fname = field
--> 498 finfo = self.ds._get_field_info("unknown", fname)
499 if finfo.particle_type:
500 ftype = self._current_particle_type
/var/lib/jenkins/jobs/yt-docs-3.0/workspace/yt/data_objects/static_output.pyc in _get_field_info(self, ftype, fname)
485 self._last_finfo = self.field_info[(ftype, fname)]
486 return self._last_finfo
--> 487 raise YTFieldNotFound((ftype, fname), self)
488
489 def _setup_classes(self):
YTFieldNotFound: Could not find field '('all', 'virial_radius')' in RedshiftOutput0006.
```
Ping @brittonsmith @hegan.
More information about the yt-dev
mailing list