This is how bin data are saved every second time level

       program analiz

       parameter (n=64,m=64,l=64)
       real*8 u(n,m,l),v(n,m,l),w(n,m,l)
       real*8  th(n,m,l),qv(n,m,l)
       parameter(ncd=40)
       real*8 dist(n,m,l,ncd)
       dimension dists(n,m,l,ncd)
       parameter(nt=80) 

           do itap=6,6
c standard bin resolution
       if(itap.eq.1) open(10,file='../fort.9.n130.nf',
     1        form='unformatted',status='old')
       if(itap.eq.2) open(10,file='../N26/fort.9.n26.nf',
     1        form='unformatted',status='old')
       if(itap.eq.3) open(10,file='../N650/fort.9.n650.nf',
     1        form='unformatted',status='old')
       if(itap.eq.4) open(10,file='../HTKE/fort.9.hTKE.N26',
     1        form='unformatted',status='old')
       if(itap.eq.5) open(10,file='../HTKE/fort.9.hTKE',
     1     form='unformatted',status='old')
       if(itap.eq.6) open(10,file='../HTKE/fort.9.hTKE.N650',
     1     form='unformatted',status='old')
c
           do it=1,nt

          read(10) !  (((u(i,j,k),i=1,n),j=1,m),k=1,l)
          read(10) !  (((v(i,j,k),i=1,n),j=1,m),k=1,l)
          read(10) !  (((w(i,j,k),i=1,n),j=1,m),k=1,l)
          read(10) !  (((th(i,j,k),i=1,n),j=1,m),k=1,l)
          read(10)   (((qv(i,j,k),i=1,n),j=1,m),k=1,l)
          read(10)    dist

        do i=1,n
        do j=1,m
        do k=1,l
        do ib=1,ncd
        dists(i,j,k,ib)=dist(i,j,k,ib)
        enddo
        enddo
        enddo
        enddo

          amx=-1.e20
          amn= 1.e20
        do i=1,n
        do j=1,m
        do k=1,l
        do ib=1,ncd
        amx=max(amx,dists(i,j,k,ib))
        amn=min(amn,dists(i,j,k,ib))
        enddo
        enddo
        enddo
        enddo
         print*,'-it,min,max: ',it,amn,amx
           
cc write here:
         if(it/2*2.eq.it) then
         write(17) dists
         endif

          enddo

          enddo

          stop
          end