readme_020320

Files with particles' position & velocity

 cr123_127.tar    contains 1000 snapshots, 0 - 200 Chariklo rotation periods
                  snapshots are stored as unformatted data 
		  =269 MB


examples:

unformatted:
  cr123_127_elon02_eta01_s150_n10_soft_possave500.dat     after 100 periods
  cr123_127_elon02_eta01_s150_n10_soft_possave1000.dat    after 200 periods

same examples as formatted ascii files:
   cr123_127_elon02_eta01_s150_n10_soft_possave500_ascii.dat 
   cr123_127_elon02_eta01_s150_n10_soft_possave1000_ascii.dat


format of the files:
(note that the type of variables has to be specified before reading)

a,b,c              scalars, double
time,omeb          scalars, double
n                  scalar, 4-byte integer
x                  float array with n elements
y                   -"-
z                   -"-
vx                  -"-
vy                  -"-
vz                  -"-
EJ                  -"-

x,y,z,vx,vy,vz  are positions and velocities in inertial system
EJ jacobi energy/mass unit

to convert to rotating one (intially Chariklo major axis = x-axis)

 sino=sin(time*omeb)
 coso=cos(time*omeb)
 xrot= x*coso  +y*sino
 yrot=-x*sino  +y*coso
     

Simulation parameters:

a=1.  GM=1.
b=0.98
c=1.
chariklo Ome=0.7600627 --> period=8.26667

refrad=00.993198
f=-0.010036

in these units corotation r=1.19904    Lz=1.09274
               2:3        r=1.57294    Lz 1.25266