Table of Contents

Name

snaprect - diagonalize moment-of-inertia of snapshot

Synopsis

snaprect in=snap_file out=snap_file [parameter=value] ...

Description

snaprect transforms the input N-body snapshot to a coordinate system which diagonalizes the moment-of-inertia tensor of a specified subset of particles (see keywords weight, rcut and shell below). The centroid of the subset is translated to the origin. As with snapkinem(1NEMO) , an arbitrary weighting may be specified for each particle.

Iterative determinations of the axis ratio, using the moments of inertia, have been implemented on an experimental basis.

Parameters

The following parameters are recognized.
in=snap_file
Input data is read from snap_file, which must be in snapshot(5NEMO) format. Default: required.
out=snap_file
Output data is written to snap_file, in snapshot(5NEMO) format. Default: required.
weight=body_weight
Expression used to compute the weight of each body. See bodytrans(1NEMO) for the form of body_weight. Default: 1.0.
rcut=cutoff_radius
If cutoff_radius is greater than zero, only particles within cutoff_radius of the body_weighted centroid will be used to compute the moment-of-inertia. Default: 0.0.
times=time_range
Only frames within time_range will be analyzed. Default: all.
shell=r_min,r_max
Two radii are expected here. They define an ellipsoidal shell within which the (modified) inertia tensor is computed and iterated. Default: all particles.
iters=max_iters
Maximum number of iterations allowed in computing the inertia tensor within a changing volume/shell. Default: 1.
tol=min_tol
Minimum tolerance required to stop iterations. [Default: 0.0001]
modified=t|f
Use modified inertia tensor? If set to true, the modified (normalized w.r.t. the squared ellipsoidal radius from a previous iteraton) inertia tensor is used instead of the regular one. Note that the body_weight is also applied. It only makes sense to use this option if max_iters > 1. [Default: f].
test=snap_file
If a filename is given, it is expected to be a snapshot. The first snapshot in this file is used to rectify all selected snapshots from the inputfile. Default: none.

Bugs

During iterations, the printed output contains momentary eigenvectors and angles, not the cumulative ones.

See Also

snapkinem(1NEMO) , snapinert(1NEMO) , snapplot(1NEMO) , snapshot(5NEMO)

Aguilar & Meritt (1990) - ApJ 354, 33.    Using a fraction of the most bound particles

Katz (1991) ApJ 368, 325. Volume iterative tensor.

Dubinski & Carlberg (1991) ApJ 378, 496. Volume iterative modified tensor.

Kuijken & Dubinksi (1994) MN 269, 13. Shell iterative modified tensor.

Author

Joshua E. Barnes, Peter Teuben

History


5-mar-89    original version?    JEB
4-jul-94    finally implemented test=, also output theta/phi of p.axes    PJT
24-jul-94    V1.4 experimental iters=, shell=    PJT
15-mar-02    V1.4d addded log= keyword    WD