I think it is likely that you have set the gravitational softening for the new particle type to 0. This will give infinity for the gravity force - as you get here.
Volker
On May 31, 2011, at 9:55 PM, Luciano del Valle wrote:
> Dear Gadget list,
>
> I'm trying to implement the change of types of particles in the way suggested by Volker Springel. My problem is that although for all the "for (i = 0; i <N_gas; i + +)" i added an "if (P [i]. Type == 0)" Gadget2 answer me with problems due to bad "timsteps" for the particles that change the type.
>
> "Error : A timestep of size zero was assigned on the integer timeline!
> We better stop.
> Task=0 Part-ID=7061 Part_Indice=617 dt=nan dt_courant=0 dt_displacement=0.01 All.Timebase_interval=4.47035e-06 tibase=4.47035e-06 ti_step=-2147483648 ac=nan xyz=(-3.6471|-7.8643|-0.00861445) tree=(nan|nannan)
>
> task 0: endrun called with an error level of 818"
>
>
> In fact if I see the aceleration of these particles is a "nan".
> Someone could help me find where I could be find the problem. .I've searched the places that makes explicit the use of particles of type 0 but can not find anything.
>
>
> Cheers!
> Luciano
>
>
>
>
> -----------------------------------------------------------
>
> If you wish to unsubscribe from this mailing, send mail to
> minimalist_at_MPA-Garching.MPG.de with a subject of: unsubscribe gadget-list
> A web-archive of this mailing list is available here:
> http://www.mpa-garching.mpg.de/gadget/gadget-list
Received on 2011-06-01 08:49:55