all tests now use generic random number functions

Author: Daniel Ruprecht

test/src/advection_boundaries_test.f90

@@ -34,10 +34,6 @@ PROGRAM advection_boundaries_test
 CALL RANDOM_NUMBER(random_real)
 Nz = N_min + FLOOR( (N_max+1-N_min)*random_real )
 
-write(*,*), Nx
-write(*,*), Ny
-write(*,*), Nz
-
 Nthreads = 8
 
 ALLOCATE(Q(   -2:Nx+3,-2:Ny+3,-2:Nz+3,0:Nthreads-1))

test/src/boundaries_test.f90

@@ -8,16 +8,25 @@ PROGRAM boundaries_test
 
 DOUBLE PRECISION, ALLOCATABLE, DIMENSION(:,:,:) :: Q, Qall
 DOUBLE PRECISION :: dx, dy, dz, x, y, z
-INTEGER :: Nx, Ny, Nz, seed, time1(8), i, j, k
-
-CALL DATE_AND_TIME(values=time1)
-seed = 1000*time1(7)+time1(8)
-CALL SRAND(seed)
+INTEGER :: Nx, Ny, Nz, seed_size, i, j, k, clock
+INTEGER, PARAMETER :: N_min = 25, N_max = 256
+INTEGER, ALLOCATABLE, DIMENSION(:) :: seed
+REAL random_real
+
+CALL RANDOM_SEED(size = seed_size)
+ALLOCATE(seed(seed_size))
+CALL SYSTEM_CLOCK(count = clock)
+seed = clock +  37 * (/ (i - 1, i = 1, seed_size) /)
+CALL RANDOM_SEED(put = seed)
+DEALLOCATE(seed)
 
 ! Create random integer values for Nx, Ny, Nz
-Nx = 25 + INT(RAND()*256)
-Ny = 25 + INT(RAND()*256)
-Nz = 25 + INT(RAND()*256)
+CALL RANDOM_NUMBER(random_real)
+Nx = N_min + FLOOR( (N_max+1-N_min)*random_real )
+CALL RANDOM_NUMBER(random_real)
+Ny = N_min + FLOOR( (N_max+1-N_min)*random_real )
+CALL RANDOM_NUMBER(random_real)
+Nz = N_min + FLOOR( (N_max+1-N_min)*random_real )
 
 dx = 1.0/DBLE(Nx)
 dy = 1.0/DBLE(Ny)
@@ -106,4 +115,4 @@ PROGRAM boundaries_test
 
 PRINT*, '\x1B[32m[0] -- Successful: Ghost cell values generated by boundaries module match directly inserted values.\x1B[0m'
 
-END PROGRAM boundaries_test
+END PROGRAM boundaries_test