[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