[yt-users] TotalEnergy

[Elizabeth Tasker]

Thanks Matt -- it is indeed really total energy.

Odd. I'll try this at a higher resolution with PPM....

Elizabeth

Matthew Turk wrote:
> 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
>>
>>     
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>