c     a very trivial program to calculate the running 
c     average of any observable of a MD run.
c     reads column of data from file data.dat and writes
c     the running average to "av.dat"
      program average

      parameter (imax=100000)

      integer i,max
      real*8 data(imax)
      real*8 sum

      write(*,*) 'How many values ?'
      read(*,*) max

      sum=0.d0

      open (33,file='data.dat')
      open (34,file='av.dat')

      do i=1,max
      read(33,*) data(i)
      sum = sum + data(i)
      write(34,*) sum/i
      end do

      close (33)
      close (34)

      end