Page 436 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 435

C and for hydrogens
lx=40+nint(xw(i,2)*2.D0)
ly=13-nint(yw(i,2))
mx=40+nint(xw(i,1)*2.D0)
my=13-nint(yw(i,1))
if(nx.lt.1.or.nx.gt.80.or.ny.lt.1.or.ny.gt.25)goto 120
if(lx.lt.1.or.lx.gt.80.or.ly.lt.1.or.ly.gt.25)goto 120
if(mx.lt.1.or.mx.gt.80.or.my.lt.1.or.my.gt.25)goto 120
screen(lx,ly)='.'
screen(mx,my)='.'
screen(nx,ny)='o'
endif
120 continue
screen(2,1)='1'
screen(3,1)='2'
screen(2,5)='8'
screen(2,9)='4'
screen(2,13)='0'
screen(2,17)='-'
screen(3,17)='4'
screen(2,21)='-'
screen(3,21)='8'
screen(2,25)='-'
screen(3,25)='1'
screen(4,25)='2'
screen(7,25)='-'
screen(8,25)='1'
screen(9,25)='6'
screen(15,25)='-'
screen(16,25)='1'
screen(17,25)='2'
screen(23,25)='-'
screen(24,25)='8'
screen(31,25)='-'
screen(32,25)='4'
screen(40,25)='0'
screen(48,25)='4'
screen(56,25)='8'
screen(64,25)='1'
screen(65,25)='2'
screen(72,25)='1'
screen(73,25)='6'
endif

C now print the constructed surface