Dear All,I am trying to create a HaloCatalog with the following code.import ytfrom 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.indexFile "/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 _facttempo = 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 _magfeedmag = (0.05*4.0*pie*data["density"]*(data["TurbVel"]**2))**0.5
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