Table of Contents

Models are not stable of course, but fun to play with or to study collapse
and violent relaxtion. More serious homogeneous models can be found in *mkhom(1NEMO)*
.
Friedman & Polyachenko (1986?, Part I) discuss how to make equilibrium models.

**out=***out_file*- The data are written to this
file in
*SnapShot(5NEMO)*format [no default]. **nbody=***num_bodies*- Number of
stars in sphere [default:
**2048**] **rmin=***min_radius*- Inner cutoff radius of
the sphere. [Default:
**0.0**]. **rmax=***max_radius*- Outer cutoff radius of the sphere.
The default values of
**rmin=0**and**rmax=1.2**give a virial radius of 1.0 [Default:**1.2**]. **2t/w=**- Virial ratio which measures deviation from virial equilibrium
(two times the ratio of kinetic to potential energy,
**2t/w**must be positive). Only accurate for a homogenous sphere (**power=0**and**rmin=0**). See also**vmax=**for an alternate means to specify scaling velocities. [default:**0**]. **vmax=***vmax*- Maximum of isotropic velocities. No scaling as function of radius is applied.
[default:
**0**]. **power=***p*- Power law exponent of the density as function of
radius, Density
*rho = r^{-p}*, thus centrally concentrated models have*p > 0*. Values of*p > 3*are illegal. [default:**0**]. **seed=***random_seed*- Use random
number seed. Non-positive have special meaning. [default:
**0**]. **zerocm=****t|f**- Center
coordinates on center of mass after generation? [default:
**t**]. **headline=***message*- Text headline for output file [default: not used].

20-mar-89V1.0 createdPJT 28-mar-89V1.1 keyword ’rmin’ addedPJT 5-apr-89V1.2 keyword 2t/w implemented when rmin=0PJT 6-apr-89V1.3 keyword zerocm addedPJT 20-may-97V1.4 added vmax=, power=PJT