transx.f

Go to the documentation of this file.

      SUBROUTINE transx
*
*
*       Transformations from KS to chain variables.
*     ---------------------------------------------
*
      include 'commonc.h'
      include 'common2.h'
      REAL*8 q0(3)
*
*
*   Transform to chain coordinates.
      DO i=1,n-1
      l1=3*(i-1)+1
      ks1=4*(i-1)+1
      CALL ksphys(q(ks1),p(ks1),xc(l1),wc(l1))
      END DO
*
*   Obtain physical variables from chain quantities.
      l=3*(n-2)
      DO k=1,3
      pi(k)=-wc(k)
      pi(l+k+3)=wc(l+k)
      END DO
      DO i=2,n-1
      l=3*(i-1)
      DO k=1,3
      pi(l+k)=wc(l+k-3)-wc(l+k)
      END DO
      END DO
*
      DO k=1,3
      xi(k)=0.0
      q0(k)=0.0
      END DO
*
      DO i=1,n-1
      l=3*(i-1)
      DO k=1,3
      xi(l+3+k)=xi(l+k)+xc(l+k)
      END DO
      END DO
*
      DO i=1,n
      l=3*(i-1)
      DO k=1,3
      q0(k)=q0(k)+xi(l+k)*mc(i)/mass
      END DO
      END DO
*
*   Rearrange according to INAME(I) and add CM (if required).
      DO i=1,n
      l=3*(i-1)
      lf=3*(iname(i)-1)
      DO k=1,3
      x(lf+k)=xi(l+k)-q0(k)
*    &       +CMX(K)
      v(lf+k)=pi(l+k)/mc(i)
*    &       +CMV(K)
      END DO
      END DO
*
      RETURN
      END