[yt-users] Halo catalog error.
Reju Sam John
rejusamjohn at gmail.com
Thu Feb 23 02:02:05 PST 2017
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
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170223/5f3a0b63/attachment.htm>
More information about the yt-users
mailing list