NBodyLab Adaptation of Dr. Sverre Aarseth's NBODY4, running with a GRAPE-6a

See overview (sample runs, timings, etc) and explanation of input and output quantities.

NBODY4 input file
  
KSTART TCOMP GPID
N NFIX NCRIT NRAND NRUN
ETA DTADJ DELTAT TCRIT QE RBAR ZMBAR
KZ1   KZ2   KZ3 KZ4   KZ5   KZ6   KZ7   KZ8   KZ9   KZ10
KZ11 KZ12 KZ13 KZ14 KZ15 KZ16 KZ17 KZ18 KZ19 KZ20
KZ21 KZ22 KZ23 KZ24 KZ25 KZ26 KZ27 KZ28 KZ29 KZ30
KZ31 KZ32 KZ33 KZ34 KZ35 KZ36 KZ37 KZ38 KZ39 KZ40
DTMIN RMIN ETAU ECLOSE GMIN GMAX
ALPHA BODY1 BODYN NBIN0 ZMET EPOCH0 DTPLOT
Additional lines depend on options - see explanations.

Time scale for diagnostic plots   
Axes length for N-body plots, or 'autoscale'



Run id: 14291        Program version: 050823

Running from TIME = 0 to TIME = 20.0 ... INTEGRATION BINARY SIZE TIME Myr <R> pc N NSTEPS DE/E E EB/E MEMBERS AMIN __ 0.0 0.0 0.00 1001 0.00E+00 0.00E+00 0.007590 0.0000 0 0 0.0E+00 __ 1.0 0.0 0.00 1001 9.21E+04 -4.54E-08 0.007590 0.0000 0 0 0.0E+00 __ 2.0 0.0 0.00 1001 1.77E+05 -1.25E-07 0.007590 0.0000 0 0 0.0E+00 __ 3.0 0.0 0.00 1001 2.58E+05 1.43E-07 0.007590 0.0000 0 0 0.0E+00 __ 4.0 0.0 0.00 1001 3.37E+05 7.27E-08 0.007590 0.0000 0 0 0.0E+00 __ 5.0 0.0 0.00 1001 4.12E+05 5.67E-09 0.007590 0.0000 0 0 0.0E+00 __ 6.0 0.0 0.00 1001 4.85E+05 -1.57E-07 0.007590 0.0000 0 0 0.0E+00 __ 7.0 0.0 0.00 1001 5.53E+05 6.86E-08 0.007590 0.0000 0 0 0.0E+00 __ 8.0 0.0 0.00 1001 6.19E+05 1.14E-07 0.007590 0.0000 0 0 0.0E+00 __ 9.0 0.0 0.00 1001 6.82E+05 3.20E-08 0.007590 0.0000 0 0 0.0E+00 __ 10.0 0.0 0.00 1001 7.44E+05 -1.27E-07 0.007590 0.0000 0 0 0.0E+00 __ 11.0 0.0 0.00 1001 8.07E+05 -5.93E-08 0.007590 0.0000 0 0 0.0E+00 __ 12.0 0.0 0.00 1001 8.68E+05 5.69E-09 0.007590 0.0000 0 0 0.0E+00 __ 13.0 0.0 0.00 1001 9.26E+05 -2.04E-09 0.007590 0.0000 0 0 0.0E+00 __ 14.0 0.0 0.00 1001 9.83E+05 1.76E-07 0.007590 0.0000 0 0 0.0E+00 __ 15.0 0.0 0.00 1001 1.04E+06 -1.27E-07 0.007590 0.0000 0 0 0.0E+00 __ 16.0 0.0 0.00 1001 1.10E+06 2.13E-08 0.007590 0.0000 0 0 0.0E+00 __ 17.0 0.0 0.00 1001 1.16E+06 -4.19E-08 0.007590 0.0000 0 0 0.0E+00 __ 18.0 0.0 0.00 1001 1.22E+06 1.99E-08 0.007590 0.0000 0 0 0.0E+00 __ 19.0 0.0 0.00 1001 1.28E+06 -9.52E-09 0.007590 0.0000 0 0 0.0E+00 __ 20.0 0.0 0.00 1001 1.33E+06 2.52E-08 0.007590 0.0000 0 0 0.0E+00 __ 21.0 0.0 0.00 1001 1.39E+06 2.42E-07 0.007590 0.0000 0 0 0.0E+00

Elapsed real time (secs): 12


Carefully review the parameter assignments and results in the full NBODY4 run output
NBODY4 run warnings and errors in stderr




Eccentricity dispersion (rms value) DISP and total energy EDISK of the disk

File used for above diagnostic plots
Contents: TTOT ERROR DETOT N NSTEP TSCALE*TTOT RC RSCALE RC*RBAR RTIDE EB/E AMIN DISP EDISK DMIN
























File (ASCII) of N-body evolution.
Rows are m x y z vx vy vz. Used to generate above N-body plots.


Click here for 3D animation of N-body evolution, using Java applet