[Dart-dev] [4113] DART/trunk: Extended the 3D version of obs_diag to be able to read a list of

nancy at ucar.edu nancy at ucar.edu
Thu Oct 22 14:33:36 MDT 2009


Revision: 4113
Author:   thoar
Date:     2009-10-22 14:33:36 -0600 (Thu, 22 Oct 2009)
Log Message:
-----------
Extended the 3D version of obs_diag to be able to read a list of 
filenames as well as the 'status quo' of a single filename and
a filename that is part of a sequence. There is a new namelist
parameter 'obs_sequence_list' that is the name of a file containing
the list of filenames to process. 'obs_sequence_file' still exists,
and defaults to 'obs_seq.final'. One of 'obs_sequence_file' or
'obs_sequence_list' MUST be an empty string.

Changed the documentation and the default input.nml settings
for those models that use the threed_ version of obs_diag.

Modified Paths:
--------------
    DART/trunk/diagnostics/threed_sphere/obs_diag.f90
    DART/trunk/diagnostics/threed_sphere/obs_diag.html
    DART/trunk/diagnostics/threed_sphere/obs_diag.nml
    DART/trunk/models/MITgcm_annulus/work/input.nml
    DART/trunk/models/MITgcm_ocean/work/input.nml
    DART/trunk/models/PBL_1d/work/input.nml
    DART/trunk/models/POP/work/input.nml
    DART/trunk/models/am2/work/input.nml
    DART/trunk/models/bgrid_solo/work/input.nml
    DART/trunk/models/cam/work/input.nml
    DART/trunk/models/pe2lyr/work/input.nml
    DART/trunk/models/rose/work/input.nml
    DART/trunk/models/wrf/regression/CONUS-V2/input.nml
    DART/trunk/models/wrf/regression/Radar/input.nml
    DART/trunk/models/wrf/work/input.nml
    DART/trunk/observations/utilities/threed_sphere/input.nml

-------------- next part --------------
Modified: DART/trunk/diagnostics/threed_sphere/obs_diag.f90
===================================================================
--- DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2009-10-22 20:33:36 UTC (rev 4113)
@@ -51,11 +51,12 @@
                              operator(>), operator(<), operator(/), &
                              operator(/=), operator(<=)
 use    utilities_mod, only : open_file, close_file, register_module, &
-                             file_exist, error_handler, E_ERR, E_WARN, E_MSG, &
-                             initialize_utilities, logfileunit, nmlfileunit,  &
-                             find_namelist_in_file, check_namelist_read,      &
-                             nc_check, do_nml_file, do_nml_term, timestamp,   &
-                             next_file, find_textfile_dims, file_to_text
+                             file_exist, error_handler, E_ERR, E_WARN, E_MSG,  &
+                             initialize_utilities, logfileunit, nmlfileunit,   &
+                             find_namelist_in_file, check_namelist_read,       &
+                             nc_check, do_nml_file, do_nml_term, timestamp,    &
+                             next_file, get_next_filename, find_textfile_dims, &
+                             file_to_text
 use         sort_mod, only : sort
 
 use typeSizes
@@ -160,6 +161,7 @@
 !-----------------------------------------------------------------------
 
 character(len = 129) :: obs_sequence_name = 'obs_seq.final'
+character(len = 129) :: obs_sequence_list = ''
 integer, dimension(6) :: first_bin_center = (/ 2003, 1, 1, 0, 0, 0 /)
 integer, dimension(6) :: last_bin_center  = (/ 2003, 1, 2, 0, 0, 0 /)
 integer, dimension(6) :: bin_separation   = (/    0, 0, 0, 6, 0, 0 /)
@@ -183,7 +185,8 @@
 logical :: print_obs_locations   = .false.
 logical :: verbose               = .false.
 
-namelist /obs_diag_nml/ obs_sequence_name, first_bin_center, last_bin_center, &
+namelist /obs_diag_nml/ obs_sequence_name, obs_sequence_list, &
+                       first_bin_center, last_bin_center, &
                        bin_separation, bin_width, time_to_skip, max_num_bins, &
                        plevel, hlevel, mlevel, rat_cri, input_qc_threshold, &
                        Nregions, lonlim1, lonlim2, latlim1, latlim2, &
