allareout=1
!C change suggested by Ana -end
do i=1,nres-1
+ mnum=molnum(i)
if ((itype(i,mnum).eq.ntyp1_molec(mnum))&
.and.(itype(i+1,mnum).eq.ntyp1_molec(mnum))) then
chain_end=i
c(1,i)=dmod(c(1,i),boxxsize)
c(2,i)=dmod(c(2,i),boxysize)
c(3,i)=dmod(c(3,i),boxzsize)
+ c(1,i+nres)=dmod(c(1,i+nres),boxxsize)
+ c(2,i+nres)=dmod(c(2,i+nres),boxysize)
+ c(3,i+nres)=dmod(c(3,i+nres),boxzsize)
endif
enddo
return