Page 405 - Handbook Of Multiphase Flow Assurance

404 10. Research methods in flow assurance

RETURN
ENDIF
ENDIF
ENDIF
C ENDIF
C
C ... Reject proposed move
Erun = Erun + Eold
CALL ZSAMPLE
C
RETURN
END
C
C-------------------------------------------------------------------
C
SUBROUTINE ACCEPT(Ener, Iflag)
IMPLICIT REAL*8(A-H,O-Z)
COMMON /ESAVE / Eold,Epold,Epoly(99,99),Ep(99,99),Esurf(99),Es(99)
COMMON /ZSAVE / zrun(100),epolymer,Erun
COMMON /types / iatomtype(1000,2),nbondfound,natom,ngroups,ncall
COMMON /chain / monodata(100,100),monocount(99),nbackbone
COMMON /NATOMS/ Norgn, Nwatrs, M1, M2, M3, M4
COMMON /NCC / Kacct, Kaccr, Kaccp
COMMON /NWPOS / xn(1000), yn(1000), zn(1000)
COMMON /SPACE / xms(1000), yms(1000), zms(1000)
C
C assign new values for polymer segments-surface interactions
do 10 i=1,nbackbone
10 Esurf(i)=Es(i)

c assign new values for the intramolecular polymer segments interactions
if(Iflag.eq.2)then
Epold=0.D0
do 20 i=1,nbackbone-2
do 20 j=i+2,nbackbone
Epoly(i,j)=Ep(i,j)
20 Epold=Epold+Epoly(i,j)
endif
C assign new coordinates to polymer
DO 30 I = 1, natom
xms(I) = xn(I)
yms(I) = yn(I)
zms(I) = zn(I)
30 CONTINUE
C