sumodl_min=1.0d10
do i=0,3**nchain-1
ii=i
- do ichain=1,nchain
+ lshift(1)=0
+ do ichain=2,nchain
lshift(ichain)=mod(ii,3)-1
ii=ii/3
enddo
enddo
#endif
do ichain=1,nchain
+c write (iout,*) "shift_coord",shift_coord(:,ichain)
do i=chain_border1(1,ichain),chain_border1(2,ichain)
do j=1,3
c(j,i)=c(j,i)+shift_coord(j,ichain)