program tmp
  
integer i, j,k  
real*8 data(10)
character*8 :: stringa

stringa='loffio'
  
open(1, file="to_be_read.dat", form="unformatted")  
rewind(1)
do  i = 1, 10  
  data(i) = ((i - 1) * 5) + (i + 1)  
  write(6,*) stringa,i,i,data(i)  
  write(1) stringa,i,i,data(i)  
enddo
end