Page 432 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 431

npts=npts*4
side=side*2.D0
xmaxh = xmax
ymaxh = ymax
xmax = xmax + xmax
ymax = ymax + ymax
endif

if(ap.eq.3.and.a.eq.1) then
C for [110] (ap=3)
do 462 ix=0, 1
do 462 i=1, npts
x1(i+ix*npts)=x1(i)+ix
y1(i+ix*npts)=y1(i)
z1(i+ix*npts)=z1(i)
x2(i+ix*npts)=x2(i)+ix
y2(i+ix*npts)=y2(i)
z2(i+ix*npts)=z2(i)
x3(i+ix*npts)=x3(i)+ix
y3(i+ix*npts)=y3(i)
z3(i+ix*npts)=z3(i)
462 continue
npts=npts + npts
xmaxh= xmax
xmax = xmax + xmax
endif

if((a.eq.1.and.ap.eq.2).or.
& (a.eq.2.and.ap.eq.2.and.ngroups.gt.8)) then
C for [111]
do 463 ix=0, 1
do 463 iy=0, 1
do 463 i=1, npts
x1(i+ix*2*npts+iy*npts)=x1(i)+(ix*xmax+iy*xmaxh)/unitside
y1(i+ix*2*npts+iy*npts)=y1(i)+iy*ymax/unitside
z1(i+ix*2*npts+iy*npts)=z1(i)
x2(i+ix*2*npts+iy*npts)=x2(i)+(ix*xmax+iy*xmaxh)/unitside
y2(i+ix*2*npts+iy*npts)=y2(i)+iy*ymax/unitside
z2(i+ix*2*npts+iy*npts)=z2(i)
x3(i+ix*2*npts+iy*npts)=x3(i)+(ix*xmax+iy*xmaxh)/unitside
y3(i+ix*2*npts+iy*npts)=y3(i)+iy*ymax/unitside
z3(i+ix*2*npts+iy*npts)=z3(i)
463 continue
npts=npts*4
C adjust parameters for [111]