findj.f

Go to the documentation of this file.

      SUBROUTINE findj(I,J,IMERGE)
*
*
*       Find merger ghost.
*       ------------------
*
      include 'common6.h'
      common/binary/  cm(4,mmax),xrel(3,mmax),vrel(3,mmax),
     &                hm(mmax),um(4,mmax),umdot(4,mmax),tmdis(mmax),
     &                namem(mmax),nameg(mmax),kstarm(mmax),iflag(mmax)
*
*
*       Check whether the search index is a c.m. body or KS component.
      j = -1
      imerge = 0
      IF (i.GE.ifirst) THEN
          icm = i
          IF (i.LE.n.OR.name(icm).GE.0) THEN
              WRITE (6,1)  i, name(i), ifirst, n
    1         FORMAT (' WARNING!    FINDJ    I NAM I* N ',4i6)
              go to 30
          END IF
      ELSE
*       Set KS pair and c.m. index (exit if no merger).
          jp = kvec(i)
          icm = n + jp
          IF (name(icm).GE.0) go to 30
      END IF
*
*       Locate current position in the merger table.
      imerge = 0
      DO 5 k = 1,nmerge
          IF (namem(k).EQ.name(icm)) THEN
              imerge = k
          END IF
    5 CONTINUE
      IF (imerge.EQ.0) go to 30
*
*       Find index of ghost particle with specified name (search NTOT).
      DO 10 k = 1,ntot
          IF (name(k).EQ.nameg(imerge)) THEN
              j = k
              go to 20
          END IF
   10 CONTINUE
*
   20 IF (j.GT.0.AND.body(j).GT.0.0d0) THEN
          WRITE (6,25)  i, jp, name(icm), j, body(j)
   25     FORMAT (' WARNING!    FINDJ    I JP NAM J MJ',2i5,2i6,1p,e9.1)
          j = -2
      END IF
*
*       Assign nominal value of merger index if not identified.
   30 IF (imerge.EQ.0) imerge = 1
*
      RETURN
*
      END