subroutine analbl
      implicit none
      include 'pimc.par'
      integer k,i
      real*8 sa,vold,vnew,ai
      real*8 sang
      external sang
      include 'pimc.cm'

      sa=sang(ndim)/ndim
c compute the g(r) of the single block
      if (np.gt.1) then
       vold=0.0
       k=igr-1
       do i=1,lptable
        k=k+1
        vnew=sa*(dr*i)**ndim
        ai=2./(vnew-vold)/dfloat(np*nslices)
        vold=vnew
        avp(k)=ai*avp(k)
       enddo
      endif

c update global averages
      call cumul(avp(ikin),av(ikin),anorm(ikin),1)
      call cumul(avp(iev),av(iev),anorm(iev),1)
      call cumul(avp(iepot),av(iepot),anorm(iepot),5)
      call cumul(avp(ietot),av(ietot),anorm(ietot),3)
      call cumul(avp(irg),av(irg),anorm(irg),2)
      call cumul(avp(ieff),av(ieff),anorm(ieff),nslices/2)
      call cumul(avp(igr),av(igr),anorm(igr),lptable)
      return
      end