[yt-users] TotalEnergy

Matthew Turk matthewturk at gmail.com
Wed Nov 10 10:17:00 PST 2010


Hi Elizabeth,

If yt finds a field in the pf.h.field_list list, it should be reading
it directly off disk.  You can, if you like, examine the data directly
from a grid:

d1 = pf.h.io._read_data_set(some_grid, "TotalEnergy")

and compare it against:

d2 = some_grid["TotalEnergy"]

The former will not apply any unit conversions, so the ratio should be
constant if they are the same:

print (na.unique( d1 / d2 ) )

That should give you a single value back, or maybe a few if there's
some jitter from floating point roundoff.

-Matt

On Wed, Nov 10, 2010 at 10:07 AM, Elizabeth Tasker <taskere at mcmaster.ca> wrote:
> Okay, new question:
>
> When I call:
>
> data.quantities["TotalQuantity"]("TotalEnergy")[0]
>
>
> Is that truly total energy or actually the second baryon field? i.e. when I
> run yt on a zeus-hydro run will it return the internal energy or the total
> energy? (I've just run this for a Sedov blast and the total energy seems to
> be increasing, whereas I was expecting it to drop, since Zeus doesn't
> perform this test well and looses energy in the first few time-steps).
>
> Elizabeth
> _______________________________________________
> 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