Starlab tools

From Modesta

Table of contents

SeBa

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:59:32
   source file     $STARLAB_INSTALL_PATH/src/star/dstar/evolve/SeBa.C

SeBa: Binary evolution program SeBa.

            computes the evolution of a binary given any 
            initial conditions (M, m, a, e).
           

Output: in the form of the following files:

            -init.dat           contains selected initial conditons.
            -SeBa.data           contains binary evolution histories
            -binev.data          contains remnant formation information
            in the form of standard output (cerr).
            Initialized parameters include:
            - mass of the most massice component (primary star)
            - mass of its binary companion (secondary star)
            - semi-major axis of the binary system
            - orbital eccentricity
            routines incuded can be found in double_star.h.
            The mass function routines are adoped from mkmass.C
            and are defined in starbase.h
            externally visible routines are:
            -get_random_mass_ratio
            -get_random_semi_major_axis
            -get_random_eccentricity
            The two utilities for the various parameter are:
            -extract_...._distribution_type_string(....)
            and
            -type_string(char*)
            The executable takes initial conditions (see Options)
            and returns randomized binary parameters.                 
               

Options: -M upper primary mass limit [100[Msun]]

          -m    lower limit to primary mass [0.1[Msun]]
          -x    mass function exponent in case of power law [-2.35]
          -F/f  mass function option: 0) Equal mass
                                      1) Power-law [default]
                                      2) Miller & Scalo
                                      3) Scalo
                                      4) Kroupa
          Option -F requires one of the following strings:
                    (mf_Power_Law, Miller_Scalo, Scalo, Kroupa)
                 -f requires the appropriate interger (see mkmass.C)
           -A   maximum semi-major axis limit [1000000[Rsun]]   
           -a   minimum semi-major axis limit [0] 
           -y   exponent for a power-law distribution [0] (flat in log)
           -G/g Semi major axis option: 0) Equal_sma
                                        1) Power Law [default]
                                        2) Duquennoy & Mayor (1987)
          Option -G requires one of the following strings:
                    (Equal_sma, sma_Power_Law, Duquennoy_Mayor)
                 -g requires appropriate interger (see double_star.h)
           -E   maximum eccentricity [1] 
           -e   minimum eccentricity [0] 
           -v   exponent for a power-law distribution 
           -U/u eccentricity option: 0) Equal eccentricity
                                     1) Power Law 
                                     2) Thermal distribution [default]
          Option -U requires one of the following strings:
                    (Equal_ecc, ecc_Power_Law, Thermal_Distribution)
                 -u requires appropriate interger (see double_star.h)
           -Q   maximum mass ratio [1]
           -q   minimum mass ratio [<-m option/selected primary mass>]
           -w   exponent for a power-law distribution  
           -P/p mass ratio option: 0) constant mass ratio
                                     1) Flat_q
                                     2) Power Law 
                                     3) Hogeveen (1992)
          Option -P requires one of the following strings:
                    (Equal_q, Flat_q, qf_Power_Law, Hogeveen)
                 -p requires appropriate interger (see double_star.h)
          -I select input file for reading initial conditions.
             -uses: double_star::dump as input format.  [no default]
          -R select random initial conditions    [false]
             with parameters as discribed above.   
          -n number of binaries to be simulated.  [1]
             Options: -I all binaries in input file are computed.
                      -R the number of binaries indicated.
                      oterwise one binary is simulated with
                      -M, -m, -a, -e as initial conditions.
          -T or -t  binary end time. [35] Myr


_dyn_io

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:16
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/_dyn_/_dyn_io.C

_dyn_io: Starlab _dyn_ class I/O functions.

        Define scan_dyn_story and print_dyn_story for
        the _dyn_ class.

Options: none

aclocal

Usage: aclocal [OPTIONS] ...


Generate `aclocal.m4' by scanning `configure.ac' or `configure.in'

 --acdir=DIR           directory holding config files
 --help                print this help, then exit
 -I DIR                add directory to search list for .m4 files
 --force               always update output file
 --output=FILE         put output in FILE (default aclocal.m4)
 --print-ac-dir        print name of directory holding m4 files
 --verbose             don't be silent
 --version             print version number, then exit

Report bugs to <bug-automake@gnu.org>.

aclocal-1.8

Usage: aclocal [OPTIONS] ...


Generate `aclocal.m4' by scanning `configure.ac' or `configure.in'

 --acdir=DIR           directory holding config files
 --help                print this help, then exit
 -I DIR                add directory to search list for .m4 files
 --force               always update output file
 --output=FILE         put output in FILE (default aclocal.m4)
 --print-ac-dir        print name of directory holding m4 files
 --verbose             don't be silent
 --version             print version number, then exit

Report bugs to <bug-automake@gnu.org>.

add_daughter_node

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:39
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/add_daughter_node.C

Add one extra node to the tree structure in the input snapshot, under a specified node.

Usage: add_daughter_node [OPTIONS] < input > output

Options:

                 -c     add a comment to snapshot [false]
            -e     echo tree structure [false]
            -i     specify index of node to add to [root]
            -j     specify index for new node [none]
            -m     specify node mass [1]
            -s     specify random seed [take from system clock]
            -r     specify node radial position (angle random) [0]
            -v     specify node speed (direction random) [0]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

add_plummer

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:19
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/add_plummer.C

Add Plummer parameters to an input snapshot and write it out again. Do not change the input data. See add_power_law.

Usage: add_plummer [OPTIONS] < input > output

Options:

                  -c    add comment [none]
             -C    specify center [(0,0,0)]
             -f    turn on/off friction [false]
             -m/M  specify mass [1]
             -a/R  specify scale [1]
             -n    force interpretation of parameters in N-body units [no]

Written by Steve McMillan

Report bugs to starlab@sns.ias.edu.

add_power_law

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:23
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/add_power_law.C

Add power-law parameters to an input snapshot and write it out again. Do not change the input data. External field is basically that of a power-law mass distribution

               M(r)  =  A r^x

except that the density is constant for r < a. Parameters are interpreted in physical units if they have been specified with add_star, are converted to N-body units if necessary, and are written to the root log story for subsequent use by set_com, scale, kira, etc. Note that -e 0 is currently implemented as a Plummer field, and is flagged as such.

Usage: add_power_law [OPTIONS] < input > output

Options:

             -A/M  specify coefficient [1]
             -c    add comment [none]
             -C    specify center [(0,0,0)]
             -e/E/x/X  specify exponent [0]
             -a/R  specify scale [1]
             -G    select parameters (physical units) appropriate
                   for the Galactic center (Mezger et al. 1999) [no]
             -n    force interpretation of all parameters in
                   N-body units [no]

Written by Steve McMillan.

Report bug to starlab@sns.ias.edu.

add_star

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:47
   source file     $STARLAB_INSTALL_PATH/src/star/sstar/init/add_star.C

addstar: Add physical stellar data to existing node structure.

Options: -c comment to put in the starbase log structure [none]

           -M    mass scaling (system mass unit, in solar units)
           -Q    virial ratio (used if time scaling is omitted) [0.5]
           -R    dynamical size scaling for stars
                     (system length unit, in parsecs)
           -s    initial types of stars [main_sequence]
           -T    dynamical time scaling  (system time unit, in Myr)
           -t    stellar age (in millions of years) [0]

Version 1.0: Apr 1993 Piet Hut, Steve McMillan, Jun Makino Version 2.0: May 1993 Simon Portegies Zwart

add_tidal

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:27
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/add_tidal.C

Add tidal parameters to an input snapshot and write it out again. Do not change the input data. Use before scaling.

Usage: add_tidal [OPTIONS] < input > output

Options:

             -c    add comment [none]
             -F    specify tidal_field_type [1]
             -J    specify Jacobi radius scaling factor [no default]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

addchydro

