endif
endif
do j=1,3
- chomo(j,i,k)=c(j,i)
cref(j,i,cou)=c(j,i)
cref(j,i+nres,cou)=c(j,i+nres)
if (i.le.nres) then
endif
enddo
enddo
+ do i=1,2*nres
+ do j=1,3
+ chomo(j,i,k)=c(j,i)
+ enddo
+ enddo
+
write (iout,*) chain_length
if (chain_length.eq.0) chain_length=nres
do j=1,3