Per #1020, add new test in unit_grid_stat.xml to exercise the set_att…

…rs functionality.

test/config/GridStatConfig_GRIB_set_attrs

@@ -0,0 +1,240 @@
+///////////////////////////////////////////////////////////////////////////////
+//
+// Grid-Stat configuration file.
+//
+// For additional information, see the MET_BASE/config/README file.
+//
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Output model name to be written
+//
+model = "GALWEM";
+
+//
+// Output description to be written
+// May be set separately in each "obs.field" entry
+//
+desc = "NA";
+
+//
+// Output observation type to be written
+//
+obtype = "GPP";
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Verification grid
+//
+regrid = {
+   to_grid    = OBS;
+   method     = BILIN;
+   width      = 2;
+   vld_thresh = 0.5;
+   shape      = SQUARE;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// May be set separately in each "field" entry
+//
+censor_thresh       = [];
+censor_val          = [];
+cat_thresh          = [];
+cnt_thresh          = [ NA ];
+cnt_logic           = UNION;
+wind_thresh         = [ NA ];
+wind_logic          = UNION;
+eclv_points         = 0.05;
+nc_pairs_var_name   = "";
+nc_pairs_var_suffix = "";
+rank_corr_flag      = FALSE;
+
+//
+// Forecast and observation fields to be verified
+//
+
+set_attrs = {
+  set_init  = "20170807_000000";
+  set_valid = "20170807_120000";
+  set_lead  = "12";
+}
+
+fcst = {
+
+   set_attrs = {
+     set_name  = "Temperature";
+     set_units = "Kelvin";
+   }
+
+   field = [
+      {
+        //tropopause
+        name = "TMP";
+        GRIB_lvl_typ = 7;
+        GRIB_lvl_val1 = 0;
+        set_attrs = {
+          set_level     = "TROPO";
+          set_long_name = "Temperature at the tropopause";
+        }
+      },
+      {
+        //max wind
+        name = "TMP";
+        GRIB_lvl_typ = 6;
+        GRIB_lvl_val1 = 0;
+        set_attrs = {
+          set_level     = "MAXWIND";
+          set_long_name = "Temperature at the maximum wind level";
+        }
+      }
+   ];
+
+}
+obs = fcst;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Verification masking regions
+//
+mask = {
+   grid = [ "FULL" ];
+   poly = [];
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Confidence interval settings
+//
+ci_alpha  = [ 0.05 ];
+
+boot = {
+   interval = PCTILE;
+   rep_prop = 1.0;
+   n_rep    = 0;
+   rng      = "mt19937";
+   seed     = "";
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Data smoothing methods
+//
+interp = {
+   field      = BOTH;
+   vld_thresh = 1.0;
+   shape      = SQUARE;
+
+   type = [
+      {
+         method = NEAREST;
+         width  = 1;
+      }
+   ];
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Neighborhood methods
+//
+nbrhd = {
+   width      = [ 1 ];
+   cov_thresh = [ >=0.5 ];
+   vld_thresh = 1.0;
+   shape      = SQUARE;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Fourier decomposition
+//
+fourier = {
+   wave_1d_beg = [];
+   wave_1d_end = [];
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Gradient statistics
+// May be set separately in each "obs.field" entry
+//
+gradient = {
+   dx = [ 1 ];
+   dy = [ 1 ];
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Distance Map statistics
+// May be set separately in each "obs.field" entry
+//
+distance_map = {
+   baddeley_p        = 2;
+   baddeley_max_dist = NA;
+   fom_alpha         = 0.1;
+   zhu_weight        = 0.5;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Statistical output types
+//
+output_flag = {
+   fho    = NONE;
+   ctc    = NONE;
+   cts    = NONE;
+   mctc   = NONE;
+   mcts   = NONE;
+   cnt    = NONE;
+   sl1l2  = STAT;
+   sal1l2 = NONE;
+   vl1l2  = NONE;
+   val1l2 = NONE;
+   vcnt   = NONE;
+   pct    = NONE;
+   pstd   = NONE;
+   pjc    = NONE;
+   prc    = NONE;
+   eclv   = NONE;
+   nbrctc = NONE;
+   nbrcts = NONE;
+   nbrcnt = NONE;
+   grad   = NONE;
+   dmap   = NONE;
+};
+
+//
+// NetCDF matched pairs output file
+//
+nc_pairs_flag = {
+   latlon       = FALSE;
+   raw          = TRUE;
+   diff         = FALSE;
+   climo        = FALSE;
+   climo_cdp    = FALSE;
+   weight       = FALSE;
+   nbrhd        = FALSE;
+   fourier      = FALSE;
+   gradient     = FALSE;
+   distance_map = FALSE;
+   apply_mask   = FALSE;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+grid_weight_flag = NONE;
+tmp_dir          = "/tmp";
+output_prefix    = "GRIB_set_attrs";
+version          = "V9.1";
+
+////////////////////////////////////////////////////////////////////////////////

test/xml/unit_grid_stat.xml

@@ -33,6 +33,23 @@
     </output>
   </test>
 
+  <test name="grid_stat_GRIB_set_attrs">
+    <exec>&MET_BIN;/grid_stat</exec>
+    <env>
+      <pair><name>OUTPUT_PREFIX</name>  <value>GRIB_set_attrs</value></pair>
+    </env>
+    <param> \
+      &DATA_DIR_MODEL;/grib2/galwem/GALWEM_20170610_CY00_FH18.GR2\
+      &DATA_DIR_MODEL;/grib2/galwem-gpp/PS.557WW_SC.U_DI.F_GP.GALWEM-GD_GR.0P25DEG_AR.GLOBAL_DD.20170610_CY.00_FH.018_DF.GR2 \
+      &CONFIG_DIR;/GridStatConfig_GRIB_set_attrs \
+      -outdir &OUTPUT_DIR;/grid_stat -v 1
+    </param>
+    <output>
+      <stat>&OUTPUT_DIR;/grid_stat/grid_stat_GRIB_set_attrs_120000L_20170807_120000V.stat</stat>
+      <grid_nc>&OUTPUT_DIR;/grid_stat/grid_stat_GRIB_set_attrs_120000L_20170807_120000V_pairs.nc</grid_nc>
+    </output>
+  </test>
+
   <test name="grid_stat_GRIB2_NAM_RTMA">
     <exec>&MET_BIN;/grid_stat</exec>
     <env>