[Dart-dev] DART/branches Revision: 12605

[dart at ucar.edu]

hendric at ucar.edu
2018-05-29 11:07:59 -0600 (Tue, 29 May 2018)
30

RMSE now working with MPI. 




Modified: DART/branches/rma_closest_member_tool/assimilation_code/programs/closest_member_tool/closest_member_tool.f90
===================================================================
--- DART/branches/rma_closest_member_tool/assimilation_code/programs/closest_member_tool/closest_member_tool.f90	2018-05-29 14:32:14 UTC (rev 12604)
+++ DART/branches/rma_closest_member_tool/assimilation_code/programs/closest_member_tool/closest_member_tool.f90	2018-05-29 17:07:59 UTC (rev 12605)
@@ -19,8 +19,8 @@
 
 use types_mod,            only : r8, i8, obstypelength, MAX_NUM_DOMS, MAX_FILES
 
-use time_manager_mod,     only : time_type, set_time_missing,               &
-                                 operator(/=), print_time
+use time_manager_mod,     only : time_type, set_time_missing, operator(/=), &
+                                 print_time
  
 use utilities_mod,        only : register_module, find_namelist_in_file,        &
                                  error_handler, nmlfileunit, E_MSG, E_ERR,      &
@@ -34,21 +34,21 @@
 
 use  sort_mod,            only : index_sort
 
-use assim_model_mod,      only : static_init_assim_model, get_model_size,   &
+use assim_model_mod,      only : static_init_assim_model, get_model_size, &
                                  get_state_meta_data
 
 use state_vector_io_mod,  only : read_state, write_state
 
-use io_filenames_mod,     only : file_info_type, io_filenames_init, &
-                                 set_io_copy_flag, set_file_metadata, &
+use io_filenames_mod,     only : file_info_type, io_filenames_init,        &
+                                 set_io_copy_flag, set_file_metadata,      &
                                  set_member_file_metadata, file_info_dump, &
-                                 stage_metadata_type, get_stage_metadata, &
+                                 stage_metadata_type, get_stage_metadata,  &
                                  get_restart_filename, READ_COPY
 
 use state_structure_mod,  only : get_num_domains
 
-use mpi_utilities_mod,    only : initialize_mpi_utilities, task_count,     &
-                                 finalize_mpi_utilities, my_task_id, &
+use mpi_utilities_mod,    only : initialize_mpi_utilities, task_count, &
+                                 finalize_mpi_utilities, my_task_id,   &
                                  send_sum_to
 
 use ensemble_manager_mod, only : ensemble_type, init_ensemble_manager, compute_copy_mean, &
@@ -65,15 +65,15 @@
 character(len=128), parameter :: revdate  = "$Date$"
 
 integer               :: iunit, io, ens, i, j, qtyindex
+integer               :: num_qtys, stype
 integer(i8)           :: ii, model_size
 integer, allocatable  :: index_list(:)
+integer, parameter    :: max_list_len = 500
 character(len = 128)  :: msgstring, msgstring1
 logical               :: allqtys, done
 logical, allocatable  :: useqty(:), useindex(:)
 type(location_type)   :: loc
-integer               :: num_qtys, stype
 type(time_type)       :: mean_time, member_time
-integer, parameter    :: max_list_len = 500
 type(file_info_type)  :: ens_file_info
 
 
@@ -86,13 +86,13 @@
 !----------------------------------------------------------------
 ! These variables are namelist-controllable.
 !
-integer              :: ens_size               = 1
-integer              :: difference_method      = 4
-logical              :: single_restart_file_in = .false.
+integer                        :: ens_size               = 1
+integer                        :: difference_method      = 4
+logical                        :: single_restart_file_in = .false.
 character(len = obstypelength) :: use_only_qtys(max_list_len) = ''
-character(len=256) :: input_restart_file_list(MAX_NUM_DOMS)  = 'null'
-character(len=256) :: mean_restart_file_list(MAX_NUM_DOMS)  = 'null'
-character(len=256) :: output_file_name
+character(len=256)             :: input_restart_file_list(MAX_NUM_DOMS) = 'null'
+character(len=256)             :: mean_restart_file_list(MAX_NUM_DOMS)  = 'null'
+character(len=256)             :: output_file_name
 
 !----------------------------------------------------------------
 ! different methods to compute 'distance' from mean:
@@ -112,23 +112,22 @@
    single_restart_file_in,       &
    difference_method,            &
    use_only_qtys         
+
 type(ensemble_type)             :: ens_handle
+type(stage_metadata_type)       :: input_restart_files
 character(len=256), allocatable :: file_array_input(:,:)
 character(len=256)              :: input_state_files(MAX_FILES) = '' 
+character(len=256)              :: my_base, my_desc
+integer(i8), allocatable        :: vars_array(:)
+integer(i8)                     :: owners_index
 integer                         :: num_domains, idom, imem
-type(stage_metadata_type)       :: input_restart_files
-character(len=256)              :: my_base, my_desc
 integer                         :: ENS_MEAN_COPY 
+integer                         :: copies, my_num_vars
 
-real(r8),    allocatable :: adiff(:), local_diffs(:), total_diff(:)
-integer(i8), allocatable :: vars_array(:)