Page 422 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 421

C use sII - prepare variables
unitside=17.1D0
side=unitside
npts=136
endif
C
if(a.eq.1) then
C use sI - set variables
unitside=12.D0
side=unitside
npts=46
endif
C
if(ap.eq.2) then
C using [111] surface - prepare variables
side=side*Dsqrt(2.D0)
diag=unitside*Dsqrt(3.D0)
ietc=111
nx=1
ny=1
nz=2
xmax = unitside * Dsqrt(2.D0)
ymax = unitside * Dsqrt(1.5D0)
xmaxh = xmax/2.D0
ymaxh = ymax/2.D0
endif
C
if (ap.eq.1) then
C using [100] surface - prepare variables]
diag=unitside
ietc=100
xmax = unitside
ymax = unitside
xmaxh = xmax/2.D0
ymaxh = ymax/2.D0
nx=0
ny=0
nz=0
endif
C
if (ap.eq.3) then
C [110] surface
diag=unitside * Dsqrt(2.d0)
ietc=110
xmax=unitside