Program md
  use mdmod
  implicit none

  integer :: i,icount

  call initialyze

  open(11,file='md.axfs', form='formatted')

  write(11,*) 'ANIMSTEPS', timesteps
  write(11,*) 'CRYSTAL'
  write(11,*) 'PRIMVEC'
  write(11,'(3f15.7)') latticevec(:,1)
  write(11,'(3f15.7)') latticevec(:,2)
  write(11,'(3f15.7)') latticevec(:,3)

  icount=1
  call writecoords(icount,11)

  do i=2,timesteps
     icount=icount+1
     call moveparticles
     call writecoords(icount,11)
  enddo

 end program