Entropy IO

From: Yves Revaz <yves.revaz_at_epfl.ch>
Date: Wed, 21 Jan 2009 10:25:33 +0100

Dear Gadget List,

In Gadget-2, when a snapshot is written,
velocities have to be synchronized with the positions.
This is done in io.c by first computing dt_gravkick and
then extrapolating using the accelerations.

        fp[k] = P[pindex].Vel[k] + P[pindex].GravAccel[k] * dt_gravkick;
        if(P[pindex].Type == 0)
          fp[k] += SphP[pindex].HydroAccel[k] * dt_hydrokick;

I expected to find the same synchronization for the entropy, something like
       
        SphP[i].Entropy + SphP[i].DtEntropy * dt_entr;


However this is not present. Maybe I'm missing something.
Thanks for your replies,

yves

-- 
                                                 (o o)
--------------------------------------------oOO--(_)--OOo-------
  Yves Revaz
  Laboratory of Astrophysics EPFL
  Observatoire de Sauverny     Tel : ++ 41 22 379 24 28
  51. Ch. des Maillettes       Fax : ++ 41 22 379 22 05
  1290 Sauverny             e-mail : Yves.Revaz_at_epfl.ch
  SWITZERLAND                  Web : http://www.lunix.ch/revaz/
----------------------------------------------------------------

Received on 2009-01-21 11:25:34