Page 435 - Handbook Of Multiphase Flow Assurance

434 10. Research methods in flow assurance

my=25-nint(yw(i,1))
if(nx.lt.1.or.nx.gt.80.or.ny.lt.1.or.ny.gt.25)goto 115
if(lx.lt.1.or.lx.gt.80.or.ly.lt.1.or.ly.gt.25)goto 115
if(mx.lt.1.or.mx.gt.80.or.my.lt.1.or.my.gt.25)goto 115
screen(lx,ly)='.'
screen(mx,my)='.'
screen(nx,ny)='o'
endif
115 continue
screen(1,1)='2'
screen(2,1)='4'
screen(1,5)='2'
screen(2,5)='0'
screen(1,9)='1'
screen(2,9)='6'
screen(1,13)='1'
screen(2,13)='2'
screen(1,17)='8'
screen(1,21)='4'
screen(1,25)='0'
screen(8,25)='4'
screen(16,25)='8'
screen(24,25)='1'
screen(25,25)='2'
screen(32,25)='1'
screen(33,25)='6'
screen(40,25)='2'
screen(41,25)='0'
screen(48,25)='2'
screen(49,25)='4'
screen(56,25)='2'
screen(57,25)='8'
screen(64,25)='3'
screen(65,25)='2'
screen(72,25)='3'
screen(73,25)='6'
endif

if(ap.eq.2) then
C construct a [111] surface
do 120 i=1,Nwatrs
if(zw(i,3).ge.-5.D0) then
C for oxygens
nx=40+nint(xw(i,3)*2.D0)
ny=13-nint(yw(i,3))