# HRRR-E PYTHON GRIB GENERATION
# CONFIGURATION FILE

## CREATE DICTIONARIES ##
grib_discipline = {}
grib_category = {}
grib_number = {}
grib_surface = {}
grib_level = {}

##################################
#  GENERAL ENSEMBLE SETTINGS     #
##################################

# total number of members in HRRR-E
mnum = 9
# minimum number of ens members to process probabilities
mmin = 9
# model resolution (km)
dx = 3

# wgrib2 executable
wgrib2 = '/home/rtrr/RAP/exec/UPP/wgrib2'


#################################
#  ENS MEAN/PMM/SPREAD SETTINGS #
#################################

# parms to process
mean_parms = []
std_parms = []
max_parms = ['1hqpf']
min_parms = []
pmm_parms = []

# neighborhood (points) for probability-matched mean
pmm_neighborhood = 15

##############################
#  MULTI-MODEL POST SETTINGS #
##############################

# number of members to use from each system
hrrre_mnum = 5
ncar_mnum = 5

# parms to process
multi_mean_parms = []
multi_std_parms = []
multi_max_parms = []
multi_min_parms = []
multi_pmm_parms = []

###################
#  GRIB TABLES    #
###################

grib_discipline['t2m'] = 0
grib_category['t2m'] = 0
grib_number['t2m'] = 0
grib_surface['t2m'] = 'sfc'
grib_level['t2m'] = 2

grib_discipline['td2m'] = 0
grib_category['td2m'] = 0
grib_number['td2m'] = 6
grib_surface['td2m'] = 'sfc'
grib_level['td2m'] = 2

grib_discipline['spd10'] = 0
grib_category['spd10'] = 2
grib_number['spd10'] = [2,3,1] # u, v, speed
grib_surface['spd10'] = 'sfc'
grib_level['spd10'] = 10

grib_discipline['capesfc'] = 0
grib_category['capesfc'] = 7
grib_number['capesfc'] = 6
grib_surface['capesfc'] = 'sfc'
grib_level['capesfc'] = 0

grib_discipline['qpf'] = 0
grib_category['qpf'] = 1
grib_number['qpf'] = 8
grib_surface['qpf'] = 'sfc'
grib_level['qpf'] = 0

grib_discipline['cref'] = 0
grib_category['cref'] = 16
grib_number['cref'] = 196
grib_surface['cref'] = 10
grib_level['cref'] = 0

grib_discipline['ref1'] = 0
grib_category['ref1'] = 16
grib_number['ref1'] = 195
grib_surface['ref1'] = 'sfc'
grib_level['ref1'] = 1000

grib_discipline['snow'] = 0
grib_category['snow'] = 1
grib_number['snow'] = 29
grib_surface['snow'] = 'sfc'
grib_level['snow'] = 0


##################################
# SEVERE WX PROBABILITY SETTINGS #
##################################

# create an output grid every X hours
severe_proc_interval = 1
# forecast interval length (hours)
severe_interval = 4
# list of severe products
severe_types = ['wind','sigwind','hail','sighail','tornado']

# Gaussian filter width (km) for all severe probabilities
severe_filter = 120
# neighborhood size (km)
severe_neighborhood = 40

# severe wind thresholds (m/s)
wind_thresh = 25.7
sigwind_thresh = 33.4
wind_muli_thresh = 0.0 # maximum most unstable lifted index

# size thresholds for severe hail (meters)
hail_thresh = 0.0254 # 1 inch
sighail_thresh = 0.0508 # 2 inches

# tornado thresholds
tor_uh_thresh = 75
tor_lcl_thresh = 1500 # m
tor_shear01_thresh = 7.5 # m/s
tor_caperatio_thresh = 0.75

# smoothing of environmental parameters (Gaussian filter width in km)
muli_smooth = 120 
lcl_smooth = 120
shear01_smooth = 120
caperatio_smooth = 120


#################
# PQPF SETTINGS #
#################

