nbhist.f

Go to the documentation of this file.

      SUBROUTINE nbhist
*
*
*       Neighbour list histogram.
*       -------------------------
*
      include 'common6.h'
      INTEGER  ihist(32)
*
*
*       Initialize histogram counters.
      DO 10 j = 1,32
          ihist(j) = 0
   10 CONTINUE
*
*       Loop over all single particles & c.m.
      jmax = 0
      fac = 1.0/log(1.9999999)
      DO 20 i = ifirst,ntot
          nnb = list(1,i)
          zz = nnb + 1
          j = 1 + log(zz)*fac
          ihist(j) = ihist(j) + 1
          jmax = max(j,jmax)
   20 CONTINUE
*
*       Print histogram of neighbour list membership.
      jmax = min(jmax,12)
      WRITE (6,30)  (ihist(j),j=1,jmax)
   30 FORMAT (' #8  NBLIST    ',12i6)
*
      RETURN
*
      END