1 c     a very trivial program to calculate the running 

2 c     average of any observable of a MD run.

3 c     reads column of data from file data.dat and writes

4 c     the running average to "av.dat"

5       program average

6

7       parameter (imax=100000)

8

9       integer i,max

10       real*8 data(imax)

11       real*8 sum

12

13       write(*,*) 'How many values ?'

14       read(*,*) max

15

16       sum=0.d0

17

18       open (33,file='data.dat')

19       open (34,file='av.dat')

20

21       do i=1,max

22       read(33,*) data(i)

23       sum = sum + data(i)

24       write(34,*) sum/i

25       end do

26

27       close (33)

28       close (34)

29

30       end