Currently a table is produced with the
following entries:

time, -2T/W, T+W, T, W, W_1, W_2, Masswhere

Wacc = sum_(m r.f)(sum over all masses (

*W_1* and W_2 in the table are the other two
potentials not selected as the prime potential, the original order of the
potentials is **acc, phi, exact**. The choosen one comes first, the others ones
follow in the order left from this array.

**in=***InFile*- Input file name (snapshot). No default.
**times=***t1:t2,t3:t4,...*- Times to select snapshot. [Default: all].
**wmode=***W_mode*- Mode how to compute the prime potential: choice
is between
**acc**(forces from the snapshot, hence Clausius virial),**phi**(potentials from the snapshot), and**exact**(computed exact newtonian potential m/r [*expensive!*]). The remaining unchoosen ones are listed as W_1 resp. W_2 in the output table. [Default:**acc**] **newton=t|f**- Do an N^2 newtonian calculation too? If true,
an N^2 algorithm will be used and [Default:
**f**] **eps=**- Standard gravitational
softening for
**mode=exact**. [Default:**0.05**]

Marmo, C.; Secco, L. *Clausius’ virial vs. total potential
energy in the dynamics of a two-component system* (**2003NewA....8..629M**)

05-nov-91V0.0 CreatedSMF 27-Mar-92V0.1 added clausius - changed columns and namesPJT 3-apr-92V0.2 added wmode= keywordPJT 13-mar-97V0.4 added total mass to the outputPJT 30-jul-97V0.5 added eps=PJT