ghost.f

Go to the documentation of this file.

      SUBROUTINE ghost(J)
*
*
*       Initialization of ghost particle.
*       ---------------------------------
*
      include 'common6.h'
      common/chainc/  xc(3,ncmax),uc(3,ncmax),bodyc(ncmax),ich,
     &                listc(lmax)
*
*
*       Form ghost and initialize integration variables for body #J.
      body(j) = 0.0d0
      t0(j) = 1.0e+06
      DO 10 k = 1,3
          x0dot(k,j) = 0.0d0
          xdot(k,j) = 0.0d0
          f(k,j) = 0.0d0
          fdot(k,j) = 0.0d0
          d0(k,j) = 0.0d0
          d1(k,j) = 0.0d0
          d2(k,j) = 0.0d0
          d3(k,j) = 0.0d0
          d0r(k,j) = 0.0d0
          d1r(k,j) = 0.0d0
          d2r(k,j) = 0.0d0
          d3r(k,j) = 0.0d0
   10 CONTINUE
*
*       Set large X0 & X to avoid perturber selection (no escape).
      x0(1,j) = 1.0e+06
      x(1,j) = 1.0e+06
*
*       Copy neighbour list for routine NBREM.
      nnb = list(1,j)
      DO 20 l = 2,nnb+1
          jpert(l-1) = list(l,j)
   20 CONTINUE
*
*       Remove ghost from lists of neighbours containing body #ICH.
      jlist(1) = j
      CALL nbrem(ich,1,nnb)
*
*       Specify zero neighbour membership.
      list(1,j) = 0
*
      RETURN
*
      END