[Dart-dev] [6144] DART/branches/development/models/lorenz_63/work/input.workshop.nml: (should have been part of commit 6143)

nancy at ucar.edu nancy at ucar.edu
Thu May 16 16:51:48 MDT 2013


Revision: 6144
Author:   thoar
Date:     2013-05-16 16:51:48 -0600 (Thu, 16 May 2013)
Log Message:
-----------
(should have been part of commit 6143)

The oned obs_diag program can now skip the spinup.
I removed the useless 'rat_cri' and 'input_qc_threshold' because
the DART QC values are available to do the same thing.

The input namelists also had some other variables removed,
like 'iskip_days' (replaced by init_skip_days, augmented by init_skip_seconds)
and 'obs_select' ... which was horrible to begin with. 'obs_select' required
you to specify the INTEGER depicting the obs TYPE you wanted to explore.

Modified Paths:
--------------
    DART/branches/development/models/lorenz_63/work/input.workshop.nml

-------------- next part --------------
Modified: DART/branches/development/models/lorenz_63/work/input.workshop.nml
===================================================================
--- DART/branches/development/models/lorenz_63/work/input.workshop.nml	2013-05-16 22:24:25 UTC (rev 6143)
+++ DART/branches/development/models/lorenz_63/work/input.workshop.nml	2013-05-16 22:51:48 UTC (rev 6144)
@@ -19,7 +19,7 @@
    output_forward_op_errors = .false.,
    print_every_nth_obs   = -1,
    silence               = .false.,
-  /
+   /
 
 &filter_nml
    async                    = 0,
@@ -38,7 +38,7 @@
    last_obs_days            = -1,
    last_obs_seconds         = -1,
    num_output_state_members = 20,
-   num_output_obs_members   = 0,
+   num_output_obs_members   = 20,
    output_interval          = 1,
    num_groups               = 1,
    input_qc_threshold       =  3.0,
@@ -49,7 +49,7 @@
    trace_execution          = .false.,
    silence                  = .false.,
 
-   inf_flavor                  = 2,                       0,
+   inf_flavor                  = 0,                       0,
    inf_initial_from_restart    = .false.,                 .false.,
    inf_sd_initial_from_restart = .false.,                 .false.,
    inf_output_restart          = .true.,                  .true.,
@@ -63,50 +63,57 @@
    inf_lower_bound             = 1.0,                     1.0,
    inf_upper_bound             = 100.0,                   1000000.0,
    inf_sd_lower_bound          = 0.6,                     0.0
-/
+   /
 
 &smoother_nml
    num_lags              = 0,
    start_from_restart    = .false.,
    output_restart        = .false.,
    restart_in_file_name  = 'smoother_ics',
-   restart_out_file_name = 'smoother_restart'  /
+   restart_out_file_name = 'smoother_restart'
+   /
 
 &ensemble_manager_nml
    single_restart_file_in  = .true.,
    single_restart_file_out = .true.,
-   perturbation_amplitude  = 0.2  /
+   perturbation_amplitude  = 0.2
+   /
 
 &assim_tools_nml
    filter_kind                     = 1,
-   cutoff                          = 0.001,
-   sort_obs_inc                    = .false.,
+   cutoff                          = 0.00001,
+   sort_obs_inc                    = .true.,
    spread_restoration              = .false.,
    sampling_error_correction       = .false.,
    adaptive_localization_threshold = -1,
    output_localization_diagnostics = .false.,
    localization_diagnostics_file   = 'localization_diagnostics',
-   print_every_nth_obs             = 0  /
+   print_every_nth_obs             = 0
+   /
 
 &cov_cutoff_nml
-   select_localization = 1  /
+   select_localization = 1
+   /
 
 &reg_factor_nml
    select_regression    = 1,
    input_reg_file       = "time_mean_reg",
    save_reg_diagnostics = .false.,
-   reg_diagnostics_file = "reg_diagnostics"  /
+   reg_diagnostics_file = "reg_diagnostics"
+   /
 
 &obs_sequence_nml
-   write_binary_obs_sequence = .false.  /
+   write_binary_obs_sequence = .false.
+   /
 
 &obs_kind_nml
-   assimilate_these_obs_types = 'RAW_STATE_VARIABLE'  /
+   assimilate_these_obs_types = 'RAW_STATE_VARIABLE'
+   /
 
 &assim_model_nml
    write_binary_restart_files = .false.,
    netCDF_large_file_support  = .false.
-  /
+   /
 
 &model_nml
    sigma  = 10.0, 
@@ -114,23 +121,25 @@
    b      = 2.6666666666667,
    deltat = 0.01,
    time_step_days = 0,
-   time_step_seconds = 3600  /
+   time_step_seconds = 3600
+   /
 
 &utilities_nml
    TERMLEVEL = 1,
    module_details = .false.,
    logfilename = 'dart_log.out',
    nmlfilename = 'dart_log.nml',
-   write_nml   = 'terminal'  /
+   write_nml   = 'terminal'
+   /
 
 &preprocess_nml
     input_obs_def_mod_file = '../../../obs_def/DEFAULT_obs_def_mod.F90', 
    output_obs_def_mod_file = '../../../obs_def/obs_def_mod.f90', 
    input_obs_kind_mod_file = '../../../obs_kind/DEFAULT_obs_kind_mod.F90', 
   output_obs_kind_mod_file = '../../../obs_kind/obs_kind_mod.f90', 
-               input_files = '../../../obs_def/obs_def_1d_state_mod.f90'  /
+               input_files = '../../../obs_def/obs_def_1d_state_mod.f90'
+   /
 
-
 &obs_sequence_tool_nml
    num_input_files   = 2,
    filename_seq      = 'obs_seq.one', 'obs_seq.two',
@@ -165,9 +174,7 @@
 #   copy_type          = 'RADIOSONDE_TEMPERATURE',
 #   min_copy           = 250.0,
 #   max_copy           = 300.0,
-#
 
-
 &restart_file_tool_nml
    input_file_name              = "filter_restart",
    output_file_name             = "filter_updated_restart",
@@ -184,17 +191,20 @@
    new_advance_days             = -1,
    new_advance_secs             = -1,
    gregorian_cal                = .false.
-/
+   /
 
 &obs_diag_nml
-   obs_sequence_name  = 'obs_seq.final',
-   iskip_days         = 0,
-   obs_select         = 1,
-   rat_cri            = 4.0,
-   input_qc_threshold = 3.0,
-   bin_width_seconds = 0,
-   lonlim1   = 0.0, 0.0, 0.5, -1.0,
-   lonlim2   = 1.0, 0.5, 1.5, -1.0,
-   reg_names = 'whole', 'yin', 'yang', 'bogus',
-   verbose   = .false.  /
+   obs_sequence_name     = 'obs_seq.final',
+   bin_width_days        = -1,
+   bin_width_seconds     = -1,
+   init_skip_days        = 0,
+   init_skip_seconds     = 0,
+   Nregions              = 3,
+   lonlim1               = 0.00, 0.00, 0.50, -1.0,
+   lonlim2               = 1.01, 0.50, 1.01, -1.0,
+   reg_names             = 'whole', 'yin', 'yang', 'bogus',
+   create_rank_histogram = .true.,
+   outliers_in_histogram = .true.,
+   verbose               = .false.
+   /
 


More information about the Dart-dev mailing list