+ imodelclust=0
+ do i=1,nclust
+ do j=1,ninclust(i)
+ imodelclust(icluster(j,i))=1
+ enddo
+ enddo
+ print *,"Models in clusters"
+ nmodel_out = 0
+ do i=1,nmodel
+ if (imodelclust(i).eq.1) then
+ print '(a32)',model(i)
+ nmodel_out = nmodel_out + 1
+ endif
+ enddo
+ print *,"Models excluded"
+ do i=1,nmodel
+ if (imodelclust(i).eq.0) print '(a32)',model(i)
+ enddo
+c Write model information
+ write(2,'(i5)') nmodel_out,nclust
+ do i=1,nmodel
+ if (imodelclust(i).eq.1) write(2,'(a32)') model(i)
+ enddo
+c Write cluster information
+ do i=1,nclust
+ write(2,'(2i5)') ninclust(i),nlenclust(i)
+ write(2,'(16i5)') (icluster(k,i),k=1,ninclust(i))
+ write(2,'(16i5)') (ifrag(k,i),k=1,nlenclust(i))
+ enddo
+ end
+c---------------------------------------------------------------------------------
+ character function charm(i)
+ integer i
+ if (i.eq.0) then
+ charm = "."
+ else
+ charm = char(ichar("A")+i-1)
+ endif
+ return