[Dart-dev] DART/branches Revision: 12446
dart at ucar.edu
dart at ucar.edu
Wed Mar 14 15:19:16 MDT 2018
nancy at ucar.edu
2018-03-14 15:19:16 -0600 (Wed, 14 Mar 2018)
101
fix syntax and variable names in the (certainly unused)
pert_model_copies() routine so it compiles.
Modified: DART/branches/rma_netcdf_utils/models/mpas_atm/model_mod.f90
===================================================================
--- DART/branches/rma_netcdf_utils/models/mpas_atm/model_mod.f90 2018-03-14 19:36:43 UTC (rev 12445)
+++ DART/branches/rma_netcdf_utils/models/mpas_atm/model_mod.f90 2018-03-14 21:19:16 UTC (rev 12446)
@@ -1685,7 +1685,7 @@
logical, allocatable :: within_range(:)
real(r8), allocatable :: min_var(:), max_var(:)
integer :: start_ind, end_ind
-real(r8) :: pert_ampl
+real(r8) :: pert_val, range
integer :: copy
integer :: num_variables
integer :: i, j
@@ -1744,12 +1744,12 @@
! make the perturbation amplitude a fraction of the
! entire variable range.
range = max_var(i) - min_var(i)
- pert_ampl = pert_amount * range
+ pert_val = model_perturbation_amplitude * range ! this is a namelist item
do j=1, ens_handle%my_num_vars
if (ens_handle%my_vars(j) >= start_ind .and. ens_handle%my_vars(j) <= end_ind) then
do copy = 1, ens_size
- ens_handle%copies(copy, j) = random_gaussian(random_seq, ens_handle%copies(copy, j), pert_ampl)
+ ens_handle%copies(copy, j) = random_gaussian(random_seq, ens_handle%copies(copy, j), pert_val)
enddo
! keep variable from exceeding the original range
More information about the Dart-dev
mailing list