@@ -319,6 +322,13 @@
 if (do_nml_file()) write(nmlfileunit, nml=obs_diag_nml)
 if (do_nml_term()) write(    *      , nml=obs_diag_nml)
 
+if ((obs_sequence_name /= '') .and. (obs_sequence_list /= '')) then
+   write(msgstring,*)'specify "obs_sequence_name" or "obs_sequence_list"'
+   call error_handler(E_MSG, 'obs_diag', msgstring, source, revision, revdate)
+   write(msgstring,*)'set other to an empty string ... i.e. ""'
+   call error_handler(E_ERR, 'obs_diag', msgstring, source, revision, revdate)
+endif
+
 !----------------------------------------------------------------------
 ! Now that we have input, do some checking and setup
 !----------------------------------------------------------------------
@@ -509,11 +519,20 @@
 ! If the sequence is completely before the time period of interest, we skip.
 ! If the sequence is completely past the time period of interest, we stop.
 
-ObsFileLoop : do ifile=1, Nepochs*4
+ObsFileLoop : do ifile=1, 1000
 !-----------------------------------------------------------------------
 
-   obs_seq_in_file_name = next_file(obs_sequence_name,ifile)
+   if (obs_sequence_list == '') then ! try to increment filename
 
+      obs_seq_in_file_name = next_file(obs_sequence_name,ifile)
+
+   else
+
+      obs_seq_in_file_name = get_next_filename(obs_sequence_list,ifile)
+      if (obs_seq_in_file_name == '') exit ObsFileLoop
+
+   endif
+
    if ( file_exist(trim(obs_seq_in_file_name)) ) then
       write(msgstring,*)'opening ', trim(obs_seq_in_file_name)
       call error_handler(E_MSG,'obs_diag',msgstring,source,revision,revdate)
@@ -2723,6 +2742,8 @@
              'WriteNetCDF', 'mlevel:long_name')
    call nc_check(nf90_put_att(ncid, MlevelVarID, 'units',     'nondimensional'), &
              'WriteNetCDF', 'mlevel:units')
+   call nc_check(nf90_put_att(ncid, MlevelVarID, 'axis',     'Z'), &
+             'WriteNetCDF', 'mlevel:axis')
    call nc_check(nf90_put_att(ncid, MlevelVarID, 'valid_range', (/1,Nmlevels/)), &
              'WriteNetCDF', 'mlevel:valid_range')
 
@@ -2734,6 +2755,8 @@
              'WriteNetCDF', 'plevel:long_name')
    call nc_check(nf90_put_att(ncid, PlevelVarID, 'units',     'hPa'), &
              'WriteNetCDF', 'plevel:units')
+   call nc_check(nf90_put_att(ncid, PlevelVarID, 'axis',     'Z'), &
+             'WriteNetCDF', 'plevel:axis')
    call nc_check(nf90_put_att(ncid, PlevelVarID, 'valid_range', &
       (/ minval(plevel(1:Nplevels)), maxval(plevel(1:Nplevels)) /)), &
              'WriteNetCDF', 'plevel:valid_range')
@@ -2746,6 +2769,8 @@
              'WriteNetCDF', 'plevel_edges:long_name')
    call nc_check(nf90_put_att(ncid, PlevIntVarID, 'units',     'hPa'), &
              'WriteNetCDF', 'plevel_edges:units')
+   call nc_check(nf90_put_att(ncid, PlevIntVarID, 'axis',     'Z'), &
+             'WriteNetCDF', 'plevel_edges:axis')
    call nc_check(nf90_put_att(ncid, PlevIntVarID, 'valid_range', &
       (/ minval(plevel_edges(1:Nplevels+1)), maxval(plevel_edges(1:Nplevels+1)) /)), &
              'WriteNetCDF', 'plevel_edges:valid_range')
@@ -2758,6 +2783,8 @@
              'WriteNetCDF', 'hlevel:long_name')
    call nc_check(nf90_put_att(ncid, HlevelVarID, 'units',     'm'), &
              'WriteNetCDF', 'hlevel:units')
