Merge pull request #636 from NCAR/WRF_Tutorial_Derecho

Update WRF scripting and tutorial for derecho

Author: mjs2369

.github/ISSUE_TEMPLATE/bug_report.md

@@ -34,5 +34,5 @@ If your code changes are available on GitHub, please provide the repository.
 
 ### Build information
 Please describe:  
- 1. The machine you are running on (e.g. windows laptop, NCAR supercomputer Cheyenne).   
+ 1. The machine you are running on (e.g. windows laptop, NSF NCAR supercomputer Derecho).   
  2. The compiler you are using (e.g. gnu, intel).  

CHANGELOG.rst

@@ -22,6 +22,13 @@ individual files.
 
 The changes are now listed with the most recent at the top.
 
+**March 13 2024 :: Update WRF-DART scripts and bug template to Derecho; remove no-op routines in ensemble manager. Tag v11.3.1**
+
+- Updated the csh scripting templates used to run WRF-DART and WRF-DART tutorial from Cheyenne to Derecho
+- Updated bug report template to use Derecho instead of Cheyenne
+- Removed the following no-op routines from ensemble manager: prepare_to_write_to_vars, prepare_to_write_to_copies,
+  prepare_to_read_from_vars, prepare_to_read_from_copies, prepare_to_update_vars, prepare_to_update_copies
+
 **March 12 2024 :: MITgcm/N-BLING with Compressed Staggered Grids. Tag v11.3.0**
 
 - The DART-MITgcm code now supports compressed grids, especially suited for areas like 

assimilation_code/modules/assimilation/assim_tools_mod.pf.f90

@@ -62,7 +62,7 @@ module assim_tools_mod
 
 use ensemble_manager_mod, only : ensemble_type, get_my_num_vars, get_my_vars,             &
                                  compute_copy_mean_var, get_var_owner_index,              &
-                                 prepare_to_update_copies, map_pe_to_task
+                                 map_pe_to_task
 
 use mpi_utilities_mod,    only : my_task_id, broadcast_send, broadcast_recv,              &
                                  sum_across_tasks, task_count, start_mpi_timer,           &
@@ -496,10 +496,6 @@ subroutine filter_assim(ens_handle, obs_ens_handle, obs_seq, keys,           &
          my_state_loc(         ens_handle%my_num_vars))
 ! end alloc
 
-! we are going to read/write the copies array
-call prepare_to_update_copies(ens_handle)
-call prepare_to_update_copies(obs_ens_handle)
-
 ! Initialize assim_tools_module if needed
 if (.not. module_initialized) call assim_tools_init()
 

assimilation_code/modules/assimilation/filter_mod.dopplerfold.f90

@@ -46,11 +46,10 @@ module filter_mod
                                   compute_copy_mean, compute_copy_mean_sd,                    &
                                   compute_copy_mean_var, duplicate_ens, get_copy_owner_index, &
                                   get_ensemble_time, set_ensemble_time, broadcast_copy,       &
-                                  map_pe_to_task, prepare_to_update_copies,  &
-                                  copies_in_window, set_num_extra_copies, get_allow_transpose, &
-                                  all_copies_to_all_vars, allocate_single_copy, allocate_vars, &
-                                  get_single_copy, put_single_copy, deallocate_single_copy,   &
-                                  print_ens_handle
+                                  map_pe_to_task, copies_in_window, set_num_extra_copies,     &
+                                  get_allow_transpose, all_copies_to_all_vars,                &
+                                  allocate_single_copy, allocate_vars, get_single_copy,       &
+                                  put_single_copy, deallocate_single_copy, print_ens_handle
 
 use adaptive_inflate_mod,  only : do_ss_inflate, mean_from_restart, sd_from_restart,  &
                                   inflate_ens, adaptive_inflate_init,                 &
@@ -806,7 +805,6 @@ subroutine filter_main()
       call trace_message('Before prior inflation damping and prep')
 
       if (inf_damping(PRIOR_INF) /= 1.0_r8) then
-         call prepare_to_update_copies(state_ens_handle)
          state_ens_handle%copies(PRIOR_INF_COPY, :) = 1.0_r8 + &
             inf_damping(PRIOR_INF) * (state_ens_handle%copies(PRIOR_INF_COPY, :) - 1.0_r8)
       endif
@@ -907,7 +905,6 @@ subroutine filter_main()
       call trace_message('Before posterior inflation damping')
 
       if (inf_damping(POSTERIOR_INF) /= 1.0_r8) then
-         call prepare_to_update_copies(state_ens_handle)
          state_ens_handle%copies(POST_INF_COPY, :) = 1.0_r8 + &
             inf_damping(POSTERIOR_INF) * (state_ens_handle%copies(POST_INF_COPY, :) - 1.0_r8)
       endif
