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
write (iout,'(8f10.5)') ((c(j,i),j=1,3),i=1,nres)
write (iout,'(8f10.5)') ((c(j,i+nres),j=1,3),i=nnt,nct)
write (iout,*) "Internal coordinates"
+ call intout
write (iout,'(8f10.4)') (vbld(k),k=nnt+1,nct)
write (iout,'(8f10.4)') (vbld(k),k=nres+nnt,nres+nct)
write (iout,'(8f10.4)') (rad2deg*theta(k),k=3,nres)