chdata.f

Go to the documentation of this file.

      SUBROUTINE chdata(XJ,VJ,BODYJ)
*
*
*       Basic data for chain members.
*       -----------------------------
*
      include 'common6.h'
      parameter(nmx=10,nmx3=3*nmx,nmxm=nmx*(nmx-1)/2)
      REAL*8  m,mass,mc,mij,mkk
      common/chain1/  xch(nmx3),vch(nmx3),m(nmx),
     &                zz(nmx3),wc(nmx3),mc(nmx),
     &                xi(nmx3),pi(nmx3),mass,rinv(nmxm),rsum,mkk(nmx),
     &                mij(nmx,nmx),tkk(nmx),tk1(nmx),iname(nmx),nn
      common/chainc/  xc(3,ncmax),uc(3,ncmax),bodyc(ncmax),ich,
     &                listc(lmax)
      common/chreg/  timec,tmax,rmaxc,cm(10),namec(6),nstep1,kz27,kz30
      REAL*4  xj(3,6),vj(3,6),bodyj(6)
*
*
*       Save global address of all chain members in common array JLIST.
      DO 2 l = 1,nch
          DO 1 j = ifirst,n
*       Identify each member sequentially from NAMEC (note NAME(ICH) = 0).
              IF (name(j).EQ.namec(l).OR.name(j).EQ.0) THEN
                  jlist(l) = j
              END IF
    1     CONTINUE
    2 CONTINUE
      
*       Obtain global coordinates & velocities from current chain & c.m.
      lk = 0
      DO 10 i = 1,nch
          bodyj(i) = bodyc(i)
          DO 5 k = 1,3
              lk = lk + 1
              xj(k,i) = xch(lk) + x(k,ich)
              vj(k,i) = vch(lk) + xdot(k,ich)
    5     CONTINUE
   10 CONTINUE
*
      RETURN
*
      END