[yt-users] cannot access pressure field

[Britton Smith]

Hi JC,

This was my suspicion.  Matt and others have been working hard to fix up the
field system in YT to deal problems like this.  Since different versions of
enzo use TotalEnergy and Total_Energy, it has been difficult to make a
lasting solution to this.  I think something is probably coming very soon.

In the mean time, I can tell you the brute force way to fix this, although
others should please chime in if there is something simpler that can be
done.
Go into your yt installation, and edit the file: yt/frontends/enzo/fields.py
and change the definition of ThermalEnergy (roughly line 128) to use the
Total_Energy field instead of TotalEnergy.  That will do it.

I am normally not in favor of solutions that require editing the source.
I'll reiterate that people are working to provide a permanent solution to
this problem and that it shouldn't be too far off.  If someone else has a
more elegant solution, then please speak up.

Britton

On Mon, Jun 27, 2011 at 8:14 PM, Jean-Claude Passy <jcpassy at gmail.com>wrote:

> **
> Hi Britton,
>
> That's actually very helpful to know you're using Zeus.  With Zeus, the
> Gas_Energy field will not exist.  There should only be the ThermalEnergy
> field.
>
> Oh that's right, I forgot about that.
>
> An issue that keeps popping up is that some Enzo simulations use
> Total_Energy instead of TotalEnergy, or the other way around.  The
> definition of ThermalEnergy calls for using TotalEnergy when the simulation
> used Zeus.  Can you try accessing the TotalEnergy field, explicitly without
> the underscore in the name?
>
> You are right, I can't access 'TotalEnergy'.
>
> JC
>
>
> On Mon, Jun 27, 2011 at 6:35 PM, Jean-Claude Passy <jcpassy at gmail.com>wrote:
>
>>  Hi Britton,
>>
>> 'Total_Energy' works fine but 'Gas_Energy' and 'ThermalEnergy' do not and
>> lead to the same error.
>> I don't know if it can be of any use but I am running with the Zeus
>> solver.
>>
>> Thanks,
>>
>> JC
>>
>>
>> On 27/06/11 12:23, Britton Smith wrote:
>>
>> Hi JC,
>>
>> It looks like the problem is coming up when the Pressure field asks for
>> the ThermalEnergy field.  If you try region['ThermalEnergy'], does that
>> work?  If not, please also try region['TotalEnergy'] and
>> region['GasEnergy'].
>>
>> Britton
>>
>> On Mon, Jun 27, 2011 at 6:19 PM, Jean-Claude Passy <jcpassy at gmail.com>wrote:
>>
>>>  Hi all,
>>>
>>> I must be doing something silly because I get a weird error when trying
>>> to access the 'Pressure' field of a region:
>>>
>>> --------------------------------- TypeError: unsupported operand type(s)
>>> for *=: 'NoneType' and 'float' ---------------------------------
>>>
>>> It works fine with 'Density', 'x-velocity', 'x', ... I have pasted the
>>> complete output error here
>>> http://paste.enzotools.org/show/7up5l2Q3Rc8nnLZA3VKu
>>> Has anyone ever encountered that ?
>>>
>>> Thanks a lot,
>>>
>>> JC
>>>
>>>
>>> PS: I really like the new home page of the yt website. Congrats !
>>>
>>> _______________________________________________
>>> 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/20110627/5f20b7ed/attachment.html>