////////////////////////////////////////////////////////////////////////////////
//
// Default point_stat configuration file
//
////////////////////////////////////////////////////////////////////////////////

//
// Specify a name to designate the model being verified.  This name will be
// written to the second column of the ASCII output generated.
//
model = "WRF";

//
// Specify a comma-separated list of fields to be verified.  Each field is
// specified as a grib code or corresponding grib code abbreviation followed
// by an accumulation or vertical level indicator.
//
// Each verification field is specified as one of the following:
//    GC/ANNN for accumulation interval NNN
//    GC/ZNNN for vertical level NNN (may only be set to 0, 2, or 10)
//    GC/PNNN for pressure level NNN in hPa
//    GC/PNNN-NNN for a range of pressure levels in hPa
//    Where GC is the number of or abbreviation for the grib code
//    to be verified.
//
//    NOTE: To verify winds as vectors rather than scalars,
//          specify UGRD (or 33) followd by VGRD (or 34) with the
//          same level values.
//
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
// e.g. vx_grib_code[] = [ "SPFH/P500", "TMP/P500", "HGT/P500", "UGRD/P500", "VGRD/P500" ];
//
vx_grib_code[] = [ "61/A0" ];

//
// Specify a comma-separated list of groups of thresholds to be applied to the
// verification fields listed above.  At least one threshold must be provided
// for each verification field listed above.  The lengths of the "vx_grib_code"
// array and the "thresholds" array must be the same.  To apply  multiple
// thresholds to a verification field, separate the threshold values with a
// space.
//
// Each threshold must be preceded by a two letter indicator for the type of
// thresholding to be performed:
//    'lt' for less than     'le' for less than or equal to
//    'eq' for equal to      'ne' for not equal to
//    'gt' for greater than  'ge' for greater than or equal to
//
// e.g. thresholds[] = [ "gt80", "gt0", "gt300", "gt5", "gt5" ];
//
thresholds[] = [ "ge0 ge0.5 ge1 ge5 ge10 ge20 ge30" ];

//
// Specify a comma-separated list of PrepBufr message types with which
// to perform the verification.  Statistics will be computed separately
// for each message type specified.  At least one PrepBufr message type
// must be provided.
// List of valid message types:
//    ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
//    MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
//    SFCSHP SPSSMI SYNDAT VADWND
//    ANYAIR (= AIRCAR, AIRCFT)
//    ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
//    ONLYSF (= ADPSFC, SFCSHP)
// http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
//
// e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
//
message_type[] = [ "ADPSFC" ];

//
// Specify a comma-separated list of grids to be used in masking the data over
// which to perform scoring.  An empty list indicates that no masking grid
// should be performed.  The standard NCEP grids are named "GNNN" where NNN
// indicates the three digit grid number.  Enter "FULL" to score over the
// entire domain.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
//
// e.g. mask_grids[] = [ "FULL" ];
//
mask_grids[] = [ "FULL" ];

//
// Specify a comma-separated list of ASCII files containing lat/lon polygons to
// be used in masking the data over which to perform scoring.  An empty list
// indicates that no polygon mask should be used.
//
// Latitude values are given in degrees north and longitude values are given in
// degrees east.  By default, the first and last points defined for the polygon
// are connected.
//
// The entries in the list should be files containing a name for the polygon
// followed by a space-separated list of lat/lon points defining the polygon:
//    "poly_name lat1 lon1 lat2 lon2... latn lonn"
// Where n is the number of polygon points.
//
// MET_BASE may be used in the path for the lat/lon polygon files.
//
// e.g. mask_polys[] = [ "MET_BASE/data/poly/EAST.poly",
//                       "MET_BASE/data/poly/WEST.poly" ];
//
mask_polys[] = [];

