[yt-users] Pressure in Enzo methods

Matthew Turk matthewturk at gmail.com
Tue May 5 10:54:41 PDT 2015


You can append ("gas", "pressure") to _show_field_errors on
FieldInfoContainer to raise whatever fails.

On Tue, May 5, 2015, 12:49 PM Nathan Goldbaum <nathan12343 at gmail.com> wrote:

> Hi Dave,
>
> The pressure field is probably coming from the definition in
> yt/fields/fluid_fields.py.  If it doesn't show up for your dataset, it's
> likely failing field detection.
>
> I'm not sure whether we have a test dataset on yt-project.org/data that
> uses HydroMethod 4 or 6.  Do you happen to have a small one I could play
> with?
>
> -Nathan
>
> On Tue, May 5, 2015 at 10:40 AM, David Collins <dcollins4096 at gmail.com>
> wrote:
>
>> Hi, Matt--
>>
>> Thanks for your input.  In frontends/enzo/feilds.py, I see where the
>> thermal_energy is computed based on HydroMethod=4,6, but there's no
>> pressure field.  Which in principle is fine, since the computation in
>> yt/fields/fluid_fields.py (in setup_fluid_fields) is fine.  Does that
>> instance not get called?
>>
>> On Tue, May 5, 2015 at 12:26 PM, Matthew Turk <matthewturk at gmail.com>
>> wrote:
>>
>>> Hi Dave,
>>>
>>> You can see all the logic inside yt/frontends/enzo/fields.py --
>>> probably not too hard to add a new one, and then PR it!
>>>
>>> -Matt
>>>
>>> On Tue, May 5, 2015 at 11:02 AM, David Collins <dcollins4096 at gmail.com>
>>> wrote:
>>> > Hi, Everybody!
>>> >
>>> > I'm having a little trouble with getting pressure when using
>>> HydroMethod=6.
>>> > For instance,
>>> >>>> ds0.all_data()['pressure']
>>> >
>>> > Traceback (most recent call last):
>>> >   File "<stdin>", line 1, in <module>
>>> >   File
>>> >
>>> "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py",
>>> > line 233, in __getitem__
>>> >     f = self._determine_fields([key])[0]
>>> >   File
>>> >
>>> "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/data_containers.py",
>>> > line 489, in _determine_fields
>>> >     finfo = self.ds._get_field_info("unknown", fname)
>>> >   File
>>> >
>>> "/Users/dcollins/local-2014-09-24-yt3/src/yt-hg/yt/data_objects/static_output.py",
>>> > line 482, in _get_field_info
>>> >     raise YTFieldNotFound((ftype, fname), self)
>>> > yt.utilities.exceptions.YTFieldNotFound: Could not find field '('all',
>>> > 'pressure')' in data0000.
>>> >
>>> >
>>> >
>>> > while this works for a run with HydroMethod = 0.
>>> > Could someone give me the nickel tour of hydro method impacts setting
>>> of
>>> > variables, and where I should look to make sure pressure is defined?
>>> >
>>> > Thanks!
>>> > d.
>>> >
>>> > --
>>> > -- Sent from a computer.
>>> >
>>> > _______________________________________________
>>> > yt-users mailing list
>>> > yt-users at lists.spacepope.org
>>> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>> >
>>> _______________________________________________
>>> yt-users mailing list
>>> yt-users at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>>
>>
>>
>>
>> --
>> -- Sent from a computer.
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150505/91d84fdd/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org


More information about the yt-users mailing list