dtcenter · JohnHalleyGotway · May 12, 2026 · May 12, 2026 · May 12, 2026 · May 13, 2026
diff --git a/.github/workflows/compilation_options.yml b/.github/workflows/compilation_options.yml
@@ -4,14 +4,7 @@ name: Compilation Options
 
 on:
 
-  # Note that scheduled cron events are run on the default branch.
-  # Enable this schedule when main_v12.0 becomes the default branch.
-
-  #  schedule:
-  #  - cron: '0 7 * * 0''
-
   push:
-
     tags:
       - '**' 
 

diff --git a/data/config/ConfigConstants b/data/config/ConfigConstants
@@ -161,6 +161,7 @@ AREA           = 3;
 // Point weight flag settings
 NONE           = 1;
 SID            = 2;
+KDE            = 3;
 
 // Duplicate flag settings
 NONE           = 1;

diff --git a/data/config/EnsembleStatConfig_default b/data/config/EnsembleStatConfig_default
@@ -323,6 +323,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "";
 version       = "V13.0.0";

diff --git a/data/config/PairStatConfig_default b/data/config/PairStatConfig_default
@@ -213,6 +213,10 @@ output_flag = {
 ////////////////////////////////////////////////////////////////////////////////
 
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 tmp_dir = "/tmp";
 version = "V13.0.0";

diff --git a/data/config/PointStatConfig_default b/data/config/PointStatConfig_default
@@ -325,7 +325,11 @@ output_flag = {
 ////////////////////////////////////////////////////////////////////////////////
 
 point_weight_flag = NONE;
-
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
+
 tmp_dir       = "/tmp";
 output_prefix = "";
 version       = "V13.0.0";

diff --git a/data/config/STATAnalysisConfig_CBS_Index b/data/config/STATAnalysisConfig_CBS_Index
@@ -93,7 +93,7 @@ alpha = [];
 line_type = [ "SL1L2" ];
 
 column    = [ "RMSE" ];
-            
+
 weight    = [ 6.4, 6.4, 6.4, 6.4, 6.4,
               2.4, 2.4, 2.4, 2.4, 2.4,
               2.4, 2.4, 2.4, 2.4, 2.4,

diff --git a/data/config/STATAnalysisConfig_GO_Index b/data/config/STATAnalysisConfig_GO_Index
@@ -97,7 +97,7 @@ alpha = [];
 line_type = [ "SL1L2" ];
 
 column    = [ "RMSE" ];
-            
+
 weight    = [ 4.0, 3.0, 2.0, 1.0,
               4.0, 3.0, 2.0, 1.0,
               4.0, 3.0, 2.0, 1.0,

diff --git a/docs/Users_Guide/config_options.rst b/docs/Users_Guide/config_options.rst
@@ -2803,10 +2803,49 @@ methods are planned for future versions:
 * SID to use the weights defined by the station ID masking configuration option,
   "mask.sid".
 
+* KDE to apply a kernel density estimator to compute weights based on observation
+  density to avoid the impact of oversampling from data rich regions. Note that
+  KDE weights are automatically rescaled a standard range of 1 to 100.
+
 .. code-block:: none
 
   point_weight_flag = NONE;
 
+kde_ref_angle
+-------------
+
+The "kde_ref_angle" entry defines the reference angle used when computing the
+weights for "point_weight_flag = KDE". This reference angle is defined in degrees
+with a default value of 0.75. As described in :ref:`Haiden et al., 2012 <Haiden-2012>`,
+an increase of this parameter increases the weights of data-sparse regions in global
+or continental-scale averages.
+
+MET checks for the divide-by-zero condition to avoid computing an infinite KDE weights.
+When that occurs, MET errors out and instructs the user to adjust the reference angle.
+
+.. code-block:: none
+
+  kde_ref_angle = 0.75;
+
+write_weights
+-------------
+
+Point weights are either pre-defined ("point_weight_flag = SID") or computed based
+on the station location density ("point_weight_flag = KDE). The "write_weights"
+entry is a boolean that can be set to TRUE or FALSE (default). If TRUE, the weights
+are written to an output ASCII file ending in "_point_weights.txt". That file can
+be used to define the "mask.sid" station ID masking file in future runs.
+
+Note that "mask.sid" can be specified separately for each verification task, and,
+technically, the same station could be assigned different weights for different
+verification tasks. If the same station name appears in multiple "mask.sid"
+lists, only the first weight encountered for each station is written to the
+output.
+
+.. code-block:: none
+
+  write_weights = FALSE;
+
 hss_ec_value
 ------------
 

diff --git a/docs/Users_Guide/ensemble-stat.rst b/docs/Users_Guide/ensemble-stat.rst
@@ -187,6 +187,8 @@ ____________________
   obtype_as_group_val_flag = FALSE;
   grid_weight_flag         = NONE;
   point_weight_flag        = NONE;
+  kde_ref_angle            = 0.75;
+  write_weights            = FALSE;
   output_prefix            = "";
   version                  = "VN.N";
 

diff --git a/docs/Users_Guide/pair-stat.rst b/docs/Users_Guide/pair-stat.rst
@@ -111,6 +111,8 @@ ________________________
   model             = "FCST";
   desc              = "NA";
   point_weight_flag = NONE;
+  kde_ref_angle     = 0.75;
+  write_weights     = FALSE;
   tmp_dir           = "/tmp";
   version           = "VN.N";
 

diff --git a/docs/Users_Guide/point-stat.rst b/docs/Users_Guide/point-stat.rst
@@ -411,6 +411,8 @@ ________________________
   message_type_group_map   = [ ... ];
   obtype_as_group_val_flag = FALSE;
   point_weight_flag        = NONE;
+  kde_ref_angle            = 0.75;
+  write_weights            = FALSE;
   tmp_dir                  = "/tmp";
   output_prefix            = "";
   version                  = "VN.N";

diff --git a/docs/Users_Guide/release-notes.rst b/docs/Users_Guide/release-notes.rst
@@ -168,6 +168,16 @@ MET Version 13.0.0 Upgrade Instructions
      However, if no spatial masking region is requested in the configuration file, then the
      "FULL" grid will automatically be added, as noted in :numref:`config_options-mask`.
 
+   * Point-Stat, Ensemble-Stat, and Pair-Stat configuration files
+
+     * The existing "point_weight_flag" entry can now be set to "KDE" to request kernel density
+       estimation weighting.
+
+     * The "kde_ref_angle" numeric entry defines a reference angle when computing KDE weights.
+
+     * The "write_weights" boolean entry instructs the tools to write an output text file
+       containing the point weight values applied based on the "point_weight_flag" setting.
+
    * Grid-Diag configuration file
 
      * The "mask.grid" and "mask.poly" entries have changed from strings to arrays of strings
@@ -190,9 +200,7 @@ MET Version 13.0.0 Upgrade Instructions
        for a description of wind rotation and derivation and :numref:`config_wind_field_names` for the
        corresponding configuration options.
 
-.. dropdown:: Output format changes - NONE
-
-.. dropdown:: Output format changes - NONE
+.. dropdown:: Output format changes
 
    MET version 13.0.0 adds or modifies the following output file formats:
 
@@ -210,7 +218,7 @@ MET Version 13.0.0 Upgrade Instructions
      * Adds new output variables with names include "EAS" and "EAS_WIDTH" for the Ensemble
        Agreement Scale algorithm.
 
-.. dropdown:: Output data changes - NONE
+.. dropdown:: Output data changes
 
    MET version 13.0.0 modifies existing output data values in the following ways:
 

diff --git a/internal/test_unit/config/EnsembleStatConfig b/internal/test_unit/config/EnsembleStatConfig
@@ -228,6 +228,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_LAND_TOPO_MASK b/internal/test_unit/config/EnsembleStatConfig_LAND_TOPO_MASK
@@ -290,6 +290,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_LAPSE_RATE_MSL_AGL b/internal/test_unit/config/EnsembleStatConfig_LAPSE_RATE_MSL_AGL
@@ -307,6 +307,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_MASK_SID b/internal/test_unit/config/EnsembleStatConfig_MASK_SID
@@ -220,6 +220,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_climo b/internal/test_unit/config/EnsembleStatConfig_climo
@@ -250,6 +250,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_grid_weight b/internal/test_unit/config/EnsembleStatConfig_grid_weight
@@ -242,6 +242,10 @@ rng = {
 
 grid_weight_flag  = ${GRID_WEIGHT};
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_one_cdf_bin b/internal/test_unit/config/EnsembleStatConfig_one_cdf_bin
@@ -234,6 +234,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_point_weight b/internal/test_unit/config/EnsembleStatConfig_point_weight
@@ -133,9 +133,9 @@ obs_window = {
 // Verification masking regions
 //
 mask = {
-   grid  = [];
+   grid  = [ ${MASK_GRID} ];
    poly  = [];
-   sid   = [ "${CONFIG_DIR}/SID_CONUS_ADPSFC_ELEV.txt" ];
+   sid   = [ ${MASK_SID} ];
    llpnt = [];
 }
 
@@ -207,6 +207,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = ${POINT_WEIGHT};
+kde_ref_angle     = 0.75;
+write_weights     = TRUE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_python b/internal/test_unit/config/EnsembleStatConfig_python
@@ -231,6 +231,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_qty_inc_exc b/internal/test_unit/config/EnsembleStatConfig_qty_inc_exc
@@ -88,7 +88,7 @@ fcst = {
          obs_quality_exc = [];
          desc            = "QTY_INC";
       },
-      
+
       {
          name            = "TMP";
          level           = [ "Z2" ];
@@ -237,6 +237,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_rps_climo_bin_prob b/internal/test_unit/config/EnsembleStatConfig_rps_climo_bin_prob
@@ -218,6 +218,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_single_file_grib b/internal/test_unit/config/EnsembleStatConfig_single_file_grib
@@ -280,6 +280,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/EnsembleStatConfig_single_file_nc b/internal/test_unit/config/EnsembleStatConfig_single_file_nc
@@ -285,6 +285,10 @@ rng = {
 
 grid_weight_flag  = NONE;
 point_weight_flag = NONE;
+kde_ref_angle     = 0.75;
+write_weights     = FALSE;
+
+////////////////////////////////////////////////////////////////////////////////
 
 output_prefix = "${OUTPUT_PREFIX}";
 

diff --git a/internal/test_unit/config/GridStatConfig_fbias_perc_thresh b/internal/test_unit/config/GridStatConfig_fbias_perc_thresh
@@ -63,7 +63,7 @@ fcst = {
    field = [
       { name = "TMP";   level = [ "P500" ];  cat_thresh = [ ==FBIAS0.9, ==FBIAS1.0, ==FBIAS1.1, ==FBIAS0.9, ==FBIAS1.0, ==FBIAS1.1, <243, <243, <243, >253, >253, >253 ]; } 
    ];
-   
+
 }
 obs = {
 

diff --git a/internal/test_unit/config/GridStatConfig_ugrid_mpas b/internal/test_unit/config/GridStatConfig_ugrid_mpas
@@ -245,8 +245,8 @@ nc_pairs_flag = {
 seeps_p1_thresh = NA;
 
 ////////////////////////////////////////////////////////////////////////////////
-// Unstructured Grid (optional ecept ugrid_dataset)
 
+// Unstructured Grid (optional except ugrid_dataset)
 ugrid_dataset = "mpas";
 //ugrid_max_distance_km = 30;
 //ugrid_map_config = "";

diff --git a/internal/test_unit/config/GridStatConfig_ugrid_mpas_diag b/internal/test_unit/config/GridStatConfig_ugrid_mpas_diag
@@ -245,8 +245,8 @@ nc_pairs_flag = {
 seeps_p1_thresh = NA;
 
 ////////////////////////////////////////////////////////////////////////////////
-// Unstructured Grid (optional ecept ugrid_dataset)
 
+// Unstructured Grid (optional except ugrid_dataset)
 ugrid_dataset = "mpas";
 //ugrid_max_distance_km = 30;
 //ugrid_map_config = "";

diff --git a/internal/test_unit/config/MODEConfig_multivar b/internal/test_unit/config/MODEConfig_multivar
@@ -75,10 +75,10 @@ fcst = {
    //${METPLUS_FCST_VLD_THRESH}
    //${METPLUS_FCST_FILTER_ATTR_NAME}
    filter_attr_name = ["AREA"];
-   
+
    //${METPLUS_FCST_FILTER_ATTR_THRESH}
    filter_attr_thresh = [>=25];
-   
+
    //${METPLUS_FCST_MERGE_THRESH}
    //${METPLUS_FCST_MERGE_FLAG}
    //${METPLUS_FCST_FILE_TYPE}
-Original file line number
+Diff line change
@@ Expand Up / @@ -63,7 +63,7 @@ fcst = { @@
        field = [
           { name = "TMP";   level = [ "P500" ];  cat_thresh = [ ==FBIAS0.9, ==FBIAS1.0, ==FBIAS1.1, ==FBIAS0.9, ==FBIAS1.0, ==FBIAS1.1, <243, <243, <243, >253, >253, >253 ]; }
        ];
     }
     obs = {
@@ Expand Down @@