Nbody6
 All Files Functions Variables
tides3.f
Go to the documentation of this file.
1  SUBROUTINE tides3(QPERI,M1,M2,VSTAR,HI,ECC,DE)
2 *
3 *
4 * GR tidal energy loss for interacting stars.
5 * -------------------------------------------
6 *
7  IMPLICIT REAL*8 (a-h,m,o-z)
8  REAL*8 de(2)
9 *
10 *
11 * Specify constants in N-body units.
12  gm = 1.0
13  c = 3.0d+05/vstar
14  fac = (gm/c**2)**3.5
15 *
16 * FE = 1.0 + 73.0/24.0*ECC2 + 37.0/96.0*ECC**4
17 * FE2 = SQRT(ECC**2 - 1.0D0)*(301.0/ ....
18 * Adopt simplified eccentricity factor near ECC = 1.
19  ge = 425.0*3.1415/(32.0*sqrt(2.0d0))
20 *
21 * Set mass and pericentre expression in N-body units.
22  fm = sqrt(m1 + m2)*m1**2*m2**2/qperi**3.5
23 * Form energy change in N-body units (total in DE(1)).
24  de(1) = 8.0/15.0*fac*fm*c**2*ge
25  de(2) = 0.0
26 *
27 * Obtain specific energy change for diagnostic purposes.
28  mu = m1*m2/(m1 + m2)
29  dh = de(1)/mu
30 *
31 * Print diagnostics for capture of hyperbolic orbit.
32  IF (hi.GT.0.0.AND.dh.GT.hi) THEN
33  WRITE (6,5) m1, m2, hi, dh, de(1), qperi
34  5 FORMAT (' CAPTURE CHECK! M1 M2 H DH DE QP ',
35  & 1p,6e9.1)
36  END IF
37 *
38  RETURN
39 *
40  END