Page 402 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 401

iy=int(xn(ibonddata(i,j))*scale)+nyoff
write(1,'(2I4)')ix,iy
end do
end do
iix=int(yn(monodata(100,1))*scale)+nxoff
iiy=int(xn(monodata(100,1))*scale)+nyoff
write(1,'(3I4)')iix,iiy,nbackbone
do i=2,nbackbone
ix=int(yn(monodata(100,i))*scale)+nxoff
iy=int(xn(monodata(100,i))*scale)+nyoff
write(1,'(2I4)')ix,iy
end do

C data for YX surface
ir=int(scale/3.d0)
nxoff=60
nyoff=240
if(ap.eq.2)then
nxoff=160
nyoff=240
endif
write(1,'(2I4)')ir,nwatrs
do i=1,nwatrs
iix=int(yw(i,2)*scale)+nxoff
iiy=int(xw(i,2)*scale)+nyoff
write(1,'(2I4)')iix,iiy
ix=int(yw(i,3)*scale)+nxoff
iy=int(xw(i,3)*scale)+nyoff
ixw1=int(yw(i,1)*scale)+nxoff
ixw3=int(yw(i,3)*scale)+nxoff
iyw1=int(xw(i,1)*scale)+nyoff
iyw3=int(xw(i,3)*scale)+nyoff
write(1,'(6I4)')ix,iy,ixw3,iyw3,ixw1,iyw1
end do

C data for YZ polymer
nxoff=320
nzoff=-240
if(ap.eq.2)then
nxoff=480
nzoff=-240
endif
do i=1,natom
ir=iatomtype(i,1)*scale/15
iix=int(yn(i)*scale)+nxoff