https://github.com/florentrenaud/nbody6tt
Raw File
Tip revision: 8aea50c213fd132d500c415511ae1e27eeabab80 authored by florent on 14 February 2015, 16:38:53 UTC
corrected mb typo in ttgalaxy
Tip revision: 8aea50c
rpmax.f
      REAL*8 FUNCTION RPMAX(M1,M2,VSTAR,DV,PERI)
*
*
*       Maximum periastron distance for GR capture.
*       -------------------------------------------
*       Quinlan & Shapiro, Ap.J. 343, 725, 1989.
*
*
      IMPLICIT REAL*8 (A-H,M,O-Z)
      DATA FAC /2.6790/
*
*
*       Note M1, M2 and DV are in N-body units.
      MB = M1 + M2
      MFAC = (M1*M2)**(2.0/7.0)*MB**(3.0/7.0)
      CLIGHT = 3.0D+05/VSTAR
      VFAC = CLIGHT**(10.0/7.0)*DV**(4.0/7.0)
      RPMAX = FAC*MFAC/VFAC
*
*       Include diagnostics for testing (suppressed).
*     WRITE (6,5)  M1, M2, DV, RPMAX, PERI
*   5 FORMAT (' RPMAX    M1 M2 DV RPX PERI ',1P,5E10.2)
*     CALL FLUSH(6)
      RETURN
      END
back to top