Feature 1184 dryline

met/Make-include

@@ -33,6 +33,7 @@ MET_CPPFLAGS = -I${top_builddir}/src/basic/vx_cal \
 	     -I${top_builddir}/src/libcode/vx_solar \
 	     -I${top_builddir}/src/libcode/vx_statistics \
 	     -I${top_builddir}/src/libcode/vx_stat_out \
+	     -I${top_builddir}/src/libcode/vx_bool_calc \
 	     -I${top_builddir}/src/libcode/vx_summary \
 	     -I${top_builddir}/src/libcode/vx_time_series \
 	     -I${top_builddir}/src/libcode/vx_series_data \
@@ -71,6 +72,7 @@ MET_LDFLAGS = -L${top_builddir}/src/basic/vx_cal \
 	     -L${top_builddir}/src/libcode/vx_pxm \
 	     -L${top_builddir}/src/libcode/vx_render \
 	     -L${top_builddir}/src/libcode/vx_regrid \
+	     -L${top_builddir}/src/libcode/vx_bool_calc \
 	     -L${top_builddir}/src/libcode/vx_shapedata \
 	     -L${top_builddir}/src/libcode/vx_solar \
 	     -L${top_builddir}/src/libcode/vx_statistics \

met/configure.ac

@@ -1209,6 +1209,7 @@ AC_CONFIG_FILES([Makefile
                  src/libcode/vx_summary/Makefile
                  src/libcode/vx_python3_utils/Makefile
                  src/libcode/vx_data2d_python/Makefile
+                 src/libcode/vx_bool_calc/Makefile
                  src/tools/Makefile
                  src/tools/core/Makefile
                  src/tools/core/ensemble_stat/Makefile

met/data/config/MODEMultivarConfig_default

@@ -0,0 +1,266 @@
+////////////////////////////////////////////////////////////////////////////////
+//
+// MODE configuration file.
+//
+// For additional information, please see the MET User's Guide.
+//
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Output model name to be written
+//
+model = "WRF";
+
+//
+// Output description to be written
+//
+desc = "NA";
+
+//
+// MODE Multivar boolean combination logic
+//
+multivar_field_logic = "#1 && #2 && #3";
+
+//
+// Output observation type to be written
+//
+obtype = "PRECIP";
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Verification grid
+//
+regrid = {
+   to_grid = NONE;
+   method  = NEAREST;
+   width   = 1;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Approximate grid resolution (km)
+//
+grid_res = 20;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Run all permutations of radius and threshold
+//
+quilt = FALSE;
+
+//
+// Forecast and observation fields to be verified
+//
+fcst = {
+
+   // field = {
+   //    name  = "GAMMA";
+   //    level = "(*,*)";
+   // }
+
+   field = [
+
+      {
+         name  = "ALPHA";
+         level = "(*,*)";
+      }, 
+
+      {
+         name  = "BETA";
+         level = "(*,*)";
+      }, 
+
+      {
+         name  = "GAMMA";
+         level = "(*,*)";
+      }
+
+   ];
+
+   conv_radius        = 2; // in grid squares
+   conv_thresh        = >=5.0;
+   vld_thresh         = 0.5;
+   filter_attr_name   = [];
+   filter_attr_thresh = [];
+   merge_thresh       = >=3.5;
+   merge_flag         = NONE;
+}
+obs = fcst;
+// obs = {
+//    field = {
+//       name  = "APCP_";
+//       level = "A06";
+//    }
+// }
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Handle missing data
+//
+mask_missing_flag = NONE;
+
+//
+// Match objects between the forecast and observation fields
+//
+match_flag = NONE;
+
+//
+// Maximum centroid distance for objects to be compared
+//
+max_centroid_dist = 800.0/grid_res;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Verification masking regions
+//
+mask = {
+   grid      = "";
+   grid_flag = NONE; // Apply to NONE, FCST, OBS, or BOTH
+   poly      = "";
+   poly_flag = NONE; // Apply to NONE, FCST, OBS, or BOTH
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Fuzzy engine weights
+//
+weight = {
+   centroid_dist    = 2.0;
+   boundary_dist    = 4.0;
+   convex_hull_dist = 0.0;
+   angle_diff       = 1.0;
+   aspect_diff      = 0.0;
+   area_ratio       = 1.0;
+   int_area_ratio   = 2.0;
+   curvature_ratio  = 0.0;
+   complexity_ratio = 0.0;
+   inten_perc_ratio = 0.0;
+   inten_perc_value = 50;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Fuzzy engine interest functions
+//
+interest_function = {
+
+   centroid_dist = (
+      (            0.0, 1.0 )
+      (  60.0/grid_res, 1.0 )
+      ( 600.0/grid_res, 0.0 )
+   );
+
+   boundary_dist = (
+      (            0.0, 1.0 )
+      ( 400.0/grid_res, 0.0 )
+   );
+
+   convex_hull_dist = (
+      (            0.0, 1.0 )
+      ( 400.0/grid_res, 0.0 )
+   );
+
+   angle_diff = (
+      (  0.0, 1.0 )
+      ( 30.0, 1.0 )
+      ( 90.0, 0.0 )
+   );
+
+   aspect_diff = (
+      (  0.00, 1.0 )
+      (  0.10, 1.0 )
+      (  0.75, 0.0 )
+   );
+
+   corner   = 0.8;
+   ratio_if = (
+      (    0.0, 0.0 )
+      ( corner, 1.0 )
+      (    1.0, 1.0 )
+   );
+
+   area_ratio = ratio_if;
+
+   int_area_ratio = (
+      ( 0.00, 0.00 )
+      ( 0.10, 0.50 )
+      ( 0.25, 1.00 )
+      ( 1.00, 1.00 )
+   );
+
+   curvature_ratio = ratio_if;
+
+   complexity_ratio = ratio_if;
+
+   inten_perc_ratio = ratio_if;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Total interest threshold for determining matches
+//
+total_interest_thresh = 0.7;
+
+//
+// Interest threshold for printing output pair information
+//
+print_interest_thresh = 0.0;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Plotting information
+//
+met_data_dir = "MET_BASE";
+
+fcst_raw_plot = {
+   color_table      = "MET_BASE/colortables/met_default.ctable";
+   plot_min         = 0.0;
+   plot_max         = 25.4;
+   colorbar_spacing = 1;
+}
+
+obs_raw_plot = {
+   color_table      = "MET_BASE/colortables/met_default.ctable";
+   plot_min         = 0.0;
+   plot_max         = 25.4;
+   colorbar_spacing = 1;
+}
+
+object_plot = {
+   color_table      = "MET_BASE/colortables/mode_obj.ctable";
+}
+
+//
+// Boolean for plotting on the region of valid data within the domain
+//
+plot_valid_flag = FALSE;
+
+//
+// Plot polyline edges using great circle arcs instead of straight lines
+//
+plot_gcarc_flag = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// NetCDF matched pairs, PostScript, and contingency table output files
+//
+ps_plot_flag    = TRUE;
+nc_pairs_flag   = TRUE;
+ct_stats_flag   = TRUE;
+
+////////////////////////////////////////////////////////////////////////////////
+
+output_prefix  = "r_";
+version        = "V9.0";
+
+////////////////////////////////////////////////////////////////////////////////

met/data/config/Makefile.am

@@ -29,6 +29,7 @@ config_DATA = \
 	Madis2NcConfig_default \
 	MODEAnalysisConfig_default \
 	MODEConfig_default \
+	MODEMultivarConfig_default \
 	MTDConfig_default \
 	PB2NCConfig_default \
 	PointStatConfig_default \