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:       
  135    continue
          write(*,*) "Number of frames processed",nf
  
         open(12,file="gr.dat",
       $            form='formatted',status='unknown')
  
  
negative.code:        double precision, dimension(:,:), allocatable :: rd
         character*5, dimension(:), allocatable :: tag
         integer, dimension(:), allocatable :: gr
         double precision, dimension(:), allocatable :: grs
         character*30 buf
         character*180 bigbuf
         double precision d,dr
  