The galaxy is scaled based on the disk with mass=1 and (exponential) scalelength=1. In addition, a selection of the solar radius (2.43 scale lenghts is the current default) will fix the rotation speed at the IAU value of 220 km/s, though the velocity is also normalized to this value.
This NEMO frontend only controls a number of parameters, to match the functionality of mkkd95(1NEMO) , however there are quite a few more options in the original program, e.g.. adding a gas disk, satellited, individual smoothing lengths. These will be implemented in a future release. It is easy for the user to manually create such models those, by setting cleanup=f and editing the input parameter file magalie.in in the temporary run directory, and running the make-it shell script. Careful not to change the number of particles, or modify the corresponding keyword in this script.
543216789 ! Random generator seed (ran1) n ! Output particles smooth length? 20000 ! Number of disc particles (M_d=1) 3. ! Mass of disk (reset = 1 for default) 1. ! Disc scalelength 2.43 ! Solar radius 1.5 ! Toomre Q parameter @ Sun ( was 3/2 -> close to 1) 2. ! z-max (= 10.*scale height) 10. ! Max disc part radius n ! Gas in disc? 0 ! Number of disc gas particles 0. ! mass of disc gas particles 0. ! temperature of gas in disc 0. ! scale height of gas in disc z0gas 0. ! max z value of gas in disc (= 10 z0gas) 0. ! max cylindrical radius of gas in disc 0. ! min " " " " " " n ! include gas self-gravity (in disc)? [y/n] n ! Add bulge? 0.75 ! Mass of bulge 2. ! Bulge scale length (hernquist ’a’) y ! bulge self-gravity? 5000 ! N part in bulge 10. ! max radius for bulge particles n ! Non-spherical bulge? 1. ! Z-max value for bulge particles 5 ! Number of Simpson integration steps n ! Include rotation ? (flip z-momentum component) n ! Include a halo? y ! halo self-gravity ? 20. ! max radius of halo 40000 ! Halo particle number 0.01 ! halo particle softening length LH ! Halo type (LH or IS) 10. ! Halo mass 2. ! Halo core radius (need two lengths if IS halo) 2. ! Second length : dummy (LH) or truncation radius (IS) n ! Non-spherical halo? n ! Galactic satellite 1. ! Satellite scale length (LH model only) 5.,0.,0. ! xyz position of the satellite in Galaxy frame of reference 5. ! maximum satellite size (compare with sat length scale) y ! Include satellite self-gravity ? (yes = live model only) 0 ! number of satellite particles 0.1 ! smoothing length for satelllite particles b ! data format for output (a)scii or (b)inary nbody ! formatted as LH, nbody roundbulge ! name of dataset (note - truncated at cr character)
Does not work on the linux/intel compiler yet.
mkkd95(1NEMO), tabtos(1NEMO), unfio(1NEMO) https://ui.adsabs.harvard.edu/abs/2001NewA....6...27B/ (Boily et al. 2001; magalie) https://ui.adsabs.harvard.edu/abs/1993ApJS...86..389H/abstract (Hernquist 1993) https://ui.adsabs.harvard.edu/abs/2015MNRAS.452.2367Y/abstract (Yurin & Springel 2015; GalIC) Files NEMO/src/nbody/init/magalie.c NEMO/usr/boily/Magalie/
Christian Boily, Pavel Kroupa, Jorge Penarrubia-Garrido (fortran code) Peter Teuben (NEMO interface) -
dark ages V1.0 See their NewA paper Boily et al. 21-Mar-04 V1.1 Created at the Nbody school PJT 24-mar-04 V1.2 added most primary keywords, at 37,000ft PJT 25-jan-05 V1.2a fixed bulge mass encoding problem 24-mar-06 V1.2b fixed bulge size encoding problem,updated docs PJT 19-jul-06 V1.3 add header= for 64 bit compilers PJT 5-aug-06 merged CVS doc versions PJT feb-2020 added refs PJT
Table of Contents