if (TUBElog.gt.0) then
call reada(controlcard,"XTUBE",tubecenter(1),0.0d0)
call reada(controlcard,"YTUBE",tubecenter(2),0.0d0)
+ call reada(controlcard,"ZTUBE",tubecenter(3),0.0d0)
call reada(controlcard,"RTUBE",tubeR0,0.0d0)
call reada(controlcard,"TUBETOP",bordtubetop,boxzsize)
call reada(controlcard,"TUBEBOT",bordtubebot,0.0d0)
iread=index(rekord,lancuch(:ilen(lancuch))//"=")
if (iread.eq.0) return
iread=iread+ilen(lancuch)+1
+
read (rekord(iread:),*,end=10,err=10) (tablica(i),i=1,dim)
10 return
end
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