parameter (nb=36,kgit=48,nk=56,igit=64)
      real an0(nb,kgit,igit,3),an1(nb,kgit,igit,3),
     $     an2(nb,kgit,igit,3),
     $     diff3(nb,kgit,igit,3),difflb(nb,igit),fphi_0(nb,igit),
     $     philb0(nb,igit),philb1(nb,igit)
c----- an0(j,k,i)
      real an0_64(nb,nk,igit,3),an1_64(nb,nk,igit,3),
     $     an2_64(nb,nk,igit,3),
     $     diff3_64(nb,nk,igit,3)
      real an0_64_av(nb,nk,3),an1_64_av(nb,nk,3),
     $     an2_64_av(nb,nk,3),
     $     diff3_64_av(nb,nk,3)
      integer ken1(64) 
      open ( 2,status='old',form='unformatted',file='f2_0_Mar_19')
      open (12,status='new',form='unformatted',file='f2_0_Mar_56')
      read (2) an0
      read (2) an1
      read (2) an2
      read (2) diff3
      read (2) philb0
      read (2) philb1
      read (2) difflb
      read (2) fphi_0
      read (2) ken1 
      close (2)
      rigit=real(igit)
      do m=1,3
        do k=1,kgit
         do j=1,nb
           an0_64_av(j,k,m)=0.
           an1_64_av(j,k,m)=0.
           an2_64_av(j,k,m)=0.
           diff3_64_av(j,k,m)=0.
       end do
       end do
       end do
       do m=1,3
        do k=1,kgit
         do j=1,nb
       do i=1,igit
           an0_64_av(j,k,m)=an0_64_av(j,k,m)+an0(j,k,i,m)/rigit
           an1_64_av(j,k,m)=an1_64_av(j,k,m)+an1(j,k,i,m)/rigit
           an2_64_av(j,k,m)=an2_64_av(j,k,m)+an2(j,k,i,m)/rigit
c           diff3_64_av(j,k,m)=diff3_64_av(j,k,m)+diff3(j,k,i,m)/rigit
       end do
      do i=1,igit
           an0_64(j,k,i,m)=an0_64_av(j,k,m)
           an1_64(j,k,i,m)=an1_64_av(j,k,m)
           an2_64(j,k,i,m)=an2_64_av(j,k,m)
c           diff3_64(j,k,i,m)=diff3_64_av(j,k,m)
      end do
         end do
        end do
      end do
      do m=1,3
       do i=1,igit
         do j=1,nb
        an0_z=an0_64(j,kgit,i,m)-an0_64(j,kgit-1,i,m)
        an1_z=an1_64(j,kgit,i,m)-an1_64(j,kgit-1,i,m)
        an2_z=an2_64(j,kgit,i,m)-an2_64(j,kgit-1,i,m)
c        diff3_z=diff3_64(j,kgit,i,m)-diff3_64(j,kgit-1,i,m)
        do k=kgit+1,nk
      an0_64(j,k,i,m)=an0_64(j,k-1,i,m)
c     $                an0_z*float(nk-k)/float(nk-kgit)
      an1_64(j,k,i,m)=an1_64(j,k-1,i,m)
c     $                an1_z*float(nk-k)/float(nk-kgit)
      an2_64(j,k,i,m)=an2_64(j,k-1,i,m)
c     $                an2_z*float(nk-k)/float(nk-kgit)
      diff3_64(j,k,i,m)=diff3_64(j,k-1,i,m)
c     $                diff3_z*float(nk-k)/float(nk-kgit)
         end do
        end do
       end do
      end do
      ken1(3)=0 
      ken1(4)=0 
      ken1(5)=0 
      ken1(16)=60
      write (12) an0_64
      write (12) an1_64
      write (12) an2_64
c      write (12) diff3_64
      write (12) philb0
      write (12) philb1
c      write (12) difflb
      write (12) fphi_0
      write (12) ken1 
      close (12)
      stop 
      end