////////////////////////////////////////////////////////////////////////////////
//
// Grid-Stat configuration file.
//
// For additional information, please see the MET User's Guide.
//
////////////////////////////////////////////////////////////////////////////////

//
// Output model name to be written
//
model = "GFS";

//
// Output description to be written
// May be set separately in each "obs.field" entry
//
desc = "NA";

//
// Output observation type to be written
//
obtype = "GFSANL";

////////////////////////////////////////////////////////////////////////////////

//
// Verification grid
//
regrid = {
   to_grid    = "${CLIMO_DIR}/mslp_mean.grib";
   method     = BILIN;
   width      = 2;
   vld_thresh = 0.5;
   shape      = SQUARE;
}

////////////////////////////////////////////////////////////////////////////////

//
// May be set separately in each "field" entry
//
censor_thresh       = [];
censor_val          = [];
mpr_column          = [];
mpr_thresh          = [];
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 = "";
hss_ec_value        = NA;
rank_corr_flag      = FALSE;

//
// Forecast and observation fields to be verified
//
fcst = {
   field = [
      { name = "TMP";  level = [ "P500" ]; },
      { name = "UGRD"; level = [ "P500" ]; },
      { name = "VGRD"; level = [ "P500" ]; },
      { name = "WIND"; level = [ "P500" ]; },
      { name = "TMP";  level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "UGRD"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "VGRD"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "WIND"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; }
   ];
}
obs = {

   cnt_thresh = [ NA, <CDP50, >=CDP50, ==CDP25 ];

   field = [
      { name = "TMP";  level = [ "P500" ]; cnt_thresh = [ NA ]; wind_thresh = [ NA ]; },
      { name = "UGRD"; level = [ "P500" ]; cnt_thresh = [ NA ]; wind_thresh = [ NA ]; },
      { name = "VGRD"; level = [ "P500" ]; cnt_thresh = [ NA ]; wind_thresh = [ NA ]; },
      { name = "WIND"; level = [ "P500" ]; cnt_thresh = [ NA ]; wind_thresh = [ NA ]; },
      { name = "TMP";  level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "UGRD"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "VGRD"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; },
      { name = "WIND"; level = [ "P850" ]; cat_thresh = [ >CDP90 ]; }
   ];
}

////////////////////////////////////////////////////////////////////////////////

//
// Climatology data
//
climo_mean = fcst;
climo_mean = {

   file_name = [ "${CLIMO_DIR}/t500hPa_mean.grib",
                 "${CLIMO_DIR}/t850hPa_mean.grib",
                 "${CLIMO_DIR}/u500hPa_mean.grib",
                 "${CLIMO_DIR}/u850hPa_mean.grib",
                 "${CLIMO_DIR}/v500hPa_mean.grib",
                 "${CLIMO_DIR}/v850hPa_mean.grib",
                 "${CLIMO_DIR}/ws500hPa_mean.grib",
                 "${CLIMO_DIR}/ws850hPa_mean.grib" ];

   regrid = {
      method     = NONE;
      width      = 1;
      vld_thresh = 0.5;
   }

   time_interp_method = DW_MEAN;
   day_interval       = 1;
   hour_interval      = 12;
}

climo_stdev = climo_mean;
climo_stdev = {

   file_name = [ "${CLIMO_DIR}/t850hPa_stdev.grib",
                 "${CLIMO_DIR}/u850hPa_stdev.grib",
                 "${CLIMO_DIR}/v850hPa_stdev.grib",
                 "${CLIMO_DIR}/ws850hPa_stdev.grib" ];
}

//
// May be set separately in each "obs.field" entry
//
climo_cdf = {
   cdf_bins    = 3;
   center_bins = TRUE;
   write_bins  = FALSE;
   direct_prob = FALSE;
}

////////////////////////////////////////////////////////////////////////////////

//
// 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, 2, 5 ];
   dy = [ 1, 3, 5 ];
}

////////////////////////////////////////////////////////////////////////////////

//
// 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;
   beta_value(n)     = n * n / 2.0;
}

////////////////////////////////////////////////////////////////////////////////

//
// Statistical output types
//
output_flag = {
   fho    = NONE;
   ctc    = STAT;
   cts    = NONE;
   mctc   = NONE;
   mcts   = NONE;
   cnt    = NONE;
   sl1l2  = STAT;
   sal1l2 = STAT;
   vl1l2  = STAT;
   val1l2 = STAT;
   vcnt   = STAT;
   pct    = NONE;
   pstd   = NONE;
   pjc    = NONE;
   prc    = NONE;
   eclv   = NONE;
   nbrctc = NONE;
   nbrcts = NONE;
   nbrcnt = NONE;
   grad   = STAT;
   dmap   = NONE;
   seeps  = NONE;
}

//
// NetCDF matched pairs output file
//
nc_pairs_flag = {
   latlon       = TRUE;
   raw          = TRUE;
   diff         = TRUE;
   climo        = TRUE;
   climo_cdp    = TRUE;
   weight       = TRUE;
   nbrhd        = FALSE;
   fourier      = FALSE;
   gradient     = TRUE;
   distance_map = FALSE;
   apply_mask   = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
// Threshold for SEEPS p1 (Probability of being dry)

seeps_p1_thresh = NA;

////////////////////////////////////////////////////////////////////////////////

grid_weight_flag = COS_LAT;
tmp_dir          = "/tmp";
output_prefix    = "${OUTPUT_PREFIX}";
version          = "V11.1.0";

////////////////////////////////////////////////////////////////////////////////
