!     This routine saves trajectories for computing vel.-vel.
!     autocorrelation functions
!     MSE 6270, Leonid Zhigilei
      SUBROUTINE Corlfn()
      INCLUDE 'common.h'
      DIMENSION XVEL(50),YVEL(50),ZVEL(50),RVEL(50)
      NTRA=50
      T=Time*1000.0d0
      DO I = 1, NTRA
          XVEL(I) = Q1D(1,I)/DELTA
          YVEL(I) = Q1D(2,I)/DELTA
          ZVEL(I) = Q1D(3,I)/DELTA
      ENDDO
!     DATA FOR THE VELOCITY AUTO-CORRELATION FUNCTION
      WRITE(23,708) T,(XVEL(I),I=1,ntra)
      WRITE(24,707) (YVEL(I),I=1,ntra)
      WRITE(25,707) (ZVEL(I),I=1,ntra)
!
 708  FORMAT(F12.4,1x,50F9.6)
 707  FORMAT(50F9.6)
      RETURN
      END