diff --git a/met/data/config/TCPairsConfig_default b/met/data/config/TCPairsConfig_default
index 8c62d1dce3..aa4afb51b6 100644
--- a/met/data/config/TCPairsConfig_default
+++ b/met/data/config/TCPairsConfig_default
@@ -50,10 +50,17 @@ init_inc = [];
 init_exc = [];
 
 //
-// Valid model time window
+// Valid model time windows to include or exclude
 //
 valid_beg = "";
 valid_end = "";
+valid_inc = [];
+valid_exc = [];
+
+//
+// Valid times for which output should be written
+//
+write_valid = [];
 
 //
 // Model initialization hours
diff --git a/met/docs/Users_Guide/config_options_tc.rst b/met/docs/Users_Guide/config_options_tc.rst
index f3fcc75450..b9542a878d 100644
--- a/met/docs/Users_Guide/config_options_tc.rst
+++ b/met/docs/Users_Guide/config_options_tc.rst
@@ -92,6 +92,8 @@ or exclude (exc). Tracks whose initial time meets the specified
 criteria will be used. An empty string indicates that all times
 should be used.
 
+In TC-Stat, the **-init_beg**, **-init_end**, **init_inc** and **-int_exc** job command options can be used to further refine these selections.
+
 For example:
 
 | init_beg = "20100101";
@@ -107,26 +109,34 @@ For example:
   init_inc = [];
   init_exc = [];
 
+.. _valid_beg end inc exc_1:
 
-.. _valid_beg, valid_end_1:
-
-:ref:`valid_beg, valid_end <valid_beg, valid_end_1>`
-
-Specify a model valid time window in YYYYMMDD[_HH[MMSS]] format.
-Tracks for which all valid times fall within the time window will be used.
-An empty string indicates that all times should be used.
+ ref:`valid_beg, valid_end, valid_inc, valid_exc <valid_beg end inc exc_1>`
+  
+Specify a model valid time window YYYYMMDD[_HH[MMSS]] format or provide a
+list of specific valid times to include (inc) or exclude (exc). If a time
+window is specified, only tracks for which all points are contained within
+the window will be used. If valid times to include or exclude are specified,
+tracks will be subset down to the points which meet that criteria. Empty
+begin/end time strings and empty include/exclude lists indicate that all
+valid times should be used.
 
+In TC-Stat, the **-valid_beg**, **-valid_end**, **valid_inc** and **-valid_exc** job command options can be used to further refine these selections.
 
 For example:
-		
+
 | valid_beg = "20100101";
-| valid_end = "20101231";
-| 
+| valid_end = "20101231_12";
+| valid_inc = [ "20101231_06" ];
+| valid_exc = [ "20101231_00" ];
+|
 
 .. code-block:: none
-		
+
   valid_beg = "";
   valid_end = "";
+  valid_inc = [];
+  valid_exc = [];
 
 .. _init_hour_1:
 
@@ -422,7 +432,6 @@ by basin or sub-basin. Note that if your model data and best track do not
 use the same basin identifier conventions, using an empty list for this
 parameter will result in missed matches.
 