+   call nc_check(nf90_put_att(ncid, HlevelVarID, 'axis',     'Z'), &
+             'WriteNetCDF', 'hlevel:axis')
    call nc_check(nf90_put_att(ncid, HlevelVarID, 'valid_range', &
       (/ minval(hlevel(1:Nhlevels)), maxval(hlevel(1:Nhlevels)) /)), &
              'WriteNetCDF', 'hlevel:valid_range')
@@ -2770,6 +2797,8 @@
              'WriteNetCDF', 'hlevel_edges:long_name')
    call nc_check(nf90_put_att(ncid, HlevIntVarID, 'units',     'm'), &
              'WriteNetCDF', 'hlevel_edges:units')
+   call nc_check(nf90_put_att(ncid, HlevIntVarID, 'axis',     'Z'), &
+             'WriteNetCDF', 'hlevel_edges:axis')
    call nc_check(nf90_put_att(ncid, HlevIntVarID, 'valid_range', &
       (/ minval(hlevel_edges(1:Nhlevels+1)), maxval(hlevel_edges(1:Nhlevels+1)) /)), &
              'WriteNetCDF', 'hlevel_edges:valid_range')

Modified: DART/trunk/diagnostics/threed_sphere/obs_diag.html
===================================================================
--- DART/trunk/diagnostics/threed_sphere/obs_diag.html	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/diagnostics/threed_sphere/obs_diag.html	2009-10-22 20:33:36 UTC (rev 4113)
@@ -281,7 +281,7 @@
    ampersand '&#38;' and terminating with a slash '/'.
    <div class=namelist><pre>
    <em class=call>namelist / obs_diag_nml / </em> &#38;
-      obs_sequence_name, first_bin_center, last_bin_center, &#38;
+      obs_sequence_name, obs_sequence_list, first_bin_center, last_bin_center, &#38;
       bin_separation, bin_width, time_to_skip, max_num_bins, &#38;
       plevel, hlevel, mlevel, <em class=removed>obs_select</em>, rat_cri, input_qc_threshold, &#38;
       Nregions, lonlim1, lonlim2, latlim1, latlim2, &#38;
@@ -296,6 +296,13 @@
    (YYYY,&nbsp;MO,&nbsp;DA,&nbsp;HR,&nbsp;MIN,&nbsp;SEC). <BR>
    The allowable ranges for the region boundaries are: latitude [-90.,90], 
    longitude [0.,Inf.]
+<br />
+<br />
+   You can only specify <strong>either</strong> 
+   <em class="code">obs_sequence_name</em>
+   <strong>or</strong> 
+   <em class="code">obs_sequence_list</em> -- not both.
+   One of them has to be an empty string ... i.e. <em class="code">''</em>.
 </P>
 
 
@@ -321,6 +328,18 @@
 	    'obsdir_002/obs_seq.final', and so on.<BR>
 	    Default 'obs_seq.final'</TD></TR>  
 
+<TR><!--contents--><TD valign=top>   obs_sequence_list </TD>
+    <!--  type  --><TD valign=top>   character         </TD>
+    <!--descript--><TD>Name of a file containing a list of observation 
+            sequence files. <BR>
+            You can only specify <strong>either</strong> 
+            <em class="code">obs_sequence_name</em>
+            <strong>or</strong> 
+            <em class="code">obs_sequence_list</em> -- not both.
+            One of them has to be an empty string ... i.e. <em class="code">''</em>.
+	    <BR>
+	    Default '' - an empty string.</TD></TR>  
+
 <TR><!--contents--><TD valign=top>   first_bin_center </TD>
     <!--  type  --><TD valign=top>  integer, dimension(6) </TD>
     <!--descript--><TD>first timeslot of the first obs_seq.final file to process.

