open(irest2,file=rest2name,status='unknown')
read(irest2,*) totT,EK,potE,totE,t_bath
totTafm=totT
- do i=1,2*nres
+ do i=0,2*nres
read(irest2,'(3e15.5)') (d_t(j,i),j=1,3)
enddo
- do i=1,2*nres
+C WARING IN OUTER FIELD CRITICAL CORRECTION TO CHANGE i=0
+ do i=0,2*nres
read(irest2,'(3e15.5)') (dc(j,i),j=1,3)
enddo
if(usampl) then