https://github.com/nbodyx/Nbody6
Raw File
Tip revision: dff3b5c68673d4d4c4c9f54c8ba110b8416098f2 authored by nitadori on 31 May 2020, 13:04:00 UTC
Avoided touching uninitialized SEMIX
Tip revision: dff3b5c
xtrnlp.f
      SUBROUTINE XTRNLP(Q1,Q3,FP)
*
*
*       External KS perturbation.
*       -------------------------
*
      INCLUDE 'common6.h'
      REAL*8  FP(3)
*
*
*       See whether to include the galactic tidal force.
      IF (TIDAL(1).GT.0.0) THEN
          FP(1) = FP(1) + TIDAL(1)*Q1
          FP(3) = FP(3) + TIDAL(3)*Q3
*       Omit Coriolis terms which do not affect the binding energy.
      END IF
*
*       NB! Second call from routine KSPOLY uses velocity argument.
*
      RETURN
*
      END
back to top