nbrest.f

Go to the documentation of this file.

      SUBROUTINE nbrest(ICM,NSYS,NP)
*
*
*       Restore ghosts in 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
*
*       Add other members of the subsystem subject to limit.
   30     DO 60 k = 1,nsys
              j = jlist(k)
              IF (j.EQ.icm) go to 60
*
*       Skip addition if body #J is already a member.
              DO 40 l = 2,nnb1
                  IF (list(l,i).EQ.j) go to 60
   40         CONTINUE
*
*       Move members down until correct location identified.
              DO 50 l = nnb1,1,-1
                  IF (list(l,i).GT.j.AND.l.GT.1) THEN
                      list(l+1,i) = list(l,i)
                  ELSE
*       Place body #J in vacated location and increase membership.
                      list(l+1,i) = j
                      list(1,i) = list(1,i) + 1
                      nnb1 = nnb1 + 1
*       See whether neighbour list has reached maximum size.
                      IF (nnb1.GT.nnbmax) go to 70
                      go to 60
                  END IF
   50         CONTINUE
   60     CONTINUE
*
*       Also check any KS perturber list.
   70     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