usage: addchydro [-R #] [-r #] [-m #] [-c ".."]

 for effective radius, core radius, core mass

adddouble

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:59:10
   source file     $STARLAB_INSTALL_PATH/src/star/dstar/stardyn/adddouble.C

adddouble: Adds double star class to excisting node structure.

Options: -A maximum initial semi-major axis [in solar units]

           -a    minimum initial semi-major axis [in solar units]
           -c    comment to put in the starbase log structure.
           -E    maximum initial eccentricity
           -e    minimum initial eccentricity
           -M    Mass scaling [total cluster mass in solar units].
           -Q    Virial ratio (used if time scaling is omitted)
           -R    Dynamical size scaling for the star
                 [in units of the virial radius].
           -S    use internal scaling?
           -s    Random seed.
           -T    Dynamical time scaling
                 [in units of the NBODY time].
           -t    start time [in million years].

Version 1.0: Feb 1995 Simon Portegies Zwart

addhydro

usage: addhydro [-R #] [-c ".."]

apply_plummer

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:09
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/apply_plummer.C

Apply a Plummer-model spatial and velocity distribution to an existing snapshot, without scaling. Masses are left unchanged.

Usage: apply_plummer [OPTIONS] < input > output

Options:

       -s    specify random seed [random from system clock]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

autoconf

Usage: ./autoconf [OPTION] ... [TEMPLATE-FILE]

Generate a configuration script from a TEMPLATE-FILE if given, or `configure.ac' if present, or else `configure.in'. Output is sent to the standard output if TEMPLATE-FILE is given, else into `configure'.

Operation modes:

 -h, --help                print this help, then exit
 -V, --version             print version number, then exit
 -v, --verbose             verbosely report processing
 -d, --debug               don't remove temporary files
 -f, --force               consider all files obsolete
 -o, --output=FILE         save output in FILE (stdout is the default)
 -W, --warnings=CATEGORY   report the warnings falling in CATEGORY [syntax]

Warning categories include:

 `cross'         cross compilation issues
 `obsolete'      obsolete constructs
 `syntax'        dubious syntactic constructs
 `all'           all the warnings
 `no-CATEGORY'   turn off the warnings on CATEGORY
 `none'          turn off all the warnings
 `error'         warnings are error

The environment variable `WARNINGS' is honored.

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

Tracing:

 -t, --trace=MACRO     report the list of calls to MACRO
 -i, --initialization  also trace Autoconf's initialization process

In tracing mode, no configuration script is created.

Report bugs to <bug-autoconf@gnu.org>.

autoheader

Usage: ./autoheader [OPTION] ... [TEMPLATE-FILE]

Create a template file of C `#define' statements for `configure' to use. To this end, scan TEMPLATE-FILE, or `configure.ac' if present, or else `configure.in'.

 -h, --help               print this help, then exit
 -V, --version            print version number, then exit
 -v, --verbose            verbosely report processing
 -d, --debug              don't remove temporary files
 -f, --force              consider all files obsolete
 -W, --warnings=CATEGORY  report the warnings falling in CATEGORY

Warning categories include:

 `cross'         cross compilation issues
 `gnu'           GNU coding standards (default in gnu and gnits modes)
 `obsolete'      obsolete features or constructions
 `override'      user redefinitions of Automake rules or variables
 `portability'   portability issues
 `syntax'        dubious syntactic constructs (default)
 `unsupported'   unsupported or incomplete features (default)
 `all'           all the warnings
 `no-CATEGORY'   turn off warnings in CATEGORY
 `none'          turn off all the warnings
 `error'         treat warnings as errors

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

Report bugs to <bug-autoconf@gnu.org>.

autom4te

Usage: ./autom4te [OPTION] ... [FILES]

Run GNU M4 on the FILES, avoiding useless runs. Output the traces if tracing, the frozen file if freezing, otherwise the expansion of the FILES.

If some of the FILES are named `FILE.m4f' they are considered to be M4 frozen files of all the previous files (which are therefore not loaded). If `FILE.m4f' is not found, then `FILE.m4' will be used, together with all the previous files.

Some files may be optional, i.e., will only be processed if found in the include path, but then must end in `.m4?'; the question mark is not part of the actual file name.

Operation modes:

 -h, --help               print this help, then exit
 -V, --version            print version number, then exit
 -v, --verbose            verbosely report processing
 -d, --debug              don't remove temporary files
 -o, --output=FILE        save output in FILE (defaults to `-', stdout)
 -f, --force              don't rely on cached values
 -W, --warnings=CATEGORY  report the warnings falling in CATEGORY
 -l, --language=LANG      specify the set of M4 macros to use
 -C, --cache=DIRECTORY    preserve results for future runs in DIRECTORY
     --no-cache           disable the cache
 -m, --mode=OCTAL         change the non trace output file mode (0666)
 -M, --melt               don't use M4 frozen files

Languages include:

 `Autoconf'   create Autoconf configure scripts
 `Autotest'   create Autotest test suites
 `M4sh'       create M4sh shell scripts
 `M4sugar'    create M4sugar output

Warning categories include:

 `cross'         cross compilation issues
 `gnu'           GNU coding standards (default in gnu and gnits modes)
 `obsolete'      obsolete features or constructions
 `override'      user redefinitions of Automake rules or variables
 `portability'   portability issues
 `syntax'        dubious syntactic constructs (default)
 `unsupported'   unsupported or incomplete features (default)
 `all'           all the warnings
 `no-CATEGORY'   turn off warnings in CATEGORY
 `none'          turn off all the warnings
 `error'         treat warnings as errors

The environment variable `WARNINGS' is honored.

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

Tracing:

 -t, --trace=MACRO      report the MACRO invocations
 -p, --preselect=MACRO  prepare to trace MACRO in a future run

Freezing:

 -F, --freeze   produce an M4 frozen state file for FILES

Report bugs to <bug-autoconf@gnu.org>.

automake

Usage: ./automake [OPTION] ... [Makefile]...

Generate Makefile.in for configure from Makefile.am.

Operation modes:

     --help               print this help, then exit
     --version            print version number, then exit
 -v, --verbose            verbosely list files processed
     --no-force           only update Makefile.in's that are out of date
 -W, --warnings=CATEGORY  report the warnings falling in CATEGORY

Dependency tracking:

 -i, --ignore-deps      disable dependency tracking code
     --include-deps     enable dependency tracking code

Flavors:

     --cygnus           assume program is part of Cygnus-style tree
     --foreign          set strictness to foreign
     --gnits            set strictness to gnits
     --gnu              set strictness to gnu

Library files:

 -a, --add-missing      add missing standard files to package
     --libdir=DIR       directory storing library files
 -c, --copy             with -a, copy missing files (default is symlink)
 -f, --force-missing    force update of standard files

Warning categories include:

 `gnu'           GNU coding standards (default in gnu and gnits modes)
 `obsolete'      obsolete features or constructions
 `override'      user redefinitions of Automake rules or variables
 `portability'   portability issues
 `syntax'        dubious syntactic constructs (default)
 `unsupported'   unsupported or incomplete features (default)
 `all'           all the warnings
 `no-CATEGORY'   turn off warnings in CATEGORY
 `none'          turn off all the warnings
 `error'         treat warnings as errors

Files which are automatically distributed, if found:

 ABOUT-GNU           README              config.rpath        ltcf-gcj.sh
 ABOUT-NLS           THANKS              config.sub          ltconfig
 AUTHORS             TODO                configure           ltmain.sh
 BACKLOG             acconfig.h          configure.ac        mdate-sh
 COPYING             aclocal.m4          configure.in        missing
 COPYING.DOC         ansi2knr.1          depcomp             mkinstalldirs
 COPYING.LESSER      ansi2knr.c          elisp-comp          py-compile
 COPYING.LIB         compile             install-sh          stamp-vti
 ChangeLog           config.guess        libversion.in       texinfo.tex
 INSTALL             config.h.bot        ltcf-c.sh           ylwrap
 NEWS                config.h.top        ltcf-cxx.sh

Report bugs to <bug-automake@gnu.org>.

automake-1.8

Usage: ./automake-1.8 [OPTION] ... [Makefile]...

Generate Makefile.in for configure from Makefile.am.

Operation modes:

     --help               print this help, then exit
     --version            print version number, then exit
 -v, --verbose            verbosely list files processed
     --no-force           only update Makefile.in's that are out of date
 -W, --warnings=CATEGORY  report the warnings falling in CATEGORY

Dependency tracking:

 -i, --ignore-deps      disable dependency tracking code
     --include-deps     enable dependency tracking code

Flavors:

     --cygnus           assume program is part of Cygnus-style tree
     --foreign          set strictness to foreign
     --gnits            set strictness to gnits
     --gnu              set strictness to gnu

Library files:

 -a, --add-missing      add missing standard files to package
     --libdir=DIR       directory storing library files
 -c, --copy             with -a, copy missing files (default is symlink)
 -f, --force-missing    force update of standard files

Warning categories include:

 `gnu'           GNU coding standards (default in gnu and gnits modes)
 `obsolete'      obsolete features or constructions
 `override'      user redefinitions of Automake rules or variables
 `portability'   portability issues
 `syntax'        dubious syntactic constructs (default)
 `unsupported'   unsupported or incomplete features (default)
 `all'           all the warnings
 `no-CATEGORY'   turn off warnings in CATEGORY
 `none'          turn off all the warnings
 `error'         treat warnings as errors

Files which are automatically distributed, if found:

 ABOUT-GNU           README              config.rpath        ltcf-gcj.sh
 ABOUT-NLS           THANKS              config.sub          ltconfig
 AUTHORS             TODO                configure           ltmain.sh
 BACKLOG             acconfig.h          configure.ac        mdate-sh
 COPYING             aclocal.m4          configure.in        missing
 COPYING.DOC         ansi2knr.1          depcomp             mkinstalldirs
 COPYING.LESSER      ansi2knr.c          elisp-comp          py-compile
 COPYING.LIB         compile             install-sh          stamp-vti
 ChangeLog           config.guess        libversion.in       texinfo.tex
 INSTALL             config.h.bot        ltcf-c.sh           ylwrap
 NEWS                config.h.top        ltcf-cxx.sh

Report bugs to <bug-automake@gnu.org>.

autoreconf

Usage: ./autoreconf [OPTION] ... [CONFIGURE-AC or DIRECTORY] ...

Run `autoconf' (and `autoheader', `aclocal', `automake', `autopoint' (formerly `gettextize'), and `libtoolize' where appropriate) repeatedly to remake the GNU Build System files in the DIRECTORIES or the directory trees driven by CONFIGURE-AC (defaulting to `.').

By default, it only remakes those files that are older than their predecessors. If you install new versions of the GNU Build System, running `autoreconf' remakes all of the files by giving it the `--force' option.

Operation modes:

 -h, --help               print this help, then exit
 -V, --version            print version number, then exit
 -v, --verbose            verbosely report processing
 -d, --debug              don't remove temporary files
 -f, --force              consider all files obsolete
 -i, --install            copy missing auxiliary files
 -s, --symlink            with -i, install symbolic links instead of copies
 -m, --make               when applicable, re-run ./configure && make
 -W, --warnings=CATEGORY  report the warnings falling in CATEGORY [syntax]

Warning categories include:

 `cross'         cross compilation issues
 `gnu'           GNU coding standards (default in gnu and gnits modes)
 `obsolete'      obsolete features or constructions
 `override'      user redefinitions of Automake rules or variables
 `portability'   portability issues
 `syntax'        dubious syntactic constructs (default)
 `unsupported'   unsupported or incomplete features (default)
 `all'           all the warnings
 `no-CATEGORY'   turn off warnings in CATEGORY
 `none'          turn off all the warnings
 `error'         treat warnings as errors

The environment variable `WARNINGS' is honored. Some subtools might support other warning types, using `all' is encouraged.

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

The environment variables AUTOCONF, AUTOHEADER, AUTOMAKE, ACLOCAL, AUTOPOINT, LIBTOOLIZE are honored.

Report bugs to <bug-autoconf@gnu.org>.

autoscan

Usage: ./autoscan [OPTION] ... [SRCDIR]

Examine source files in the directory tree rooted at SRCDIR, or the current directory if none is given. Search the source files for common portability problems, check for incompleteness of `configure.ac', and create a file `configure.scan' which is a preliminary `configure.ac' for that package.

 -h, --help          print this help, then exit
 -V, --version       print version number, then exit
 -v, --verbose       verbosely report processing
 -d, --debug         don't remove temporary files

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

Report bugs to <bug-autoconf@gnu.org>.

autoscan.log

tcsh: ./autoscan.log: Permission denied.

autoupdate

Usage: ./autoupdate [OPTION] ... [TEMPLATE-FILE...]

Update the TEMPLATE-FILE... if given, or `configure.ac' if present, or else `configure.in', to the syntax of the current version of Autoconf. The original files are backed up.

Operation modes:

 -h, --help                 print this help, then exit
 -V, --version              print version number, then exit
 -v, --verbose              verbosely report processing
 -d, --debug                don't remove temporary files
 -f, --force                consider all files obsolete

Library directories:

 -B, --prepend-include=DIR  prepend directory DIR to search path
 -I, --include=DIR          append directory DIR to search path

Report bugs to <bug-autoconf@gnu.org>.

chydro_leapfrog

usage: chydro_leapfrog -t # -a # -e # -D # -d # [-c ".."] for t (time span), a (time step length), d (output interval) D (snapshot output interval) and e (softening length)

coagulation

Usage: ./coagulation [-N #] [-n #] [-M] [-m] [-e] [-o] [-T #] [-v]

compute_com

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:12
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/compute_com.C

Determine the center of mass position and velocity of the input N-body system. Center of mass position and velocity are written to the dyn story of the top-level node; they are also optionally returned as function arguments in the library version. Note: The computed center of mass is defined in absolute terms, and so includes the pos and vel of the parent node.

Usage: compute_com [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

compute_density

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:21
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/compute_density.C

Compute local densities around particles in the input N-body system, based on k-th nearest neighbors. Save results in the particle dyn stories.

Usage: compute_density [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -k    specify k [12]
            -v    verbose mode [off]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

compute_max_cod

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:25
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/compute_max_cod.C

Determine the max density center position and velocity for the input N-body system. The max density center is defined as the position of the particle with the highest local density. Densities are not computed here -- run compute_density before invoking compute_max_cod. Density center position and velocity are written to the dyn story of the top-level node; they are also optionally returned as function arguments in the library version. Note: The computed density center is defined in absolute terms, .and so includes the pos and vel of the parent node.

Usage: compute_max_cod [OPTIONS] < input > output

Options:

                  -c    add a comment to the output snapshot [false]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

compute_mcom

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:17
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/compute_mcom.C

Determine the "modified center of mass" position and velocity of the input N-body system, recursively defined as the center of mass of the f [default 90%] particles closest to the center of mass. Center of mass position and velocity are written to the dyn story of the top-level node; they are also optionally returned as function arguments in the library version. Note: The computed center of mass is defined in absolute terms, and so includes the pos and vel of the parent node.

Usage: compute_mcom [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -f    specify the fraction of particles to retain [0.9]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

compute_mean_cod

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:29
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/compute_mean_cod.C

Determine the density center position and velocity for the input N-body system. The density center is defined analogously to the center of mass, but instead of a mass weighting, the weighting used is proportional to the square of the density (the original suggestion in the paper below(*), to use a weighting factor linear in the density, may not converge well). Densities are not computed here -- run compute_density before invoking compute_mean_cod. Density center position and velocity are written to the dyn story of the top-level node; they are also optionally returned as function arguments in the library version.

(*) Stefano Casertano and Piet Hut: Astroph.J. 298,80 (1985). To use their recipe, k >= 2 is required.

Note: The computed density center is defined in absolute terms, and so includes the pos and vel of the parent node.

Usage: compute_mean_cod [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]

Written by Piet Hut, Steve McMilland, and Jun Makino.

Report bugs to starlab@sns.ias.edu.


cutrandomsample

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:02
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/cutrandomsample.C

Create a linked list of equal-mass nodes.

Usage: cutrandomsample [OPTIONS]

Options:

     -m       specify total mass [1]
     -n       specify number of nodes [1]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

display_subtrees

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:13
   source file     $STARLAB_INSTALL_PATH/src/node/util/display_subtrees.C

Output all subtrees in input snapshot(s), neglecting the root and top-level nodes.

Usage: mass_dist [OPTIONS]

Options:

       -c    add a comment to snapshot [false]

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

dstar_state

Hello

dstar_stats

Usage: ./dstar_stats [-b] [-n] [-e] [-o] [-s] [-t]

dumbp

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:33
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/dumbp.C

Dump out an N-body system in a dumb format suitable for digestion by NBODY1-5 and starcluster. This is the inverse function to readp.

Usage: dumbp [OPTIONS] < input > output

Options:

                 -p    specify precision of output [6 sig. fig.]
            -t    include time in output [no]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

dyn_io

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:37
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/dyn_io.C

Test Starlab dyn class I/O functions. Define scan_dyn_story and print_dyn_story for the dyn class.

Usage: dyn_io [OPTIONS] < input > output

Options: None.

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

dyndiff

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:41
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/dyndiff.C

Compute the rms distance in phase space between two (flat) N-body systems.

Usage: dyndiff [OPTIONS] < input > output

Options:

                 -r    use spatial differences only.

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

dyndist

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:46
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/dyndist.C

Get statistics over N 6-dimensional particle pair separations. Compute statistics for the separations between corresponding particles in two N-body systems, using both configuration space and velocity data. The distribution of particle separations are reported as n-tiles.

Usage: dyndist [OPTIONS] < input > output

Options:

                 -n   specify number of n-tiles [4]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

encountervelocity

Usage: ./encountervelocity [-n #] [-M] [-m] [-e] [-o] [-r #] [-v]

energy

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:50
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/energy.C

Print out the energy of an N-body system. Does not include the effects of any external tidal field.

Usage: energy [OPTIONS] < input > output

Options:

                 -e    specify softening parameter [0]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

extract_snap

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:54
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/extract_snap.C

Find and print the last snapshot, or the first snapshot following a specified time. If -n is set to 0 and -v is set, then all we do is count snapshots.

Usage: extract_snap [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -t    specify time [last snapshot]
            -n    specify number of snapshots to extract [1]
            -v    set verbose mode [off]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

flat_hermite

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:11
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/_dyn_/flat_hermite.C


flatten

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:59
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/flatten.C

Flatten a dyn tree to a single-level linked list under the root node.

Usage: flatten [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -C    force col output [take from input format]
            -v    print diagnostic info [no info]

Written by Piet Hut

Report bugs to starlab@sns.ias.edu.

freeze

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:03
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/freeze.C

Reduce all top-level velocities, leaving positions unchanged.

Usage: freeze [OPTIONS] < input > output

Options:

            -c    add a comment to the output snapshot [false]
            -f    specify freeze factor [0]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

freezebin

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:07
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/freezebin.C

Reduce all top-level binary CM velocities, leaving positions unchanged.

Usage: freezebin [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -f    specify freeze factor [0]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

get_densities

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:57:02
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/util/get_densities.C

Compute densities of top-level nodes. Densities are saved in node dyn stories. The program only actually does the computation if a GRAPE is attached.

Usage: get_densities

Options: None.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

get_std_center

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:15
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/get_std_center.C

Determine the "standard" center of the input N-body system, determined as follows: (1) if densities are computed and up to date, use the density center (but *don't* recompute the densities), as determined by compute_max_cod(), (2) otherwise, use the modified center of mass, as returned by compute_mcom(). Center position and velocity are written to the dyn story of the top-level node; they are also optionally returned as function arguments in the library version. Note: The computed center is defined in absolute terms, and so includes the pos and vel of the parent node.

Usage: get_std_center [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

gethist

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:25
   source file     $STARLAB_INSTALL_PATH/src/std/gethist.C

Check Starlab function for saving command line and date of execution.

Options: None.

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

getsnap

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:43
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/getsnap.C

Select a specific snapshot from the input stream.

Usage: getsnap [OPTIONS] < input > output

Options:

                 -c    add comment to snapshots [false]
            -t    specify time of snapshot [0]

Written by Simon Portieges Zwart.

Report bugs to starlab@sns.ias.edu.

grow_black_hole

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:52
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/grow_black_hole.C

Replace the specified inner mass in stars with a single black hole of the same mass.

Usage: grow_black_hole [OPTIONS]

Options:

      -M    select black hole mass (if <1 mass is read as a fraction)

Written by Simon Portegies Zwart

Report bugs to starlab@sns.ias.edu.

hdyn_pp3

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:32
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/util/hdyn_pp3.C

Starlab hdyn-specific debugging tools. Perform pp3_tree on all nodes, descending from the root node.

Usage: hdyn_pp3 [OPTIONS]

Options: None.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

hop

Usage: ./hop [-i] [-n #] [-v]

Make sure you run compute_density before running hop. Otherwise, it won't work.

hscale

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:55
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/util/hscale.C

(Re)scale an N-body system to specified M, Q (=T/U), R, or E. This hdyn version is identical to the dyn "scale", except that GRAPE-4 or GRAPE-6 is used if available.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

hsys_stats

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:46
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/util/hsys_stats.C

This is the "hdyn" (kira output) version of sys_stats. It contains all of the "dyn" sys_stats functionality (and will work with standard dyn data files), but includes additional functions specific to the hdyn class.

Print out various diagnostic statistics on the input system. These include:

       system time, number, mass, mass distribution;
       relaxation time;
       system energy;
       core parameters;
       lagrangian radii for quartiles [default], for ten-percentiles,
       and for "special" choice of Lagrangian masses, currently 0.005,
       0.01, 0.02, 0.05, 0.1, 0.25, 0.5, 0.75, 0.9;
       mass distribution by lagrangian zone;
       anisotropy by lagrangian zone;
       binary parameters.

In addition, Lagrangian radii are written to the root dyn story. If N^2_ops is selected, then core parameters and particle densities are also written to the root and particle dyn stories.

If N^2_ops is selected, Lagrangian radii are computed relative to the density center, as are binary radial coordinates. Otherwise, the modified center of mass (center of mass with outliers excluded) is used.

Options:

            -b    specify level of binary statistics [2]
                  0: none
                  1 (or no arg): short binary output
                  2: full binary output
            -e    recalculate the total energy (even if -n is no) [yes]
            -l    specify percentile choice [2]
                  0: quartiles (1-3)
                  1: 10-percentiles (10-90)
                  2 (or no arg): nonlinear Lagrangian masses
                  (0.5, 1, 2, 5, 10, 25, 50, 75, 90%)
            -n    perform/don't perform actions requiring O(N^2)
                  operations (e.g. computation of energy and core
                  radius; see -e) [no]
            -o    pipe system to cout [no]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

hxstarplot

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:57:44
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/xutil/hxstarplot.C

hxstarplot: Plot an N-body system in an X environment.

           This is the "hdyn" (kira output) version of xstarplot.
           It contains all of the "dyn" xstarplot functionality
           (and will work with standard dyn data files), but
           includes additional functions specific to the hdyn class.

Options:

           -a    specify viewing axis [1/2/3 = x/y/z]
           -b    allow backward steps (up to 512 frames retained) [no]
           -C    center each frame on the standard system center [no]
           -d    dimensionality of plot [2]
           -D    delay time between frames, in ms [16]
           -e    color by energy [no]
           -f    solid-color stars [true]
           -l/L  specify limits [take from initial data]
           -m    highlight multiple systems [no]
           -o    pipe cin to cout [no]
           -p    point scale mode [0]
           -P    specify point size (relative to axis/30) [1, min = 0.5]
           -r    "reverse video" (black background) [true]
           -s    window size, relative to "standard" 400 pixels [1]
           -t    show tree structure [no]
           -u    highlight unperturbed systems [no]
           -v    verbose mode (echo snap times) [no]

Note: There appears to be a problem with the X interface on some

      systems that causes xstarplot to crash if too many run-time
      user commands are received in a short space of time (e.g.
      try holding down the ">" key in 3D mode).  Presently, the
      only option is to avoid excessive keystrokes in the graphics
      window...

hydro_leapfrog

usage: hydro_leapfrog -t # -a # -e # -D # -d # [-c ".."] for t (time span), a (time step length), d (output interval) D (snapshot output interval) and e (softening length)

ifnames

Usage: ./ifnames [OPTION] ... [FILE] ...

Scan all of the C source FILES (or the standard input, if none are given) and write to the standard output a sorted list of all the identifiers that appear in those files in `#if', `#elif', `#ifdef', or `#ifndef' directives. Print each identifier on a line, followed by a space-separated list of the files in which that identifier occurs.

 -h, --help      print this help, then exit
 -V, --version   print version number, then exit

Report bugs to <bug-autoconf@gnu.org>.

imf_evolve

Usage: ./imf_evolve [-N #] [-n #] [-M #] [-m #] [-T #] [-x #] [-r #] [-v #]

jiggle

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:07
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/jiggle.C

Add a small, energy-preserving perturbation to all top-level velocities, leaving positions unchanged.

Usage: jiggle [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -f    specify (fractional) scale of perturbation [0.01]
            -s    specify random seed [take from system clock]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

kepler

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:35
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/kepler/kepler.C

kepler: Starlab kepler structure functions. Initialize an orbit

        from mass, pos, and vel, or mass, semi-major axis and
        eccentricity, and allow the user to manipulate the
        resulting structure interactively, or animate the orbit.

Command-line arguments:

        -a    specify semi-major axis [1]
        -A    animate to apocenter [no]
        -d    specify animation timestep [1/64 orbit]
        -e    specify eccentricity [0]
        -m    specify total mass or primary mass (animation only) [1]
        -m2   specify mass of second component (animation only) [primary]
        -M    specify mean anomaly [0]
        -O    turn on animation [no]
        -P    animate to pericenter [no]
        -r(R) specify sum of radii (to stop animation) [0]
        -t    specify time [0]
        -T    specify end time of animation [0]
        -x    specify separation vector [(1,0,0)]
        -v    specify relative velocity vector [(0,1,0)]

kingfit

Usage: ./kingfit [-L #] [-l #] [-R #] [-r #]

kira

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:57:09
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/evolve/kira.C

Hermite N-body integrator with evolving hierarchical tree structure, stellar and binary evolution, and an arbitrary external field. The program reads a snapshot from standard input and writes snapshot(s) to standard output. Periodic log output is sent to standard error.

Basic options are listed below. Default values for a new simulation are indicated in square brackets. For restart (continuation of a previous kira calculation), the defaults for many options [*] are determined from the input snapshot, making it possible to continue a run without having to re-specify all the command-line parameters previously used. (These new defaults may still be overridden on the command line, of course.) Some options typically having to do with initial setup may be overridden by data from the input snapshot (if present), as noted below. Kira may also turn *on* some options (the B, G, S, and u settings) if they were turned on in the previous run, but are not specified on the current command line. To prevent this, use the "-o" switch.

The first page of log output gives detailed information on all parameter settings adopted and any modifications made during initialization. In case of doubt, read the log file!

Usage: kira [OPTIONS] < infile > output

Options:

            -0    force non-GRAPE operation, if relevant
            -1    suppress density calculation [compute with GRAPE]
            -2    enable DMA GRAPE access (experimental) [no DMA]
            -3    enable special treatment of isolated multiples [no]
            -a    specify accuracy parameter [0.1][*]
            -A    enable "alternate" output [off]
            -b    specify frequency of full binary output, in (integer)
                  units of the log output interval
                  [10; no argument or 0 ==> no output]
            -B    turn on binary evolution [off][*]
            -c    include comment [none]
            -C    specify GRAPE release interval, in seconds [15]
            -d    specify log output interval [1][*]
            -D    specify snapshot interval [end of run]. 
                  Special values:
                              xN: formatted full dump, frequency N;
                              XN: unformatted full dump, frequency N;
                              full/all: same as x1;
                              b: track binary changes only (formatted);
                              B: track binary changes (unformatted);
            -e    specify softening length [0][*]
            -E    use exact calculation [false]
            -f    turn on/off internal dynamical friction on stars [0][*]
            -F    turn on external dynamical friction on the cluster,
                  and optionally specify a scaling coefficient [none]
            -g    specify hysteresis factor [2.5][*]
            -G    specify initial stripping radius [none][*]
            -h    specify stellar-evolution time step [0.015625 = 1/64][*]
            -i    ignore all internal forces (i.e. external only) [false]
            -I    specify (re)initialization timescale [1][*]
            -k    specify perturbation factor [1.e-7][*]
            -K    specify log2(maximum slowdown factor) (integer): [0][*]
            -l    specify close-encounter distance [0.25 --> 0.25/N][*]
                  [option renamed from -f, 7/04]
            -L    specify CPU time limit, in seconds [none]
            -n    stop at specified number of particles [5]
            -N    specify frequency of CPU check output [50000]
            -o    prevent kira from overriding some settings (BGSu)
                  based on input snapshot data [allow]
            -O    save (and overwrite) extra snapshot at each output [no]
            -q    specify initial virial ratio [0.5]
            -r    specify initial virial radius (may not be
                  specified in the input snap) [no]
            -R    specify snapshot file for (re)start [none: use stdin]
            -s    specify random seed [take from system clock]
            -S    turn on stellar evolution [off][*]
            -t    specify time span of calculation [10]
            -T    enable experimental threading and specify n_threads [0]
            -u    toggle unperturbed multiple motion [disabled][*]
            -U    toggle all unperturbed motion [enabled][*]
            -v    toggle "verbose" mode [on]
            -W    specify full-dump (worldline) timescale [1]
            -x    toggle output of extended-precision time [on]
            -X    specify escaper removal timescale [reinit][*]
            -y    specify stellar encounter criterion
                  [0 N-body units or solar radii][*]
            -z    specify stellar merger criterion [0 stellar radii][*]
            -Z    specify stellar tidal dissipation criterion
                  [0 stellar radii][*]

As a convenient shorthand, any "dt" interval specified less than zero is interpreted as a power of 2, i.e. "-d -3" sets dt_log = 0.125. In the case of dt_snap, this will also cause snapshot data to be written immediately on restart (usually we wait until time dt_snap).

Written by J. Makino, S. McMillan, S. Portegies Zwart, and P. Hut, .

Report bugs to starlab@sns.ias.edu.

kira_smallN

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:57:18
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/evolve/kira_smallN.C

Self-contained few-body integrator, using a fourth-order Hermite scheme incorporating a modified unperturbed treatment of close approaches and time symmetrization. The program reads a snapshot from standard input and optionally writes snapshot(s) to standard output. Optional periodic log output is sent to standard error. This version runs as a standalone program, or incorporated into kira to handle multiple systems containing perturbed hard binaries.

Usage: kira_smallN [OPTIONS] < input > ouptut

Options:

       -a    set accuracy parameter [0.03]
       -d    set log output interval [0 --> no output]
       -D    set snap output interval [0 --> no output]
       -E    set energy output interval [0 --> no output]
       -g    set unperturbed limit [1.e-5]
       -n    set number of symmetrization iterations [1]
       -r    set termination radius [infinite]
       -t    set termination time [200]

Authors: Fan-Chi Lin and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

lagrad

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:11
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/lagrad.C

Compute Lagrangian (mass) radii for an N-body system. The radii are stored in the system (root) dyn story. If a current density center is found in the root dyn story, the Lagrangian radii are calculated relative to it. Otherwise, if a valid center of mass is found, it is used. If neither center is found, the geometric center is used.

Usage: lagrad [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -n    specify number of Lagrangian zones (linear in mass) [4]
            -s    use "special" nonlinear binning:
                      0.005, 0.01, 0.02, 0.05,
                      0.1, 0.25, 0.5, 0.75, 0.9
            -t    same as -n 10

Written by Piet Hut, Steve McMillan, and Jun Makino.

Report bugs to starlab@sns.ias.edu.

lagradplot

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:16
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/lagradplot.C

Compute and plot Lagrangian radii for input N-body system(s), based on the coordinate center. (Probably only marginally useful to use this center -- better to use the density center or modified center of mass, as in lagrad...)

Usage: lagradplot [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -t    use ten-percentiles rather than quartiles [quartiles]
            -d    use double width output, 158 rather than 79 columns
            -q    use quadruple width output, 158 rather than 79 columns

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

lagradplotsmooth

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:11
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/lagradplotsmooth.C

Plot Lagrangian radii for input N-body system(s), based on the geometric center; takes output from "lagradplot -q" and smooths over w output lines. (Probably only marginally useful to use this center -- better to use the density center or modified center of mass, as in lagrad...) Note: This really is a terrible kluge to smooth Lagrangian radii, and in addition it is programmed in an opaque way, just to get some quick and dirty results; it works is all I can say. We really should clean this all up soon(ish) -- hope springs eternal.

Usage: lagradplotsmooth [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -w    number of output lines to smooth over

Examples:

mkplummer -n 128 -i -s 123 | ( kira -t 1000 -d 10 -D 0.2 -n 0 | compute_density | to_cod | lagradplot -q > run128q.out ) >& run128q.log

lagradplotsmooth -w 30 < run128q.out

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

leapfrog

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:06
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/evolve/leapfrog.C

leapfrog: Leapfrog integrator for flat N-body dyn systems.

Options: -a specify timestep [0.02]

          -c    add a comment to the output snapshot [false]
          -d    specify log output time interval [0.25]
          -D    specify snapshot output time interval [none]
          -e    specify softening parameter [0.05]
          -q    echo parameters to cerr [true]
          -t    specify time span of integration [1]
          -x    terminate precisely at specified end time [no]

libtool

Usage: libtool [OPTION]... [MODE-ARG]...

Provide generalized library-building support services.

   --config          show all configuration variables
   --debug           enable verbose shell tracing

-n, --dry-run display commands without modifying any files

   --features        display basic configuration information and exit
   --finish          same as `--mode=finish'
   --help            display this help message and exit
   --mode=MODE       use operation mode MODE [default=inferred from MODE-ARGS]
   --quiet           same as `--silent'
   --silent          don't print informational messages
   --tag=TAG         use configuration variables from tag TAG
   --version         print version information

MODE must be one of the following:

     clean           remove files from the build directory
     compile         compile a source file into a libtool object
     execute         automatically set library path, then run a program
     finish          complete the installation of libtool libraries
     install         install libraries or executables
     link            create a library or an executable
     uninstall       remove libraries from an installed directory

MODE-ARGS vary depending on the MODE. Try `libtool --help --mode=MODE' for a more detailed description of MODE.

Report bugs to <bug-libtool@gnu.org>.

libtoolize

Usage: libtoolize [OPTION]...

Prepare a package to use libtool.

   --automake        work silently, and assume that Automake is in use

-c, --copy copy files rather than symlinking them

   --debug           enable verbose shell tracing

-n, --dry-run print commands rather than running them -f, --force replace existing files

   --help            display this message and exit
   --ltdl            install libltdl in a subdirectory
   --ltdl-tar        install the libltdl tarball
   --version         print version information and exit

You must `cd' to the top directory of your package before you run `libtoolize'.

Report bugs to <bug-libtool@gnu.org>.

list_snap

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:19
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/list_snap.C

Print times of all snapshots in the input stream.

Usage: list_snap [OPTIONS] < input > output

Options: None.

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

lmt_to_ubvri

   Starlab version 4.4.2
   program created on Jul 25 2005 at 18:00:39
   source file     $STARLAB_INSTALL_PATH/src/star/util/lmt_to_ubvri.C
tl_to_ubv: conversion from effective temperature and luminosity to
           Johnston UBVRI.
           For single stars and for binaries.

Options: -T Temperature [Kelvin] of primary star [no default]

            -t    Temperature [Kelvin] of secondary star [no default]
            -L    Luminosity [Lsun] of primary star [no default]
            -l    Luminosity [Lsun] of secondary star [no default]
            -M    Mass [Msun] of primary star [no default]
            -m    Mass [Msun] of secondary star [no default]


m_bar

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:54
   source file     $STARLAB_INSTALL_PATH/src/node/util/m_bar.C

Cdetermine the mean mass of a power-law mass distribution.

Usage: m_bar [OPTIONS]

Options:

          -f        multiply result by this factor [1]
          -e/x      exponent [-2.35 (Salpeter)]
          -l/L      lower mass limit [1]
          -u/U      upper mass limit [1]
           -v       set verbose mode [false]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

make_image

Usage: make_image [-c] [-C comment] [-f] [-i] [-F fmt] [-g] [-p]

      [-m color-map] [-s sizex [sizey]]

make_single_node

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:25
   source file     $STARLAB_INSTALL_PATH/src/node/util/make_single_node.C

Create a single node, to which other nodes can be added later, using `add_daughter_node'.

Usage: make_single_node [OPTIONS]

Options:

           -c  comment            optional comment
           -m  mass               optional mass [1]
           -i                     label nodes

Written by Piet Hut and Peter Teuben.

Report bugs to starlab@sns.ias.edu.

makebinary

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:13
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makebinary.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: makebinary [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
          -e    maximum eccentricity [1]
          -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) [unit: Rsun]
                1: semi-major axis,
                2: semi-major axis, detached and hard
                (note: -u gives an upper limit on sma),
                3: -l gives peri, -u gives apo, and detached
                (-f 3) [N-body units]
                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@sns.ias.edu.

makeblack_hole

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:48
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makeblack_hole.C

Replace the star closest to com by a black hole.

Usage: makeblack_hole [OPTIONS]

Options:

    -c    add comment to snapshot
    -M    select black hole mass; if <1 mass is read as a fraction
    -r    select black hole position
    -v    select black hole velocity [relative to circular speed]
    -i    select black hole identity

Written by Simon Portegies Zwart.

Report bugs to starlab@sns.ias.edu.

makecube

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:40
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makecube.C

Construct a simple homogeneous cube, with (default)

            M = 1, T/U = -1/2, E = -1/4.

If the "-u" flag is set, the particles are left unscaled, with masses 1/n, positions uniformly distributed in [-L, L], and velocities uniformly distributed in a range giving approximate virial equilibrium.

Usage: makecube [OPTIONS]

Options:

            -c    add a comment to the output snapshot [false]
            -C    output data in 'col' format [no]
            -i    number the particles sequentially [don't number]
            -l    write cube size to dyn story [don't write]
            -L    specify cube size (+/-L) [1]
            -n    specify number of particles [no default]
            -o    echo value of random seed [don't echo]
            -s    specify random seed [random from system clock]
            -u    leave unscaled [scale to E=-1/4, M = 1, R = 1]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makedisk

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:45
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makedisk.C

Construct a near-keplerian disk, with N low-mass objects orbiting a single massive object on almost circular paths.

Units: masses are measured in millions of solar masses, lengths are measured in parsecs. Hence (with G = 1) the time unitis 1.49e4 yr and the velocity unit is 65.7 km/s.

Usage: makedisk [OPTIONS]

Options:

            -c    add a comment to the output snapshot [false]
            -C    output data in 'col' format [no]
            -i    number the particles sequentially [don't number]
            -l    specify particle radius [0]
            -m    specify the total mass of the disk [no default]
            -M    specify the mass of the central object [no default]
            -n    specify number of disk particles [no default]
            -o    echo value of random seed [don't echo]
            -r    specify inner disk radius [no default]
            -R    specify outer disk radius [no default]
            -s    specify random seed [random from system clock]
            -v    specify 3-D disk velocity dispersion [0]
            -V    verbose mode [off]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makedyn_from_input

   Starlab version 4.4.2
   program created on Jul 25 2005 at 18:00:46
   source file     $STARLAB_INSTALL_PATH/src/star/util/makedyn_from_input.C

mkdyn_from_input: create a linked list of dyns and read pos etc.

                  from input

Options: -F Specify input filename with the following format:

                   number of stars
                   id x m v t j
        where:     id: identity
                   x:  position vector
                   m:  mass
                   v:  velocity vector
                   t:  time


makeheavystar

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:29
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makeheavystar.C

Double the mass of one or more stars in a snapshot to approximate the presence of a binary.

Usage: makeheavystar [OPTIONS]

Options:

         -f    specify fraction of stars to double [0.1]
         -l    specify lower mass limit (unused) [1]
         -s    specify random seed [random from system clock]

Written by Kimberly Engle.

Report bugs to starlab@sns.ias.edu.

makeking

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:23
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makeking.C

Construct a King model.

Usage: makeking [OPTIONS]

Options:

            -b    specify Steve's rescaling parameter (< 1) [0]
                  [models with b > 0 are just rescaled King models;
                  models with b < 0 approach isothermal spheres
                  as b --> -infinity]
            -c    add a comment to the output snapshot [false]
            -C    output data in 'col' format [no]
            -i    number the particles sequentially [don't number]
            -n    specify number of particles [no default]
                  if unspecified assumes an input snapshot with masses.
            -o    echo value of random seed [don't echo]
            -s    specify random seed [random from system clock]
            -T    test options (print to cerr) [0]
                    1: print King model, with unit central
                    density and core radius, and exit;
                    2: print King model, scaled to unit mass
                    and virial radius, and exit;
                    3: test realized velocity distribution
            -u    leave final N-body system unscaled
                  [scale to E=-1/4, M = 1, R = 1]
            -w    specify King dimensionless depth [no default]

Written by Steve McMillan and Kimberly Engle.

Report bugs to starlab@sns.ias.edu.

makemass

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:08
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/makemass.C

Add a mass spectrum to an input snapshot. Existing node masses are overwritten. Dyn version is copied from the node version.

Usage: makemass [OPTIONS] < input > output

Options:

         -C        force output data to be in 'col' format [no]
         -e/E/x/X  exponent [-2.35 (Salpeter)]
         -F/f      mass function option: 1) Power-law [default]
                                         2) Miller & Scalo
                                         3) Scalo
                                         4) Kroupa
                                         5) GdeMarchi
                   Option -F requires one of the following strings:
                   (Power_Law, Miller_Scalo, Scalo, Kroupa, GdeMarchi).
                   Option -f requires the appropriate integer.
         -i        (re)number stellar index from highest to lowest mass.
         -l/L      lower mass limit [1]
         -m/M      scale to specified total mass [don't scale]
         -u/U      upper mass limit [1]
         -s        random seed

Written by Steve McMillan and Simon Portegies Zwart.

Report bugs to starlab@sns.ias.edu.

makenode

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:39
   source file     $STARLAB_INSTALL_PATH/src/node/util/makenode.C

Create a linked list of equal-mass nodes.

Usage: makenode [OPTIONS]

Options:

       -m       specify total mass [1]
       -n       specify number of nodes [1]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makeplanetary

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:43
   source file     $STARLAB_INSTALL_PATH/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@sns.ias.edu.

makeplummer

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:50
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makeplummer.C

Construct a Plummer model, with a spatial or mass cut-off to ensure finite radius. The new model system is written to standard output. The model system is shifted to place its center of mass at rest at the origin of coordinates. Unscaled systems will be in approximate virial equilibrium, based on the continuum limit.

Usage: makeplummer [OPTIONS]

Options:

        -c    add a comment to the output snapshot [false]
        -C    output data in 'col' format [no]
        -i    number the particles sequentially [don't number]
        -m    specify mass cutoff (for finite radius) [0.999]
        -n    specify number of particles [no default]
        -o    echo value of random seed [don't echo]
        -r    specify radius cutoff [22.804 for default mass cutoff]
        -R    toggle reshuffle of particles to remove correlation
              between index and distance from cluster center [true]
        -s    specify random seed [random from system clock]
        -u    leave unscaled [scale to E=-1/4, M = 1, R = 1]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makepowerlaw

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:56
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makepowerlaw.C

Construct a system of test particles drawn from a power-law mass distribution, in virial equilibrium with the field due to that distribution. The mass distribution is of the form:

             M(r)  =  A r^x            (0 < x < 3)

except that r is modified to include a softening scale parameter:

                   r  -->  sqrt(r^2 + a^2)

The power-law particle distribution is created between some minimum and maximum radii. The inner radius is for numerical convenience, to limit the number of stars on very tight orbits. In addition, we allow the possibility that the power-law potential is cut off within some radius b and that there a point mass at the center. The default mass for the central point is the power-law mass within the cutoff. The point-mass potential is softened to avoid numerical problems. If b is nonzero, a is set to zero.

We assume that x > 0. Particle parameters will be chosen so that the total mass of the N-body system is 1, independent of the actual mass of the background distribution. For now, we are interested only in test particles. The particles are currently always drawn from the power-law distribution, even within the cutoff radius.

The output snapshot will have the external field already enabled, and will contain a flag to disable internal interactions.

Usage: makepowerlaw [OPTIONS]

Options:

            -A    specify the coefficient A [1]
            -a/R  specify scale [1]
            -b    cutoff radius [0]
            -c    add a comment to the output snapshot [false]
            -C    output data in 'col' format [no]
            -d    make orbits disk-like within the cutoff [isotropic]
            -e    softening parameter [b/100]
            -i    number the particles sequentially [don't number]
            -l    minimum radius of the N-body system [0.1]
            -m/M  central mass [mass at cutoff]
            -n    specify number of particles [no default]
            -o    echo value of random seed [don't echo]
            -u    maximum radius of the N-body system [10]
            -s    specify random seed [random from system clock]
            -x    specify exponent [1]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makepyth

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:55
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makepyth.C

Set up a 3-body system corresponding to the Pythagorean problem.

Usage: makepyth [OPTIONS]

Options:

           -c    add a comment to the output snapshot [false]
           -C    output data in 'col' format [no]
           -i    number the particles sequentially [don't number]

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

makesecondary

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:12
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/makesecondary.C

Create binary secondary components for randomly selected stars in the input snapshot, placing the results in a binary tree. Only secondary masses are set here; orbital parameters are set by dyn::makebinary. No attempt is made to maintain any existing mass or energy scaling. Use scale after this function is called, but before makebinary, if it is desired to specify energies in kT units. Dyn version is copied from the node version.

Usage: makesecondary [OPTIONS] < input > output

Options:

            -C    force output data to be in 'col' format [no]
                  (note: loses binary data)
            -f    specify binary fraction [0.1]
                  for stars with mass >= spacified with -m
            -i    use (a,b) as component indices [false]
            -I    don't limit masses to primary mass range [false]
            -l    specify lower limit on mass ratio or 
                  secondary mass [0]
            -M    specify upper limit for primaries to be binaries [inf]
            -m    specify lower limit for primaries to be binaries [0]
            -q    select choice of minimum mass ratio [false]
                  If true, secondary mass ratio is chosen
                  uniformly on [lower_limit, upper_limit].
                  If false, secondary mass is chosen uniformly
                  on [mmin, primary_mass], where mmin and
                  mmax are specified on the command line
            -S    split primary star [false]
            -s    specify random seed [random from system clock]
            -u    specify upper limit on mass ratio or 
                  secondary mass [1 or m_primary]

Written by Steve McMillan and Simon Portegies Zwart.

Report bugs to starlab@sns.ias.edu.

makesphere

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:00
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makesphere.C

Construct a simple homogeneous sphere, with (default)

            M = 1, T/U = -1/2, E = -1/4.

If the "-u" flag is set, the particles are left unscaled, with masses 1/n, positions uniformly distributed in a sphere of radius R, and velocities uniformly distributed in a range giving approximate virial equilibrium.

Usage: makesphere [OPTIONS]

Options:

            -c    add a comment to the output snapshot [false]
            -C    output data in 'col' format [no]
            -i    number the particles sequentially [don't number]
            -l    write sphere radius to dyn story [don't write]
            -n    specify number of particles [no default]
            -o    echo value of random seed [don't echo]
            -R    specify sphere radius [1]
            -s    specify random seed [random from system clock]
            -u    leave unscaled [scale to E=-1/4, M = 1, R = 1]
            -U    leave unscaled and don't place in center of mass
                  frame [scale and center]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

makewrite

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:05
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/makewrite.C

Turn a text file into a snapshot. Used in some demonstrations only. File may be read top to bottom or left to right.

Usage: makewrite [OPTIONS]

Options:

         -c    add a comment to the output snapshot              [false]
         -C    output data in 'col' format                          [no]
         -d    toggle direction: top-bottom/left-right      [left-right]
         -s    random seed                          [from system cloock]
         -v    random velocity scale                                 [0]

Examples: banner -w 40 STARLAB | makewrite -c "STARLAB"

mass_dist

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:51
   source file     $STARLAB_INSTALL_PATH/src/node/util/mass_dist.C

Compute a mass histogram from the input snapshot(s).

Usage: mass_dist [OPTIONS]

Options:

           -b    specify number of bins [25; 0 is OK]
           -v    toggle verbose mode [off]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

merge_low

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:23
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/merge_low.C

Merge all low-level nodes.

Usage: merge_low [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -f    merge the specified fraction [1]
            -o    specify how fraction f is chosen. [1]
                  Option -1: random.
                  Option  0: every 1/f-th binary.
                  Option  n > 0: n-th block of size f.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

merge_snaps

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:17
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/merge_snaps.C

Sequentially merge all snapshots in the input stream. Do *not* rescale masses, positions, or velocities, but always set the system center of mass to 0. Root node information is retained for all snapshots read.

Usage: merge_snaps [OPTIONS] < input > output

Options: None.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

molecules

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:30
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/molecules.C

Print a hierarchical decomposition of a (small!) system.

Usage: molecules [OPTIONS] < input > output

Options: None.

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

nn_distribution

   Starlab version 4.4.2
   program created on Jul 25 2005 at 18:00:51
   source file     $STARLAB_INSTALL_PATH/src/star/util/nn_distribution.C

nn: find nearest neighbor.

        bladiebla

Options: -c add a comment to the output snapshot [false]

no_xreal

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:21
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/_dyn_/no_xreal.C

no_xreal: Rewrite the input with real (not xreal) data.

Options: none

node_io

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:09
   source file     $STARLAB_INSTALL_PATH/src/node/util/node_io.C

Test/check Starlab node I/O functions.

Options: None.

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

nstarev

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:59:04
   source file     $STARLAB_INSTALL_PATH/src/star/sstar/evolve/nstarev.C

nstarev: evolve a cluster of single stars.

        the single stars should be provided in the input stream.

Options: -c comment to put in the starbase log structure.

           -M    Mass scaling [total cluster mass in solar units].
           -n    number of output timesteps (timesteps are taken
                 with constant time intervals) 
           -R    Dynamical size scaling for the star
                 [in units of the virial radius].
           -S    Random seed.
           -s    Initial stellar type [default is main_sequence].
           -T    Dynamical time scaling
                 [in units of the NBODY time].
           -t    end time of the stellar evolution [in Million year].

Latest version (SPZ:1.0) February 1993. &&&&&&&&&&&&&&&&&&&&&&

pgetopt

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:37
   source file     $STARLAB_INSTALL_PATH/src/std/pgetopt.C

Portable and simplified version of getopt() in UNIX, system V. This is a simple driver/test program for the library function. The length of an option is one character, and the argument has to follow the corresponding option, separated by a space (or tab). Options must start with a minus sign, but more than one option can be combined after the same minus sign. Thus, the next four command lines all have the same effect, while the following three are illegal and will give error messages.

            good -a -b 10 -c
                 good -c -a -b 10
                 good -b 10 -ca
                 good -acb 10
                 bad -a -b10 -c
                 bad -ab10c 
                 bad -a -c -b

Usage: pgetopt

Options:

        -a        test: no argument expected
        -b        test: one argument expected
        -c        test: two arguments expected
        -d        test: three arguments expected
        -e        test: gour arguments expected
        -f        test: no argument expected
        -g        test: one optional argument
        -h        test: no argument expected

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

plot_data

Usage: plot_data [-c xcol ycol1 ycol2 ...] [-c[c][xy]] [-C color1...]

                  [-e]  [-h header]  [-i]  [-l xmin xmax ymin ymax] 
                  [-L scale]  [-N nmax]  [-o]  [-O xo yo]  [-p[p]] 
                  [-P point_size]  [-q]  [-s xs ys]  [-S skip]  [-t ntrail] 
                  [-w[w][xy]]  [-W]  [-x x-label]  [-X]  [-y y-label] 
                  [-z zcol] [[-]-help] 

Options:

       -c xcol ycol1...  plot data in xcol horizontally, ycol vertically [1 2]
       -c[c][xy]       crop [x or y] data to plot limits
       -C color1...    specify line/point colors [all black]
       -e              echo current settings
       -h header       specify overall label for plot [none]
       -i              ignore inline commands [don't ignore]
       -l xmin xmax ymin ymax  specify limits for plot [get from data]
       -ll             force plot lines only [plot lines]
       -L scale        specify limits to be +/- scale for both axes
       -N nmax         specify maximum number of points to store [50000]
       -o              echo stdin to stdout [do not echo]
       -O xo yo        specify top left corner of plotting box [150, 50]
       -p              toggle plot points only [plot lines]
       -pp             force plot points only [plot lines]
       -P size         specify point size, in x-axis units [0 ==> pixel]
       -q              suppress most output [don't suppress]
       -Q              quit
       -s xs ys        specify box size [500, 500]
       -S skip         skip leading lines [0]
       -t ntrail       specify number of trailing points [infinite]
       -w[w][xy]       wrap [x or y] data to plot limits
       -W              wait for keyboard input [inline only]
       -x xlabel       specify label for x-axis ["column 'xcol'"]
       -X              clear the display and redraw axes [inline only]
       -y ylabel       specify label for y-axis ["column 'ycol1...'"]
       -z zcol         specify column for color data [none]
       [-]-help                print this help message

plot_stars

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:43
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/plot_stars.C


precision

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:51
   source file     $STARLAB_INSTALL_PATH/src/std/precision.C

Set Starlab output precision using environment variable STARLAB_PRECISION, if available.

Options: None.

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

pretty_print_tree

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:32
   source file     $STARLAB_INSTALL_PATH/src/node/util/pretty_print_tree.C

Print out the tree structure of the input snapshot(s).

Options: None.

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

print_id

unique id of "--help" is 0 (id>>20) = 0, id-((id>>20)<<20) = 0 id_n_clump(id) = 0

print_normal

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:36
   source file     $STARLAB_INSTALL_PATH/src/node/util/print_normal.C

Print out the tree structure of the input snapshot(s) in normal form.

Options: None.

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

print_radii

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:21
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/print_radii.C

Print out indices (if any), masses, and radii (distances from the origin) for an N-body system.

Usage: print_radii [OPTIONS] < input > output

Options:

                 -p    specify precision of output [6 sig. fig.]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

print_time

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:48
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/print_time.C


proj_lagr

   Starlab version 4.4.2
   program created on Jul 25 2005 at 18:00:34
   source file     $STARLAB_INSTALL_PATH/src/star/util/proj_lagr.C

lagrad: Compute Lagrangian (mass) radii for an N-body system.

        The radii are stored in the system (root) dyn story.
        If a current density center is found in the root dyn story,
        the Lagrangian radii are calculated relative to it.
        Otherwise, if a valid center of mass is found, it is used.
        If neither center is found, the geometric center is used.

Options: -c add a comment to the output snapshot [false]

            -n    specify number of Lagrangian zones (linear in mass) [4]
            -s    use "special" nonlinear binning:
                      0.005, 0.01, 0.02, 0.05,
                      0.1, 0.25, 0.5, 0.75, 0.9
            -t    same as -n 10

radial_density

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:51
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/radial_density.C

Compute the 1-dimensional radial density profile of an N-body system. The function fills in the density array corresponding to the provided radius array. The tool prints out radius and density 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_density [OPTIONS] < input > output

Options:

                 -n    specify number of radial bins [100]
            -r    specify maximum radius [take from data]

Written by Steve McMillan and Ernest Mamikonyan.

Report bugs to starlab@sns.ias.edu.

radial_numdensity

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:55
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/radial_numdensity.C

Compute the 1-dimensional radial density profile of an N-body system. The function fills in the density array corresponding to the provided radius array. The tool prints out radius and density 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_numdensity [OPTIONS] < input > output

Options:

                 -n    specify number of radial bins [100]
            -r    specify maximum radius [take from data]

Written by Steve McMillan and Ernest Mamikonyan.

Report bugs to starlab@sns.ias.edu.

radial_profile

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:54:04
   source file     $STARLAB_INSTALL_PATH/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@sns.ias.edu.

radial_vdisp

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:59
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/radial_vdisp.C

Compute the 1-dimensional radial velocity dispersion profile of an N-body system. The function fills in the vdisp array corresponding to the provided radius array. The tool prints out radius and velocity dispersion 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_vdisp [OPTIONS] < input > output

Options:

                 -n    specify number of radial bins [100]
            -r    specify maximum radius [take from data]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

randinter

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:43
   source file     $STARLAB_INSTALL_PATH/src/std/randinter.C

Check the Starlab random number generators. Print out a specified number (N) of random numbers uniformly distributed on [0,1).

Usage: randinter N seed

Options:

   first argument = number of "throws of the dice," N  [no default]
   second argument = random seed  [take from system clock]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

rcore2

Usage: ./rcore2 [-n #] [-M] [-m] [-e] [-o] [-v]

readp

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:33
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/readp.C

Convert ASCII "dumbp" format:

                  (mass1, pos1, vel1,
                   mass2, pos2, vel2,
                   mass3, pos3, vel3,
                  etc.)

data into a Starlab snapshot (flat tree). This is the inverse function to dumbp.

Usage: readp [OPTIONS]

Options:

            -c    add a comment to the output snapshot [false]
            -i    number the particles sequentially [don't number]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

readstoa

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:38
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/readstoa.C

Convert NEMO "stoa" format:

         (N, ndim, time,
          mass[i], i = 1,...,N,
          pos[i],  i = 1,...,N,
          vel[i],  i = 1,...,N)

data into a Starlab snapshot.

Usage: readstoa [OPTIONS]

Options:

         -i    number the particles sequentially [don't number]
         -w    write stoa insted of reading it

Written by Jun Makino.

Report bugs to starlab@sns.ias.edu.

rearrange_mf

   Starlab version 4.4.2
   program created on Jul 25 2005 at 18:00:06
   source file     $STARLAB_INSTALL_PATH/src/star/util/rearrange_mf.C
rearrange_mf.C: sort stars by stellar mass.
Input options:
   -B:    Binning option of the Lagrangian radii.
   -C: Cut off criterion, mass, luminosity, number of stars.
   -l:    Lower luminosity limit (see below).
   -n: Number of Lagrangian radii bins.
               Might for some choises be forced.
   -O:    Output option, what information should be studied.
   -o: Output the story with data written in root.
   -S:    Sort option, on what parameter must information be sorted.
Concerning the Lower luminosity limit option (-l, see above).
Currently this option is only applied to the sorting functions.
This means that the binning (using Lagrangian radii) is applied
on all cluster members.
whether or not this is realistic depends on the application.
Consequently, the luminosity cut-off does not affect the
Lagrangian radii.

red_star_cluster

Usage: ./red_star_cluster [-a #] [-B #] [-C #] [-c #] [-l #] [-N #] [-n] [-o] [-S #] [-s] [-v] Parameters usage for: red_stellar_system. options:

    -B:  Binning option of the Lagrangian radii.
         0= lineair, 1= predefined.
    -C:  Cut off creterium, mass, luminosity, number of stars.
         0= by number, 1= by mass, 2= by luminosity.
    -l:  Lower luminosity limit.
    -N:  Number of Lagrangian radii bins.
              Might for some choices be forced.
    -n:  Allow for n-squared operations.
    -o:  Output the story with data written in root.
    -S:  Sort option, on which parameter should be sorted.
         0= distance to com, 1= stellar mass.
    -s:  Scale to astrophysical units.
    -v:  give extended output

red_stellar_system

Usage: ./red_stellar_system [-B #] [-C #] [-c #] [-l #] [-N #] [-n] [-o] [-O #] [-S #] [-s] [-v] Parameters usage for: red_stellar_system. options:

    -B:  Binning option of the Lagrangian radii.
         0= lineair, 1= predefined.
    -C:  Cut off creterium, mass, luminosity, number of stars.
         0= by number, 1= by mass, 2= by luminosity.
    -l:  Lower luminosity limit.
    -N:  Number of Lagrangian radii bins.
              Might for some choices be forced.
    -n:  Allow for n-squared operations.
    -O:  Output option, what information should be studied.
         0= velocity anisotropy, 1= 3d dispersion velocity, 2= number of stars,
         3= number of stellar types, 4= total mass, 5= mass to light ratio, 
         6= mean mass, 7= number density, 8= mass density.
    -o:  Output the story with data written in root.
    -S:  Sort option, on which parameter should be sorted.
         0= distance to com, 1= stellar mass.
    -s:  Scale to astrophysical units.
    -v:  give extended output

reflect_velocities

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:34
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/reflect_velocities.C

Multiply all velocity components by -1, thereby effectively reversing the direction of time in the input N-body system(s).

Usage: reflect_velocities [OPTIONS] < input > output

Options:

                 -c     add a comment to the output snapshot [false]

Written by Piet Hut.

Reports bugs to starlab@sns.ias.edu.

renumber

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:58
   source file     $STARLAB_INSTALL_PATH/src/node/util/renumber.C

Renumber stars in a specific order

Usage: reunmber [OPTIONS]

Options:

          -c c      add a comment to the output snapshot [false]
          -I/i      start numbering number               [1]
          -M        Renumber the stars in order of mass
                    (highest mass=I/lowest mass=i)       [false]
          -N        Name the stars                       [false]
          -S/s      Single fixed number for all stars    [false]
                    Except if a star was already numbered

Written by Piet Hut and Simon Portegies Zwart.

Report bugs to starlab@sns.ias.edu. Options:

reset_counters

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:40
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/util/reset_counters.C

Clear all time/force counters in the input snapshot(s).

Usage: reset_counters [OPTIONS]

Options:

       -c    add a comment to the output snapshot [false]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

rmq

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:17
   source file     $STARLAB_INSTALL_PATH/src/node/util/rmq.C

Remove a quantity from the dyn (or star or hydro) story of the input snapshot(s) (FLAT TREE ONLY).

Usage rmq [OPTIONS]

Options:

           -c     add a comment to snapshot [false]
           -q     specify the quantity to remove [no default]
           -h     remove from hydro story [false]
           -s     remove from star story [false]

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

runtime_help

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:48
   source file     $STARLAB_INSTALL_PATH/src/std/runtime_help.C

Test Starlab run-time help function

Options:

     --help     print this message

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

scale

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:38
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/scale.C

(Re)scale an N-body system to specified M, Q (=T/U), and E. Note that if the required Q implies E > 0, no additional energy scaling is applied.

Note: Only top-level nodes are considered in scale_virial and scale_energy. Also, only internal positions and velocities are modified -- the root pos and vel are left untouched. Option "-s" is equivalent to "-M 1 -R 1 -Q 0.5". For systems in which stellar data have been specified, the physical mass, length, and time scales are already defined. Normally those will correspond to the cluster mass, length, and time units, so the "-s" option is probably most appropriate in most cases. This tool completes the connection between properties and physical units. In all cases, command-line parameters are specified in N-body units. As of 7/01, systems with embedded tidal or other external fields are also properly scaled (assuming that the Jacobi radius scales with the virial radius). Note that the virial radius is defined in terms of the

  • internal* potential energy only. For now, only allow the energy

to be specified if there are no external (non-tidal) fields.

Usage: scale [OPTIONS] < input > output

Options:

            -c    zero the center of mass position and velocity [no]
            -d    debug mode [off]
            -e    specify softening parameter [0]
            -E    specify total energy [don't scale]
            -m/M  specify total mass [don't scale]
            -q/Q  specify virial ratio [don't scale]
            -r/R  specify virial radius [don't scale]
            -s/S  scale to "standard" units (-m 1 -r 1 -q 0.5) [not set]

Scale writes the following quantities to the root log story if system_time = 0:

    initial_mass
    initial_total_energy
    initial_rvirial
    total_energy

If no arguments are provided and system_time = 0, scale will compute whatever quantities are needed to set these values.

Written by Steve McMillan.

Report bugs by starlab@sns.ias.edu.

scanbinary

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:01
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/tdyn/scanbinary.C

scanbinary: read and print binary statistics on a series of worbundles.

Options: -F input file [stdin]

                 *** file input is much faster (why?) ***

Program reads a collection of worldbundles, then repeatedly prompts for the name of a particle or node and prints out the interaction history of that object, in lines of the form

               time: top-level node name

End with exit, quit, or ^D.

Created by Steve, 9/2001.

scanbundle

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:05
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/tdyn/scanbundle.C

scanbundle: read and print statistics on a series of worldbundles.

Options: -F specify input file [stdin]

                 *** file input is much faster (why?) ***
           -v    more (too?) detailed output                  [no]

semi_to_period

Usage: ./semi_to_period [-a #] [-P #] [-M #] [-m #] [-q #]

Comprehensive tool for transforming orbital separations (-a in solar radii) to orbital periods (-P in days) given the primary (-M in solar masses) and secondary (-m in solar masses) masses or the mass ratio (-q).

If the semi-major axis (-a) is identified the code returns the orbital period (in days), if the orbital period (-P) is given the code returns the semi-major axis (in solar radii).

set_com

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:59
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/set_com.C

Modify positions and velocities to set the center-of-mass position and velocity. Write com_pos and com_vel to the root dyn story. If an external field has been specified, the velocity is taken to be in units of the circular orbit speed at the specified location. Otherwise, the velocity is taken as is. Positions and velocities are in physical units (parsecs and km/s, if relevant) if physical stellar parameters are known, and N-body units otherwise (or if -n is specified).

Usage: set_com [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]
            -n    force interpretation of r and v in N-body units [no]
            -r    specify center of mass position [(0,0,0)]
            -v    specify center of mass velocity [(0,0,0)]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

set_ignore_internal

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:35
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/set_ignore_internal.C

Add a flag specifying that internal interactions should be suppressed i.e. create a system of test particles.

Usage: set_ignore_internal [OPTIONS] < input > output

Options: None.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

set_radius

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:56:26
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/_dyn_/set_radius.C

setradius: set radiuses of specified particle(s) in input snapshot

         to specified values.

Usage: setradius -l l1 -m radius1 -l l2 -m radius2 ...

  or:  setradius -A radius

Options: -A set selected radius for all stars

           -l    specify label of next particle to modify [no default]
           -m    specify new radius for particle [no default]

setmass

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:01
   source file     $STARLAB_INSTALL_PATH/src/node/util/setmass.C

Set masses of specified particle(s) in input snapshot to specified values.

Usage: setmass -l l1 -m mass1 -l l2 -m mass2 ...

Options:

       -l    specify label of next particle to modify [no default]
       -m    specify new mass for particle [no default]

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

simple

Xlib: connection to ":0.0" refused by server Xlib: No protocol specified

Cannot connect to X server :0

snap2speck

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:50
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/snap2speck.C

Dump out an N-body system in a dumb format suitable for digestion by partiview.

Usage: snap2speck [OPTIONS] < input > output

Options:

                 -p    specify precision of output [6 sig. fig.]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

snap_to_image

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:57:30
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/hdyn/xutil/snap_to_image.C

snap_to_image: Construct images of a series of snapshots.

Options:

         -1           combine all frames in a simgle image          [yes]
         -a           produce a series of frames for animation       [no]
         -c           compress the image file(s) using gzip          [no]
         -C colormap  specify a colormap file name                   [no]
         -d           delete frames once the animation is made       [no]
         -f filename  specify root name of image files   ["-" --> stdout]
         -F format    specify image file format
                          (0 = PNG, 1 = SUN, 2 = GIF)                 [2]
         -g           write GIF files -- same as "-F 2"             [yes]
         -G           toggle forcing particles to grid (nicer single
                          frames, but jerkier movies)              [true]
         -H           toggle Herzsprung-Russel diagram or positional
                          plot                                 [position]
         -i index     specify (real) color index (0-1) for all stars
                                                     [use internal index]
         -l scale     specify width of field of view (+/- scale)      [3]
         -L loop      specify number of loops in animation           [10]
         -m           use mass to determine star color and/or size   [no]
         -n nmax      specify maximum number of images to produce   [Inf]
         -N nbody     color using a (small-N) colormap               [no]
         -o filename  same as -f (more standard name)
         -O option    specify how to choose the plot center           [0]
                          0:  as is (don't adjust)
                          1:  initial center of mass
                          2:  modified center of mass of each frame
                          (x,y,z):  specify coordinates (single argument)
                              - "(" can also be "[" or "{", but note
                                that all must be escaped or quoted,
                                since they have special meaning to
                                the shell
                              - separate components by commas or
                                spaces (again quoted or escaped, so
                                the vector is a single argument)
                              - the "{" version specifies an incremental
                                offest between successive frames.
         -p psize     specify (maximum) star radius/scale, in pixels
                      (psize < 0 ==> lower limit on pixel size = 0,
                      otherwise, limit = 1)
                                        [0 (single image), 1 (animation)]
         -P axis      specify projection axis                         [z]
         -q           toggle suppression of diagnostic output
                                                         [don't suppress]
         -r           use stellar radius to set point size           [no]
         -R           animate in reverse                             [no]
         -s nx ny     specify image size, in pixels                 [256]
         -S nskip     specify snaps to skip between images            [0]
         -t           test the color map [don't test]
         -T           specify precedence scheme for points in the image
                          (c = color, r = radius, z = depth)          [z]
         -x           specify right (log effective temparature) edge
                          of HRD (-H only)                            [3]
         -X           specify left (log effective temparature) edge
                          of HRD (-H only)                            [5]
         -y           specify minimum (log luminosity/Lsun) limit
                          of HRD (-H only)                           [-3]
         -Y           specify maximum (log luminosity/Lsun) limit
                          of HRD (-H only)                            [3]

In the case of GIF output, the command line responsible for creation of the image is encoded in the text segmant of the output file. Note that the default output format has been changed to GIF now that the LZW patent issues seem to have gone away.

Notes: 1. If animations are specified, an MNG or animated GIF file will

         be created.  The individual frames will be retained unless
         the "-d" option is set.
      2. If PNG output is requested and the PNG libraries are unavailable,
         then GIF output is produced instead.
      3. The "T" option determines which attribute is favored in
         displaying a particle in the image.  The default is depth,
         but we can also raise particles based on color or radius.

snapprune

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:21
   source file     $STARLAB_INSTALL_PATH/src/node/util/snapprune.C

Pass through only one out of every k snapshots from the input stream.

Usage snapprune [OPTIONS]

Options:

          -c    add comment to snapshots [false]
           -k    specify k [2]
           -s    specify random seed [random from system clock]
           -x    specify output snapshot after which to exit [at end]

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

split_particles

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:55:18
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/init/split_particles.C

Split the specified particles in input snapshot into binaries having specified orbital properties.

Naming convention: particle X is split into components X1 and X2. Splitting may be applied recursively to produce multiple systems, e.g.

      split_particles -l 1 -a 0.1 -e 0 -l 1b -a 0.01 -e 0.5

creates a triple system (11, (121, 122)).

Unlike mkscat, structure is never implicit, but must be created from the top down.

Usage: split_particles [-s #] -l l1 -a a1 -e e1 -E E1 -q q1 -l l2 (etc.)

Options:

           -s    specify random seed [random from system clock]
           -i    use "a" and "b" instead of "1" and "2" [false]
           -l    specify label of next particle to split [no default]
           -a    specify semi-major axis of current binary [0.1]
           -e    specify eccentricity of current binary [random]
           -E    specify energy of current binary [0]
           -q    specify mass ratio (secondary[b]/primary[a]) of
                 the current binary [1]

The "-E" option takes precedence over the "-a" option.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

sstar_stats

Usage: ./sstar_stats [-b] [-n] [-e] [-o] [-s] [-t]

starev

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:59
   source file     $STARLAB_INSTALL_PATH/src/star/sstar/evolve/starev.C

starev: evolve a single star.

       creates a single star and evolves it in time.

Options: -c comment to put in the starbase log structure.

           -M    initial mass of the star [in solar units].
           -n    number of output timesteps (timesteps are taken
                 with constant time intervals) 
           -R    Dynamical size scaling for the star
                 [in units of the virial radius].
           -S    Random seed.
           -s    Initial stellar type [default is main_sequence].
           -T or -t end time of the stellar evolution [in Million year].

Latest version (SPZ:1.0) February 1993.

starplot

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:47
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/starplot.C

Plot input N-body system(s) on a non-graphics screen.

Usage: starplot [OPTIONS] < input > output

Options:

                 -a    view axis [z]
            -c    clear the screen before each new plot [scroll screen]
            -l    specify plot limits [get from first snapshot]
            -n    number of lines to use [100]

Written by Piet Hut.

Report bugs to starlab@sns.ias.edu.

std_rename_nodes

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:51:07
   source file     $STARLAB_INSTALL_PATH/src/node/util/std_rename_nodes.C

Rename all CM nodes bottom up to reflect the standard Starlab naming scheme (a,b).

Options: None.

Written by Steve McMillan.

Report bugs to starlab@sns.ias.edu.

story

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:49:54
   source file     $STARLAB_INSTALL_PATH/src/std/story.C

Test Starlab story manipulation functions.

Options: None.

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

sys_stats

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:52:54
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/sys_stats.C

Print out various diagnostic statistics on the input system. These include:

       system time, number, mass, mass distribution;
       relaxation time;
       system energy;
       core parameters;
       lagrangian radii for quartiles [default], for ten-percentiles,
       and for "special" choice of Lagrangian masses, currently 0.005,
       0.01, 0.02, 0.05, 0.1, 0.25, 0.5, 0.75, 0.9;
       mass distribution by lagrangian zone;
       anisotropy by lagrangian zone;
       binary parameters.

In addition, Lagrangian radii are written to the root dyn story. If N^2_ops is selected, then core parameters and particle densities are also written to the root and particle dyn stories.

If N^2_ops is selected, Lagrangian radii are computed relative to the density center, as are binary radial coordinates. Otherwise, the modified center of mass (center of mass with outliers excluded) is used.

If GRAPE is available, the hdyn version of this function will compute the energies even if the "-n" flag is set false.

Usage: sys_stats [OPTIONS] < input > output

Options:

            -b    specify level of binary statistics [2]
                  0: none
                  1 (or no arg): short binary output
                  2: full binary output
            -e    recalculate the total energy (even if -n is no) [yes]
            -l    specify percentile choice [2]
                  0: quartiles (1-3)
                  1: 10-percentiles (10-90)
                  2 (or no arg): nonlinear Lagrangian masses
                  (0.5, 1, 2, 5, 10, 25, 50, 75, 90%)
            -n    perform/don't perform actions requiring O(N^2)
                  operations (e.g. computation of energy and core
                  radius; see -e) [no]
            -o    pipe system to cout [no]

Written by the Starlab development group.

Report bugs to starlab@sns.ias.edu.

tdyn_io

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:09
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/tdyn/tdyn_io.C

tdyn_io: Starlab pdyn and tdyn class I/O functions.

         Define scan_star_story, scan_dyn_story and print_dyn_story
         for the pdyn and tdyn classes.

Options: none

to_com

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:53:03
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/util/to_com.C

Bring all positions and velocities to center-of-mass frame. Forces all daughter nodes to have a CM at rest at the origin, and zeroes the root pos and vel. Uses compute_com and writes to the root dyn story. Does not correct the virial radius in the case of a tidal field.

Usage: to_com [OPTIONS] < input > output

Options:

                 -c    add a comment to the output snapshot [false]

Written by Piet Hut and Steve McMillan.

Report bugs to starlab@sns.ias.edu.

worldlines

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:14
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/tdyn/worldlines.C
  worldlines:  Create output snapshots from input worldbundles.
               Produce a snapshot or diagnostic starting at time t,
               optionally with time step dt.
  options:     -d    output interval (negative integer ==> 2^[-d])   [0]
               -D    diagnostic: no snapshsot                       [no]
               -F    specify input file                          [stdin]
                     *** file input is much faster (why?) ***
               -t    initial output time                 [start of data]
               -v    verbosity level                                 [1]
  Snapshots are generated as soon as the data are available
  -- i.e. we don't read in the entire dataset first.
  Only one snapshot is produced if the time step is 0 (default).
  Note: Output is in pdyn format, which currently cannot be used
  directly as input to kira.
  Created by  SPZ at MIT in Dec 2000
  Modified by SLWM at DU in Dec 2003

xhrdplot

Usage: ./xhrdplot [-a #] [-d #] [-D #] [-e] [-f] [-l #] [-o] [-p #] [-P #] [-r] [-s #] [-t]

ximage

ximage: [nx] [ny] [-c colormapfile] [-o xo [yo]] [-s nx [ny]]

xreal

xreal data type is real

xstarplot

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:41
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/xutil/xstarplot.C

xstarplot: plot an N-body system in an X environment.

Options:

           -a    specify viewing axis [3 = z]
           -c    don't display in the cluster CM frame [work in CM frame]
           -d    dimensionality of plot [2]
           -D    delay time between plots, in ms [0]
           -e    color by energy [no]
           -f    solid-color stars [true]
           -l    specify limits [take from initial data]
           -p    pipe cin to cout [no]
           -P    specify point size (relative to axis/30) [1]
           -r    "reverse video" (black background) [true]
           -s    window size, relative to "standard" 400 pixels [1]
           -t    show tree structure [no]

xstarplot22

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:58:22
   source file     $STARLAB_INSTALL_PATH/src/node/dyn/tdyn/xstarplot22.C

xstarplot2: plot an N-body system in an X environment.

            This version uses the NEW interpolation routines.

Options: -a specify viewing axis [1/2/3 = x/y/z]

           -b    allow backward steps [no]
           -d    dimensionality of plot [2]
           -D    time interval between frames [0.015625 = 1/64]
           -e    color by energy [no]
           -E    toggle display of stars flagged as escapers [true]
           -f    solid-color stars [true]
           -F    input file [run.out]
           -k    maximum size of CM indicator is k/N [2]
           -l    show links and nodes [no]
           -L    specify limits [take from initial data]
           -m    highlight multiple systems [no]
           -o    pipe cin to cout [no]
           -p    point scale mode [0]
           -P    specify point size (relative to axis/30) [1, min = 0.5]
           -r    "reverse video" (black background) [true]
           -s    window size, relative to "standard" 400 pixels [1]
           -t    show tree structure [no]
           -u    highlight unperturbed systems [no]
           -v    verbose mode (worldbundle details) [0]


Note: There appears to be a problem with the X interface on some

      systems that causes xstarplot to crash if too many run-time
      user commands are received in a short space of time (e.g.
      try holding down the ">" key in 3D mode).  Presently, the
      only fix is to avoid excessive keystrokes in the graphics
      window at run time...

xtest

   Starlab version 4.4.2
   program created on Jul 25 2005 at 17:50:02
   source file     $STARLAB_INSTALL_PATH/src/std/xtest.C