Page 380 - Handbook Of Multiphase Flow Assurance

Computer code (Makogon, 1994, 1997) 379

C assume that equilibration was not done yet
Neq=0
C randomize number generator
Print *, ' '
Print *, ' Enter a seed for random number'
read *,randseed
print *,'input was', randseed
call srand(randseed+0.)
ccc CALL SEED(randseed+0.)
C
if (ncrash.GT.0) then
do NNRn=1,ncrash-1
ra=roulet()
end do
endif
C
Print *, ' '
Print *, ' Enter number of equilibration cycles'
Read *, Ndead
print *,'input was', ndead
C
Print *, ' '
Print *, ' Enter number of run cycles'
Read *, Ncycle
print *,'input was', ncycle
C
Print *, ' '
Print *, ' Enter print interval, in cycles'
Read *, Kprint
print *,'input was', kprint
C
Print *, ' '
Print *, ' Enter interval for graphical output'
Read *, Kdprnt
print *,'input was', kdprnt
C
Print *, ' '
RETURN
END
C
C-------------------------------------------------------------------
C
REAL*8 FUNCTION ROULET()
IMPLICIT REAL*8(A-H,O-Z)
REAL*4 aa
COMMON /CRASH / ncrash, ncalls