diff --git a/assimilation_code/modules/utilities/random_seq_mod.f90 b/assimilation_code/modules/utilities/random_seq_mod.f90 index 22fa38fa5d..b750176593 100644 --- a/assimilation_code/modules/utilities/random_seq_mod.f90 +++ b/assimilation_code/modules/utilities/random_seq_mod.f90 @@ -25,7 +25,8 @@ module random_seq_mod twod_gaussians, & random_gamma, & random_inverse_gamma, & - random_exponential + random_exponential, & + ran_unif, init_ran character(len=*), parameter :: source = 'random_seq_mod.f90' @@ -415,7 +416,7 @@ function ran_unif(s) ! at this point we have an integer value for k ! this routine returns 0.0 <= real < 1.0, so do ! the divide here. return range: [0,1). - +print*, 'k=', k ran_unif = real(real(k, digits12) / 4294967296.0_digits12, r8) end function ran_unif diff --git a/developer_tests/random_seq/work/quickbuild.sh b/developer_tests/random_seq/work/quickbuild.sh index 3dc912094d..99ea98c5c3 100755 --- a/developer_tests/random_seq/work/quickbuild.sh +++ b/developer_tests/random_seq/work/quickbuild.sh @@ -26,6 +26,7 @@ test_hist test_inv_gamma test_random test_reseed +test_ran_unif )