[Dart-dev] DART/branches Revision: 12441
dart at ucar.edu
dart at ucar.edu
Mon Mar 12 16:40:40 MDT 2018
nancy at ucar.edu
2018-03-12 16:40:40 -0600 (Mon, 12 Mar 2018)
167
tim renamed random_number() to random_array() to avoid a builtin
in one of the fortran compilers. keep the doppler fold version in
sync with the main filter module.
Modified: DART/branches/rma_netcdf_utils/assimilation_code/modules/assimilation/filter_mod.dopplerfold.f90
===================================================================
--- DART/branches/rma_netcdf_utils/assimilation_code/modules/assimilation/filter_mod.dopplerfold.f90 2018-03-12 15:24:12 UTC (rev 12440)
+++ DART/branches/rma_netcdf_utils/assimilation_code/modules/assimilation/filter_mod.dopplerfold.f90 2018-03-12 22:40:40 UTC (rev 12441)
@@ -1989,7 +1989,7 @@
integer :: i, j ! loop variables
type(random_seq_type) :: r(ens_size)
-real(r8) :: random_number(ens_size) ! array of random numbers
+real(r8) :: random_array(ens_size) ! array of random numbers
integer :: local_index
! Need ens_size random number sequences.
@@ -2005,11 +2005,11 @@
do j = 1, ens_size
! Can use %copies here because the random number
! is only relevant to the task than owns element i.
- random_number(j) = random_gaussian(r(j), ens_handle%copies(j, local_index), perturbation_amplitude)
+ random_array(j) = random_gaussian(r(j), ens_handle%copies(j, local_index), perturbation_amplitude)
enddo
if (ens_handle%my_vars(local_index) == i) then
- ens_handle%copies(1:ens_size, local_index) = random_number(:)
+ ens_handle%copies(1:ens_size, local_index) = random_array(:)
local_index = local_index + 1 ! task is ready for the next random number
local_index = min(local_index, ens_handle%my_num_vars)
endif
More information about the Dart-dev
mailing list