if (rstcount.eq.1000.or.itime.eq.n_timestep) then
open(irest2,file=rest2name,status='unknown')
write(irest2,*) totT,EK,potE,totE,t_bath
if (rstcount.eq.1000.or.itime.eq.n_timestep) then
open(irest2,file=rest2name,status='unknown')
write(irest2,*) totT,EK,potE,totE,t_bath
write (irest2,'(3e15.5)') (d_t(j,i),j=1,3)
enddo
write (irest2,'(3e15.5)') (d_t(j,i),j=1,3)
enddo
write (irest2,'(3e15.5)') (dc(j,i),j=1,3)
enddo
close(irest2)
write (irest2,'(3e15.5)') (dc(j,i),j=1,3)
enddo
close(irest2)