c(j,i+nres+nnt-1)=xoord(j,i+nres)
enddo
enddo
+c write (iout,*) "Before boxshift"
+c call flush(iout)
c Box shift
call oligomer
+c write (iout,*) "After oligomer"
+c call flush(iout)
do i=1,nres
do j=1,3
xoord(j,i)=c(j,i)
enddo
enddo
c end change
-
+c write (iout,*) "Before islice"
+c call flush(iout)
if (islice.gt.0 .and. islice.le.nslice .and. (.not.separate_parset
& .or. iset.eq.myparm)) then
ii=ii+1