c Program to measure time to initialize a nxn matrix to a value.
c  a(i,j) = 1 for i=1,...,n, and j=1,...,n.
      program main
      implicit none
      integer i, j, k, n
      integer, allocatable :: a(:,:)
      integer*8 time0, time1, dtime
      real*8 time
c
      write(*,*)"Matrix size    Elapsed time [sec]"
      do k=1,30
c Matrix size
        n = k * 1000

c Allocation
        allocate(a(n,n))

c Start time
        call system_clock(time0)

c Main calculation: initialization
        do i=1,n
          do j=1,n
            a(i,j) = 1
          end do
        end do

c Finish time
        call system_clock(time1, dtime)

c Output time
        time = 1d0*(time1-time0)/dtime
        write(*,"(i12,f16.7)")n, time

c Deallocation
        deallocate(a)
      end do

      end program