This is how superdroplet 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)
       real*8  sum,sum1,sum11,sum12
      parameter(nppg=40)
      parameter(npp=(n-1)*(m-1)*(l-1)*nppg)
c superdroplet positions
      real*8 xp(npp),yp(npp),zp(npp)
c superdroplet cloud water, multiplicity, radius
      real*8 qcp(npp),plic(npp),radp(npp)
c as above but in single precision to be save to fort.17
      dimension xps(npp),yps(npp),zps(npp)
      dimension qcps(npp),plics(npp),radps(npp)

       parameter(nt=80)

           do itap=7,7
        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='../150SDS/fort.9.150.nf',
     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.N26',
     1   form='unformatted',status='old')
        if(itap.eq.7) 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)    xp,yp,zp,qcp,plic,radp 

c check:
          print*,it
c          print*,' u,v,w: ',u(2,2,2),v(2,2,2),w(2,2,2)
c          print*,' th,qv: ',th(2,2,2),qv(2,2,2)
          print*,' xp,yp,zp: ',xp(3),yp(3),zp(3)
          print*,' qcp,plic,radp: ',qcp(3),plic(3),radp(3)

          do i=1,npp
      xps(i)=xp(i)
      yps(i)=yp(i)
      zps(i)=zp(i)
      qcps(i)=qcp(i)
      plics(i)=plic(i)
      radps(i)=radp(i)
      enddo
          print*,' xp,yp,zp: ',xps(3),yps(3),zps(3)
cc write here:
         if(it/2*2.eq.it) then
         write(17) xps,yps,zps,qcps,plics,radps
         endif

          enddo

          enddo

          stop
          end