Table of Contents

name

mkhalo - initial conditions from quite general spherical model (public version)

synopsis

mkhalo [parameter=value]

requirement

requires NEMO version 3.0.13 or higher, see http://www.astro.umd.edu/nemo.

description

mkhalo creates N-body initial conditions with spherical density 
                            sech (r/r_t)
    rho(r) = C ------------------------------------,
                inner     eta     (outer-inner)/eta
               x      * (x    + 1) 
with

                   2    2
        x  = sqrt(r  + r_c)/r_s
Default parameters will construct a Dehnen McLaughlin (2005) model.
parameters out=file The initial conditions are written in NEMO snapshot
format to file. If file=-, output is made to stdout. Default: required. nbody=num
The initial conditions will contain num body masses, positions, and velocities.
Default: required. inner=val The inner logarithmic slope (parameter inner
above) is set to val. Default: 7/9. outer=val The outer logarithmic slope
(parameter outer above) is set to val. Default: 31/9. eta=val The parameter
eta, controlling the transitions in power-law slopes, is set to val. Default:
4/9. mass=val The total mass of the model is set to val; used to determined
the constant C. Default: 1. r_s=val The scale radius of the Dehnen sphere
is set to val. Default: 1. r_2=val If given, r_s is ignored and the scale
radius set such that the negative logarithmic density slope equals 2 at
radius equal  val The scale radius of the Dehnen sphere is set to val. No
Default. r_c=val The core radius of the model set to val. Default: 0. r_t=val
The truncation radius of the model is set to val. Note if val=0, the truncation
radius is set to inifinity. Default: 0. seed=num We will use num as seed
for the (pseudo-) random number generator. If num=0, a unique new seed is
created from the current value of the time, otherwise, the seed actually
provided is used. Default: 0. q_ran=t|f If true, we use quasi-random rather
than pseudo-random numbers for drawing position and angles in velocity space.
Pseudo-random numbers are still used in the rejection-method to draw the
total velocity from the DF. Default: f. eps=val If given, all individual
softening lengths are set to val. No Default. f_pos=val The fraction of bodies
with positive sense of rotation about the z-axis is set to val. Default:
0.5. giveF=t|f If true, the phase-space densities of the model are written
with the snapshot in the auxiliary data fields. Default: f. accname=name
If given, the potential name is taken as external potential (in addition
to the potential generated by the density model sampled) and taken into
account when computing the equilibrium. No Default. accpars=vals If accname
is given, the parameters vals are provided as parameters for the external
potential. No Default. accfile=file If accname is given, the file name is
provided as data file for the external potential. No Default.  see alsomkdehnen(1falcON),
mkplum(1falcON), mkking(1falcON) 
McMillan & Dehnen, 2007, MNRAS, 378, 541

files
falcON/doc/user_guide.pdf                         User Guide for falcON
falcON/doc/mkgalaxy_user_guide.pdf                User Guide for mkhalo

author


Walter Dehnen                              walter.dehnen@astro.le.ac.uk
update history of this man page
17-Sep-2004    first version of manual  WD

Table of Contents