Revision dd8daa98413988b11c594e5a3bff8f5659d7ecf7 authored by Florent Renaud on 26 June 2014, 09:41:19 UTC, committed by Florent Renaud on 26 June 2014, 09:41:19 UTC
Originally released on 18 March 2013 Based on Nbody6 downloaded on 29 January 2013
1 parent edab1b7
peters.f
IMPLICIT REAL*8 (A-H,O-Z)
REAL M1,M2
*
M1 = 3.0D-04
M2 = 3.0D-04
ECCGR = 0.999
SEMIGR = 2.0D-05
CVEL = 18000.0
E2 = ECCGR**2
ADOT = 64.0/5.0*M1*M2*(M1+M2)/(CVEL**5*SEMIGR**3*(1.0-E2)**3.5)
ADOT = ADOT*(1.0 + 73.0/24.0*E2 + 37.0/96.0*E2**2)
EDOT = 304.0/15.0*ECCGR*M1*M2*(M1+M2)/(CVEL**5*SEMIGR**4)
EDOT = EDOT/(1.0-E2)**2.5*(1.0 + 121.0/304.0*E2)
TGR = SEMIGR/ADOT
EGR = ECCGR/EDOT
*
SMU = 6.0D+04
ASUN = SEMIGR*2.0D+05*215.0
* Include Banerjee formula in yrs (SEMI in solar radii & SMU in Msun).
TM = 1.5D+08*ASUN**4/(M1*SMU)**3*(1.0-E2)**3.5
THETA = 3.0*6.2830*(M1 + M2)/(SEMIGR*CVEL**2*(1.0-E2))
WRITE (6,10) ECCGR, SEMIGR, TGR, EGR, THETA, TM
10 FORMAT (' E A TGR EGR THETA TM ',F10.5,1P,5E10.2)
END
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...