Click here to download a gzipped tar file of the Clustin code.
The code is a piece of a larger piece of software described in Sigurdsson & Phinney, ApJS v 99 p 609 (1995) The code creates a multi-mass King model, by recursively solving for the density profile of N separate mass groups. Typically the code generates a grid of N mass groups on a radial grid of 100-200 points. The code generates an N-body realisation of the density as well as files giving the intrinsic and projected distributions. The IMF is set by "fremnts.in" which is an input file to subroutine fremnts.f. The subroutine maps some ZAMS IMF to an evolved mass function - current choices have linear WD-mass MS-mass relation, and set critical mass for type II SN at 8 Msun and 25 Msun - above 25 Msun stars become black holes Input parameters are: fret - fraction of neutron stars retained, likely value ~ 0.1 or less x - the Salpeter IMF slope - canonical value 1.35 eps - tolerance parameter for IMF calculation Turnoff mass - highest current main sequence mass fmass - determines if there is a break in the IMF (x=0 for mass < fmass in current formulation) lowm, highm - upper, lower ZAMS IMF masses nmgrp - number of mass bins - recommend ~ 10 lower - lower boundary of lowest mass mass group, usually same as lowm Sequence of nmgrp monotonically increasing numbers showing mass group boundaries ----------------------- The input file for the Cluster structure parameters is "runtri.d" Most of the parameters in runtri.d should not be adjusted! The adjustable parameters are: sigsq - this is the 1-D mean velocity dispersion in the core in (cm/s)^2 dens0 - this is the mean NUMBER DENSITY of stars in the center of the core w0 - this is the King concentration parameter - note that it is not the same as the single mass King parameter, it measures the ratio of the dynamical core of a mean ensemble of stars to the tidal radius - for multiple mass models the visible projected core can signifcantly differ from the dynamical mean core iseed - a flag for the random number generator initialisation
Last updated 10/00
Back to the top of my home page.