!     THIS SUBROUTINE CLEANS UP AFTER THE MD LOOP
!     MSE 6270, Leonid Zhigilei

      SUBROUTINE SetEnd()
      INCLUDE 'common.h'

!     Since we multiplied velocities by timestep DELTA
!     at the beginning of the run (in SetInit.f) we have
!     to divide by DELTA here to get real velocities
      DO I=1,NAN3
        Q1(I)=Q1(I)/DELTA
      ENDDO

      RETURN
      END