[FLASH-USERS] Dark Matter Density

[Anshu Dubey]

On Mon, Jul 7, 2008 at 9:12 AM, M.A. Latife <latife at astro.rug.nl> wrote:
> Hi
> Dear All,
> I am including particles in my simulations to include the effect of dark
> matter. In user guide it is mentioned that pden is mesh mapped particle
> density but it seems me that it is sum of particle density + gas density. Is
> it true?

pden is as stated in the users guide outside of the gravity unit. In the routine
Grid_potentialListOfBlocks, pden briefly becomes a combination of partlcle and
gas density to calculate the potential. Once potential is calculated,
the value of
gas density is deleted from pden before returning from the routine.

> My second question is ,I want to create  density profile  like
> rho=rho0*Radius^(-alpha)  for  dark matter . To create a dark matter density
> profile   should i initialize  more particles where i want to create profile
> (i mean changing Particle  positions file )or FLASH will do it by putting
> profile of dark matter in simulation _init Block.f90  as it works for gas.
> Thanks in advance
> Best regards
> Latif
>cd
You will have to write a customized "Particles_initPositions" (see
PoisParticles setup for an example)
or pt_initPositions (see Orbit setup for example)  in your setup directory.

>