Thankyou Nathan.. thankyou very much...<div><br></div><div><br></div><div><br><br><div class="gmail_quote"><div dir="ltr">On Thu 23 Feb, 2017, 9:13 PM Nathan Goldbaum, <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="gmail_msg"><div class="gmail_extra gmail_msg"><div class="gmail_quote gmail_msg">On Thu, Feb 23, 2017 at 4:02 AM, Reju Sam John <span dir="ltr" class="gmail_msg"><<a href="mailto:rejusamjohn@gmail.com" class="gmail_msg" target="_blank">rejusamjohn@gmail.com</a>></span> wrote:<br class="gmail_msg"><blockquote class="gmail_quote gmail_msg" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="gmail_msg">Dear All,<div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">I am trying to create a HaloCatalog with the following code.</div><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg"><br clear="all" class="gmail_msg"></font><div class="gmail_msg"><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">import yt</font></div><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">from yt.analysis_modules.halo_analysis.api import HaloCatalog</font></div></div><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">data_ds = yt.load('../RD0057/RedshiftOutput0057')<br class="gmail_msg"></font></div><div class="gmail_msg"><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">hc = HaloCatalog(data_ds=data_ds, finder_method='hop', finder_kwargs={"threshold": 500.0, "dm_only": False, "ptype": "all", "padding": 0.02})</font></div><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">hc.add_filter("quantity_value", "particle_mass", ">", 5e12, "Msun")</font></div><div class="gmail_msg"><font face="monospace, monospace" color="#ff00ff" class="gmail_msg">hc.create()</font></div></div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">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. </div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg"><div class="gmail_msg"> <font face="monospace, monospace" class="gmail_msg">File "tra_flux_on_earth.py", line 161, in <module></font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> hc = HaloCatalog(halos_ds=halos)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/analysis_modules/halo_analysis/halo_catalog.py", line 126, in __init__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> halos_ds.index</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 424, in index</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> self.create_field_info()</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 481, in create_field_info</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> self.field_info.load_all_plugins()</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 279, in load_all_plugins</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> self.find_dependencies(loaded)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 292, in find_dependencies</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> deps, unavailable = self.check_derived_fields(loaded)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_info_container.py", line 362, in check_derived_fields</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> fd = fi.get_dependencies(ds = self.ds)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 178, in get_dependencies</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> e[<a href="http://self.name" class="gmail_msg" target="_blank">self.name</a>]</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 99, in __missing__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> vv = finfo(self)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 204, in __call__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> dd = self._function(self, data)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "tra_flux_on_earth.py", line 72, in _fact</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> tempo = Obs_freq/(1.6e6*data["MagField"]/for_mag)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 99, in __missing__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> vv = finfo(self)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/derived_field.py", line 204, in __call__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> dd = self._function(self, data)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "tra_flux_on_earth.py", line 36, in _magfeed</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> mag = (0.05*4.0*pie*data["density"]*(data["TurbVel"]**2))**0.5</font></div></div></div></div></blockquote><div class="gmail_msg"><br class="gmail_msg"></div></div></div></div><div dir="ltr" class="gmail_msg"><div class="gmail_extra gmail_msg"><div class="gmail_quote gmail_msg"><div class="gmail_msg">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.</div><div class="gmail_msg"> </div><blockquote class="gmail_quote gmail_msg" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"></blockquote></div></div></div><div dir="ltr" class="gmail_msg"><div class="gmail_extra gmail_msg"><div class="gmail_quote gmail_msg"><blockquote class="gmail_quote gmail_msg" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="gmail_msg"><div class="gmail_msg"><div class="gmail_msg"><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/fields/field_detector.py", line 89, in __missing__</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> finfo = self.ds._get_field_info(*field)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> File "/home/john/anaconda3/lib/python3.6/site-packages/yt/data_objects/static_output.py", line 666, in _get_field_info</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg"> raise YTFieldNotFound((ftype, fname), self)</font></div><div class="gmail_msg"><font face="monospace, monospace" class="gmail_msg">yt.utilities.exceptions.YTFieldNotFound: Could not find field '('all', 'density')' in catalog.0.h5.</font></div></div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">What could be the problem?</div><span class="m_-765490673758994392HOEnZb gmail_msg"><font color="#888888" class="gmail_msg"><div class="gmail_msg"> <br class="gmail_msg"></div>-- <br class="gmail_msg"><div class="m_-765490673758994392m_7982120193208649507gmail_signature gmail_msg"><div dir="ltr" class="gmail_msg">Reju Sam John</div></div>
</font></span></div></div>
<br class="gmail_msg"></blockquote></div></div></div><div dir="ltr" class="gmail_msg"><div class="gmail_extra gmail_msg"><div class="gmail_quote gmail_msg"><blockquote class="gmail_quote gmail_msg" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">_______________________________________________<br class="gmail_msg">
yt-users mailing list<br class="gmail_msg">
<a href="mailto:yt-users@lists.spacepope.org" class="gmail_msg" target="_blank">yt-users@lists.spacepope.org</a><br class="gmail_msg">
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" class="gmail_msg" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br class="gmail_msg">
<br class="gmail_msg"></blockquote></div><br class="gmail_msg"></div></div>
_______________________________________________<br class="gmail_msg">
yt-users mailing list<br class="gmail_msg">
<a href="mailto:yt-users@lists.spacepope.org" class="gmail_msg" target="_blank">yt-users@lists.spacepope.org</a><br class="gmail_msg">
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" class="gmail_msg" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br class="gmail_msg">
</blockquote></div></div>