Modified: DART/trunk/diagnostics/threed_sphere/obs_diag.nml
===================================================================
--- DART/trunk/diagnostics/threed_sphere/obs_diag.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/diagnostics/threed_sphere/obs_diag.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -10,6 +10,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/MITgcm_annulus/work/input.nml
===================================================================
--- DART/trunk/models/MITgcm_annulus/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/MITgcm_annulus/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -159,6 +159,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 6, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/MITgcm_ocean/work/input.nml
===================================================================
--- DART/trunk/models/MITgcm_ocean/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/MITgcm_ocean/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -240,6 +240,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  1996, 1, 1,12, 0, 0 ,
    last_bin_center  =  1996, 1,14,12, 0, 0 ,
    bin_separation   =     0, 0, 1, 0, 0, 0 ,

Modified: DART/trunk/models/PBL_1d/work/input.nml
===================================================================
--- DART/trunk/models/PBL_1d/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/PBL_1d/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -217,6 +217,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 5, 6, 0, 0, 0,
    last_bin_center  =  2003, 5, 7, 0, 0, 0,
    bin_separation   =     0, 0, 0, 1, 0, 0,

Modified: DART/trunk/models/POP/work/input.nml
===================================================================
--- DART/trunk/models/POP/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/POP/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -9,6 +9,7 @@
    last_obs_days         = -1,
    last_obs_seconds      = -1,
    output_interval       = 1,
+   trace_execution       = .false.,
    restart_in_file_name  = "perfect_ics",
    restart_out_file_name = "perfect_restart",
    obs_seq_in_file_name  = "obs_seq.in",
@@ -36,12 +37,12 @@
    num_output_obs_members   = 10,
    output_interval          = 1,
    num_groups               = 1,
-   input_qc_threshold       =  4.0,
+   input_qc_threshold       =  1.0,
    outlier_threshold        =  3.0,
    output_forward_op_errors = .false.,
    output_timestamps        = .false.,
    output_inflation         = .true.,
-   trace_execution          = .true.,
+   trace_execution          = .false.,
 
    inf_flavor                  = 0,                       0,
    inf_initial_from_restart    = .false.,                 .false.,
@@ -101,44 +102,24 @@
 
 &obs_kind_nml
    assimilate_these_obs_types = 'FLOAT_TEMPERATURE',
-     evaluate_these_obs_types = 'SEA_SURFACE_HEIGHT',
                                 'FLOAT_SALINITY',
+                                'BOTTLE_SALINITY',
+                                'BOTTLE_TEMPERATURE',
+                                'DRIFTER_SALINITY',
+                                'DRIFTER_TEMPERATURE',
+                                'MOORING_SALINITY',
+                                'MOORING_TEMPERATURE',
+                                'CTD_SALINITY',
+                                'CTD_TEMPERATURE',
+                                'XCTD_SALINITY',
+                                'XCTD_TEMPERATURE',
+                                'XBT_TEMPERATURE',
+                                'MBT_TEMPERATURE',
+     evaluate_these_obs_types = 'SEA_SURFACE_HEIGHT',
                                 'U_CURRENT_COMPONENT',
                                 'V_CURRENT_COMPONENT'
   /
 
-#          xxx_these_obs_types = 'SALINITY',
-#                                'TEMPERATURE',
-#                                'U_CURRENT_COMPONENT',
-#                                'V_CURRENT_COMPONENT',
-#                                'SEA_SURFACE_HEIGHT',
-#                                'ARGO_U_CURRENT_COMPONENT',
-#                                'ARGO_V_CURRENT_COMPONENT',
-#                                'ARGO_SALINITY',
-#                                'ARGO_TEMPERATURE',
-#                                'ADCP_U_CURRENT_COMPONENT',
-#                                'ADCP_V_CURRENT_COMPONENT',
-#                                'ADCP_SALINITY',
-#                                'ADCP_TEMPERATURE',
-#                                'FLOAT_SALINITY',
-#                                'FLOAT_TEMPERATURE',
-#                                'DRIFTER_U_CURRENT_COMPONENT',
-#                                'DRIFTER_V_CURRENT_COMPONENT',
-#                                'DRIFTER_SALINITY',
-#                                'DRIFTER_TEMPERATURE',
-#                                'GLIDER_U_CURRENT_COMPONENT',
-#                                'GLIDER_V_CURRENT_COMPONENT',
-#                                'GLIDER_SALINITY',
-#                                'GLIDER_TEMPERATURE',
-#                                'MOORING_U_CURRENT_COMPONENT',
-#                                'MOORING_V_CURRENT_COMPONENT',
-#                                'MOORING_SALINITY',
-#                                'MOORING_TEMPERATURE',
-#                                'SATELLITE_MICROWAVE_SST',
-#                                'SATELLITE_INFRARED_SST',
-#                                'SATELLITE_SSH',
-#                                'SATELLITE_SSS',
-
 &preprocess_nml
     input_obs_kind_mod_file = '../../../obs_kind/DEFAULT_obs_kind_mod.F90',
    output_obs_kind_mod_file = '../../../obs_kind/obs_kind_mod.f90',
