Nbody6
 All Files Functions Variables
grrad.f
Go to the documentation of this file.
1 ***
2  SUBROUTINE grrad(m1,m2,sep,ecc,jorb,djgr,de)
3  implicit none
4  real*8 m1,m2,sep,ecc,jorb,djgr,de,acrit
5  real*8 ecc2,sqome2,sqome5,f1
6 *
7 * Include gravitational radiation for very close or massive systems.
8 *
9  djgr = 0.d0
10  de = 0.d0
11  acrit = 10.0
12  IF (m1 + m2.GT.10.0) acrit = 30.0
13  if(sep.le.acrit)then
14  djgr = 8.315d-10*m1*m2*(m1+m2)/(sep*sep*sep*sep)
15  ecc2 = ecc*ecc
16  sqome2 = sqrt(1.d0-ecc2)
17  sqome5 = sqome2**5
18  f1 = (19.d0/6.d0) + (121.d0/96.d0)*ecc2
19  de = djgr*ecc*f1/sqome5
20  djgr = djgr*jorb*(1.d0+0.875d0*ecc2)/sqome5
21  endif
22 *
23  RETURN
24  END
25 ***