[Dart-dev] [3311] DART/trunk: Add a new namelist item 'inf_damping'
which allows you to
nancy at subversion.ucar.edu
nancy at subversion.ucar.edu
Mon Apr 14 15:43:46 MDT 2008
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080414/db774b1b/attachment-0001.html
-------------- next part --------------
Modified: DART/trunk/filter/filter.dopplerfold.f90
===================================================================
--- DART/trunk/filter/filter.dopplerfold.f90 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/filter/filter.dopplerfold.f90 2008-04-14 21:43:46 UTC (rev 3311)
@@ -21,7 +21,8 @@
static_init_obs_sequence, destroy_obs, read_obs_seq_header, &
set_qc_meta_data, get_expected_obs, get_first_obs, &
get_obs_time_range, delete_obs_from_seq, delete_seq_head, &
- delete_seq_tail, replace_obs_values, replace_qc
+ delete_seq_tail, replace_obs_values, replace_qc, &
+ destroy_obs_sequence
use obs_def_mod, only : obs_def_type, get_obs_def_error_variance, get_obs_def_time
use time_manager_mod, only : time_type, get_time, set_time, operator(/=), operator(>), &
operator(-)
@@ -118,6 +119,7 @@
inf_diag_file_name(2) = 'not_initialized'
real(r8) :: inf_initial(2) = 1.0_r8
real(r8) :: inf_sd_initial(2) = 0.0_r8
+real(r8) :: inf_damping(2) = 1.0_r8
real(r8) :: inf_lower_bound(2) = 1.0_r8
real(r8) :: inf_upper_bound(2) = 1000000.0_r8
real(r8) :: inf_sd_lower_bound(2) = 0.0_r8
@@ -131,7 +133,7 @@
num_output_obs_members, output_interval, num_groups, outlier_threshold, &
input_qc_threshold, output_forward_op_errors, output_timestamps, &
inf_flavor, inf_initial_from_restart, inf_sd_initial_from_restart, &
- inf_output_restart, inf_deterministic, inf_in_file_name, &
+ inf_output_restart, inf_deterministic, inf_in_file_name, inf_damping, &
inf_out_file_name, inf_diag_file_name, inf_initial, inf_sd_initial, &
inf_lower_bound, inf_upper_bound, inf_sd_lower_bound, output_inflation
@@ -209,6 +211,10 @@
write(msgstring, *) 'inf_flavor=', inf_flavor(i), ' Must be 0, 1, 2, 3 '
call error_handler(E_ERR,'filter_main', msgstring, source, revision, revdate)
endif
+ if(inf_damping(i) < 0.0_r8 .or. inf_damping(i) > 1.0_r8) then
+ write(msgstring, *) 'inf_damping=', inf_damping(i), ' Must be 0.0 <= d <= 1.0'
+ call error_handler(E_ERR,'filter_main', msgstring, source, revision, revdate)
+ endif
end do
! Observation space inflation for posterior not currently supported
@@ -281,17 +287,17 @@
! Initialize the output sequences and state files and set their meta data
if(my_task_id() == 0) then
call filter_generate_copy_meta_data(seq, prior_inflate, &
- PriorStateUnit, PosteriorStateUnit, in_obs_copy, output_state_mean_index, &
- output_state_spread_index, prior_obs_mean_index, posterior_obs_mean_index, &
- prior_obs_spread_index, posterior_obs_spread_index)
+ PriorStateUnit, PosteriorStateUnit, in_obs_copy, output_state_mean_index, &
+ output_state_spread_index, prior_obs_mean_index, posterior_obs_mean_index, &
+ prior_obs_spread_index, posterior_obs_spread_index)
if(ds) call smoother_gen_copy_meta_data(num_output_state_members, output_inflation)
else
- output_state_mean_index = 0
- output_state_spread_index = 0
- prior_obs_mean_index = 0
- posterior_obs_mean_index = 0
- prior_obs_spread_index = 0
- posterior_obs_spread_index = 0
+ output_state_mean_index = 0
+ output_state_spread_index = 0
+ prior_obs_mean_index = 0
+ posterior_obs_mean_index = 0
+ prior_obs_spread_index = 0
+ posterior_obs_spread_index = 0
endif
! Need to find first obs with appropriate time, delete all earlier ones
@@ -381,6 +387,11 @@
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
if(do_single_ss_inflate(prior_inflate) .or. do_varying_ss_inflate(prior_inflate)) then
+ if (inf_damping(1) /= 1.0_r8) then
+ ens_handle%copies(PRIOR_INF_COPY, :) = 1.0_r8 + &
+ inf_damping(1) * (ens_handle%copies(PRIOR_INF_COPY, :) - 1.0_r8)
+ endif
+
call filter_ensemble_inflate(ens_handle, PRIOR_INF_COPY, prior_inflate, ENS_MEAN_COPY)
! Recompute the the mean and spread as required for diagnostics
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
@@ -455,6 +466,11 @@
!-------- Test of posterior inflate ----------------
if(do_single_ss_inflate(post_inflate) .or. do_varying_ss_inflate(post_inflate)) then
+ if (inf_damping(2) /= 1.0_r8) then
+ ens_handle%copies(POST_INF_COPY, :) = 1.0_r8 + &
+ inf_damping(2) * (ens_handle%copies(POST_INF_COPY, :) - 1.0_r8)
+ endif
+
call filter_ensemble_inflate(ens_handle, POST_INF_COPY, post_inflate, ENS_MEAN_COPY)
! Recompute the mean or the mean and spread as required for diagnostics
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
@@ -504,34 +520,34 @@
! sd < 0, then bypass this entire code block altogether for speed.
if ((inf_sd_initial(2) >= 0.0_r8) .or. inf_sd_initial_from_restart(2)) then
- ! Ship the ensemble mean to the model; some models need this for computing distances
- ! Who stores the ensemble mean copy
- call get_copy_owner_index(ENS_MEAN_COPY, mean_owner, mean_owners_index)
- ! Broadcast it to everybody else
- if(my_task_id() == mean_owner) then
- ens_mean = ens_handle%vars(:, mean_owners_index)
- call broadcast_send(mean_owner, ens_mean)
- else
- call broadcast_recv(mean_owner, ens_mean)
- endif
-
- ! Now send the mean to the model in case it's needed
- call ens_mean_for_model(ens_mean)
+ ! Ship the ensemble mean to the model; some models need this for computing distances
+ ! Who stores the ensemble mean copy
+ call get_copy_owner_index(ENS_MEAN_COPY, mean_owner, mean_owners_index)
+ ! Broadcast it to everybody else
+ if(my_task_id() == mean_owner) then
+ ens_mean = ens_handle%vars(:, mean_owners_index)
+ call broadcast_send(mean_owner, ens_mean)
+ else
+ call broadcast_recv(mean_owner, ens_mean)
+ endif
- ! Need obs to be copy complete for assimilation: IS NEXT LINE REQUIRED???
- call all_vars_to_all_copies(obs_ens_handle)
- if (output_timestamps) then
- if (do_output()) call timestamp("Before posterior inflation", pos='debug')
- endif
- call filter_assim(ens_handle, obs_ens_handle, seq, keys, ens_size, num_groups, &
- obs_val_index, post_inflate, ENS_MEAN_COPY, ENS_SD_COPY, &
- POST_INF_COPY, POST_INF_SD_COPY, OBS_KEY_COPY, OBS_GLOBAL_QC_COPY, &
- OBS_PRIOR_MEAN_START, OBS_PRIOR_MEAN_END, OBS_PRIOR_VAR_START, &
- OBS_PRIOR_VAR_END, inflate_only = .true.)
- if (output_timestamps) then
- if (do_output()) call timestamp("After posterior inflation", pos='debug')
- endif
- call all_copies_to_all_vars(ens_handle)
+ ! Now send the mean to the model in case it's needed
+ call ens_mean_for_model(ens_mean)
+
+ ! Need obs to be copy complete for assimilation: IS NEXT LINE REQUIRED???
+ call all_vars_to_all_copies(obs_ens_handle)
+ if (output_timestamps) then
+ if (do_output()) call timestamp("Before posterior inflation", pos='debug')
+ endif
+ call filter_assim(ens_handle, obs_ens_handle, seq, keys, ens_size, num_groups, &
+ obs_val_index, post_inflate, ENS_MEAN_COPY, ENS_SD_COPY, &
+ POST_INF_COPY, POST_INF_SD_COPY, OBS_KEY_COPY, OBS_GLOBAL_QC_COPY, &
+ OBS_PRIOR_MEAN_START, OBS_PRIOR_MEAN_END, OBS_PRIOR_VAR_START, &
+ OBS_PRIOR_VAR_END, inflate_only = .true.)
+ if (output_timestamps) then
+ if (do_output()) call timestamp("After posterior inflation", pos='debug')
+ endif
+ call all_copies_to_all_vars(ens_handle)
endif ! sd >= 0 or sd from restart file
endif ! if doing state space posterior inflate
@@ -582,8 +598,9 @@
! Master task must close the log file
if(my_task_id() == 0) call timestamp(source,revision,revdate,'end')
-! Free up the observation kind
+! Free up the observation kind and obs sequence
call destroy_obs(observation)
+call destroy_obs_sequence(seq)
if(ds) call smoother_end()
Modified: DART/trunk/filter/filter.f90
===================================================================
--- DART/trunk/filter/filter.f90 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/filter/filter.f90 2008-04-14 21:43:46 UTC (rev 3311)
@@ -119,6 +119,7 @@
inf_diag_file_name(2) = 'not_initialized'
real(r8) :: inf_initial(2) = 1.0_r8
real(r8) :: inf_sd_initial(2) = 0.0_r8
+real(r8) :: inf_damping(2) = 1.0_r8
real(r8) :: inf_lower_bound(2) = 1.0_r8
real(r8) :: inf_upper_bound(2) = 1000000.0_r8
real(r8) :: inf_sd_lower_bound(2) = 0.0_r8
@@ -132,7 +133,7 @@
num_output_obs_members, output_interval, num_groups, outlier_threshold, &
input_qc_threshold, output_forward_op_errors, output_timestamps, &
inf_flavor, inf_initial_from_restart, inf_sd_initial_from_restart, &
- inf_output_restart, inf_deterministic, inf_in_file_name, &
+ inf_output_restart, inf_deterministic, inf_in_file_name, inf_damping, &
inf_out_file_name, inf_diag_file_name, inf_initial, inf_sd_initial, &
inf_lower_bound, inf_upper_bound, inf_sd_lower_bound, output_inflation
@@ -206,6 +207,10 @@
write(msgstring, *) 'inf_flavor=', inf_flavor(i), ' Must be 0, 1, 2, 3 '
call error_handler(E_ERR,'filter_main', msgstring, source, revision, revdate)
endif
+ if(inf_damping(i) < 0.0_r8 .or. inf_damping(i) > 1.0_r8) then
+ write(msgstring, *) 'inf_damping=', inf_damping(i), ' Must be 0.0 <= d <= 1.0'
+ call error_handler(E_ERR,'filter_main', msgstring, source, revision, revdate)
+ endif
end do
! Observation space inflation for posterior not currently supported
@@ -378,6 +383,11 @@
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
if(do_single_ss_inflate(prior_inflate) .or. do_varying_ss_inflate(prior_inflate)) then
+ if (inf_damping(1) /= 1.0_r8) then
+ ens_handle%copies(PRIOR_INF_COPY, :) = 1.0_r8 + &
+ inf_damping(1) * (ens_handle%copies(PRIOR_INF_COPY, :) - 1.0_r8)
+ endif
+
call filter_ensemble_inflate(ens_handle, PRIOR_INF_COPY, prior_inflate, ENS_MEAN_COPY)
! Recompute the the mean and spread as required for diagnostics
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
@@ -451,6 +461,11 @@
!-------- Test of posterior inflate ----------------
if(do_single_ss_inflate(post_inflate) .or. do_varying_ss_inflate(post_inflate)) then
+ if (inf_damping(2) /= 1.0_r8) then
+ ens_handle%copies(POST_INF_COPY, :) = 1.0_r8 + &
+ inf_damping(2) * (ens_handle%copies(POST_INF_COPY, :) - 1.0_r8)
+ endif
+
call filter_ensemble_inflate(ens_handle, POST_INF_COPY, post_inflate, ENS_MEAN_COPY)
! Recompute the mean or the mean and spread as required for diagnostics
call compute_copy_mean_sd(ens_handle, 1, ens_size, ENS_MEAN_COPY, ENS_SD_COPY)
Modified: DART/trunk/filter/filter.html
===================================================================
--- DART/trunk/filter/filter.html 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/filter/filter.html 2008-04-14 21:43:46 UTC (rev 3311)
@@ -175,6 +175,7 @@
inf_diag_file_name,
inf_initial,
inf_sd_initial,
+ inf_damping,
inf_lower_bound,
inf_upper_bound,
inf_sd_lower_bound
@@ -376,6 +377,17 @@
values will adapt through time, so they are time-varying.
Default: (/0.0_r8, 0.0_r8/)</TD></TR>
+ <TR><!--contents--><TD valign=top>inf_damping</TD>
+ <!-- type --><TD valign=top>real(r8) (len=2)</TD>
+ <!--descript--><TD>Damping factor for inflation mean values.
+ The difference between the current inflation value and 1.0
+ is multiplied by this factor before the next assimilation cycle.
+ The value should be between 0.0 and 1.0.
+ 0.0 turns all inflation off (by clamping the inflation value to 1.0).
+ 1.0 turns inflation damping off (leaving the original inflation
+ value unchanged).
+ Default: (/1.0_r8, 1.0_r8/)</TD></TR>
+
<TR><!--contents--><TD valign=top>inf_lower_bound</TD>
<!-- type --><TD valign=top>real(r8) (len=2)</TD>
<!--descript--><TD>Lower bound for inflation value.
Modified: DART/trunk/filter/filter.nml
===================================================================
--- DART/trunk/filter/filter.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/filter/filter.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -34,6 +34,7 @@
inf_diag_file_name = 'not_initialized', 'not_initialized',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/9var/work/input.nml
===================================================================
--- DART/trunk/models/9var/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/9var/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/MITgcm_annulus/work/input.nml
===================================================================
--- DART/trunk/models/MITgcm_annulus/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/MITgcm_annulus/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/MITgcm_ocean/work/input.nml
===================================================================
--- DART/trunk/models/MITgcm_ocean/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/MITgcm_ocean/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/PBL_1d/work/input.nml
===================================================================
--- DART/trunk/models/PBL_1d/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/PBL_1d/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/am2/work/input.nml
===================================================================
--- DART/trunk/models/am2/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/am2/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -49,6 +49,7 @@
inf_initial = 1.00, 1.00,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
+ inf_damping = 1.0, 1.0,
inf_sd_initial = 0.1, 0.1,
inf_sd_lower_bound = 0.0, 0.0
/
Modified: DART/trunk/models/bgrid_solo/work/input.nml
===================================================================
--- DART/trunk/models/bgrid_solo/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/bgrid_solo/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/cam/work/input.nml
===================================================================
--- DART/trunk/models/cam/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/cam/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -32,6 +32,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.1, 0.1,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/forced_lorenz_96/work/input.nml
===================================================================
--- DART/trunk/models/forced_lorenz_96/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/forced_lorenz_96/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/ikeda/work/input.nml
===================================================================
--- DART/trunk/models/ikeda/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/ikeda/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/lorenz_04/work/input.nml
===================================================================
--- DART/trunk/models/lorenz_04/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/lorenz_04/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/lorenz_63/work/input.nml
===================================================================
--- DART/trunk/models/lorenz_63/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/lorenz_63/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/lorenz_84/work/input.nml
===================================================================
--- DART/trunk/models/lorenz_84/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/lorenz_84/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/lorenz_96/work/input.nml
===================================================================
--- DART/trunk/models/lorenz_96/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/lorenz_96/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/lorenz_96_2scale/work/input.nml
===================================================================
--- DART/trunk/models/lorenz_96_2scale/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/lorenz_96_2scale/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/null_model/work/input.nml
===================================================================
--- DART/trunk/models/null_model/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/null_model/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/pe2lyr/work/input.nml
===================================================================
--- DART/trunk/models/pe2lyr/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/pe2lyr/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/rose/work/input.nml
===================================================================
--- DART/trunk/models/rose/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/rose/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/simple_advection/work/input.nml
===================================================================
--- DART/trunk/models/simple_advection/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/simple_advection/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.2,
inf_sd_initial = 0.2, 0.5,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.2, 0.5
Modified: DART/trunk/models/template/work/input.nml
===================================================================
--- DART/trunk/models/template/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/template/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.0, 1.0,
inf_sd_initial = 0.0, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
Modified: DART/trunk/models/wrf/work/input.nml
===================================================================
--- DART/trunk/models/wrf/work/input.nml 2008-04-14 16:03:08 UTC (rev 3310)
+++ DART/trunk/models/wrf/work/input.nml 2008-04-14 21:43:46 UTC (rev 3311)
@@ -51,6 +51,7 @@
inf_diag_file_name = 'prior_inflate_diag', 'post_inflate_diag',
inf_initial = 1.02, 1.00,
inf_sd_initial = 0.1, 0.0,
+ inf_damping = 1.0, 1.0,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
inf_sd_lower_bound = 0.0, 0.0
More information about the Dart-dev
mailing list