Page 424 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 423

if (bflag.eq.0) then
write(6,34)B
else
write(6,35)B
endif
C end of check for proper slice position input
else
C input slice position for sI hydrate
write(*,*)'Using structure I hydrate',ietc,' slice'
write(*,*)'Slice size for sI will be increased to 2x2x1'
write(*,*)'Reading the unit cell data from the file pm3n.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='pm3n.hyd', status='old')
30 FORMAT(F6.2,'|____________')
31 FORMAT(1X,'Choose from 0 to',F8.4,' Angstroms')
32 FORMAT(1X,'Choose from 0 to',F7.3,' Angstroms')
read(*,*) B
C set flag that B is within limits
BFlag=0

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