program read2

  integer, dimension(12)  :: head

  real, dimension(15,64)  :: var
  real                     :: sigi(65),sigl(64)
  character (len=60) :: filename


  write(*,'("Enter input file name: ")',advance='no')
  read(*,'(a60)') filename

  open (1,file=trim(filename),form='unformatted')
  open (7,file=trim(filename)//'.asc',form='formatted')

1000 continue

  !print *, 'New date read in!'

  read(1,end=999) head

  write(*,'(1x,a10,12i5)') "Header :  ",head ; write(7,'(1x,a10,12i5)') "Header :  ",head

  read(1) sigi,sigl

  !  print *, 'number of stations:' , head(9)

  num_times = (head(11)-head(10))/head(12)+1 ! ; print *, num_times


  do l=1,num_times
     write(7,'(a12,i5)') ' TIMESTEP = ',(l-1)*head(12) 
     read(1) var(3,1:head(5)) ; write(7,'(1x,a7)') 'REC = 3' ;write (7,'(8e12.4)') var(3,1:head(5))
     read(1) var(4,1:head(6)) ; write(7,'(1x,a7)') 'REC = 4' ;write (7,'(8e12.4)') var(4,1:head(6)) 
     if (l == 1) write(*,'(1x,a12,2f9.2)') 'Location :  ',var(3,1),var(3,2)
!     if (l == 1) write(7,'(1x,a12,2f9.2)') 'Position :  ',var(3,1),var(3,2)
     do m=5,4+head(7)
        write(7,'(1x,a4,i6)') 'REC = ',m
        read(1,end=999) var(m,1:head(8)) ; write (7,'(8e12.4)') var(m,1:head(8))       
     end do

     write(*,'(11e12.4)') (var(j,1),j=5,15)

  enddo

  goto 1000

999 print *, m,kkk
  close(1)

end program read2