[yt-dev] Issue #949: Halo analysis cookbook notebook fails (yt_analysis/yt)

[Nathan Goldbaum]

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.