Page 403 - Handbook Of Multiphase Flow Assurance

402 10. Research methods in flow assurance

iiz=int(zn(i)*scale)+nzoff
write(1,'(5I4)')ir,iix,-iiz,nbonds(i),iatomtype(i,1)
do j=1,nbonds(i)
ix=int(yn(ibonddata(i,j))*scale)+nxoff
iz=int(zn(ibonddata(i,j))*scale)+nzoff
write(1,'(2I4)')ix,-iz
end do
end do
iix=int(yn(monodata(100,1))*scale)+nxoff
iiz=int(zn(monodata(100,1))*scale)+nzoff
write(1,'(3I4)')iix,-iiz,nbackbone
do i=2,nbackbone
ix=int(yn(monodata(100,i))*scale)+nxoff
iz=int(zn(monodata(100,i))*scale)+nzoff
write(1,'(2I4)')ix,-iz
end do
C data for YZ surface
ir=int(scale/3.D0)
nxoff=320
nzoff=-240
if(ap.eq.2)then
nxoff=480
nzoff=-240
endif
write(1,'(I4)')ir
do i=1,nwatrs
iix=int(yw(i,2)*scale)+nxoff
iiz=int(zw(i,2)*scale)+nzoff
write(1,'(2I4)')iix,-iiz
ix=int(yw(i,3)*scale)+nxoff
iz=int(zw(i,3)*scale)+nzoff
ixw1=int(yw(i,1)*scale)+nxoff
ixw3=int(yw(i,3)*scale)+nxoff
izw1=int(zw(i,1)*scale)+nzoff
izw3=int(zw(i,3)*scale)+nzoff
write(1,'(6I4)')ix,-iz,ixw3,-izw3,ixw1,-izw1
end do
C current simulation information
Deno = 3.D0*DFLOAT(Kcyc)
if(kcyc.eq.0)deno=1.D0
Eave = Erun/Deno
Epol = Epolymer/Deno
nbb=nbackbone
nhbs=0
do 10 i=1,nbackbone