+ enddo
+ else
+ write (iout,
+ & '("Free energies and probabilities of clusters at",
+ & f6.1," K")') 1.0d0/(1.987d-3*beta_h(ib))
+ write(iout,'(a5,2x,a6,3a10)')"clust","efree","prob","sumprob"
+ do i=1,ngr
+ sumprob=sumprob+prob(i)
+ write (iout,'(i3,2x,f8.1,2f10.4)')
+ & i,totfree_gr(i)/beta_h(ib),prob(i),sumprob
+ enddo
+ endif