[yt-users] Pressure in Enzo methods

David Collins dcollins4096 at gmail.com
Tue May 5 11:21:22 PDT 2015


Ok.  Here's the set:

http://labs.nationaldataservice.org/upload/5b8d380a

d.

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

> I often debug these issues by inserting a breakpoint at line 315 in
> yt/fields/field_info_container.py:
>
>
> https://bitbucket.org/yt_analysis/yt/src/d683286da1623d41e858dc37409910b704f25589/yt/fields/field_info_container.py?at=yt#cl-315
>
> This is where we handle any exceptions that a derived field raises during
> field detection.  If a field raises an error here, the field fails field
> detection and is removed from consideration.  If other fields depend on
> this one, then they won't show up either.
>
> On Tue, May 5, 2015 at 11:03 AM, Matthew Turk <matthewturk at gmail.com>
> wrote:
>
>> Hi Dave,
>>
>> Sorry for the abruptness!  If you're already editing the source code,
>> just toss it on to the field_info_container.py definition.  Stick it
>> in the __init__ method like so:
>>
>> self._show_field_errors = [("gas","pressure")]
>>
>> On Tue, May 5, 2015 at 12:58 PM, David Collins <dcollins4096 at gmail.com>
>> wrote:
>> > What instance of FieldInfoContainer should I attach that to?
>> >
>> > On Tue, May 5, 2015 at 1:54 PM, Matthew Turk <matthewturk at gmail.com>
>> wrote:
>> >>
>> >> 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
>> >>
>> >>
>> >> _______________________________________________
>> >> 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
>>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>


-- 
-- Sent from a computer.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150505/039c9e2c/attachment.html>
-------------- 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