Skip to content

Commit 5ac6a2a

Browse files
program to print out the integer k generated by ran_unif
see #495 for discussion on cce compiler
1 parent 502af78 commit 5ac6a2a

File tree

2 files changed

+4
-2
lines changed

2 files changed

+4
-2
lines changed

assimilation_code/modules/utilities/random_seq_mod.f90

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -25,7 +25,8 @@ module random_seq_mod
2525
twod_gaussians, &
2626
random_gamma, &
2727
random_inverse_gamma, &
28-
random_exponential
28+
random_exponential, &
29+
ran_unif, init_ran
2930

3031
character(len=*), parameter :: source = 'random_seq_mod.f90'
3132

@@ -415,7 +416,7 @@ function ran_unif(s)
415416
! at this point we have an integer value for k
416417
! this routine returns 0.0 <= real < 1.0, so do
417418
! the divide here. return range: [0,1).
418-
419+
print*, 'k=', k
419420
ran_unif = real(real(k, digits12) / 4294967296.0_digits12, r8)
420421

421422
end function ran_unif

developer_tests/random_seq/work/quickbuild.sh

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -26,6 +26,7 @@ test_hist
2626
test_inv_gamma
2727
test_random
2828
test_reseed
29+
test_ran_unif
2930
)
3031

3132

0 commit comments

Comments
 (0)