[Dart-dev] DART/branches Revision: 11500

[dart at ucar.edu]

hendric at ucar.edu
2017-04-17 11:44:41 -0600 (Mon, 17 Apr 2017)
65
using wrong index to get copy name for writting augmented state




Modified: DART/branches/rma_output_files/assimilation_code/modules/assimilation/filter_mod.f90
===================================================================
--- DART/branches/rma_output_files/assimilation_code/modules/assimilation/filter_mod.f90	2017-04-14 22:03:29 UTC (rev 11499)
+++ DART/branches/rma_output_files/assimilation_code/modules/assimilation/filter_mod.f90	2017-04-17 17:44:41 UTC (rev 11500)
@@ -2203,6 +2203,15 @@
       call set_io_copy_flag(file_info, STAGE_COPIES(POSTINF_SD),    READ_COPY, has_units=.false.)
 endif
 
+! This is for single file augmented state mean and sd if requested
+if(single_file_in) then
+   if (output_mean) &
+     call set_io_copy_flag(file_info, STAGE_COPIES(ENS_MEAN),  WRITE_COPY, has_units=.false.)
+   
+   if (output_sd) &
+     call set_io_copy_flag(file_info, STAGE_COPIES(ENS_SD),    WRITE_COPY, has_units=.false.)
+endif
+
 end subroutine set_input_file_info
 
 !------------------------------------------------------------------

Modified: DART/branches/rma_output_files/assimilation_code/modules/io/single_file_io_mod.f90
===================================================================
--- DART/branches/rma_output_files/assimilation_code/modules/io/single_file_io_mod.f90	2017-04-14 22:03:29 UTC (rev 11499)
+++ DART/branches/rma_output_files/assimilation_code/modules/io/single_file_io_mod.f90	2017-04-17 17:44:41 UTC (rev 11500)
@@ -724,8 +724,7 @@
 
 do copy_index = ens_size+1, ens_handle%num_copies
    if ( file_info%stage_metadata%io_flag(copy_index) == WRITE_COPY ) then
-      ! check to see if the variable in in the files
-      copyname = get_copy_name(file_info,icopy)
+      copyname = get_copy_name(file_info,copy_index)
       call get_copy(map_task_to_pe(ens_handle, 0), ens_handle, copy_index, temp_ens)
       if(my_task_id() == 0) then
          call write_extra_attributes(ncFileID, TimeDimID, copyname)