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

Running from TIME TIME Myr <R> pc __ 0.0 0.0 __ 2.0 1.2 __ 4.0 2.4 __ 6.0 3.6 __ 8.0 4.9 __ 10.0 6.1 __ 12.0 7.3 __ 14.0 8.5 __ 16.0 9.7 __ 18.0 10.9 __ 20.0 12.1 __ 22.0 13.3 __ 24.0 14.6 __ 26.0 15.8 __ 28.0 17.0 __ 30.0 18.2 __ 32.0 19.4 __ 34.0 20.6 __ 36.0 21.8 __ 38.0 23.0 __ 40.0 24.3 __ 42.0 25.5 __ 44.0 26.7 __ 46.0 27.9 __ 48.0 29.1 __ 50.0 30.3 __ 52.0 31.5 __ 54.0 32.7 __ 56.0 34.0 __ 58.0 35.2 __ 60.0 36.4 __ 62.0 37.6 __ 64.0 38.8 __ 66.0 40.0 __ 68.0 41.2 __ 70.0 42.4 __ 72.0 43.7 __ 74.0 44.9 __ 76.0 46.1 __ 78.0 47.3 __ 80.0 48.5 __ 82.0 49.7 __ 84.0 50.9 __ 86.0 52.1 __ 88.0 53.4 __ 90.0 54.6 __ 92.0 55.8 __ 94.0 57.0 __ 96.0 58.2 __ 98.0 59.4 __ 100.0 60.6 = 0 to TIME = 100.0 ... INTEGRATION BINARY SIZE N NSTEPS DE/E E EB/E MEMBERS AMIN 0.78 1000 0.00E+00 0.00E+00 -0.249563 0.0000 0 0 0.0E+00 0.83 1000 3.12E+05 -1.35E-07 -0.249564 0.0000 0 0 0.0E+00 0.76 1000 6.21E+05 6.74E-08 -0.249564 0.0000 0 0 0.0E+00 0.74 1000 9.46E+05 6.46E-08 -0.249563 0.0000 0 0 0.0E+00 0.79 1000 1.23E+06 1.05E-07 -0.249563 0.0000 0 0 0.0E+00 0.73 1000 1.52E+06 -2.31E-07 -0.249564 0.0000 0 0 0.0E+00 0.76 1000 1.81E+06 1.37E-07 -0.249563 0.0000 0 0 0.0E+00 0.79 1000 2.08E+06 2.23E-08 -0.249563 0.0000 0 0 0.0E+00 0.85 999 2.33E+06 1.46E-06 -0.250144 0.0000 0 0 0.0E+00 0.86 999 2.60E+06 -8.56E-08 -0.250144 0.0000 0 0 0.0E+00 0.90 998 2.87E+06 -1.39E-06 -0.250691 0.1059 1 23 1.8E-03 0.91 998 3.14E+06 2.31E-06 -0.250690 0.0000 0 0 0.0E+00 0.96 996 3.41E+06 7.44E-06 -0.252447 0.0000 0 0 0.0E+00 0.97 996 3.65E+06 -4.33E-07 -0.252447 0.1300 3 1 3.8E-03 1.00 996 3.88E+06 -1.30E-06 -0.252448 0.2033 1 3 2.6E-03 0.93 996 4.10E+06 -7.92E-07 -0.252448 0.2068 1 3 2.5E-03 1.03 996 4.32E+06 -2.55E-06 -0.252449 0.2348 15 1 1.1E-03 1.13 996 4.51E+06 3.66E-06 -0.252448 0.3376 1 2 1.9E-03 1.17 995 4.70E+06 4.61E-06 -0.252652 0.3992 1 2 1.7E-03 1.13 994 4.87E+06 -8.55E-06 -0.252827 1.0495 9 1 5.8E-04 1.27 991 5.02E+06 -3.18E-07 -0.254033 1.1990 9 1 5.8E-04 1.35 989 5.15E+06 -1.58E-07 -0.254232 1.1653 9 1 5.9E-04 1.37 989 5.29E+06 -3.32E-06 -0.254233 1.1283 9 1 6.0E-04 1.30 988 5.43E+06 -5.26E-06 -0.254416 0.0577 2 11 3.6E-03 1.40 985 5.56E+06 5.17E-07 -0.280915 0.0415 2 11 5.2E-03 1.48 982 5.69E+06 -6.22E-06 -0.400010 0.3254 11 2 3.1E-03 1.54 980 5.81E+06 -1.04E-07 -0.401287 0.3010 11 2 3.2E-03 1.51 979 5.92E+06 1.84E-07 -0.461368 0.1523 11 2 3.2E-03 1.53 979 6.03E+06 -8.71E-06 -0.461371 0.0000 0 0 0.0E+00 1.52 977 6.13E+06 2.99E-09 -0.462270 0.0000 0 0 0.0E+00 1.50 976 6.23E+06 1.25E-07 -0.586273 0.0000 0 0 0.0E+00 1.56 976 6.33E+06 6.79E-07 -0.586273 0.0000 0 0 0.0E+00 1.61 974 6.42E+06 -2.08E-06 -0.586591 0.0000 0 0 0.0E+00 1.62 973 6.51E+06 9.70E-08 -0.630166 0.0000 0 0 0.0E+00 1.64 971 6.60E+06 -4.49E-07 -0.631410 0.0000 0 0 0.0E+00 1.72 971 6.70E+06 -4.71E-06 -0.631411 0.0000 0 0 0.0E+00 1.68 968 6.78E+06 5.78E-08 -0.631672 0.0000 0 0 0.0E+00 1.76 967 6.87E+06 2.54E-07 -0.745841 0.0000 0 0 0.0E+00 1.75 964 6.95E+06 5.21E-08 -0.745854 0.0000 0 0 0.0E+00 1.79 963 7.03E+06 -1.23E-06 -0.746345 0.0000 0 0 0.0E+00 1.66 961 7.12E+06 -2.11E-08 -0.746726 0.0000 0 0 0.0E+00 1.68 960 7.20E+06 -7.07E-07 -0.746828 0.0000 0 0 0.0E+00 1.67 958 7.32E+06 -2.39E-05 -0.747092 0.0000 0 0 0.0E+00 1.63 958 7.40E+06 -3.12E-07 -0.747092 0.0000 0 0 0.0E+00 1.70 957 7.48E+06 6.71E-10 -0.747108 0.0000 0 0 0.0E+00 1.77 956 7.56E+06 1.01E-07 -0.747160 0.0000 0 0 0.0E+00 1.76 956 7.64E+06 -1.52E-06 -0.747160 0.0000 0 0 0.0E+00 1.72 956 7.72E+06 1.53E-08 -0.747160 0.0000 0 0 0.0E+00 1.71 954 7.80E+06 7.74E-08 -0.747155 0.0000 0 0 0.0E+00 1.84 952 7.92E+06 3.11E-05 -0.747192 0.0000 0 0 0.0E+00 1.83 951 8.01E+06 -1.37E-06 -0.747655 0.0000 0 0 0.0E+00

Elapsed real time (secs): 88

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.

File used for above density plots (time, distance, density)
File used for above velocity plots (time, distance, velocity)
File DENS, density as a function of average radii in increasing shells (option KZ7=5)
File VEL, velocity as a function of average radii in increasing shells (option KZ7=5)

HR diagram singles T=0 data file
HR diagram singles T=100. data file

HR diagram(s) for singles, log(lum) vs log(eff temp)

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