Page 423 - Handbook Of Multiphase Flow Assurance

422 10. Research methods in flow assurance

ymax=unitside * Dsqrt(2.d0)
nx=0
ny=1
nz=0
xmaxh=xmax * 0.5D0
ymaxh=ymax * 0.5D0
endif
C
C
C user dialog requesting surface dimensions
if (A.EQ.2) then
C input slice position for sII hydrate
write(*,*)'Using structure II hydrate',ietc,' slice'
write(*,*)'Reading the unit cell data from the file fd3m.dat'
write(*,*)' z ^'
write(*,*)' | unit cell'
write(6,30)diag
write(*,*)' | |'
write(*,*)' |____________|top position'
write(*,*)' |^ |'
write(*,*)' | thickness |'
write(*,*)' |______v_____|'
write(*,*)' | |'
write(*,*)' 0-|------------|---> x'
write(*,*)' 0'
write(*,*)'Enter the distance from 0 to top surface in Z direction
% - '
if(a.eq.1) then
write(6,31)diag
else
write(6,32)diag
endif
open(unit=1,file='fd3m.hyd',status='old')
C set flag that B is within limits
read(*,*) B
BFlag=0

do while (b.LT.0.)
b=b+diag
bflag=1
end do
do while (b.GT.diag)
b=b-diag
bflag=1
end do