On Thu, Feb 23, 2017 at 4:02 AM, Reju Sam John <rejusamjohn@gmail.com> wrote:
Dear All,

I am trying to create a HaloCatalog with the following code.

import yt
from yt.analysis_modules.halo_analysis.api import HaloCatalog
data_ds = yt.load('../RD0057/RedshiftOutput0057')
hc = HaloCatalog(data_ds=data_ds, finder_method='hop', finder_kwargs={"threshold": 500.0, "dm_only": False, "ptype": "all", "padding": 0.02})
hc.add_filter("quantity_value", "particle_mass", ">", 5e12, "Msun")
hc.create()

With this, I am able to create the catalog. But when we are trying to load it, I am unable to get almost all fields. Here i am copying the error message. 

 File "tra_flux_on_earth.py", line 161, in <module>
    hc = HaloCatalog(halos_ds=halos)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/analysis_modules/halo_analysis/halo_catalog.py", line 126, in __init__
    halos_ds.index
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 424, in index
    self.create_field_info()
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 481, in create_field_info
    self.field_info.load_all_plugins()
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 279, in load_all_plugins
    self.find_dependencies(loaded)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 292, in find_dependencies
    deps, unavailable = self.check_derived_fields(loaded)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 362, in check_derived_fields
    fd = fi.get_dependencies(ds = self.ds)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 178, in get_dependencies
    e[self.name]
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 99, in __missing__
    vv = finfo(self)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 204, in __call__
    dd = self._function(self, data)
  File "tra_flux_on_earth.py", line 72, in _fact
    tempo = Obs_freq/(1.6e6*data["MagField"]/for_mag)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 99, in __missing__
    vv = finfo(self)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 204, in __call__
    dd = self._function(self, data)
  File "tra_flux_on_earth.py", line 36, in _magfeed
    mag = (0.05*4.0*pie*data["density"]*(data["TurbVel"]**2))**0.5

The problematic line is here, in your script. You've defined a derived field but that derived field cannot be calculated for the halo catalog dataset. If you created this field using yt.add_field, then I think the solution is to create it using ds.add_field, but only for the dataset instance where this field makes sense to use.
 
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 89, in __missing__
    finfo = self.ds._get_field_info(*field)
  File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 666, in _get_field_info
    raise YTFieldNotFound((ftype, fname), self)
yt.utilities.exceptions.YTFieldNotFound: Could not find field '('all', 'density')' in catalog.0.h5.


What could be the problem?
   
--
Reju Sam John

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org