Re: [yt-users] cannot access pressure field

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 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 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 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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org