@@ -268,14 +249,15 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
-   first_bin_center =  2000, 1, 7, 0, 0, 0 ,
+   obs_sequence_list = '',
+   first_bin_center =  2000, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2000, 1,31, 0, 0, 0 ,
    bin_separation   =     0, 0, 1, 0, 0, 0 ,
    bin_width        =     0, 0, 1, 0, 0, 0 ,
    time_to_skip     =     0, 0, 0, 0, 0, 0 ,
    max_num_bins     = 1000,
    rat_cri            = 3.0,
-   input_qc_threshold = 4.0,
+   input_qc_threshold = 1,
    Nregions   = 3,
    lonlim1    = 30.0,  130.0, 290.0,
    lonlim2    = 130.0, 270.0, 380.0,
@@ -285,25 +267,24 @@
    print_mismatched_locs = .false.,
    print_obs_locations   = .false.,
    verbose               = .true.,
-   hlevel =  10.0, 20.0, 30.0, 40.0,
-            100.0,  200.0,  500.0, 1000.0,  2000.0,
-            3000.0, 4000.0, 8000.0, 16000.0, 32000.0,
+   hlevel =  10.0,  20.0,   30.0,   40.0,  100.0,  
+            200.0, 500.0, 1000.0, 2000.0, 3000.0
   /
 
 &schedule_nml
    calendar        = 'Gregorian',
    first_bin_start =  1999, 12, 31, 12, 0, 0 ,
-   first_bin_end   =  2000,  1, 31, 12, 0, 0 ,
-   last_bin_end    =  2000,  1, 31, 12, 0, 0 ,
-   bin_interval_days    = 32,
+   first_bin_end   =  2000,  1,  1, 12, 0, 0 ,
+   last_bin_end    =  2000,  2,  1, 12, 0, 0 ,
+   bin_interval_days    = 1,
    bin_interval_seconds = 0,
    max_num_bins         = 1000,
    print_table          = .true.
   /
 
 &obs_seq_to_netcdf_nml
-   obs_sequence_name = ''
-   obs_sequence_list = 'obs_seq_file_list',
+   obs_sequence_name = 'obs_seq.final'
+   obs_sequence_list = '',
    lonlim1    =    0.0,
    lonlim2    =  360.0,
    latlim1    =  -90.0,

Modified: DART/trunk/models/am2/work/input.nml
===================================================================
--- DART/trunk/models/am2/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/am2/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -237,6 +237,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 6, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/bgrid_solo/work/input.nml
===================================================================
--- DART/trunk/models/bgrid_solo/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/bgrid_solo/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -273,6 +273,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  1601, 1, 1, 0, 0, 0 ,
    last_bin_center  =  1601, 1, 1, 3, 0, 0 ,
    bin_separation   =     0, 0, 0, 0, 3, 0 ,

Modified: DART/trunk/models/cam/work/input.nml
===================================================================
--- DART/trunk/models/cam/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/cam/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -228,6 +228,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 6, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,
@@ -246,4 +247,21 @@
    print_obs_locations = .false.,
    verbose = .false.  /
 
