chfind.f

Go to the documentation of this file.

      SUBROUTINE chfind
*
*
*       Identification of global chain index.
*       -------------------------------------
*
      include 'common6.h'
      common/chainc/  xc(3,ncmax),uc(3,ncmax),bodyc(ncmax),ich,
     &                listc(lmax)
*
*
*       Re-determine global chain index & perturber list after switching.
      DO 10 i = ifirst,n
          IF (name(i).EQ.0) THEN
              ich = i
*       Update the chain perturber list consistently with new sequence.
              CALL chlist(ich)
*       Resolve global coordinates & velocities of subsystem.
              CALL xcpred(0)
              go to 20
          END IF
   10 CONTINUE
*
   20 RETURN
*
      END