My software


Halo - A small Fortran program to calculate properties (virial radius, circular velocity, and temperature; scaling radius, mass, density and circular velocity; half-mass radius and corresponding crossing time; concentration) of cosmological halos - as a function of virial mass, redshift, and deviation from the mean concentration. Handles both cuspy NFW halos and flat-core Burkert halos. Based on correlations of Bullock et al. (2001).  You'll need the source and const.f to compile the program. The Unix script (containing the WMAP values for the cosmological parameters) can be used to run the code.

King (updated Jan. 12, 2004) - a program to calculate radial profiles (density and 1D velocity dispersion) and generate an N-body snapshot for an isotropic King model (often used to describe structure and kinematics of globular clusters). The Fortran code can be compiled by g77 compiler on modern Unix/Linux OS (you will need also const.f file). Input parameters: number of particles N, King model parameter alpha=psi_0/sgm^2 (typically between 3 and 9 - determines the concentration of the cluster), core radius r0 in kpc, and velocity dispersion parameter sigma in km/s. Outputs: concentration c=log(r_t/r0), tidal radius r_t in kpc, and total mass in M_Sun. Two files are produced: "king.dat", which contains density (M_Sun/kpc^3) and 1D velocity dispersion (km/s) as a function of radius (kpc), and "king_snap.dat", which is an N-body snapshot of the cluster in ASCII Nemo format.

LCDM_z_time- A utility to make conversions between redshift z and time t (Gyr) for a flat LCDM Universe. The Fortran code can be used with the Unix shell script (containing the WMAP values for the cosmological parameters). To compile the Fortran code, you will need the file const.f containing the values of astronomical constants.

Shell finder 1.01 (updated Sept. 17, 2001) - a program designed to look for best spherical expanding shell candidates in HI data-cubes. Short description is in the README file. TARed and GZIPed archive (520 KB) contains the source code for Sun Solaris and x86 Linux platforms. The example FITS file can be downloaded here (8 MB). The algorithm is described in Mashchenko & St-Louis (2000).

find_field (updated Apr. 11, 2002) - A small utility for CGPS (Canadian Galactic Plane Survey) users. It allows to find the best mosaic field for given coordinates (galactic, equatorial B1950 or J2000) of an object and a square subcube size {dX,dX,dV}, where dX and dV can be either in physical (dgr,km/s), or pixel units. It also tries to extract the subcube using the program subfits (included into the package). Available as a gzipped tar archive, and includes a Fortran source and binaries for Linux (x86) and Solaris (Sun). Here is the README file.

coords - Program for galactic <-> equatorial (B1950 or J2000) celestial coordinates transformation. It is self-explantory, and is available as a Fortran source code, as a gzipped Linux x86 binary and as a gzipped Solaris Sun binary. If you are looking for a web based coordinates transformation calculator - there's a good one at NED.