NBodyLab Adaptation of NBODY4 Running with a GRAPE-6a

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: 13360 Program version: 050823

Running from TIME = 0 to TIME = 100.0 TIME Myr <R> pc N __ 0.0 0.0 3.13 1500 __ 2.0 0.9 3.07 1500 __ 4.0 1.8 3.09 1500 __ 6.0 2.7 2.99 1500 __ 8.0 3.6 3.01 1500 __ 10.0 4.5 2.97 1500 __ 12.0 5.4 2.89 1500 __ 14.0 6.3 2.74 1500 __ 16.0 7.3 2.64 1500 __ 18.0 8.2 2.54 1500 __ 20.0 9.1 2.40 1500 __ 22.0 10.0 2.24 1500 __ 24.0 10.9 2.09 1500 __ 26.0 11.8 1.84 1500 __ 28.0 12.7 1.71 1500 __ 30.0 13.6 1.66 1498 __ 32.0 14.5 1.63 1498 __ 34.0 15.4 1.64 1498 __ 36.0 16.3 1.80 1498 __ 38.0 17.2 1.93 1498 __ 40.0 18.1 2.05 1498 __ 42.0 19.0 2.16 1498 __ 44.0 19.9 2.27 1496 __ 46.0 20.9 2.40 1496 __ 48.0 21.8 2.46 1495 __ 50.0 22.7 2.48 1490 __ 52.0 23.6 2.46 1488 __ 54.0 24.5 2.50 1487 __ 56.0 25.4 2.43 1484 __ 58.0 26.3 2.41 1482 __ 60.0 27.2 2.24 1479 __ 62.0 28.1 2.33 1478 __ 64.0 29.0 2.30 1476 __ 66.0 29.9 2.17 1475 __ 68.0 30.8 2.05 1474 __ 70.0 31.7 2.04 1473 __ 72.0 32.6 2.06 1471 __ 74.0 33.5 2.12 1468 __ 76.0 34.5 2.16 1466 __ 78.0 35.4 2.20 1465 __ 80.0 36.3 2.21 1461 __ 82.0 37.2 2.31 1460 __ 84.0 38.1 2.36 1459 __ 86.0 39.0 2.36 1456 __ 88.0 39.9 2.35 1453 __ 90.0 40.8 2.42 1447 __ 92.0 41.7 2.49 1445 __ 94.0 42.6 2.47 1443 __ 96.0 43.5 2.48 1441 __ 98.0 44.4 2.47 1438 __ 100.0 45.3 2.49 1432 ... INTEGRATION BINARY SIZE NSTEPS DE/E E EB/E MEMBERS AMIN 0.00E+00 0.00E+00 -0.203295 0.0000 0 0 0.0E+00 4.75E+05 3.77E-08 -0.203295 0.0000 0 0 0.0E+00 9.70E+05 -1.20E-07 -0.203295 0.0000 0 0 0.0E+00 1.44E+06 3.25E-07 -0.203295 0.0000 0 0 0.0E+00 1.90E+06 -1.45E-07 -0.203295 0.0183 1006 1007 4.4E-03 2.32E+06 2.03E-07 -0.203295 0.0182 1006 1007 4.4E-03 2.73E+06 7.04E-07 -0.203295 0.0000 0 0 0.0E+00 3.20E+06 4.75E-07 -0.203295 0.0000 0 0 0.0E+00 3.63E+06 4.81E-07 -0.203295 0.0230 1006 1007 3.5E-03 4.04E+06 5.18E-07 -0.203295 0.0793 1003 1013 1.5E-03 4.44E+06 6.23E-07 -0.203294 0.0916 1003 1013 1.6E-03 4.80E+06 -9.16E-07 -0.203295 0.0854 1003 1013 1.4E-03 5.15E+06 5.17E-07 -0.203294 0.1312 1001 1006 2.2E-03 5.53E+06 -1.56E-06 -0.203295 0.3080 1003 1013 3.9E-04 5.89E+06 8.88E-07 -0.203295 0.3292 1003 1013 3.9E-04 6.24E+06 3.14E-06 -0.203583 0.3094 1003 1013 3.9E-04 6.57E+06 8.45E-07 -0.203582 0.2322 1003 1013 3.9E-04 6.93E+06 3.76E-06 -0.203582 0.3496 1003 1013 3.9E-04 7.28E+06 4.80E-06 -0.203581 0.5692 1003 1013 3.9E-04 7.56E+06 -1.68E-06 -0.203581 0.7158 1003 1013 3.9E-04 7.78E+06 1.25E-06 -0.203581 0.8059 1003 1013 3.9E-04 8.00E+06 6.14E-07 -0.203581 0.7460 1003 1013 3.9E-04 8.22E+06 -2.94E-06 -0.205279 1.6463 5 2 4.1E-04 8.41E+06 -6.83E-07 -0.205279 1.6449 5 2 4.1E-04 8.59E+06 -7.42E-06 -0.210285 1.5535 5 2 4.1E-04 8.75E+06 2.61E-06 -0.211498 1.5609 5 2 4.1E-04 8.90E+06 -4.54E-07 -0.211728 1.5514 5 2 4.1E-04 9.05E+06 -9.74E-07 -0.212002 18.3880 1001 1007 1.7E-05 9.20E+06 -8.44E-05 -0.214570 17.2606 1001 1007 1.8E-05 9.35E+06 -5.44E-07 -0.214918 16.9013 1001 1007 1.8E-05 9.50E+06 1.20E-06 -0.225282 14.8529 1001 1007 1.8E-05 9.64E+06 1.63E-05 -0.226408 0.6230 5 2 4.1E-04 9.77E+06 -1.20E-08 -1.645002 5.9547 5 2 4.1E-04 9.89E+06 -2.83E-08 -1.704954 1.3273 5 2 4.1E-04 1.00E+07 -1.05E-07 -1.704962 1.3120 5 2 4.1E-04 1.01E+07 1.24E-07 -1.705000 1.2967 5 2 4.1E-04 1.03E+07 -4.56E-07 -1.705105 1.2812 5 2 4.1E-04 1.04E+07 -1.53E-08 -1.705475 1.2788 5 2 4.1E-04 1.05E+07 -4.08E-07 -1.705589 0.2235 1003 1013 3.9E-04 1.06E+07 4.18E-07 -1.798966 0.2324 1003 1013 3.9E-04 1.08E+07 2.08E-08 -1.857855 0.8576 1003 1013 3.9E-04 1.09E+07 -6.78E-07 -1.857856 0.8635 1003 1013 3.9E-04 1.10E+07 2.82E-07 -1.857905 0.8513 1003 1013 3.9E-04 1.11E+07 8.47E-07 -1.892107 0.4686 1003 1013 4.0E-04 1.12E+07 2.04E-07 -1.892208 0.4688 1003 1013 4.0E-04 1.13E+07 4.09E-07 -1.892473 0.4601 1003 1013 4.0E-04 1.14E+07 -5.45E-05 -1.892580 0.0000 0 0 0.0E+00 1.16E+07 6.20E-07 -1.892584 0.0000 0 0 0.0E+00 1.17E+07 3.03E-07 -1.892606 0.0000 0 0 0.0E+00 1.18E+07 4.83E-07 -1.957105 0.0000 0 0 0.0E+00 1.19E+07 3.79E-08 -2.033176 0.0000 0 0 0.0E+00

Elapsed real time (secs): 133

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

File ESC, escaper removal (option KZ23 > 1)

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

Density and velocity as a function of average radii in increasing shells.
Plotted as log10(density) versus log10(distance), and average velocity versus central distance.