rpmax.f

Go to the documentation of this file.

      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