# QPF accumulation intervals (hours)
pqpf_acc_intervals = [1,3,6,24]
pqpf_interval = 6
# create an output grid every X hours
pqpf_proc_intervals = {}
pqpf_proc_intervals[1] = 1
pqpf_proc_intervals[3] = 3
pqpf_proc_intervals[6] = 3
pqpf_proc_intervals[24] = 3
# PQPF thresholds (inches) for each accumulation interval
pqpf_thresh = {}
pqpf_thresh[1] = [1.0,2.0]
pqpf_thresh[3] = [1.0,2.0,3.0]
pqpf_thresh[6] = [1.0,2.0,3.0]
pqpf_thresh[24] = [1.0,3.0,6.0]
# product types
pqpf_types = ['neighbor']
#pqpf_types = ['point','neighbor']
# neighborhood size (km)
pqpf_neighborhood = 40
# gaussian smoother widths (km)
pqpf_smooth = 25
pqpf_neighbor_smooth = 25

# Flash flood guidance accumulation intervals
ffg_acc_intervals = [3,6,24]
# Probability of exceeding flash flood guidance by X inches
ffg_thresh = [0]
# product types
ffg_types = ['neighbor']
#ffg_types = ['point','neighbor']
# gaussian smoother widths (km)
ffg_smooth = 25
ffg_neighbor_smooth = 25
# neighborhood size (km)
ffg_neighborhood = 40

# Atlas-14 accumulation intervals
ari_acc_intervals = [6,24]
# Probability of exceeding Atlas-14 recurrence intervals
ari_years = {}
ari_years[6] = [100]
ari_years[24] = [100]
# product types
ari_types = ['neighbor']
#ari_types = ['point','neighbor']
# gaussian smoother widths (km)
ari_smooth = 25
ari_neighbor_smooth = 25
# neighborhood size (km)
ari_neighborhood = 40


######################
# WINTER WX SETTINGS #
######################

# PROBABILISTIC SNOW
# snow accumulation intervals
psnow_acc_intervals = [1,6]
# create an output grid every X hours
psnow_proc_intervals = {}
psnow_proc_intervals[1] = 1
psnow_proc_intervals[6] = 6
# snow thresholds (inches)
psnow_thresh = {}
psnow_thresh[1] = [0.5,1.0]
psnow_thresh[6] = [1.0,3.0,6.0]
# probability types
#psnow_types = ['neighbor']
psnow_types = ['point','neighbor']
# neighborhood size (km)
psnow_neighborhood = 40
# gaussian smoother width (km)
psnow_neighbor_smooth = 25
psnow_smooth = 0

# PROBABILISTIC PTYPE
# create an output grid every X hours
ptype_proc_interval = 1
# gaussian smoother width (km)
ptype_smooth = 25
# precipitation rate threshold (in/h)
#prate_thresh = 0.01
prate_thresh = -0.0001

# PROBABILISTIC REFLECTIVITY
# create an output grid every X hours
ref_proc_interval = 1
# reflectivity probabilities
ref_thresh = [20]
# gaussian smoother width (km)
ref_smooth = 25

#####################
# AVIATION SETTINGS #
#####################

# create an output grid every X hours
aviation_proc_interval = 1
# aviation thresholds
echotop_thresh = [25,30,35,40]  #kft
vis_thresh = [1,3,5]  # mi
ceil_thresh = [500,1000,3000]  # ft
flr_cat = ['lifr','ifr','mvfr','vfr']
flr_num = {}
flr_num['lifr'] = 1
flr_num['ifr'] = 2
flr_num['mvfr'] = 3
flr_num['vfr'] = 4
flr_vis_thresh = {}
flr_vis_thresh['lifr'] = 1.0
flr_vis_thresh['ifr'] = 3.0
flr_vis_thresh['mvfr'] = 5.0
flr_ceil_thresh = {}
flr_ceil_thresh['lifr'] = 500
flr_ceil_thresh['ifr'] = 1000
flr_ceil_thresh['mvfr'] = 3000
icing_vlev_min = 10
icing_vlev_max = 25
icing_clw_thresh = 0.1 # g/kg
avn_eas_radii = [6,12,20,40,80,120]
avn_alpha = 0.1
avn_dcrit_exp = 3
avn_rad_smooth = 12
avn_p_smooth = 25
# gaussian smoother width (km)
aviation_smooth = 25