https://github.com/florentrenaud/nbody6tt
Tip revision: 8a4716382ead3ece116c48a4ae5c65f8c9534437 authored by Florent on 29 January 2015, 12:19:28 UTC
Nbody6 - 29 January 2015 (added GPU2/Build/.gitkeep)
Nbody6 - 29 January 2015 (added GPU2/Build/.gitkeep)
Tip revision: 8a47163
zero.f
SUBROUTINE ZERO
*
*
* Initialization of global scalars.
* ---------------------------------
*
INCLUDE 'common6.h'
*
*
* Initialize parameters & counters and set useful constants.
TIME = 0.0D0
TOFF = 0.0D0
TTOT = 0.0D0
TADJ = 0.0D0
TBLIST = 0.0D0
TPLOT = 0.0D0
DTB = 0.0D0
TNEXT = 0.0D0
TMDOT = 0.0D0
ZMDOT = 0.0D0
TLASTS = 0.0D0
TLASTT = 0.0D0
RMAX = 0.0
EPOCH0 = 0.0
DMIN1 = 100.0
DMIN2 = 100.0
DMIN3 = 100.0
DMIN4 = 100.0
DMINC = 100.0
CPUTOT = 0.0
WTOT = 0.0D0
ERROR = 0.0D0
ERRTOT = 0.0D0
DETOT = 0.0D0
ETCORR = 0.0
EBIN0 = 0.0D0
ESUB = 0.0D0
EMERGE = 0.0D0
EDISS = 0.0D0
ECOLL = 0.0D0
ECDOT = 0.0D0
EMDOT = 0.0D0
EKICK = 0.0
ESYNC = 0.0D0
ETIDE = 0.0D0
EPL = 0.0
ESESC = 0.0D0
EBESC = 0.0D0
EMESC = 0.0D0
EGRAV = 0.0
SBCOLL = 0.0D0
BBCOLL = 0.0D0
CHCOLL = 0.0D0
EBH = -5.0D-06 ! temporary def for KSINIT with primordials.
NPAIRS = 0
NBPREV = 0
NNTB = 0
NMERGE = 0
NEWHI = 0
NSUB = 0
NCH = 0
NBIN0 = 0
NHI0 = 0
NDISR = 0
IFIRST = 1
IPHASE = 0
ICOMP = 0
JCOMP = 0
KSPAIR = 0
NPRINT = 0
MODEL = 0
NDUMP = 0
NSTEPS = 0
NTIMER = 0
LISTR(1) = 0
LISTD(1) = 0
LISTV(1) = 0
NTAIL = 0
NTTOT = 0
NRSAVE = 0
NNBOPT = 0 ! NBODY6++ compatibility.
*
DO 10 J = 1,KMAX
KVEC(2*J-1) = J
KVEC(2*J) = J
10 CONTINUE
*
DO 20 J = 1,12
E(J) = 0.0D0
20 CONTINUE
*
DO 30 J = 1,9
TLASTB(J) = 0.0D0
30 CONTINUE
*
DO 40 J = 1,NMAX
RADIUS(J) = 0.0D0
BODY0(J) = 0.0
EPOCH(J) = 0.0
TEV(J) = 0.0
TEV0(J) = 0.0
KSTAR(J) = 0
40 CONTINUE
*
DO 50 K = 1,3
RDENS(K) = 0.0D0
TIDAL(K) = 0.0
50 CONTINUE
TIDAL(4) = 0.0
*
* Set fractional constants & two PI.
ONE3 = 1.0D0/3.0D0
ONE6 = 1.0D0/6.0D0
ONE9 = 1.0D0/9.0D0
ONE12 = 1.0D0/12.0D0
TWOPI = 8.0D0*ATAN(1.0D0)
TINY = 1.0D-14
*
RETURN
*
END