+ DO 81 IGR=1,NGR
+ LI=LICZ(IGR)
+ DO 82 I=1,LI
+ 82 IASS(NCONF(IGR,I))=IGR
+ 81 CONTINUE
+ if (lgrp) then
+ do i=1,ncon
+ iass_tot(i,icut)=iass(i)
+c write (iout,*) icut,i,iass(i),iass_tot(i,icut)
+ enddo
+ endif
+ return
+ end
+c----------------------------------------------------------------------
+ SUBROUTINE WRITE_STATS(ICUT,NCON,IB)
+ implicit real*8 (a-h,o-z)
+ include 'DIMENSIONS'
+ include 'sizesclu.dat'
+ include 'COMMON.CLUSTER'
+ include 'COMMON.FREE'
+ include 'COMMON.IOUNITS'
+ double precision prob(maxgr)
+ write (iout,
+ & '("Free energies, probabilities and rmsds of clusters at",