do i=1,nres
if (molnum(i).eq.5) then
c(1,i)=dmod(c(1,i),boxxsize)
+ if (c(1,i).lt.0) c(1,i)=c(1,i)+boxxsize
c(2,i)=dmod(c(2,i),boxysize)
+ if (c(2,i).lt.0) c(2,i)=c(2,i)+boxysize
c(3,i)=dmod(c(3,i),boxzsize)
+ if (c(3,i).lt.0) c(3,i)=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)