PROGRAM CONVERGE
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER(LD=250)
      DIMENSION G(LD),PCRIT(5),CRIT(5),WPP(3),NOVER(5)
      DATA PCRIT/.001D0,.01D0,.05D0,.10D0,.5D0/
      DATA NCRIT,NLINE/5,40/
C  Preliminaries
      CALL RDSIM0(1,LD,NITER,NPARS)
      READ(5,*)KPARS
      DKPARS=KPARS
      DO 10 ICRIT=1,NCRIT
   10 CRIT(ICRIT)=.5D0*DCHIIN(1.0D0-PCRIT(ICRIT),DKPARS)
      NG=NITER/NLINE
      DG=NG
C  First pass to find largest log data pdf
  100 DO 120 MITER=1,NLINE*NG
      CALL RDSIM(1,ITER,WPP,G)
      PDFLOG=WPP(3)
      IF(MITER.EQ.1)GO TO 110
      IF(PDFLOG.GT.GMAX)GO TO 110
      GO TO 120
  110 JTER=ITER
      GMAX=PDFLOG
  120 CONTINUE
      WRITE(6,130)GMAX,JTER
  130 FORMAT(/,' Largest log data pdf:',1PE14.6,/,
     1         '         at iteration:',I14)
C  Second pass to get summary statistics by group
  200 WRITE(6,210)KPARS,(PCRIT(ICRIT),ICRIT=1,NCRIT)
  210 FORMAT(35X,'Fraction downhill from chisquare(',I3,') at p =',/,
     1       ' Group    Iterations   Max log data pdf',5F8.3)
      write(6,215)(crit(icrit),icrit=1,ncrit)
  215 format(5f8.3)
      REWIND(1)
      READ(1,*)NDUM1,NDUM2
      DO 250 LINE=1,NLINE
      DO 220 ICRIT=1,NCRIT
  220 NOVER(ICRIT)=0
      DO 240 IG=1,NG
      CALL RDSIM(1,ITER,WPP,G)
      IF(IG.EQ.1)ITER1=ITER
      IF(IG.EQ.NG)ITER2=ITER
      PDFLOG=WPP(3)
      IF(IG.EQ.1)AMAX=PDFLOG
  230 IF(PDFLOG.GT.AMAX)AMAX=PDFLOG
      DO 240 ICRIT=1,NCRIT
  240 IF((GMAX-PDFLOG).GT.CRIT(ICRIT))NOVER(ICRIT)=NOVER(ICRIT)+1
  250 WRITE(6,260)LINE,ITER1,ITER2,AMAX,(NOVER(ICRIT)/DG,ICRIT=1,NCRIT)
  260 FORMAT(I6,I7,'-',I7,1PE18.6,0P5F8.3)
      END