//
// Specify the name of an ASCII file containing a space-separated list of
// station ID's at which to perform verification.  Each station ID specified
// is treated as an individual masking region.
//
// An empty list file name indicates that no station ID masks should be used.
//
// MET_BASE may be used in the path for the station ID mask file name.
//
// e.g. mask_stations = "MET_BASE/data/stations/CONUS.stations";
//
//mask_stations = "/home/arc/VERIFICA/list_id_spazio.txt";
mask_stations = "";

//
// Specify a comma-separated list of values for alpha to be used when computing
// confidence intervals.  Values of alpha must be between 0 and 1.
//
// e.g. ci_alpha[] = [ 0.05, 0.10 ];
//
ci_alpha[] = [ 0.025,0.05,0.10 ];

//
// Specify the method to be used for computing bootstrap confidence intervals.
// The value for this is interpreted as follows:
//    (0) Use the BCa interval method (computationally intensive)
//    (1) Use the percentile interval method
//
boot_interval = 1;

//
// Specify a proportion between 0 and 1 to define the replicate sample size
// to be used when computing percentile intervals.  The replicate sample
// size is set to boot_rep_prop * n, where n is the number of raw data points.
//
// e.g boot_rep_prop = 0.80;
//
boot_rep_prop = 1.0;

//
// Specify the number of times each set of matched pair data should be
// resampled when computing bootstrap confidence intervals.  A value of
// zero disables the computation of bootstrap condifence intervals.
//
// e.g. n_boot_rep = 1000;
//
n_boot_rep = 1000;

//
// Specify the name of the random number generator to be used.  See the MET
// Users Guide for a list of possible random number generators.
//
boot_rng = "mt19937";

//
// Specify the seed value to be used when computing bootstrap confidence
// intervals.  If left unspecified, the seed will change for each run and
// the computed bootstrap confidence intervals will not be reproducable.
//
boot_seed = "3461077832";

//
// Specify a comma-separated list of interpolation method(s) to be used
// for comparing the forecast grid to the observation points.  String values
// are interpreted as follows:
//    MIN     = Minimum in the neighborhood
//    MAX     = Maximum in the neighborhood
//    MEDIAN  = Median in the neighborhood
//    UW_MEAN = Unweighted mean in the neighborhood
//    DW_MEAN = Distance-weighted mean in the neighborhood
//    LS_FIT  = Least-squares fit in the neighborhood
//
// In all cases, vertical interpolation is performed in the natural log
// of pressure of the levels above and below the observation.
//
// e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
//
interp_method[] = [ "DW_MEAN" ];

//
// Specify a comma-separated list of box widths to be used by the
// interpolation techniques listed above.  A value of 1 indicates that
// the nearest neighbor approach should be used.  For a value of n
// greater than 1, the n*n grid points closest to the observation define
// the neighborhood.
//
// e.g. interp_width = [ 1, 3, 5 ];
//
interp_width[] = [ 3 ];

//
// When interpolating, compute a ratio of the number of valid data points
// to the total number of points in the neighborhood.  If that ratio is
// less than this threshold, do not include the observation.  This
// threshold must be between 0 and 1.  Setting this threshold to 1 will
// require that each observation be surrounded by n*n valid forecast
// points.
//
// e.g. interp_threshold = 1.0;
//
interp_threshold = 0.1;

