+ iset=1
+ nperm=1
+ print *, "AFTER read fragments"
+ write (iout,*) "iset",iset
+ if (loc_qlike) then
+ write(iout,*) "fragment, weights, q0:"
+ do i=1,nfrag_back
+ write(iout,'(2i5,3(f8.1,f8.2))') ifrag_back(1,i,iset),
+ & ifrag_back(2,i,iset),
+ & wfrag_back(1,i,iset),qin_back(1,i,iset),
+ & wfrag_back(2,i,iset),qin_back(2,i,iset),
+ & wfrag_back(3,i,iset),qin_back(3,i,iset)
+ enddo
+ else
+ write(iout,*) "fragment, weights:"
+ do i=1,nfrag_back
+ write(iout,'(2i5,3f8.1)') ifrag_back(1,i,iset),
+ & ifrag_back(2,i,iset),wfrag_back(1,i,iset),
+ & wfrag_back(2,i,iset),wfrag_back(3,i,iset)
+ enddo
+ endif
+ call read_REMDpar
+ call PMFread
+ call rescale_weights(t_bath)