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

Complete the binary formation process by adding binary orbits to
an existing binary tree.  It is assumed that the binary tree
structure and all masses have already been set (e.g. by
makesecondary), and that positions and velocities of all top-level
nodes are already known.

If possible, use the system parameters to scale the binary parameters.
If the total system energy is already known (saved in the snapshot
dyn story), then energies are in units of kT.  Otherwise, energies
are in absolute units.

Usage:  makeplanetary [OPTIONS]

Options:
           -f    function select option [3]

                 1: angular momentum per unit reduced mass
                 (L^2 = am[1-e^2]), solar units;

                 2: semi-major axis or peri/apo, solar units;

                 3: energy

           -l    lower limit on selected binary parameter [1]

           -o    specify interpretation of limits [1]

                 (-f 1)
                 1: angular momentum,
                 2: angular momentum, detached binary

                 (-f 2)
                 1: semi-major axis,
                 2: semi-major axis, detached,
                 3: -l gives peri, -u gives apo, and detached

                 (-f 3)
                 1: |binary energy|,
                 2: |binary energy| per unit reduced mass,
                 3: |binary energy| per unit binary mass

           -s    specify random seed [take from system clock]
           -u    upper limit on selected binary parameter [1]

Written by Steve McMillan and Simon Portegies Zwart..

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