HTML automatically generated with rman
Table of Contents
snapmradii - print lagrangian mass radii in a snapshot
snapscale
in=in_file [parameter=value] ...
snapmradii prints a table of the
fractional lagrangian mass-radii of a snapshot. The snapshot will be sorted
internally by increasing radius (which can be changed by using sort=),
there is no need to use snapsort(1NEMO)
preceding this call.
The
following parameters are recognized:
- in=in_file
- input file, in snapshot
format. Default: none.
- fraction=f1,f2,...
- Fractional mass radii. They must be
sorted and between 0 and 1 (inclusive). [Default: 0.1:0.9:0.1]
- tab=t|f
- Output
in the form of a table of mass fraction vs. radius? The normal output mode
is one snapshot per output row, with time, and the mass fractions [Default:
f]
- log=t|f
- Output the radii in logarithmic units [Default: f]
- sort=sort_var
- Sorting variable used to sort the particles. Any bodytrans(3NEMO)
expression
can be used. Default: r. Default
The following example shows the evolution
of 10% lagrangian mass radius and outputs the result in a graph
snapmradii snap.dat 0.05 |\
tabplot - xlab=time ylab="10% Mass radius"
The following example shows the cumulative mass distribution function for
the third snapshot in a large simulation, and plots it:
snaptrim snap.dat - times=#3 | snapmradii - 0.01:0.99:01 |\
tabplot - xlab="Mass fraction" ylab=Radius
Doesn’t handle snapshots with zero mass very well.
Doesn’t deal with snapshots
where there is no stars in one of the fraction bins requested. Since linear
interpolation is used, some crazy larger radius is reported.
snapcenter(1NEMO)
,
snapstat(1NEMO)
, snapprint(1NEMO)
W.L.Sweatman - (1993) MNRAS 261, 497.
Peter
Teuben
~/src/nbody/reduc snapmradii.c
12-May-92 V1.0: created PJT
10-aug-95 V1.2: added tab= PJT
15-aug-96 V1.3: fixed bug which assumed total mass =1 PJT
27-jul-05 V1.5: add sort= PJT
1-apr-21 V1.6: handle massless snapshots for Tjeerd PJT
Table of Contents