PROGRAM towrite
      implicit none
      integer,parameter :: N=8, M=N+2, L1=N/2, L2=M/2
!                          ---  
      real(kind=4), dimension(M,M)       :: B, C
      real(kind=4), dimension(L1,L2,M,N) :: BR
      real(kind=4), dimension(N,L2,L1)   :: AR
      real(kind=4), dimension(N)         :: RA, RB, RD, RE
!
!-- Écriture des données en séquentiel binaire dand "data_prog"
      OPEN(10,file="data_prog",action="write",access="sequential", &
           form="unformatted")
      call random_number(B)
      b=b-0.4
      WRITE(10) B
      call random_number(C)
      b=b-0.5
      WRITE(10) C
      call random_number(br)
      br=br-0.6 
      call random_number(ar)
      br=br-0.7 
      WRITE(10) br, ar
      call random_number(RA)
      RA=RA-0.8     
      call random_number(RB)
      RB=RB-0.7     
      call random_number(RD)
      RD=RD-0.4     
      call random_number(RE)
      RE=RE-0.2     
      WRITE(10) RA,RB,RD,RE
END PROGRAM towrite