levels.f

Go to the documentation of this file.

      SUBROUTINE levels
*
*
*       Histograms of block time-step levels.
*       -------------------------------------
*
      include 'common6.h'
      INTEGER  ihist(32),ihistr(32)
*
*
*       Initialize histogram counters.
      DO 10 j = 1,32
          ihist(j) = 0
          ihistr(j) = 0
   10 CONTINUE
*
*       Loop over all single particles & c.m.
      jmax = 0
      jmaxr = 0
      fac = 1.0/log(1.9999999)
      DO 20 i = ifirst,ntot
          IF (body(i).EQ.0.0d0) go to 20
          j = 1 - log(step(i))*fac
          ihist(j) = ihist(j) + 1
          jmax = max(j,jmax)
          j = 1 - log(stepr(i))*fac
          ihistr(j) = ihistr(j) + 1
          jmaxr = max(j,jmaxr)
   20 CONTINUE
*
*       Print histograms of block-steps (STEPR with KZ(33) > 1).
      jmax = min(jmax,27)
      WRITE (6,30)  (ihist(j),j=1,jmax)
   30 FORMAT (' #6  STEP   ',12i6,5i5,10i4)
*
      IF (kz(33).GT.1) THEN
          jmaxr = min(jmaxr,27)
          WRITE (6,35)  (ihistr(j),j=1,jmaxr)
   35     FORMAT (' #7  STEPR  ',10i6,5i5,10i4)
      END IF
*
      RETURN
*
      END