[Dart-dev] DART/branches Revision: 10780

[dart at ucar.edu]

hendric at ucar.edu
2016-12-01 11:22:11 -0700 (Thu, 01 Dec 2016)
70
not storing spare copies at the moment. only writting out in
stages.




Modified: DART/branches/rma_fixed_filenames/filter/filter_mod.f90
===================================================================
--- DART/branches/rma_fixed_filenames/filter/filter_mod.f90	2016-12-01 18:13:35 UTC (rev 10779)
+++ DART/branches/rma_fixed_filenames/filter/filter_mod.f90	2016-12-01 18:22:11 UTC (rev 10780)
@@ -635,9 +635,9 @@
    !>@todo JPH : have namelist option to write out things in steps or all at the end
    !> in the sake of memory.
 
-   ! If needed, store copies(mean, sd, inf_mean, inf_sd) that would have
-   ! gone in Prior_Diag.nc and write them at the end.
-   call store_prior(state_ens_handle)
+   !#! ! If needed, store copies(mean, sd, inf_mean, inf_sd) that would have
+   !#! ! gone in Prior_Diag.nc and write them at the end.
+   !#! call store_prior(state_ens_handle)
 
    if ((output_interval > 0) .and. &
        (time_step_number / output_interval * output_interval == time_step_number)) then
@@ -770,9 +770,9 @@
    call trace_message('Before posterior state space diagnostics')
    call timestamp_message('Before posterior state space diagnostics')
 
-   ! If needed store the copies (inf_mean, inf_sd) that would have
-   ! gone in Posterior_Diag.nc and write them at the end
-   call store_posterior(state_ens_handle)
+   !#! ! If needed store the copies (inf_mean, inf_sd) that would have
+   !#! ! gone in Posterior_Diag.nc and write them at the end
+   !#! call store_posterior(state_ens_handle)
 
    !>@todo FIXME JPH: option to write copies after inflation and before
    !> posterior inflation.
@@ -1844,37 +1844,30 @@
 
 !>@todo FIXME TJH,NSC ... assign strings/copy metadata at this point
 
-if (get_output_preassim()) then
-   ! Not stopping to write prior_members so keep these Prior copies
-   ! as extra copies and write them and the end.
-   SPARE_PRIOR_MEAN       = ens_size + 7
-   SPARE_PRIOR_SPREAD     = ens_size + 8
-   SPARE_PRIOR_INF_MEAN   = ens_size + 9
-   SPARE_PRIOR_INF_SPREAD = ens_size + 10
-   !#! call set_has_units(SPARE_PRIOR_MEAN,       .false.)
-   !#! call set_has_units(SPARE_PRIOR_SPREAD,     .false.)
-   !#! call set_has_units(SPARE_PRIOR_INF_MEAN,   .false.)
-   !#! call set_has_units(SPARE_PRIOR_INF_SPREAD, .false.)
-   num_extras = num_extras + 4
-   if(get_output_postassim()) then
-      ! need to store posterior inflation mean and inflation spread since
-      ! these are overwritten in filter_assim(inflate_only=.true.)
-      SPARE_POST_INF_MEAN    = ens_size + 11
-      SPARE_POST_INF_SPREAD  = ens_size + 12
-      !#! call set_has_units(SPARE_POST_INF_MEAN,   .false.)
-      !#! call set_has_units(SPARE_POST_INF_SPREAD, .false.)
-      num_extras = num_extras + 2
-   endif
-elseif (get_output_postassim()) then
-   ! Prior members and extra copies are written out prior to filter_assim.
-   ! Need to store posterior inflation mean and inflation spread since
-   ! these are overwritten in filter_assim(inflate_only=.true.)
-   SPARE_POST_INF_MEAN    = ens_size + 7
-   SPARE_POST_INF_SPREAD  = ens_size + 8
-   !#! call set_has_units(SPARE_POST_INF_MEAN,   .false.)
-   !#! call set_has_units(SPARE_POST_INF_SPREAD, .false.)
-   num_extras = num_extras + 2
-endif
+!#! !>@todo JPH : not using extra copies at the moment.
+!#! if (get_output_preassim()) then
+!#!    ! Not stopping to write prior_members so keep these Prior copies
+!#!    ! as extra copies and write them and the end.
+!#!    SPARE_PRIOR_MEAN       = ens_size + 7
+!#!    SPARE_PRIOR_SPREAD     = ens_size + 8
+!#!    SPARE_PRIOR_INF_MEAN   = ens_size + 9
+!#!    SPARE_PRIOR_INF_SPREAD = ens_size + 10
+!#!    num_extras = num_extras + 4
+!#!    if(get_output_postassim()) then
+!#!       ! need to store posterior inflation mean and inflation spread since
+!#!       ! these are overwritten in filter_assim(inflate_only=.true.)
+!#!       SPARE_POST_INF_MEAN    = ens_size + 11
+!#!       SPARE_POST_INF_SPREAD  = ens_size + 12
+!#!       num_extras = num_extras + 2
+!#!    endif
+!#! elseif (get_output_postassim()) then
+!#!    ! Prior members and extra copies are written out prior to filter_assim.
+!#!    ! Need to store posterior inflation mean and inflation spread since
+!#!    ! these are overwritten in filter_assim(inflate_only=.true.)
+!#!    SPARE_POST_INF_MEAN    = ens_size + 7
+!#!    SPARE_POST_INF_SPREAD  = ens_size + 8
+!#!    num_extras = num_extras + 2
+!#! endif
 
 end subroutine set_state_copies
 
@@ -1901,13 +1894,14 @@
 call set_file_metadata(file_info, stage_name, PRIOR_INF_SD_COPY, 'priorinf_sd', 'prior inflation standard deviation')
 call set_file_metadata(file_info, stage_name, POST_INF_COPY,     'postinf_mean', 'posterior inflation mean')
 call set_file_metadata(file_info, stage_name, POST_INF_SD_COPY,  'postinf_sd', 'posterior inflation standard deviation')
-call set_file_metadata(file_info, stage_name, SPARE_PRIOR_MEAN,       'prior_mean',     'prior ensemble mean')
-call set_file_metadata(file_info, stage_name, SPARE_PRIOR_SPREAD,     'prior_sd',       'prior ensemble standard deviation')
-!>@todo JPH : do we still need these copies??
-call set_file_metadata(file_info, stage_name, SPARE_PRIOR_INF_MEAN,   'priorinf_mean',  'prior inflation mean')
-call set_file_metadata(file_info, stage_name, SPARE_PRIOR_INF_SPREAD, 'priorinf_sd',    'prior inflation standard deviation')