Page 407 - Handbook Of Multiphase Flow Assurance

406 10. Research methods in flow assurance

SUBROUTINE REINIT
IMPLICIT REAL*8(A-H,O-Z)
COMMON /PROPTY/ Bmax, Gmax, Rdel
COMMON /ESAVE / Eold,Epold,Epoly(99,99),Ep(99,99),Esurf(99),Es(99)
COMMON /NCC / Kacct, Kaccr, Kaccp
COMMON /NHB / nhb(100),nhbpoly,nhbcurr(100)
COMMON /ZSAVE / zrun(100),epolymer,Erun
COMMON /STEPS / Deltrn, Delrot, Delpiv, Acctrn, Accrot, Accpiv
COMMON /POTPAM/ Casqof,Co,Rsqon,Rsqoff,Rcut,Pi,Tk,Temp
COMMON /NRUN / Ncycle,Ndead,Kprint,Kdprnt,Kcyc,Neq
C
Erun = 0.D0
Epolymer=0.D0

do 10 i=1,100
10 nhb(i)= 0
C
Kacct = 0
Kaccr = 0
Kaccp = 0
C ... acceptance ratios for translation, rotation, pivot
if(Neq.eq.1) then
Acctrn = 0.35D0
Accrot = 0.25D0
Accpiv = 0.35D0
endif
C
C ... output column headings
open(unit=1,POSITION='APPEND',file='output.dat',status='unknown')
ccc open(unit=1,access='APPEND',file='output.dat',status='unknown')
write(1,'(A80)')' Cyc Dtr Rtr Dro Rro Dpv Rpv <Xc> <Yc> <Zc> HB
&p HBs <R> <Epol> <Esur> <Ecur>'
write(1,'(A80)')' x1K A nrm rad nrm rad nrm Angs Angs Angs <cu
&r><cur> A kcal/m kcal/m kcal/m'
close(1)
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'
C
RETURN
END
C
C-------------------------------------------------------------------
C