Page 416 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 415

c write (1,'(1000X)')
close (1)
RETURN
END
C
C-------------------------------------------------------------------
C
SUBROUTINE RESTARTCHECK(Nstart)
IMPLICIT REAL*8(A-H,O-Z)
COMMON /CRASH / ncrash, ncalls
COMMON /ESAVE / Eold,Epold,Epoly(99,99),Ep(99,99),Esurf(99),Es(99)
COMMON /chain / monodata(100,100),monocount(99),nbackbone
COMMON /NCC / Kacct, Kaccr, Kaccp
COMMON /NHB / nhb(100),nhbpoly,nhbcurr(100)
COMMON /NRUN / Ncycle,Ndead,Kprint,Kdprnt,Kcyc,Neq
COMMON /SPACE / xms(1000), yms(1000), zms(1000)
COMMON /ZSAVE / zrun(100),epolymer,Erun
COMMON /STEPS / Deltrn, Delrot, Delpiv, Acctrn, Accrot, Accpiv
COMMON /types / iatomtype(1000,2),nbondfound,natom,ngroups,ncall
COMMON /RATIOS/ Ratiot, Ratior, Ratiop

Nstart=1
if (ncrash.GT.0) then
read (*,'(2I9)') Kcyc,Neq
Nstart=Kcyc+1
read (*,'(4F32.18)') EOLD,Epold,ERUN,epolymer
C E segments
do 40 i=1,nbackbone-2
40 read (*,'(99F23.18)') (Epoly(i,j),j=i+2,nbackbone)
read (*,'(99F23.18)') (Esurf(i),i=1,nbackbone)
read (*,'(3I16)') Kacct, Kaccr, Kaccp
read (*,'(6F24.18)') Deltrn,Delrot,Delpiv,Ratiot, Ratior, Ratiop
read (*,'(99I10)') (nhb(i),i=1,nbackbone)
read (*,'(99F27.18)') (zrun(i),i=1,100)
read (*,'(2999F22.18)') (xms(i),yms(i),zms(i),i=1,natom)

write(6,'(A80)')' Cyc Dtr Rtr Dro Rro Dpv Rpv <Xc> <Yc> <Zc> HB
&p HBs <R> <Epol> <Esur> <Ecur>'
write(6,'(A80)')' x1K A nrm rad nrm rad nrm Angs Angs Angs <cu
&r><cur> A kcal/m kcal/m kcal/m'
call prtout
c re-zero the current period counters after prtout
do i=1,100
zrun(i)=0.D0
end do