-
 .. code-block:: none
 
   basin_map = [
@@ -481,9 +490,9 @@ For example:
   valid_inc = [];
   valid_exc = [];
 
-.. _ini valid_hour lead req:
+.. _init valid_hour lead req:
 
-:ref:`ini_hour, valid_hour, lead, lead_req <ini valid_hour lead req>`
+:ref:`init_hour, valid_hour, lead, lead_req <init valid_hour lead req>`
 
 Stratify by the initialization and valid hours and lead time.
 Specify a comma-separated list of initialization hours,
diff --git a/met/docs/Users_Guide/tc-pairs.rst b/met/docs/Users_Guide/tc-pairs.rst
index f2ce129846..b32cbc15b1 100644
--- a/met/docs/Users_Guide/tc-pairs.rst
+++ b/met/docs/Users_Guide/tc-pairs.rst
@@ -99,6 +99,8 @@ ____________________
   init_exc     = [];
   valid_beg    = "";
   valid_end    = "";
+  valid_inc    = [];
+  valid_exc    = [];
   init_hour    = [];
   init_mask    = [];
   lead_req     = [];
@@ -118,6 +120,14 @@ The **model** variable contains a list of comma-separated models to be used. Eac
 
 ____________________
 
+.. code-block:: none
+
+  write_valid = [ "20101231_06" ];
+
+The **write_valid** entry specifies a comma-separated list of valid time strings in YYYYMMDD[_HH[MMSS]] format for which output should be written. An empty list indicates that data for all valid times should be written. This option may be useful when verifying track forecasts in realtime. If evaluating performance for a single valid time, this option can limit the output to that time and skip output for earlier track points.
+
+____________________
+
 .. code-block:: none
 
   check_dup = FALSE;
@@ -315,7 +325,7 @@ TC-Pairs produces output in TCST format. The default output file name can be ove
     - Valid time masking grid applied
   * - 13
     - LINE_TYPE
-    - Output line type (TCMPR or PROBRI)
+    - Output line type (TCMPR or PROBRIRW)
 
 .. _TCMPR Line Type:
 
@@ -462,9 +472,9 @@ TC-Pairs produces output in TCST format. The default output file name can be ove
     - A/BDEPTH
     - system depth, D-deep, M-medium, S-shallow, X-unknown
 
-.. _PROBRI Line Type:
+.. _PROBRIRW Line Type:
 
-.. list-table:: Format information for PROBRIRW (Probability of Rapid Intensification) output line type.
+.. list-table:: Format information for PROBRIRW (Probability of Rapid Intensification/Weakening) output line type.
   :widths: auto
   :header-rows: 2
 
@@ -475,8 +485,8 @@ TC-Pairs produces output in TCST format. The default output file name can be ove
     - Header Column Name
     - Description
   * - 13
-    - PROBRI
-    - Probability of Rapid Intensification line type
+    - PROBRIRW
+    - Probability of Rapid Intensification/Weakening line type
   * - 14
     - ALAT
     - Latitude position of edeck model
diff --git a/met/docs/Users_Guide/tc-stat.rst b/met/docs/Users_Guide/tc-stat.rst
index ee774eee0d..7052e80958 100644
--- a/met/docs/Users_Guide/tc-stat.rst
+++ b/met/docs/Users_Guide/tc-stat.rst
@@ -139,6 +139,8 @@ _________________________
   init_exc     = [];
   valid_beg    = "";
   valid_end    = "";
+  valid_inc    = [];
+  valid_exc    = [];
   init_hour    = [];
   lead_req     = [];
   init_mask    = [];
@@ -161,15 +163,6 @@ The **amodel** and **bmodel** fields stratify by the amodel and bmodel columns b
 
 _________________________
 
-.. code-block:: none
-
-  valid_inc = [];
-  valid_exc = [];
-
-The **valid_inc** and **valid_exc** fields stratify by valid times, based on a comma-separated list of specific valid times to include (inc) or exclude (exc). Time strings are defined by YYYYMMDD[_HH[MMSS]]. Using the **-valid_inc** and **-valid_exc** options within the job command lines may further refine these selections.
-
-_________________________
-
 .. code-block:: none
 
   valid_hour = [];
@@ -372,7 +365,7 @@ _________________________
         -job summary -line_type TCMPR  -column TK_ERR -dump_row ./tc_summary_job.tcst
         -job rirw    -line_type TCMPR  -rirw_time 24 -rirw_exact false -rirw_thresh ge20
         -job probrirw -line_type PROBRIRW -column_thresh RI_WINDOW ==24 \
-                      -probri_thresh 30 -probri_prob_thresh ==0.25
+                      -probrirw_thresh 30 -probrirw_prob_thresh ==0.25
 
 .. _tc_stat-output:
 
@@ -466,12 +459,12 @@ The PROBRIRW job produces probabilistic contingency table counts and statistics
 
 • The **-prob_exact bool** option is a boolean defining whether the exact or maximum BEST track intensity change over the time window should be used. If true, the values in the **BDELTA** column are used. If false, the values in the **BDELTA_MAX** column are used. The default is true.
 
-• The **-probri_bdelta_thresh** threshold option defines the BEST track intensity change event threshold. This should typically be set consistent with the probability threshold (**-prob_thresh**) chosen above. The default is greater than or equal to 30 kts.
+• The **-probrirw_bdelta_thresh** threshold option defines the BEST track intensity change event threshold. This should typically be set consistent with the probability threshold (**-prob_thresh**) chosen above. The default is greater than or equal to 30 kts.
 
-• The **-probri_prob_thresh threshold_list** option defines the probability thresholds used to create the output Nx2 contingency table. The default is probability bins of width 0.1. These probabilities may be specified as a list (>0.00,>0.25,>0.50,>0.75,>1.00) or using shorthand notation (==0.25) for bins of equal width.
+• The **-probrirw_prob_thresh threshold_list** option defines the probability thresholds used to create the output Nx2 contingency table. The default is probability bins of width 0.1. These probabilities may be specified as a list (>0.00,>0.25,>0.50,>0.75,>1.00) or using shorthand notation (==0.25) for bins of equal width.
 
 • The **-out_line_type** option defines the output data that should be written. This job can write PCT, PSTD, PJC, and PRC output line types. The default is PCT and PSTD.  Please see :numref:`table_PS_format_info_PCT` through :numref:`table_PS_format_info_PRC` for more details.
 
 Users may also specify the **-out_alpha** option to define the alpha value for the confidence intervals in the PSTD output line type. Multiple values in the **RI_WINDOW** column cannot be combined in a single PROBRIRW job since the BEST track intensity threshold should change for each. Using the **-by RI_WINDOW** option or **-column_thresh RI_WINDOW ==24** option provide convenient ways avoiding this problem.
 
-Users should note that for the PROBRIRW line type, **PROBRI_PROB** is a derived column name. The -probri_thresh option defines the probabilities of interest (e.g. **-probri_thresh 30**) and the **PROBRI_PROB** column name refers to those probability values, regardless of their column number. For example, the job command options **-probri_thresh 30 -column_thresh PROBRI_PROB >0** select 30 kt probabilities and match probability values greater than 0.
+Users should note that for the PROBRIRW line type, **PROBRI_PROB** is a derived column name. The -probrirw_thresh option defines the probabilities of interest (e.g. **-probrirw_thresh 30**) and the **PROBRI_PROB** column name refers to those probability values, regardless of their column number. For example, the job command options **-probrirw_thresh 30 -column_thresh PROBRI_PROB >0** select 30 kt probabilities and match probability values greater than 0.
diff --git a/met/src/basic/vx_config/config_constants.h b/met/src/basic/vx_config/config_constants.h
index 34fe74c978..81bd59589e 100644
--- a/met/src/basic/vx_config/config_constants.h
+++ b/met/src/basic/vx_config/config_constants.h
@@ -1012,6 +1012,7 @@ static const char conf_key_valid_end[]                = "valid_end";
 static const char conf_key_valid_inc[]                = "valid_inc";
 static const char conf_key_valid_exc[]                = "valid_exc";
 static const char conf_key_valid_hour[]               = "valid_hour";
+static const char conf_key_write_valid[]              = "write_valid";
 static const char conf_key_lead[]                     = "lead";
 static const char conf_key_lead_req[]                 = "lead_req";
 static const char conf_key_init_mask[]                = "init_mask";
diff --git a/met/src/libcode/vx_tc_util/prob_rirw_pair_info.cc b/met/src/libcode/vx_tc_util/prob_rirw_pair_info.cc
index 61b69cc125..ce38bdfe97 100644
--- a/met/src/libcode/vx_tc_util/prob_rirw_pair_info.cc
+++ b/met/src/libcode/vx_tc_util/prob_rirw_pair_info.cc
@@ -453,3 +453,20 @@ bool ProbRIRWPairInfoArray::add(const ProbRIRWInfo &p, const TrackInfo &t) {
 
 ////////////////////////////////////////////////////////////////////////
 
+void ProbRIRWPairInfoArray::subset_write_valid(const TimeArray &ta) {
+
+   // Check for no work to do
+   if(ta.n() == 0) return;
+
+   ProbRIRWPairInfoArray new_pairs;
+   for(int i=0; i<Pairs.size(); i++) {
+      if(ta.has(Pairs[i].prob_rirw().valid())) new_pairs.add(Pairs[i]);
+   }
+
+   // Save the subset
+   *this = new_pairs;
+
+   return;
+}
+
+////////////////////////////////////////////////////////////////////////
diff --git a/met/src/libcode/vx_tc_util/prob_rirw_pair_info.h b/met/src/libcode/vx_tc_util/prob_rirw_pair_info.h
index d3465e2743..69e1e29835 100644
--- a/met/src/libcode/vx_tc_util/prob_rirw_pair_info.h
+++ b/met/src/libcode/vx_tc_util/prob_rirw_pair_info.h
@@ -178,6 +178,7 @@ class ProbRIRWPairInfoArray {
 
       void add(const ProbRIRWPairInfo &);
       bool add(const ProbRIRWInfo &, const TrackInfo &);
+      void subset_write_valid(const TimeArray &);
 };
 
 ////////////////////////////////////////////////////////////////////////
diff --git a/met/src/libcode/vx_tc_util/tc_columns.cc b/met/src/libcode/vx_tc_util/tc_columns.cc
index b4b383268d..f3f9eea17e 100644
--- a/met/src/libcode/vx_tc_util/tc_columns.cc
+++ b/met/src/libcode/vx_tc_util/tc_columns.cc
@@ -160,7 +160,7 @@ void write_tc_mpr_row(TcHdrColumns &hdr, const TrackPairInfo &p,
 
       // Set the description
       if(p.n_lines() > i) {
-	hdr.set_desc((string)p.line(i)->get_item("DESC", false));
+         hdr.set_desc((string)p.line(i)->get_item("DESC", false));
       }
 
       // Write the header columns
@@ -191,7 +191,7 @@ void write_prob_rirw_row(TcHdrColumns &hdr, const ProbRIRWPairInfo &p,
 
    // Pass the description from the input line to the output
    if(p.line().n_items() > 0) {
-     hdr.set_desc((string)p.line().get_item("DESC", false));
+      hdr.set_desc((string)p.line().get_item("DESC", false));
    }
 
    // Write one line for all the probabilities
diff --git a/met/src/libcode/vx_tc_util/track_pair_info.cc b/met/src/libcode/vx_tc_util/track_pair_info.cc
index b82d0a5f0d..5a03a9c0fc 100644
--- a/met/src/libcode/vx_tc_util/track_pair_info.cc
+++ b/met/src/libcode/vx_tc_util/track_pair_info.cc
@@ -309,6 +309,14 @@ void TrackPairInfo::initialize(const TCStatLine &l) {
 
 void TrackPairInfo::set_keep(int i, int val) {
 
+   // Check range
+   if(i < 0 || i >= NPoints) {
+      mlog << Error
+           << "\nTrackPairInfo::set_keep(int, int) -> "
+           << "range check error for index value " << i << "\n\n";
+      exit(1);
+   }
+
    Keep.set(i, val);
 
    return;
@@ -780,13 +788,23 @@ TrackPairInfo TrackPairInfo::keep_subset() const {
    TrackPairInfo tpi;
 
    // Loop over the points
-   for(i=0; i<NLines; i++) {
+   for(i=0; i<NPoints; i++) {
 
       // Check the retention status
       if(Keep[i] == 0) continue;
 
-      // Add the current track pair point
-      tpi.add(Line[i]);
+      // Add from TC-Stat line data
+      if(NLines == NPoints) {
+         tpi.add(Line[i]);
+      }
+      // Otherwise, add from TC-Pairs track pair
+      else {
+         if(tpi.n_points() == 0) tpi.initialize(ADeck, BDeck);
+         tpi.add(ADeck[i], BDeck[i],
+                 ADeckDLand[i], BDeckDLand[i],
+                 TrackErr[i], XErr[i], YErr[i],
+                 AlongTrackErr[i], CrossTrackErr[i]);
+      }
    }
 
    return(tpi);
@@ -1008,4 +1026,26 @@ void TrackPairInfoArray::add_watch_warn(const ConcatString &ww_sid,
    return;
 }
 
+////////////////////////////////////////////////////////////////////////
+
+void TrackPairInfoArray::subset_write_valid(const TimeArray &ta) {
+
+   // Check for no work to do
+   if(ta.n() == 0) return;
+
+   // Check each point for requested valid times
+   int i, j, keep;
+   for(i=0; i<NPairs; i++) {
+      for(j=0; j<Pair[i].n_points(); j++) {
+         keep = (ta.has(Pair[i].valid(j)) ? 1 : 0);
+         Pair[i].set_keep(j, keep);
+      }
+
+      // Subset the track
+      Pair[i] = Pair[i].keep_subset();
+   }
+
+   return;
+}
+
 ///////////////////////////////////////////////////////////////////////
diff --git a/met/src/libcode/vx_tc_util/track_pair_info.h b/met/src/libcode/vx_tc_util/track_pair_info.h
index c2cc548b93..9bbd8fbb44 100644
--- a/met/src/libcode/vx_tc_util/track_pair_info.h
+++ b/met/src/libcode/vx_tc_util/track_pair_info.h
@@ -215,6 +215,7 @@ class TrackPairInfoArray {
 
       void add(const TrackPairInfo &);
       void add_watch_warn(const ConcatString &, WatchWarnType, unixtime);
+      void subset_write_valid(const TimeArray &);
 };
 
 ////////////////////////////////////////////////////////////////////////
diff --git a/met/src/tools/tc_utils/tc_pairs/tc_pairs.cc b/met/src/tools/tc_utils/tc_pairs/tc_pairs.cc
index 46c3a80b04..0013656a81 100644
--- a/met/src/tools/tc_utils/tc_pairs/tc_pairs.cc
+++ b/met/src/tools/tc_utils/tc_pairs/tc_pairs.cc
@@ -166,7 +166,7 @@ int main(int argc, char *argv[]) {
    process_decks();
 
    // List the output files
-   for(int i=0; i<out_files.n_elements(); i++) {
+   for(int i=0; i<out_files.n(); i++) {
       mlog << Debug(1)
            << "Output file: " << out_files[i] << "\n";
    }
@@ -201,9 +201,9 @@ void process_command_line(int argc, char **argv) {
    cline.parse();
 
    // Check for the minimum number of arguments
-   if((adeck_source.n_elements() == 0 &&
-       edeck_source.n_elements() == 0) ||
-      bdeck_source.n_elements()  == 0  ||
+   if((adeck_source.n() == 0 &&
+       edeck_source.n() == 0) ||
+      bdeck_source.n()  == 0  ||
       config_file.length()       == 0) {
       mlog << Error
            << "\nprocess_command_line(int argc, char **argv) -> "
@@ -215,25 +215,25 @@ void process_command_line(int argc, char **argv) {
    }
 
    // List the input ADECK track files
-   for(i=0; i<adeck_source.n_elements(); i++) {
+   for(i=0; i<adeck_source.n(); i++) {
       mlog << Debug(1)
-           << "[Source " << i+1 << " of " << adeck_source.n_elements()
+           << "[Source " << i+1 << " of " << adeck_source.n()
            << "] ADECK Source: " << adeck_source[i] << ", Model Suffix: "
            << adeck_model_suffix[i] << "\n";
    }
 
    // List the input EDECK track files
-   for(i=0; i<edeck_source.n_elements(); i++) {
+   for(i=0; i<edeck_source.n(); i++) {
       mlog << Debug(1)
-           << "[Source " << i+1 << " of " << edeck_source.n_elements()
+           << "[Source " << i+1 << " of " << edeck_source.n()
            << "] EDECK Source: " << edeck_source[i] << ", Model Suffix: "
            << edeck_model_suffix[i] << "\n";
    }
 
    // List the input BDECK track files
-   for(i=0; i<bdeck_source.n_elements(); i++) {
+   for(i=0; i<bdeck_source.n(); i++) {
       mlog << Debug(1)
-           << "[Source " << i+1 << " of " << bdeck_source.n_elements()
+           << "[Source " << i+1 << " of " << bdeck_source.n()
            << "] BDECK Source: " << bdeck_source[i] << ", Model Suffix: "
            << bdeck_model_suffix[i] << "\n";
    }
@@ -266,12 +266,12 @@ void process_decks() {
    process_bdecks(bdeck_tracks);
 
    // Process ADECK files
-   if(adeck_source.n_elements() > 0) {
+   if(adeck_source.n() > 0) {
       process_adecks(bdeck_tracks);
    }
 
    // Process EDECK files
-   if(edeck_source.n_elements() > 0) {
+   if(edeck_source.n() > 0) {
       process_edecks(bdeck_tracks);
    }
 
@@ -291,7 +291,7 @@ void process_bdecks(TrackInfoArray &bdeck_tracks) {
                   files, files_model_suffix);
 
    mlog << Debug(2)
-        << "Processing " << files.n_elements() << " BDECK file(s).\n";
+        << "Processing " << files.n() << " BDECK file(s).\n";
    process_track_files(files, files_model_suffix, bdeck_tracks, false,
                        (conf_info.AnlyTrack == TrackType_BDeck ||
                         conf_info.AnlyTrack == TrackType_Both));
@@ -314,7 +314,7 @@ void process_adecks(const TrackInfoArray &bdeck_tracks) {
                   files, files_model_suffix);
 
    mlog << Debug(2)
-        << "Processing " << files.n_elements() << " ADECK file(s).\n";
+        << "Processing " << files.n() << " ADECK file(s).\n";
    process_track_files(files, files_model_suffix, adeck_tracks, true,
                        (conf_info.AnlyTrack == TrackType_ADeck ||
                         conf_info.AnlyTrack == TrackType_Both));
@@ -340,7 +340,7 @@ void process_adecks(const TrackInfoArray &bdeck_tracks) {
 
    // Derive lag forecasts from the ADECK tracks
    mlog << Debug(2)
-        << "Deriving " << conf_info.LagTime.n_elements()
+        << "Deriving " << conf_info.LagTime.n()
         << " ADECK lag model(s).\n";
    i = derive_lag(adeck_tracks);
    mlog << Debug(2)
@@ -348,8 +348,8 @@ void process_adecks(const TrackInfoArray &bdeck_tracks) {
 
    // Derive CLIPER/SHIFOR forecasts from the ADECK/BDECK tracks
    mlog << Debug(2)
-        << "Deriving " << conf_info.OperBaseline.n_elements() +
-                          conf_info.BestBaseline.n_elements()
+        << "Deriving " << conf_info.OperBaseline.n() +
+                          conf_info.BestBaseline.n()
         << " CLIPER/SHIFOR baseline model(s).\n";
    i = derive_baseline(adeck_tracks, bdeck_tracks);
    mlog << Debug(2)
@@ -403,6 +403,14 @@ void process_adecks(const TrackInfoArray &bdeck_tracks) {
            << pairs.serialize_r() << "\n";
    }
 
+   // Subset pairs based on requested valid output times
+   if(conf_info.WriteValid.n() > 0) {
+      mlog << Debug(3) << "Subsetting output for "
+           << conf_info.WriteValid.n()
+           << " requested valid time(s).\n";
+      pairs.subset_write_valid(conf_info.WriteValid);
+   }
+
    // Write out the track pairs
    write_tracks(pairs);
 
@@ -423,7 +431,7 @@ void process_edecks(const TrackInfoArray &bdeck_tracks) {
                   files, files_model_suffix);
 
    mlog << Debug(2)
-        << "Processing " << files.n_elements()
+        << "Processing " << files.n()
         << " EDECK file(s).\n";
    process_prob_files(files, files_model_suffix, edeck_probs);
 
@@ -483,6 +491,14 @@ void process_edecks(const TrackInfoArray &bdeck_tracks) {
            << prob_rirw_pairs.serialize_r() << "\n";
    }
 
+   // Subset pairs based on requested valid output times
+   if(conf_info.WriteValid.n() > 0) {
+      mlog << Debug(3) << "Subsetting output for "
+           << conf_info.WriteValid.n()
+           << " requested valid time(s).\n";
+      prob_rirw_pairs.subset_write_valid(conf_info.WriteValid);
+   }
+
    // Write out the ProbRIRW pairs
    if(prob_rirw_pairs.n_pairs() > 0) write_prob_rirw(prob_rirw_pairs);
 
@@ -498,7 +514,7 @@ void get_atcf_files(const StringArray &source,
    StringArray cur_source, cur_files;
    int i, j;
 
-   if(source.n_elements() != model_suffix.n_elements()) {
+   if(source.n() != model_suffix.n()) {
       mlog << Error
            << "\nget_atcf_files() -> "
            << "the source and suffix arrays must be equal length!\n\n";
@@ -510,12 +526,12 @@ void get_atcf_files(const StringArray &source,
    files_model_suffix.clear();
 
    // Build list of files and corresponding model suffix list
-   for(i=0; i<source.n_elements(); i++) {
+   for(i=0; i<source.n(); i++) {
       cur_source.clear();
       cur_source.add(source[i]);
       cur_files = get_filenames(cur_source, NULL, atcf_suffix);
 
-      for(j=0; j<cur_files.n_elements(); j++) {
+      for(j=0; j<cur_files.n(); j++) {
          files.add(cur_files[j]);
          files_model_suffix.add(model_suffix[i]);
       }
@@ -547,7 +563,7 @@ void process_track_files(const StringArray &files,
    line.set_oper_technique(&conf_info.OperTechnique);
 
    // Process each of the input ATCF files
-   for(i=0; i<files.n_elements(); i++) {
+   for(i=0; i<files.n(); i++) {
 
       // Open the current file
       if(!f.open(files[i].c_str())) {
@@ -583,7 +599,7 @@ void process_track_files(const StringArray &files,
 
       // Dump out the current number of lines
       mlog << Debug(4)
-           << "[File " << i+1 << " of " << files.n_elements()
+           << "[File " << i+1 << " of " << files.n()
            << "] Used " << cur_add << " of " << cur_read
            << " lines read from file \"" << files[i] << "\"\n";
 
@@ -595,7 +611,7 @@ void process_track_files(const StringArray &files,
    // Dump out the total number of lines
    mlog << Debug(3)
         << "Used " << tot_add << " of " << tot_read
-        << " lines read from " << files.n_elements() << " file(s).\n";
+        << " lines read from " << files.n() << " file(s).\n";
 
    // Dump out the track information
    mlog << Debug(3)
@@ -640,7 +656,7 @@ void process_prob_files(const StringArray &files,
    line.set_oper_technique(&conf_info.OperTechnique);
 
    // Process each of the input ATCF files
-   for(i=0; i<files.n_elements(); i++) {
+   for(i=0; i<files.n(); i++) {
 
       // Open the current file
       if(!f.open(files[i].c_str())) {
@@ -676,7 +692,7 @@ void process_prob_files(const StringArray &files,
 
       // Dump out the current number of lines
       mlog << Debug(4)
-           << "[File " << i+1 << " of " << files.n_elements()
+           << "[File " << i+1 << " of " << files.n()
            << "] Used " << cur_add << " of " << cur_read
            << " lines read from file \"" << files[i] << "\"\n";
 
@@ -688,7 +704,7 @@ void process_prob_files(const StringArray &files,
    // Dump out the total number of lines
    mlog << Debug(3)
         << "Used " << tot_add << " of " << tot_read
-        << " lines read from " << files.n_elements() << " file(s).\n";
+        << " lines read from " << files.n() << " file(s).\n";
 
    // Dump out the track information
    mlog << Debug(3)
@@ -715,7 +731,7 @@ void process_prob_files(const StringArray &files,
 //
 // Check if the ATCFLineBase should be kept.  Only check those columns
 // that remain constant across the entire track:
-//    model, storm id, basin, cyclone, and init time
+//    model, storm id, basin, cyclone, and timing information
 //
 ////////////////////////////////////////////////////////////////////////
 
@@ -727,23 +743,23 @@ bool is_keeper(const ATCFLineBase * line) {
    unix_to_mdyhms(line->warning_time(), m, d, y, h, mm, s);
 
    // Check model
-   if(conf_info.Model.n_elements() > 0 &&
+   if(conf_info.Model.n() > 0 &&
       !conf_info.Model.has(line->technique()))
       keep = false;
 
    // Check storm id
-   else if(conf_info.StormId.n_elements() > 0 &&
+   else if(conf_info.StormId.n() > 0 &&
            !has_storm_id(conf_info.StormId, line->basin(),
                          line->cyclone_number(), line->warning_time()))
       keep = false;
 
    // Check basin
-   else if(conf_info.Basin.n_elements() > 0 &&
+   else if(conf_info.Basin.n() > 0 &&
            !conf_info.Basin.has(line->basin()))
       keep = false;
 
    // Check cyclone
-   else if(conf_info.Cyclone.n_elements() > 0 &&
+   else if(conf_info.Cyclone.n() > 0 &&
            !conf_info.Cyclone.has(line->cyclone_number()))
       keep = false;
 
@@ -752,17 +768,24 @@ bool is_keeper(const ATCFLineBase * line) {
             conf_info.InitBeg > line->warning_time()) ||
            (conf_info.InitEnd > 0 &&
             conf_info.InitEnd < line->warning_time()) ||
-           (conf_info.InitInc.n_elements() > 0 &&
+           (conf_info.InitInc.n() > 0 &&
             !conf_info.InitInc.has(line->warning_time())) ||
-           (conf_info.InitExc.n_elements() > 0 &&
+           (conf_info.InitExc.n() > 0 &&
             conf_info.InitExc.has(line->warning_time())))
       keep = false;
 
    // Initialization hour
-   else if(conf_info.InitHour.n_elements() > 0 &&
+   else if(conf_info.InitHour.n() > 0 &&
            !conf_info.InitHour.has(hms_to_sec(h, mm, s)))
       keep = false;
 
+   // Valid time include/exclude
+   else if((conf_info.ValidInc.n() > 0 &&
+            !conf_info.ValidInc.has(line->valid())) ||
+           (conf_info.ValidExc.n() > 0 &&
+            conf_info.ValidExc.has(line->valid())))
+      keep = false;
+
    // Return the keep status
    return(keep);
 }
@@ -781,11 +804,11 @@ void filter_tracks(TrackInfoArray &tracks) {
 
    // Loop through the tracks and determine which should be retained
    // The is_keeper() function has already filtered by model, storm id,
-   // basin, cyclone, initialization time, and initialization hour.
+   // basin, cyclone, and timing information.
    for(i=0; i<t.n(); i++) {
 
       // Check storm name
-      if(conf_info.StormName.n_elements() > 0 &&
+      if(conf_info.StormName.n() > 0 &&
          !conf_info.StormName.has(t[i].storm_name())) {
          mlog << Debug(4)
               << "Discarding track " << i+1 << " for storm name mismatch: "
@@ -812,10 +835,10 @@ void filter_tracks(TrackInfoArray &tracks) {
       // These are used in determining whether to keep or discard a track; keep a track
       // if all the required lead times are present.  If no required lead times are
       // specified in the config file, then ignore checking and proceed as usual.
-      if(conf_info.LeadReq.n_elements() > 0) {
+      if(conf_info.LeadReq.n() > 0) {
 
          // Loop over the required lead times
-         for(j=0, status=true; j<conf_info.LeadReq.n_elements(); j++) {
+         for(j=0, status=true; j<conf_info.LeadReq.n(); j++) {
 
             // If required lead time is missing, break out
             if(t[i].lead_index(conf_info.LeadReq[j]) == -1) {
@@ -1142,11 +1165,11 @@ int derive_consensus(TrackInfoArray &tracks) {
    } // end for i
 
    mlog << Debug(3)
-        << "Building consensus track(s) for " << case_list.n_elements()
+        << "Building consensus track(s) for " << case_list.n()
         << " cases.\n";
 
    // Loop through the cases and process each consensus model
-   for(i=0; i<case_list.n_elements(); i++) {
+   for(i=0; i<case_list.n(); i++) {
 
       // Break the case back out into Basin, Cyclone, InitTime
       cur_case = case_list[i];
@@ -1162,7 +1185,7 @@ int derive_consensus(TrackInfoArray &tracks) {
 
          // Loop through the consensus members
          for(k=0, skip=false;
-             k<conf_info.Consensus[j].Members.n_elements(); k++) {
+             k<conf_info.Consensus[j].Members.n(); k++) {
 
             // Add required members to the list
             if(conf_info.Consensus[j].Required[k]) {
@@ -1262,10 +1285,10 @@ int derive_lag(TrackInfoArray &tracks) {
    int n_add = 0;
 
    // If no time lags are requested, nothing to do
-   if(conf_info.LagTime.n_elements() == 0) return(0);
+   if(conf_info.LagTime.n() == 0) return(0);
 
    // Loop through the time lags to be applied
-   for(i=0; i<conf_info.LagTime.n_elements(); i++) {
+   for(i=0; i<conf_info.LagTime.n(); i++) {
 
       mlog << Debug(3)
            << "Building time-lagged track(s) for lag time \""
@@ -1335,12 +1358,12 @@ int derive_baseline(TrackInfoArray &atracks, const TrackInfoArray &btracks) {
    int n_add = 0;
 
    // If no baseline models are requested, nothing to do
-   if(conf_info.OperBaseline.n_elements() == 0 &&
-      conf_info.BestBaseline.n_elements() == 0) return(0);
+   if(conf_info.OperBaseline.n() == 0 &&
+      conf_info.BestBaseline.n() == 0) return(0);
 
    mlog << Debug(3)
         << "Building CLIPER/SHIFOR operational baseline forecasts using "
-        << conf_info.OperBaseline.n_elements() << " method(s).\n";
+        << conf_info.OperBaseline.n() << " method(s).\n";
 
    // Loop over the ADECK tracks
    for(i=0; i<atracks.n(); i++) {
@@ -1357,7 +1380,7 @@ int derive_baseline(TrackInfoArray &atracks, const TrackInfoArray &btracks) {
       if(!atracks[i].is_oper_track()) continue;
 
       // Loop over the operational baseline methods
-      for(j=0; j<conf_info.OperBaseline.n_elements(); j++) {
+      for(j=0; j<conf_info.OperBaseline.n(); j++) {
 
          // Derive the baseline model
          derive_baseline_model(conf_info.OperBaseline[j],
@@ -1368,8 +1391,8 @@ int derive_baseline(TrackInfoArray &atracks, const TrackInfoArray &btracks) {
 
    mlog << Debug(3)
         << "Building CLIPER/SHIFOR BEST track baseline forecasts using "
-        << conf_info.BestBaseline.n_elements() << " method(s) for "
-        << case_list.n_elements() << " cases.\n";
+        << conf_info.BestBaseline.n() << " method(s) for "
+        << case_list.n() << " cases.\n";
 
    // Loop over the BDECK tracks
    for(i=0; i<btracks.n(); i++) {
@@ -1389,7 +1412,7 @@ int derive_baseline(TrackInfoArray &atracks, const TrackInfoArray &btracks) {
          if(!case_list.has(cur_case)) continue;
 
          // Loop over the BEST track baseline methods
-         for(k=0; k<conf_info.BestBaseline.n_elements(); k++) {
+         for(k=0; k<conf_info.BestBaseline.n(); k++) {
 
             // Derive the baseline model
             derive_baseline_model(conf_info.BestBaseline[k],
@@ -1644,13 +1667,13 @@ void process_match(const TrackInfo &adeck, const TrackInfo &bdeck,
    valid_list.sort_array();
 
    mlog << Debug(4)
-        << "Processing " << valid_list.n_elements()
+        << "Processing " << valid_list.n()
         << " unique valid times: "
         << unix_to_yyyymmdd_hhmmss(valid_list.min()) << " to "
         << unix_to_yyyymmdd_hhmmss(valid_list.max()) << "\n";
 
    // Loop through the valid times
-   for(i=0; i<valid_list.n_elements(); i++) {
+   for(i=0; i<valid_list.n(); i++) {
 
       // Get the indices for this time
       i_adeck = adeck.valid_index(valid_list[i]);
@@ -2045,7 +2068,7 @@ void write_prob_rirw(const ProbRIRWPairInfoArray &p) {
       tchc.set_cyclone(p[i].prob_rirw().cyclone());
       tchc.set_storm_name(p[i].bdeck()->storm_name());
 
-      // Write the current TrackPairInfo object
+      // Write the current ProbRIRWPairInfo object
       write_prob_rirw_row(tchc, p[i], out_at, i_row);
    }
 
@@ -2159,11 +2182,11 @@ void set_atcf_source(const StringArray & a,
    ConcatString cs, suffix;
 
    // Check for optional suffix sub-argument
-   for(i=0; i<a.n_elements(); i++) {
+   for(i=0; i<a.n(); i++) {
       cs = a[i];
       if(cs.startswith("suffix")) {
          sa = cs.split("=");
-         if(sa.n_elements() != 2) {
+         if(sa.n() != 2) {
             mlog << Error
                  << "\nset_atcf_source() -> "
                  << "the model suffix must be specified as "
@@ -2177,7 +2200,7 @@ void set_atcf_source(const StringArray & a,
    }
 
    // Parse the remaining sources
-   for(i=0; i<a.n_elements(); i++) {
+   for(i=0; i<a.n(); i++) {
       cs = a[i];
       if(cs.startswith("suffix")) continue;
       source.add(a[i]);
diff --git a/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.cc b/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.cc
index 9bb8105282..9e038226ee 100644
--- a/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.cc
+++ b/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.cc
@@ -71,6 +71,9 @@ void TCPairsConfInfo::clear() {
    InitExc.clear();
    InitHour.clear();
    ValidBeg = ValidEnd = (unixtime) 0;
+   ValidInc.clear();
+   ValidExc.clear();
+   WriteValid.clear();
    InitMaskName.clear();
    InitPolyMask.clear();
    InitGridMask.clear();
@@ -172,6 +175,21 @@ void TCPairsConfInfo::process_config() {
    ValidBeg = Conf.lookup_unixtime(conf_key_valid_beg);
    ValidEnd = Conf.lookup_unixtime(conf_key_valid_end);
 
+   // Conf: ValidInc
+   sa = Conf.lookup_string_array(conf_key_valid_inc);
+   for(i=0; i<sa.n_elements(); i++)
+      ValidInc.add(timestring_to_unix(sa[i].c_str()));
+
+   // Conf: ValidExc
+   sa = Conf.lookup_string_array(conf_key_valid_exc);
+   for(i=0; i<sa.n_elements(); i++)
+      ValidExc.add(timestring_to_unix(sa[i].c_str()));
+
+   // Conf: WriteValid
+   sa = Conf.lookup_string_array(conf_key_write_valid);
+   for(i=0; i<sa.n_elements(); i++)
+      WriteValid.add(timestring_to_unix(sa[i].c_str()));
+
    // Conf: LeadReq
    sa = Conf.lookup_string_array(conf_key_lead_req);
    for(i=0; i<sa.n_elements(); i++){
diff --git a/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.h b/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.h
index ffbe04442a..cc2d218f9e 100644
--- a/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.h
+++ b/met/src/tools/tc_utils/tc_pairs/tc_pairs_conf_info.h
@@ -61,6 +61,9 @@ class TCPairsConfInfo {
       NumArray  InitHour;
       NumArray  LeadReq;
       unixtime  ValidBeg, ValidEnd;
+      TimeArray ValidInc;
+      TimeArray ValidExc;
+      TimeArray WriteValid;
 
       // Polyline masking regions
       ConcatString InitMaskName;
diff --git a/test/config/TCPairsConfig_ALAL2010 b/test/config/TCPairsConfig_ALAL2010
index 8468de62e7..790c667f28 100644
--- a/test/config/TCPairsConfig_ALAL2010
+++ b/test/config/TCPairsConfig_ALAL2010
@@ -50,10 +50,17 @@ init_inc = [];
 init_exc = [];
 
 //
-// Valid model time window
+// Valid model time windows to include or exclude
 //
 valid_beg = "";
 valid_end = "";
+valid_inc = [ ${VALID_INC} ];
+valid_exc = [];
+
+//
+// Valid times for which output should be written
+//
+write_valid = [ ${WRITE_VALID} ];
 
 //
 // Model initialization hours
diff --git a/test/config/TCPairsConfig_BASIN_MAP b/test/config/TCPairsConfig_BASIN_MAP
index a66e91e031..83bb491385 100644
--- a/test/config/TCPairsConfig_BASIN_MAP
+++ b/test/config/TCPairsConfig_BASIN_MAP
@@ -50,10 +50,17 @@ init_inc = [];
 init_exc = [];
 
 //
-// Valid model time window
+// Valid model time windows to include or exclude
 //
 valid_beg = "";
 valid_end = "";
+valid_inc = [];
+valid_exc = [];
+
+//
+// Valid times for which output should be written
+//
+write_valid = [];
 
 //
 // Model initialization hours
diff --git a/test/config/TCPairsConfig_INTERP12 b/test/config/TCPairsConfig_INTERP12
index 8183e6b73a..7c2301d50c 100644
--- a/test/config/TCPairsConfig_INTERP12
+++ b/test/config/TCPairsConfig_INTERP12
@@ -50,10 +50,17 @@ init_inc = [];
 init_exc = [];
 
 //
-// Valid model time window
+// Valid model time windows to include or exclude
 //
 valid_beg = "";
 valid_end = "";
+valid_inc = [];
+valid_exc = [];
+
+//
+// Valid times for which output should be written
+//
+write_valid = [];
 
 //
 // Model initialization hours
diff --git a/test/config/TCPairsConfig_PROBRIRW b/test/config/TCPairsConfig_PROBRIRW
index db7d915ee7..e75b3a3bd1 100644
--- a/test/config/TCPairsConfig_PROBRIRW
+++ b/test/config/TCPairsConfig_PROBRIRW
@@ -50,10 +50,17 @@ init_inc = [];
 init_exc = [];
 
 //
-// Valid model time window
+// Valid model time windows to include or exclude
 //
 valid_beg = "";
 valid_end = "";
+valid_inc = [];
+valid_exc = [];
+
+//
+// Valid times for which output should be written
+//
+write_valid = [ ${WRITE_VALID} ];
 
 //
 // Model initialization hours
diff --git a/test/xml/unit_tc_pairs.xml b/test/xml/unit_tc_pairs.xml
index 9ad3738de7..c0e22dcf37 100644
--- a/test/xml/unit_tc_pairs.xml
+++ b/test/xml/unit_tc_pairs.xml
@@ -16,7 +16,9 @@
 
   <test name="tc_pairs_ALAL2010">
     <env>
-      <pair><name>LEAD_REQ</name> <value></value></pair>
+      <pair><name>VALID_INC</name>   <value></value></pair>
+      <pair><name>WRITE_VALID</name> <value></value></pair>
+      <pair><name>LEAD_REQ</name>    <value></value></pair>
     </env>
     <exec>&MET_BIN;/tc_pairs</exec>
     <param> \
@@ -70,6 +72,9 @@
 
   <test name="tc_pairs_PROBRIRW">
     <exec>&MET_BIN;/tc_pairs</exec>
+    <env>
+      <pair><name>WRITE_VALID</name> <value></value></pair>
+    </env>
     <param> \
       -edeck &DATA_DIR;/edeck/*2015*.dat \
       -bdeck &DATA_DIR;/bdeck/*2015*.dat \
@@ -100,7 +105,9 @@
 
   <test name="tc_pairs_LEAD_REQ">
     <env>
-      <pair><name>LEAD_REQ</name> <value>"96","108"</value></pair>
+      <pair><name>VALID_INC</name>   <value></value></pair>
+      <pair><name>WRITE_VALID</name> <value></value></pair>
+      <pair><name>LEAD_REQ</name>    <value>"96","108"</value></pair>
     </env>
     <exec>&MET_BIN;/tc_pairs</exec>
     <param> \
@@ -116,4 +123,42 @@
     </output>
   </test>
 
+  <test name="tc_pairs_WRITE_VALID">
+    <env>
+      <pair><name>VALID_INC</name>   <value>"20100814_12", "20100815_00"</value></pair>
+      <pair><name>WRITE_VALID</name> <value>"20100815_00"</value></pair>
+      <pair><name>LEAD_REQ</name>    <value></value></pair>
+    </env>
+    <exec>&MET_BIN;/tc_pairs</exec>
+    <param> \
+      -adeck &DATA_DIR;/adeck/*2010.dat \
+      -bdeck &DATA_DIR;/bdeck/*2010.dat \
+      -config &CONFIG_DIR;/TCPairsConfig_ALAL2010 \
+      -out &OUTPUT_DIR;/tc_pairs/alal2010_WRITE_VALID \
+      -log &OUTPUT_DIR;/tc_pairs/tc_pairs_WRITE_VALID.log \
+      -v 2
+    </param>
+    <output>
+      <stat>&OUTPUT_DIR;/tc_pairs/alal2010_WRITE_VALID.tcst</stat>
+    </output>
+  </test>
+
+  <test name="tc_pairs_WRITE_VALID_PROBRIRW">
+    <exec>&MET_BIN;/tc_pairs</exec>
+    <env>
+      <pair><name>WRITE_VALID</name> <value>"20150512"</value></pair>
+    </env>
+    <param> \
+      -edeck &DATA_DIR;/edeck/*2015*.dat \
+      -bdeck &DATA_DIR;/bdeck/*2015*.dat \
+      -config &CONFIG_DIR;/TCPairsConfig_PROBRIRW \
+      -out &OUTPUT_DIR;/tc_pairs/alal2015_WRITE_VALID \
+      -log &OUTPUT_DIR;/tc_pairs/tc_pairs_WRITE_VALID_PROBRIRW.log \
+      -v 3
+    </param>
+    <output>
+      <stat>&OUTPUT_DIR;/tc_pairs/alal2015_WRITE_VALID_PROBRIRW.tcst</stat>
+    </output>
+  </test>
+
 </met_test>