[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