Environmental Engineering Reference
In-Depth Information
APPENDIX A
PROGRAM 1: TEST A RANDOM NUMBER GENERATOR IN MONTE
CARLO SIMULATION
Note, as an exercise the student may wish to insert other random number
generators or add tests to this simple program.
**************************************************************
This program is used to perform a few very simple tests of a random number gen-
erator. A congruential generator is being tested
**************************************************************
Real*8 Rnum(100,000),Rave,R2Ave,Correl,SDev
Integer Iseed,num
open(Unit=1,®le='result_testrng_02')
PMod = 2147483647.0D0
DMax = 1.0D0/PMod
*******
Input
*******
write(*,800)
800 format('enter the random number generator seed ')
read(*,921) Iseed
921 format(i5)
write(*,801) Iseed
write(1,801) Iseed
801 format(' The random number seed is ', I8)
write(*,802)
802 format('enter the number of random numbers to be generated')
read(*,921) num
write(*,803) num
write(1,803)num
803 format ('number of random numbers to be generated = ',i8)
******************************
Initialize variables, vectors
******************************
do 1 i=2,10000
Rnum(i)=0.0D0
Rave=0.D0
Correl=0.0D0
R2Ave=0.0D0
SDev=0.0D0
*************************
Search WWH ::
Custom Search