Starlab version 4.4.1
    program created on Apr 27 2005 at 22:54:28
    source file     /home/steve/starlab/src/node/dyn/util/radial_profile.C

Compute the specified 1-dimensional radial profile of an N-body system.
The function fills in a profile array corresponding to the provided
radius array, then prints out radius and profile data in a form
suitable for plotting.  If a current density center is found in the root
dyn story it is used as the system center for the density computation.
Otherwise, if a valid center of mass is found, it is used.  If neither
center is found, the geometric center is used.

Usage: radial_profile [OPTIONS] < input > output

Options:     
		  -n    specify number of radial bins [100]
             -o    choose the quantity computed [0]
                   0 ==> density of all stars,
                   1 ==> density of stars having mass >= average,
                   2 ==> velocity dispersion,
                   3 ==> radial velocity dispersion,
                   4 ==> 2, 3, and radial velocity dispersion too
             -r    specify maximum radius [4]

Written by Steve McMillan.

Report bugs to starlab (at) sns.ias.edu.