<div dir="ltr"><div>Hey Matt,</div><div><br></div><div>Thanks for looking into it - the field list is printed below </div><div><br></div><div><p style="margin:0px;font-size:11px;font-family:Monaco">(Pdb) pprint.pprint(sorted(pf.field_list))</p>
<p style="margin:0px;font-size:11px;font-family:Monaco">[('PartType0', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'Density'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'ElectronAbundance'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'InternalEnergy'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'Metallicity'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'NeutralHydrogenAbundance'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'SmoothingLength'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'StarFormationRate'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType0', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType1', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType1', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType1', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType1', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType1', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType2', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType2', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType2', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType2', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType2', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType3', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType3', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType3', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType3', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType3', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'Metallicity'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'StellarFormationTime'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType4', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'BH_Mass'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'BH_Mdot'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('PartType5', 'Velocities'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('all', 'Coordinates'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('all', 'Masses'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('all', 'ParticleIDs'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('all', 'Potential'),</p>
<p style="margin:0px;font-size:11px;font-family:Monaco"> ('all', 'Velocities')]</p></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, May 15, 2014 at 6:49 AM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Desika,<br>
<br>
Looks like a bug, yes. I think the problem is that "metal_density" is<br>
not getting aliased correctly for your dataset. What is the output<br>
of:<br>
<br>
import pprint<br>
pprint.pprint(sorted(ds.field_list))<br>
<br>
where ds is your dataset.<br>
<br>
-Matt<br>
<div><div class="h5"><br>
On Wed, May 14, 2014 at 9:05 AM, Desika Narayanan<br>
<<a href="mailto:dnarayan@haverford.edu">dnarayan@haverford.edu</a>> wrote:<br>
> Hi YT,<br>
><br>
> I have a function that loads a pf from a Gadget data set, and then smooths a<br>
> few physical quantities (gas metallicity; mass; density) onto an octree, and<br>
> then returns them:<br>
><br>
> <a href="http://paste.yt-project.org/show/4624/" target="_blank">http://paste.yt-project.org/show/4624/</a><br>
><br>
> in an older version of yt (3.0) that I'm using: 7cf523a36d10 (yt-3.0)<br>
><br>
> this works, no problem. if I run yt-update, to get to: 4d1488b8fa47<br>
><br>
> this gives the following traceback copied below. Has something important<br>
> changed that should affect the way I use smoothing functions?<br>
><br>
> -desika<br>
><br>
><br>
><br>
><br>
> File "/Users/desika/Dropbox/powderday/particle_smooth_yt.py", line 41, in<br>
> yt_smooth<br>
> saved["metallicity"] = ad["gas", "metallicity"]<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 214, in __getitem__<br>
> self.get_data(f)<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 623, in get_data<br>
> self._generate_fields(fields_to_generate)<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 640, in _generate_fields<br>
> fd = self._generate_field(field)<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 251, in _generate_field<br>
> tr = self._generate_fluid_field(field)<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 271, in _generate_fluid_field<br>
> rv = finfo(gen_obj)<br>
> File "/Users/desika/yt-x86_64/src/yt-hg/yt/fields/derived_field.py", line<br>
> 178, in __call__<br>
> dd = self._function(self, data)<br>
> File "/Users/desika/yt-x86_64/src/yt-hg/yt/fields/fluid_fields.py", line<br>
> 153, in _metallicity<br>
> tr = data[ftype, "metal_density"] / data[ftype, "density"]<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 207, in __getitem__<br>
> f = self._determine_fields([key])[0]<br>
> File<br>
> "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line<br>
> 455, in _determine_fields<br>
> finfo = self.pf._get_field_info(ftype, fname)<br>
> File "/Users/desika/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py",<br>
> line 487, in _get_field_info<br>
> raise YTFieldNotFound((ftype, fname), self)<br>
> yt.utilities.exceptions.YTFieldNotFound: Could not find field '('gas',<br>
> 'metal_density')' in snapshot_006.<br>
><br>
><br>
><br>
><br>
</div></div>> _______________________________________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br></div>