nbrem.f

Go to the documentation of this file.

      SUBROUTINE nbrem(ICM,NSYS,NP)
*
*
*       Removal of ghosts from neighbour lists.
*       ---------------------------------------
*
      include 'common6.h'
*
*
*       Examine all members inside neighbour radius of body #ICM.
      DO 100 ll = 1,np
          i = jpert(ll) 
   10     nnb1 = list(1,i) + 1
          IF (nnb1.EQ.1) go to 100
*
*       First see whether body #ICM is a neighbour.
          DO 20 l = 2,nnb1
              IF (list(l,i).EQ.icm) THEN
                  go to 30
              ELSE
                  IF (list(l,i).GT.icm) go to 100
              END IF
   20     CONTINUE
*
*       Remove any other members of the subsystem.
   30     DO 60 k = 1,nsys
              j = jlist(k)
              IF (j.EQ.icm) go to 60
*
*       Determine location of body #J.
              DO 50 l = 2,nnb1
                  IF (list(l,i).EQ.j) THEN
*       Reduce membership and move all subsequent members up by one.
                      list(1,i) = list(1,i) - 1
                      nnb1 = nnb1 - 1
                      DO 40 lj = l,nnb1
                          list(lj,i) = list(lj+1,i)
   40                 CONTINUE
                      go to 60
                  ELSE
                      IF (list(l,i).GT.j) go to 60
                  END IF
   50         CONTINUE
   60     CONTINUE
*
*       Also check any KS perturber list.
          IF (i.GT.n) THEN
              i = 2*(i - n) - 1
              IF (list(1,i).GT.0) go to 10
          END IF
  100 CONTINUE
*
      RETURN
*
      END