Hi Florian,
One possibility could be simply to extend the number of type to a value
higher than 6.
Unfortunately in the code the number of type is not coded by a constant
but by
the number 6, so you will have to change this value in a lot of places
in the code :-( .
Then of course, you will have to add new parameters in the gadget param
file,
corresponding to the new types (and read them correctly in begrun.c and
set_softenings in gravtree.c).
SofteningGas
SofteningHalo
SofteningDisk
SofteningBulge
SofteningStars
SofteningBndry
SofteningXxxx
etc...
SofteningGasMaxPhys
SofteningHaloMaxPhys
SofteningDiskMaxPhys
SofteningBulgeMaxPhys
SofteningStarsMaxPhys
SofteningBndryMaxPhys
SofteningXxxxMaxPhys
etc...
Be also carefull to the way the snapshot header will be written, as
npart (and npartTotal)
will have a larger size.
I hope it will help,
yves
Florian Beutler wrote:
> hello
> I am working with a DM halo containing particles with different mass (shell model). I would like to have different softening values for particles with different mass. The softening value is connected to the particle type therefore I am using particle types 1-5 for halo particles (I think these particle types are all collisionless and I can use them all as halo type particles... or am I wrong?). Is there any possibility to extend the particle types to more then 5?
> thanks and regards
> Florian Beutler
>
--
(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-05-26 12:11:59