<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">Hi, Matt--<div><br></div><div>No worries. That's where it looked like it should go, but there wasn't any prior so I figured it best to ask.</div><div><br></div><div>So that's useful, now it tells me the problem is with thermal_energy. Which I find awful strange, since I see (and have printf indicating) that the self.add_field('thermal_energy') bit does get called for my sim. Here's the trace. </div><div><br></div><div>Ahh, I see. That's a super useful feature, thermal_energy is complaining about MagneticEnergy now. Why is MagneticEnergy is not defined anywhere. Why is that? </div><div><br></div><div>Hm. Why is Bx not properly aliased to magnetic_field_x? </div><div><br></div><div><div>Traceback (most recent call last):</div><div> File "<stdin>", line 1, in <module></div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py", line 233, in __getitem__</div><div> f = self._determine_fields([key])[0]</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py", line 489, in _determine_fields</div><div> finfo = self.ds._get_field_info("unknown", fname)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py", line 444, in _get_field_info</div><div> self.index</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py", line 279, in index</div><div> self.create_field_info()</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py", line 334, in create_field_info</div><div> self.field_info.load_all_plugins()</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/field_info_container.py", line 228, in load_all_plugins</div><div> self.find_dependencies(loaded)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/field_info_container.py", line 241, in find_dependencies</div><div> deps, unavailable = self.check_derived_fields(loaded)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/field_info_container.py", line 310, in check_derived_fields</div><div> fd = fi.get_dependencies(ds = self.ds)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/derived_field.py", line 158, in get_dependencies</div><div> e[<a href="http://self.name">self.name</a>]</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/field_detector.py", line 113, in __missing__</div><div> vv = finfo(self)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/derived_field.py", line 184, in __call__</div><div> dd = self._function(self, data)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/fluid_fields.py", line 129, in _pressure</div><div> * (data[ftype, "density"] * data[ftype, "thermal_energy"])</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/fields/field_detector.py", line 97, in __missing__</div><div> finfo = self.ds._get_field_info(*field)</div><div> File "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py", line 482, in _get_field_info</div><div> raise YTFieldNotFound((ftype, fname), self)</div><div>yt.utilities.exceptions.YTFieldNotFound: Could not find field '('gas', 'thermal_energy')' in data0000.</div></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 5, 2015 at 2:03 PM, 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 Dave,<br>
<br>
Sorry for the abruptness! If you're already editing the source code,<br>
just toss it on to the field_info_container.py definition. Stick it<br>
in the __init__ method like so:<br>
<br>
self._show_field_errors = [("gas","pressure")]<br>
<div class="HOEnZb"><div class="h5"><br>
On Tue, May 5, 2015 at 12:58 PM, David Collins <<a href="mailto:dcollins4096@gmail.com">dcollins4096@gmail.com</a>> wrote:<br>
> What instance of FieldInfoContainer should I attach that to?<br>
><br>
> On Tue, May 5, 2015 at 1:54 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:<br>
>><br>
>> You can append ("gas", "pressure") to _show_field_errors on<br>
>> FieldInfoContainer to raise whatever fails.<br>
>><br>
>><br>
>> On Tue, May 5, 2015, 12:49 PM Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
>> wrote:<br>
>>><br>
>>> Hi Dave,<br>
>>><br>
>>> The pressure field is probably coming from the definition in<br>
>>> yt/fields/fluid_fields.py. If it doesn't show up for your dataset, it's<br>
>>> likely failing field detection.<br>
>>><br>
>>> I'm not sure whether we have a test dataset on <a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a> that<br>
>>> uses HydroMethod 4 or 6. Do you happen to have a small one I could play<br>
>>> with?<br>
>>><br>
>>> -Nathan<br>
>>><br>
>>> On Tue, May 5, 2015 at 10:40 AM, David Collins <<a href="mailto:dcollins4096@gmail.com">dcollins4096@gmail.com</a>><br>
>>> wrote:<br>
>>>><br>
>>>> Hi, Matt--<br>
>>>><br>
>>>> Thanks for your input. In frontends/enzo/feilds.py, I see where the<br>
>>>> thermal_energy is computed based on HydroMethod=4,6, but there's no pressure<br>
>>>> field. Which in principle is fine, since the computation in<br>
>>>> yt/fields/fluid_fields.py (in setup_fluid_fields) is fine. Does that<br>
>>>> instance not get called?<br>
>>>><br>
>>>> On Tue, May 5, 2015 at 12:26 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
>>>> wrote:<br>
>>>>><br>
>>>>> Hi Dave,<br>
>>>>><br>
>>>>> You can see all the logic inside yt/frontends/enzo/fields.py --<br>
>>>>> probably not too hard to add a new one, and then PR it!<br>
>>>>><br>
>>>>> -Matt<br>
>>>>><br>
>>>>> On Tue, May 5, 2015 at 11:02 AM, David Collins <<a href="mailto:dcollins4096@gmail.com">dcollins4096@gmail.com</a>><br>
>>>>> wrote:<br>
>>>>> > Hi, Everybody!<br>
>>>>> ><br>
>>>>> > I'm having a little trouble with getting pressure when using<br>
>>>>> > HydroMethod=6.<br>
>>>>> > For instance,<br>
>>>>> >>>> ds0.all_data()['pressure']<br>
>>>>> ><br>
>>>>> > Traceback (most recent call last):<br>
>>>>> > File "<stdin>", line 1, in <module><br>
>>>>> > File<br>
>>>>> ><br>
>>>>> > "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py",<br>
>>>>> > line 233, in __getitem__<br>
>>>>> > f = self._determine_fields([key])[0]<br>
>>>>> > File<br>
>>>>> ><br>
>>>>> > "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py",<br>
>>>>> > line 489, in _determine_fields<br>
>>>>> > finfo = self.ds._get_field_info("unknown", fname)<br>
>>>>> > File<br>
>>>>> ><br>
>>>>> > "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py",<br>
>>>>> > line 482, in _get_field_info<br>
>>>>> > raise YTFieldNotFound((ftype, fname), self)<br>
>>>>> > yt.utilities.exceptions.YTFieldNotFound: Could not find field<br>
>>>>> > '('all',<br>
>>>>> > 'pressure')' in data0000.<br>
>>>>> ><br>
>>>>> ><br>
>>>>> ><br>
>>>>> > while this works for a run with HydroMethod = 0.<br>
>>>>> > Could someone give me the nickel tour of hydro method impacts setting<br>
>>>>> > of<br>
>>>>> > variables, and where I should look to make sure pressure is defined?<br>
>>>>> ><br>
>>>>> > Thanks!<br>
>>>>> > d.<br>
>>>>> ><br>
>>>>> > --<br>
>>>>> > -- Sent from a computer.<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>
>>>>> ><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>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>> --<br>
>>>> -- Sent from a computer.<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>
>>>><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>
>><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>
>><br>
><br>
><br>
><br>
> --<br>
> -- Sent from a computer.<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>
><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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">-- Sent from a computer.<br></div></div>
</div>