@@ -1549,9 +1546,6 @@ subroutine filter_ensemble_inflate(ens_handle, inflate_copy, inflate, ENS_MEAN_C
 
 integer :: j, group, grp_bot, grp_top, grp_size
 
-! Assumes that the ensemble is copy complete
-call prepare_to_update_copies(ens_handle)
-
 ! Inflate each group separately;  Divide ensemble into num_groups groups
 grp_size = ens_size / num_groups
 
@@ -2827,8 +2821,6 @@ subroutine update_observations_radar(obs_ens_handle, ens_size, seq, keys, prior_
 ! for quiet execution, set it to false.
 verbose = .true.
 
-call prepare_to_update_copies(obs_ens_handle)
-
 do j = 1, obs_ens_handle%my_num_vars
    ! get the key number associated with each of my subset of obs
    ! then get the obs and extract info from it.

assimilation_code/modules/assimilation/filter_mod.f90

@@ -46,11 +46,10 @@ module filter_mod
                                   compute_copy_mean, compute_copy_mean_sd,                    &
                                   compute_copy_mean_var, duplicate_ens, get_copy_owner_index, &
                                   get_ensemble_time, set_ensemble_time, broadcast_copy,       &
-                                  map_pe_to_task, prepare_to_update_copies,  &
-                                  copies_in_window, set_num_extra_copies, get_allow_transpose, &
-                                  all_copies_to_all_vars, allocate_single_copy, allocate_vars, &
-                                  get_single_copy, put_single_copy, deallocate_single_copy,   &
-                                  print_ens_handle
+                                  map_pe_to_task, copies_in_window, set_num_extra_copies,     &
+                                  get_allow_transpose, all_copies_to_all_vars,                &
+                                  allocate_single_copy, allocate_vars, get_single_copy,       &
+                                  put_single_copy, deallocate_single_copy, print_ens_handle
 
 use adaptive_inflate_mod,  only : do_ss_inflate, mean_from_restart, sd_from_restart,  &
                                   inflate_ens, adaptive_inflate_init,                 &
@@ -809,7 +808,6 @@ subroutine filter_main()
       call trace_message('Before prior inflation damping and prep')
 
       if (inf_damping(PRIOR_INF) /= 1.0_r8) then
-         call prepare_to_update_copies(state_ens_handle)
          state_ens_handle%copies(PRIOR_INF_COPY, :) = 1.0_r8 + &
             inf_damping(PRIOR_INF) * (state_ens_handle%copies(PRIOR_INF_COPY, :) - 1.0_r8)
       endif
@@ -910,7 +908,6 @@ subroutine filter_main()
       call trace_message('Before posterior inflation damping')
 
       if (inf_damping(POSTERIOR_INF) /= 1.0_r8) then
-         call prepare_to_update_copies(state_ens_handle)
          state_ens_handle%copies(POST_INF_COPY, :) = 1.0_r8 + &
             inf_damping(POSTERIOR_INF) * (state_ens_handle%copies(POST_INF_COPY, :) - 1.0_r8)
       endif
@@ -1566,9 +1563,6 @@ subroutine filter_ensemble_inflate(ens_handle, inflate_copy, inflate, ENS_MEAN_C
 real(r8) :: lower_bound,   upper_bound
 integer  :: dist_type
 
-! Assumes that the ensemble is copy complete
-call prepare_to_update_copies(ens_handle)
-
 ! Inflate each group separately;  Divide ensemble into num_groups groups
 grp_size = ens_size / num_groups
 

assimilation_code/modules/assimilation/obs_model_mod.f90

@@ -20,8 +20,8 @@ module obs_model_mod
                                  operator(/=), operator(>), operator(-),              &
                                  operator(/), operator(+), operator(<), operator(==), &
                                  operator(<=), operator(>=)
-use ensemble_manager_mod, only : get_ensemble_time, ensemble_type, map_task_to_pe, &
-                                 prepare_to_update_vars
+use ensemble_manager_mod, only : get_ensemble_time, ensemble_type, map_task_to_pe
+
 use mpi_utilities_mod,    only : my_task_id, task_sync, block_task, &
                                  sum_across_tasks, shell_execute, my_task_id
 use io_filenames_mod,     only : file_info_type
@@ -348,8 +348,6 @@ subroutine advance_state(ens_handle, ens_size, target_time, async, adv_ens_comma
    ! Ok, this task does need to advance something. 
    need_advance = 1
 
-   call prepare_to_update_vars(ens_handle)
-
    ! Increment number of ensemble member copies I have.
    my_num_state_copies = my_num_state_copies + 1
 

assimilation_code/modules/observations/forward_operator_mod.f90

@@ -30,8 +30,6 @@ module forward_operator_mod
 use obs_kind_mod,          only : assimilate_this_type_of_obs, evaluate_this_type_of_obs
 
 use ensemble_manager_mod,  only : ensemble_type, compute_copy_mean_var, &
-                                  prepare_to_read_from_vars,            &
-                                  prepare_to_write_to_vars,             &
                                   get_my_num_copies, copies_in_window,  &
                                   get_allow_transpose, all_vars_to_all_copies, &
                                   all_copies_to_all_vars, allocate_single_copy, &
@@ -127,11 +125,6 @@ subroutine get_obs_ens_distrib_state(ens_handle, obs_fwd_op_ens_handle, &
 istatus = 999123
 expected_obs = MISSING_R8
 
-! FIXME: these no longer do anything?
-! call prepare_to_write_to_vars(obs_fwd_op_ens_handle)
-! call prepare_to_write_to_vars(qc_ens_handle)
-! call prepare_to_read_from_vars(ens_handle)
-
 ! Set up access to the state
 call create_state_window(ens_handle, obs_fwd_op_ens_handle, qc_ens_handle)