+&restart_file_utility_nml 
+   input_file_name              = "filter_restart",        
+   output_file_name             = "filter_updated_restart",
+   ens_size                     = 1,
+   single_restart_file_in       = .true.,
+   single_restart_file_out      = .true.,
+   write_binary_restart_files   = .true.,
+   overwrite_data_time          = .false.,
+   new_data_days                = -1,  
+   new_data_secs                = -1,
+   input_is_model_advance_file  = .false.,
+   output_is_model_advance_file = .false.,
+   overwrite_advance_time       = .false.,
+   new_advance_days             = -1, 
+   new_advance_secs             = -1
+/
 
+

Modified: DART/trunk/models/pe2lyr/work/input.nml
===================================================================
--- DART/trunk/models/pe2lyr/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/pe2lyr/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -196,6 +196,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2003, 1, 7, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/rose/work/input.nml
===================================================================
--- DART/trunk/models/rose/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/rose/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -163,6 +163,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/wrf/regression/CONUS-V2/input.nml
===================================================================
--- DART/trunk/models/wrf/regression/CONUS-V2/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/wrf/regression/CONUS-V2/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -129,6 +129,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2007, 4, 26, 00, 0, 0 ,
    last_bin_center  =  2007, 4, 27, 00, 0, 0 ,
    bin_separation   =     0, 0, 0, 6, 0, 0 ,

Modified: DART/trunk/models/wrf/regression/Radar/input.nml
===================================================================
--- DART/trunk/models/wrf/regression/Radar/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/wrf/regression/Radar/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -175,6 +175,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0,12, 0, 0 ,

Modified: DART/trunk/models/wrf/work/input.nml
===================================================================
--- DART/trunk/models/wrf/work/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/models/wrf/work/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -223,6 +223,7 @@
 
 &obs_diag_nml
    obs_sequence_name  = 'obs_seq.final',
+   obs_sequence_list  = '',
    first_bin_center   =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center    =  2003, 1, 2, 0, 0, 0 ,
    bin_separation     =     0, 0, 0,12, 0, 0 ,
@@ -239,8 +240,29 @@
    reg_names  = 'Northern Hemisphere', 'Southern Hemisphere', 'Tropics', 'North America',
    print_mismatched_locs = .false.,
    print_obs_locations   = .false.,
-   verbose               = .false.  /
+   verbose               = .false.
+   /
 
+&schedule_nml
+   calendar        = 'Gregorian',  
+   first_bin_start =  2008, 9, 4, 3, 0, 0 ,
+   first_bin_end   =  2008, 9, 4, 9, 0, 0 ,
+   last_bin_end    =  2008, 9, 15, 0, 0, 0 ,
+   bin_interval_days    = 0,
+   bin_interval_seconds = 21600,   
+   max_num_bins         = 1000,
+   print_table          = .true.
+   /
+   
+&obs_seq_to_netcdf_nml  
+   obs_sequence_name = 'obs_seq.final',
+   lonlim1    =    0.0,
+   lonlim2    =  360.0,
+   latlim1    =  -90.0,
+   latlim2    =   90.0,
+   verbose    = .false.   
+   /
+
 &restart_file_tool_nml 
    input_file_name              = "filter_restart",        
    output_file_name             = "filter_updated_restart",
@@ -256,9 +278,8 @@
    overwrite_advance_time       = .false.,
    new_advance_days             = -1, 
    new_advance_secs             = -1
-/
+   /
 
-
 &obs_sequence_tool_nml
    num_input_files   = 1, 
    filename_seq      = 'obs_seq.out',

Modified: DART/trunk/observations/utilities/threed_sphere/input.nml
===================================================================
--- DART/trunk/observations/utilities/threed_sphere/input.nml	2009-10-21 19:48:01 UTC (rev 4112)
+++ DART/trunk/observations/utilities/threed_sphere/input.nml	2009-10-22 20:33:36 UTC (rev 4113)
@@ -68,6 +68,7 @@
 
 &obs_diag_nml
    obs_sequence_name = 'obs_seq.final',
+   obs_sequence_list = '',
    first_bin_center =  2003, 1, 1, 0, 0, 0 ,
    last_bin_center  =  2003, 1, 2, 0, 0, 0 ,
    bin_separation   =     0, 0, 0,12, 0, 0 ,


More information about the Dart-dev mailing list