Feature 1184 dryline

.github/workflows/testing.yml

@@ -148,7 +148,7 @@ jobs:
       matrix:
         tests:
           - 'ascii2nc_indy pb2nc_indy tc_dland tc_pairs tc_stat plot_tc tc_rmw rmw_analysis tc_gen'
-          - 'met_test_scripts mode_graphics mtd regrid airnow gsi_tools netcdf modis series_analysis gen_ens_prod wwmca_regrid gen_vx_mask grid_weight interp_shape grid_diag grib_tables lidar2nc shift_data_plane trmm2nc aeronet wwmca_plot ioda2nc gaussian'
+          - 'met_test_scripts mode_multivar mode_graphics mtd regrid airnow gsi_tools netcdf modis series_analysis gen_ens_prod wwmca_regrid gen_vx_mask grid_weight interp_shape grid_diag grib_tables lidar2nc shift_data_plane trmm2nc aeronet wwmca_plot ioda2nc gaussian'
       fail-fast: false
     steps:
       - uses: actions/checkout@v2

met/Make-include

@@ -34,6 +34,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 \
@@ -73,6 +74,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

@@ -1240,6 +1240,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/libcode/vx_pointdata_python/Makefile
                  src/tools/Makefile
                  src/tools/core/Makefile

met/data/config/MODEConfig_default

@@ -25,7 +25,6 @@ obtype = "ANALYS";
 
 //
 // Verification grid
-// May be set separately in each "field" entry
 //
 regrid = {
    to_grid    = NONE;

met/data/config/MODEMultivarConfig_default

@@ -0,0 +1,259 @@
+////////////////////////////////////////////////////////////////////////////////
+//
+// 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";
+
+//
+// Output observation type to be written
+//
+obtype = "ANALYS";
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Verification grid
+//
+regrid = {
+   to_grid    = NONE;
+   method     = NEAREST;
+   width      = 1;
+   vld_thresh = 0.5;
+   shape      = SQUARE;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Approximate grid resolution (km)
+//
+grid_res = 20;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// Run all permutations of radius and threshold
+//
+quilt = FALSE;
+
+//
+// MODE Multivar boolean combination logic
+//
+multivar_logic = "#1 && #2 && #3";
+
+//
+// Forecast and observation fields to be verified
+//
+fcst = {
+
+   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;
+
+////////////////////////////////////////////////////////////////////////////////
+
+//
+// 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    = 0.0;
+}
+
+obs_raw_plot = {
+   color_table = "MET_BASE/colortables/met_default.ctable";
+   plot_min    = 0.0;
+   plot_max    = 0.0;
+}
+
+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;
+
+////////////////////////////////////////////////////////////////////////////////
+
+shift_right = 0;   //  grid squares
+
+////////////////////////////////////////////////////////////////////////////////
+
+output_prefix  = "";
+version        = "V10.1.0";
+
+////////////////////////////////////////////////////////////////////////////////

met/data/config/Makefile.am

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