HTML automatically generated with rman
Table of Contents
orbstat - tabulate some orbital statistics
orbstat in=orbit
[parameter=value]
orbstat will tabulate some statistical properties
of a series or orbits, as given in an input file. Currently time, energy
(or first integral of motion), the maximum extent in X, Y, VX and VY are
given, as well as the mean and dispersion in the angular momentum for
(XY) planar orbits are given.
The following parameters are recognized
in any order if the keyword is also given:
- in=in-file
- Input file, which
must be an orbit(5NEMO)
. [No default].
- ellipse=t|f
- Check how well the orbit
is ellipsoidal. If set to true, it will print T, E, r2_min, r2_max, v2_min,
v2_max, r2_mean, r2_sigma, v2_mean, v2_sigma, Lmin/Lmax. Default: f.
- relative=t|f
- Should energy statistics be done relative to the "known" value by subtracting
it? This is to improve the statistics if the integration is very accurate.
Default: f.
- tab=t|f
- Output a table with quantities Default: f.
The ellipse=t
option only makes sense for closed orbits, since the maximum extent in
position and velocity is used for scaling factors to check the ellipsoidability.
% mkorbit orb1 1 0 0 0 0.1 1 potname=plummer ’potpars=0,1/sqrt(2),1/sqrt(2)’
% orbint orb1 - 100 1 mode=euler | orbstat -
# T E x_max y_max u_max v_max j_mean j_sigma e_mean
e_sigma [abs]
100 -0.0723503 1 4.04463 0 0.1 0.184855 0.0101289 0.163057 0.0239529
% orbint orb1 - 200 0.5 mode=euler | orbstat -
100 -0.0723503 1 2.34815 0 0.1 0.137628 0.00395845 0.0401559 0.00887043
% orbint orb1 - 100 1 mode=rk4 | orbstat -
100 -0.0723503 1.91846 0.745393 0.872376 0.1 0.0998526 7.33814e-05 -0.0746802 0.000728993
% orbint orb1 - 200 1/2 mode=rk4 | orbstat -
100 -0.0723503 1.93045 0.710908 0.871285 0.1 0.0999952 2.51481e-06 -0.0724567 3.22187e-05
orblist(1NEMO)
Peter Teuben
~/src/orbit/misc orbstat.c
23-mar-95 V1.0 Created PJT
21-feb-03 V1.1 check how well ellipsoidal PJT
25-apr-05 V2.0 add relative= and tab= PJT
Table of Contents