//
// Specify flags to indicate the type of data to be output:
//    (1) VSDB and FHO Text Files, FHO rates:
//           Total (TOTAL),
//           Forecast Rate (F_RATE),
//           Hit Rate (H_RATE),
//           Observation Rate (O_RATE)
//
//    (2) VSDB and CTC Text Files, Contingency Table Counts:
//           Total (TOTAL),
//           Forecast Yes and Observation Yes Count (FY_OY),
//           Forecast Yes and Observation No Count (FY_ON),
//           Forecast No and Observation Yes Count (FN_OY),
//           Forecast No and Observation No Count (FN_ON)
//
//    (3) VSDB and CTP Text Files, Contingency Table Proportions:
//           Total (TOTAL),
//           FY_OY/Total (FY_OY_TP),
//           FY_ON/Total (FY_ON_TP),
//           FN_OY/Total (FN_OY_TP),
//           FN_ON/Total (FN_ON_TP),
//           FY/Total    (FY_TP),
//           FN/Total    (FN_TP),
//           OY/Total    (OY_TP),
//           ON/Total    (ON_TP)
//
//    (4) VSDB and CFP Text Files, Contingency Table Forecast Proportions:
//           Total (TOTAL),
//           FY_OY/FY (FY_OY_FP),
//           FY_ON/FY (FY_ON_FP),
//           FN_OY/FN (FN_OY_FP),
//           FN_ON/FN (FN_ON_FP),
//           FY       (FY),
//           FN       (FN)
//
//    (5) VSDB and COP Text Files, Contingency Table Observation Proportions:
//           Total (TOTAL),
//           FY_OY/OY (FY_OY_OP),
//           FY_ON/ON (FY_ON_OP),
//           FN_OY/OY (FN_OY_OP),
//           FN_ON/ON (FN_ON_OP),
//           OY       (OY),
//           ON       (ON)
//
//    (6) VSDB and CTS Text Files, Contingency Table Scores:
//           Total (TOTAL),
//           Base Rate (BASER), BASER_CL, BASER_CU,
//           Forecast Mean (FMEAN), FMEAN_CL, FMEAN_CU,
//           Accuracy (ACC), ACC_CL, ACC_CU,
//           Frequency Bias (FBIAS),
//           Probability of Detecting Yes (PODY), PODY_CL, PODY_CU,
//           Probability of Detecting No (PODN), PODN_CL, PODN_CU,
//           Probability of False Detection (POFD), POFD_CL, POFD_CU,
//           False Alarm Ratio (FAR), FAR_CL, FAR_CU,
//           Critical Success Index (CSI), CSI_CL, CSI_CU,
//           Gilbert Skill Score (GSS),
//           Hanssen and Kuipers Discriminant (HK), HK_CL, HK_CU,
//           Heidke Skill Score (HSS),
//           Odds Ratio (ODDS), ODDS_CL, ODDS_CU
//
//    (7) VSDB and CNT Text Files, Statistics of Continuous Variables:
//           Total (TOTAL),
//           Forecast Mean (FBAR), FBAR_CL, FBAR_CU,
//           Forecast Standard Deviation (FSTDEV), FSTDEV_CL, FSTDEV_CU
//           Observation Mean (OBAR), OBAR_CL, OBAR_CU,
//           Observation Standard Deviation (OSTDEV), OSTDEV_CL, OSTDEV_CU,
//           Pearson's Correlation Coefficient (PR_CORR), PR_CORR_CL, PR_CORR_CU,
//           Spearman's Rank Correlation Coefficient (SP_CORR),
//           Kendall Tau Rank Correlation Coefficient (KT_CORR),
//           Number of ranks compared (RANKS),
//           Number of tied ranks in the forecast field (FRANK_TIES),
//           Number of tied ranks in the observation field (ORANK_TIES),
//           Mean Error (ME), ME_CL, ME_CU,
//           Standard Deviation of the Error (ESTDEV), ESTDEV_CL, ESTDEV_CU,
//           Bias (BIAS = FBAR - OBAR),
//           Mean Absolute Error (MAE),
//           Mean Squared Error (MSE),
//           Bias-Corrected Mean Squared Error (BCMSE),
//           Root Mean Squared Error (RMSE),
//           Percentiles of the Error (E10, E25, E50, E75, E90)
//
//           NOTE: CL and CU values define lower and upper
//                 confidence interval limits.
//
//    (8) VSDB and SL1L2 Text Files, SL1L2
//           Total (TOTAL),
//           Forecast Mean (FBAR),
//              = mean(f)
//           Observation Mean (OBAR),
//              = mean(o)
//           Forecast*Observation Product Mean (FOBAR),
//              = mean(f*o)
//           Forecast Squared Mean (FFBAR),
//              = mean(f^2)
//           Observation Squared Mean (OOBAR)
//              = mean(o^2)
//
//    (9) VSDB and SAL1L2 Text Files, SAL1L2
//           Total (TOTAL),
//           Forecast Anomaly Mean (FABAR),
//              = mean(f-c)
//           Observation Anomaly Mean (OABAR),
//              = mean(o-c)
//           Product of Forecast and Observation Anomalies Mean (FOABAR),
//              = mean((f-c)*(o-c))
//           Forecast Anomaly Squared Mean (FFABAR),
//              = mean((f-c)^2)
//           Observation Anomaly Squared Mean (OOABAR)
//              = mean((o-c)^2)
//
//    (10) VSDB and VL1L2 Text Files, VL1L2
//           Total (TOTAL),
//           U-Forecast Mean (UFBAR),
//              = mean(uf)
//           V-Forecast Mean (VFBAR),
//              = mean(vf)
//           U-Observation Mean (UOBAR),
//              = mean(uo)
//           V-Observation Mean (VOBAR),
//              = mean(vo)
//           U-Product Plus V-Product (UVFOBAR),
//              = mean(uf*uo+vf*vo)
//           U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
//              = mean(uf^2+vf^2)
//           U-Observation Squared Plus V-Observation Squared (UVOOBAR)
//              = mean(uo^2+vo^2)
//
//    (11) VSDB and VAL1L2 Text Files, VAL1L2
//           U-Forecast Anomaly Mean (UFABAR),
//              = mean(uf-uc)
//           V-Forecast Anomaly Mean (VFABAR),
//              = mean(vf-vc)
//           U-Observation Anomaly Mean (UOABAR),
//              = mean(uo-uc)
//           V-Observation Anomaly Mean (VOABAR),
//              = mean(vo-vc)
//           U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
//              = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
//           U-Forecast Anomaly Squared Plus V-Forecast Anomaly Squared (UVFFABAR),
//              = mean((uf-uc)^2+(vf-vc)^2)
//           U-Observation Anomaly Squared Plus V-Observation Anomaly Squared (UVOOABAR)
//              = mean((uo-uc)^2+(vo-vc)^2)
//
//    (12) MPR Text File, Matched Pairs Data:
//           Total (TOTAL),
//           Index (INDEX),
//           Latitude (LAT),
//           Longitude (LON),
//           Level (LEVEL),
//           Forecast Value (FCST),
//           Observation Value (OBS),
//           Climatological Value (CLIMO),
//           Interpolation Methold (INTERP_MTHD),
//           Interpolation Points (INTERP_PNTS)
//
//   In the expressions above, f are forecast values, o are observed values,
//   and c are climatological values.
//
// Values for these flags are interpreted as follows:
//    (0) Do not generate output of this type
//    (1) Write output to a VSDB file
//    (2) Write output to a VSDB file and a text file
//
output_flag[] = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];

//
// Flag to indicate whether Kendall's Tau and Spearman's Rank Correlation
// Coefficients should be computed.  Computing them over large datasets is
// computationally intensive and slows down the runtime execution significantly.
//    (0) Do not compute these correlation coefficients
//    (1) Compute these correlation coefficients
//
rank_corr_flag = 1;

//
// Specify whether to use the NCEP/NWS Grib conventions.
// If set to zero, the following defaults will be used:
//    (1) Grib codes 128 to 254 will be named "GC_NNN" where NNN is the grib code
//    (2) Grib levels 204 to 254 will be named "GL_NNN" where NNN is the grib level
//
ncep_defaults = 1;

//
// Directory where temp files should be written by the point_stat tool
//
tmp_dir = "/home/arc/VERIFICA/tmp";

//
// Indicate a version number for the contents of this configuration file.
// The value should generally not be modified.
//
version = "V1.1";