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
! AL 4/17/17: Now writing d_t(0,:) too
do i=0,2*nres
write (irest2,'(3e15.5)') (d_t(j,i),j=1,3)
! AL 4/17/17: Now writing d_t(0,:) too
do i=0,2*nres
write (irest2,'(3e15.5)') (d_t(j,i),j=1,3)
!el common /przechowalnia/ GGinv,gdc,Cmat,nbond
!el common /przechowalnia/ nbond
integer :: max_rattle = 5
!el common /przechowalnia/ GGinv,gdc,Cmat,nbond
!el common /przechowalnia/ nbond
integer :: max_rattle = 5