Table of Contents

name

lagrange_radii - find the Lagrange radii of a stellar system

synopsis

lagrange_radii [parameter=value]

requires NEMO version 3.0.13 or higher, see http://www.astro.umd.edu/nemo.

description

lagrange_radii finds the radii of spheres centred at the origin x=(0,0,0) that contain a given fraction of the total mass. The code is significantly faster than a full sort of all body radii.

lagrange_radii may be used to create a stopfile to trigger the stopping of a running N-body simulations that is being analysed by lagrange_radii.

parameters

in=file
Input snapshots are read from file, which must be in NEMO snapshot format. If file=-, input is expected from stdin, useful for piping data from another program. Default: required
masses=list
The lagrange radii are computed for every fractional mass contained in the comma separated list (up to 100). Default: 0.01,0.03,0.1,0.3,0.5,0.7,0.9,0.97,0.99
times=time range(s)
Range(s) of simulation times, for which snapshots shall be processed. Default: all.
out=file
If given, snapshots are written in NEMO snapshot format to file with same body data as the input, unless give specifies otherwise. No default.
give=IOstring
If given, only body data indicated by IOstring (and present in the input) are written out with the snapshot output. No default.
step=val
If snapshots are written out, the outputs are made every val simulation time units. If val=0, every snapshot processed is written out. Default: 0.
tabfile=file
A table of the simulation time and Lagrange radii (in the order the fractional masses are given) is written to file, appending an existing file, if applicable. If not given, Lagrange radii are written to stdout. No default.
stopfile=file
If given, the a stopfile will be created, IF the condition defined with the next options is satisfied. See also gyrfalcON(1falcON) . No default.
stopindex=num
If stopfile is given, R_stop is the numth Lagrange radius. Default: 0.
stopvalue=val
If val>0, the stop condition is R_stop < val.

If val<0, the stop condition is R_stop > |val|.

No default.

stoprelative=t|f
If true, use R_stop/R_stop(t_initial) instead of R_stop in the stop condition. Default: f.
stopafter=val
If given, do not stop before simulation time val. No default.
stopdelay=val
Create the stopfile val simulation time units after the stop condition was first satisfied. Default: 0.

remark

Useful for isolated stellar systems only.

see also

density_centre(1falcON) , gyrfalcON(1falcON)

files


falcON/doc/user_guide.pdf                         User Guide for falcON

author


Walter Dehnen                              walter.dehnen@astro.le.ac.uk
update history of this man page
18-May-2004    first version of manual  WD


Table of Contents