Jianling Wang wrote:
> Dear Yves,
>
> Thanks for your reply.
>
> > What it seem you forget here is to change the particle type (this is
> > really what makes the difference between gas and other particles).
> > When a particle moves from gas to stars its P[i].Type become 5
> instead of 1.
> > In fact, this is not explicitly written in my point 1) of
> > http://www.mpa-garching.mpg.de/gadget/gadget-list/0025.html
> > but its what I mean when saying : "the routine that changes particle
> > types..."
> >
>
Well, it's difficult to understand what happends with so few informations.
Could you give the error message ?
> Actually I followed your step 1, and 2. I have changed the type of
> particles,
> the result is the same.
> To test this problem, I came back to the original code of gadget, and
> just changed the number of gas between every_timestep_stuff();,
> before domain_Decomposition();.
>
> N_gas = N_gas - 1
> All.TotN_gas = All.TotN_gas -1
> All.NumForcesSinceLastDomainDecomp = 1 +
> All.TreeDomainUpdateFrequency * All.TotNumPart;
>
> The result is still the same.
>
> Best,
> Jianling
>
>
>
>
>
>
>
--
(o o)
--------------------------------------------oOO--(_)--OOo-------
Dr. 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 2010-01-27 08:57:27