anchor.diff:           rl = real(k-1)*dr
            ru = rl+dr
            norm = const*(ru**3-rl**3)
  -         sum = sum + gr(k)*x
  +         sum = sum + grs(k)*x
   c         write(12,'(i5,F12.6,I5,F12.6,F12.6)') 
   c     >             k,ru,gr(k),sum,real(gr(k))*x/norm
            write(12,'(3F12.6)') 
  -     >             ru,real(gr(k))*x/norm,sum
  +     >             ru,real(grs(k))*x/norm,sum
          end do
          return
   c      error section
  
positive.code:        double precision, dimension(:,:), allocatable :: rd
         character*5, dimension(:), allocatable :: tag
         integer, dimension(:), allocatable :: gr
         character*30 buf
         character*180 bigbuf
         double precision d,dr
  
negative.code:       
  135    continue
          write(*,*) "Number of frames processed",nf
  c
  c      smoothing g(r) per Allan/Tild. p 204 (6.48)
  c      ------------------------------------------
         grs(1)=(69.0*gr(1)+4.0*gr(2)-6.0*gr(3)+
       >         4.0*gr(4)-gr(5))/70.0d0
         grs(2)=(2.0*gr(1)+27.0*gr(2)+12.0*gr(3)
       >  -8.0*gr(4)+2.0*gr(5))/35.0d0
         do i=3,nb-2
           grs(i) =(-3.0*gr(i-2)+12.0*gr(i-1)+17.0*gr(i)+
       >            12.0*gr(i+1)-3.0*gr(i+2))/35.0d0 
         end do
         grs(nb-1)=(2.0*gr(nb)+27*gr(nb-1)+12.0*gr(nb-2)-
       >            8.0*gr(nb-3)+2.0*gr(nb-4))/35.0d0
         grs(nb)=(69.0*gr(nb)+4.0*gr(nb-1)-6.0*gr(nb-2)+
       >          4.0*gr(nb-3)-gr(nb-4))/70.0d0
         open(12,file="gr.dat",
       $            form='formatted',status='unknown')
  
  