create spherical mass models from a text table
This program is a simplified version of ``mkspherical'' utility from the SMILE toolbox, intended to create spherical mass models from a text table with r, m(r) values specifying enclosed mass as a function of radius.
The input is a text file with >=5 (preferrably 10-20) lines with two columns: r, m(r); both must be increasing with line number; if m(r=0) is nonzero, it is considered to be a supermassive black hole at the galaxy center.
The output is a text file with a number of parameters computed as spline-interpolated functions of radius (potential, radial/circular period, distribution function via Eddington inversion formula, diffusion coefficients, etc.)
The methods implemented in this code are based on the following paper:
E.Vasiliev, "A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems", MNRAS, 2013 (in press)