This code allows you to model a galactic nucleus or a globular cluster with a scheme that approaches the system as a self-gravitating, conducting gas sphere.
For more information, please visit
http://astro-gr.org/modelling-galactic-nuclei-self-gravitating-conducting-gas-spheres/
Pau Amaro Seoane, Berlin 14/03/2016
- Copyright (c) 2016, Pau Amaro Seoane ([email protected])
-
and Rainer Spurzem ([email protected]) - Permission to use, copy, modify, and distribute this software for any
- purpose with or without fee is hereby granted, provided that the above
- copyright notice and this permission notice appear in all copies.
- THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
- WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
- MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
- ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
- WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
- ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
- OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
@run1 is a tcsh-script required to start the run
@run1r is the same but it can be used to restart ("r") a simulation without deleting previous data. In practise I was using this script all the time
kug1.dat contains what you want to be the initial system
To start a simulation: nohup @run1 spedi &
where "spedi" is the executable binary you get out when you compile the main code with "make"
(0) General information: Refer to star.f and params.f
(1) see Run_SglMass_Expl for a directory ready to run single-mass model (use /@run1rstar_sglmass) see Run_MultiMass_Expl for a directory ready to run multi-component model (use star_Ncomp)
To run a multicomponent model, one has to compile spedi with NCMPMX=number_of_comp in params.f
Use XMINF<0 in kug1.dat (IMF exponent) to generate automatically the mass component with a given IMF exponent. The minimum mass corresponds to the first component in kug1.dat. The max mass corresponds to the second component in the file. The other components are not used but have to be present in the file (dummies). Beware of the format of the file: strictly respect the number of characters for each field. For instance "XIMF=-2.35000+000" is correct but not "XIMF= -2.35000+000"...
(2) In Simula_BH_test/ you'll find an example to start a simulation of a star-accreting BH
(3) In plot/ you'll find scripts to process/ analyse data relative to ndcode/ (check plot/README)
(4) In kug1_examples/ you'll find "good" kug1.dat examples
(5) utils/ and smac/ are tools for plotting etc (smac contains supermongo macros and in utils you'll find scripts)