% mkplummer out=p100 nbody=100

% hackcode1 in=p100 out=p100.dat tstop=1 > p100.log

% snapprint in=p100.dat header=t> p100.posvel
x y z vx vy vz 

% head p100.posvel
100 0
-1.29186 -0.664918 -0.330791 -0.0899679 -0.138913 0.165424 
-1.16527 0.361872 -0.638704 -0.510376 0.179608 0.0785804 
-1.21846 1.53669 0.622052 0.319839 0.0801127 0.0406664 
0.217503 -0.570365 0.0600823 -0.465651 -0.439758 0.6246 
...

% wc p100.posvel 
     505    3010   29151 p100.posvel

% tabtos in=p100.posvel out=q100.dat header=nbody,time block1=pos,vel   
[reading 100 bodies at time 0.000000]
[reading 100 bodies at time 0.250000]
[reading 100 bodies at time 0.500000]
[reading 100 bodies at time 0.750000]
[reading 100 bodies at time 1.000000]

% snapplot q100.dat nxy=3,3

% snapplot q100.dat nxy=3,3 yapp=/vps

% lpr -Pastro2 pgplot.ps

% snapdiagplot p100.dat 
33 diagnostic frames read
### Warning [snapdiagplot]: Autoscaling time. MinMax=-0.05 1.05
Worst fractional energy loss dE/E = (E_t-E_0)/E_0 = -0.00659881 at T = 0.75