write (iout,'(i3,3f10.5,3x,3f10.5)') i,(d_t(j,i),j=1,3),
& (d_t(j,i+nres),j=1,3)
enddo
+ call flush(iout)
c Zeroing the total angular momentum of the system
write(iout,*) "Calling the zero-angular
& momentum subroutine"
if(me.eq.king.or..not.out1file)then
write (iout,*) "vcm right after adjustment:"
write (iout,*) (vcm(j),j=1,3)
+ call flush(iout)
endif
if (.not.rest) then
call chainbuild
c & " d_t_work_new",d_t_work_new(ii)
enddo
enddo
+ call flush(iout)
c diagnostics
c Ek1=0.0d0
c ii=0
c do i=0,nres-1
c write (iout,*) "d_t",i,(d_t(j,i),j=1,3)
c enddo
+c call flush(iout)
do i=nnt,nct
if (itype(i).ne.10) then
do j=1,3