#!/bin/bash -e

# =============================================================================
# 
# NAME: config.evs.prod.stats.cam.atmos.precip.hireswarw 
# CONTRIBUTOR(S): Marcel Caron, marcel.caron@noaa.gov, NOAA/NWS/NCEP/EMC-VPPPGB
# PURPOSE: Set up configurations to run EVS Hires Window ARW Precipitation - 
#          Statistics
# DEPENDENCIES: $HOMEevs/scripts/stats/cam/exevs_hireswarw_precip_stats.sh
# 
# =============================================================================

echo "BEGIN: $(basename ${BASH_SOURCE[0]})"

# METplus Settings
export MET_PLUS_OUT=${DATA}/${VERIF_CASE}/METplus_output
export MET_CONFIG_OVERRIDES=""

# Restart Settings
export RESTART_DIR="${COMOUTsmall}/restart/c${vhr}"
export COMPLETED_JOBS_FILE="completed_jobs.txt"

# Time Settings
export FHR_GROUP_LIST="FULL"
export FHR_END_FULL="48"
export FHR_INCR_FULL="12"
export FHR_END_SHORT="48"
export FHR_INCR_SHORT="12"

# Model/Obs Settings
export COMINfcst=${COMINhiresw}
export OBS_ACC="01"

# Additional Verification Settings
if [ "$ACC" = "01" ]; then
    export MODEL_ACC="01"
    export BUCKET_INTERVAL=$MODEL_ACC
    export FCST_LEV="A1"
    export FCST_THRESH="ge0.254,ge1.27,ge2.54,ge6.35,ge12.7,ge19.05,ge25.4"
    if [ "$NEST" = "conus" ]; then
        export OBS_LEV="A1"
    elif [ "$NEST" = "ak" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "hi" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "pr" ]; then
        export OBS_LEV="Z0"
    else
        err_exit "The provided NEST, $NEST, is not supported for $MODELNAME"
    fi
    export OBS_THRESH="ge0.254,ge1.27,ge2.54,ge6.35,ge12.7,ge19.05,ge25.4"
    export MODEL_PCP_COMBINE_METHOD="ADD"
    export MODEL_PCP_COMBINE_COMMAND=""
elif [ "$ACC" = "03" ]; then
    export MODEL_ACC="01"
    export BUCKET_INTERVAL=$MODEL_ACC
    export FCST_LEV="A3"
    export FCST_THRESH="ge0.254,ge1.27,ge2.54,ge6.35,ge12.7,ge19.05,ge25.4,ge50.8"
    if [ "$NEST" = "conus" ]; then
        export OBS_LEV="A3"
    elif [ "$NEST" = "ak" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "hi" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "pr" ]; then
        export OBS_LEV="Z0"
    else
        err_exit "The provided NEST, $NEST, is not supported for $MODELNAME"
    fi
    export OBS_THRESH="ge0.254,ge1.27,ge2.54,ge6.35,ge12.7,ge19.05,ge25.4,ge50.8"
    export MODEL_PCP_COMBINE_METHOD="ADD"
    export MODEL_PCP_COMBINE_COMMAND=""
elif [ "$ACC" = "24" ]; then
    export MODEL_ACC="01"
    export BUCKET_INTERVAL=$MODEL_ACC
    export FCST_LEV="A24"
    export FCST_THRESH="ge0.254,ge2.54,ge6.35,ge12.7,ge25.4,ge38.1,ge50.8,ge76.2,ge101.6,ge152.4"
    if [ "$NEST" = "conus" ]; then
        export OBS_LEV="A24"
    elif [ "$NEST" = "ak" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "hi" ]; then
        export OBS_LEV="Z0"
    elif [ "$NEST" = "pr" ]; then
        export OBS_LEV="Z0"
    else
        err_exit "The provided NEST, $NEST, is not supported for $MODELNAME"
    fi
    export OBS_THRESH="ge0.254,ge2.54,ge6.35,ge12.7,ge25.4,ge38.1,ge50.8,ge76.2,ge101.6,ge152.4"
    export MODEL_PCP_COMBINE_METHOD="ADD"
    export MODEL_PCP_COMBINE_COMMAND=""
else
    err_exit "The provided ACC, $ACC, is not supported for $MODELNAME"
fi
if [ "$BOOL_NBRHD" = True ]; then
    export OUTPUT_FLAG_NBRHD="STAT"
    export OUTPUT_FLAG_CATEG="NONE"
    export NBRHD_WIDTHS="1,3,5,7,9,11,13,17,21,25,31,35,39,43,47,53,57,63"
    if [ "$NEST" = "conus" ]; then
        export VERIF_TYPE="ccpa"
        export COMINobs=${EVSINccpa}
        export OBSNAME="ccpa"
        export MIN_IHOUR="00"
        export GRID="G240"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.conus.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Bukovsky_G240_CONUS.nc, ${FIXevs}/masks/Bukovsky_G240_CONUS_East.nc, ${FIXevs}/masks/Bukovsky_G240_CONUS_West.nc, ${FIXevs}/masks/Bukovsky_G240_CONUS_Central.nc, ${FIXevs}/masks/Bukovsky_G240_CONUS_South.nc"
    elif [ "$NEST" = "ak" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="06"
        export GRID="G091"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.ak.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Alaska_G091.nc"
    elif [ "$NEST" = "hi" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="00"
        export GRID="G196"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.hi.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Hawaii_G196.nc"
    elif [ "$NEST" = "pr" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="06"
        export GRID="G194"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.pr.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Puerto_Rico_G194.nc"
    else
        err_exit "The provided NEST, $NEST, is not supported for $MODELNAME"
    fi
else
    export OUTPUT_FLAG_NBRHD="NONE"
    export OUTPUT_FLAG_CATEG="STAT"
    export NBRHD_WIDTHS="1"
    if [ "$NEST" = "conus" ]; then
        export VERIF_TYPE="ccpa"
        export COMINobs=${EVSINccpa}
        export OBSNAME="ccpa"
        export MIN_IHOUR="00"
        export GRID="G212"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.conus.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Bukovsky_G212_CONUS.nc, ${FIXevs}/masks/Bukovsky_G212_CONUS_East.nc, ${FIXevs}/masks/Bukovsky_G212_CONUS_West.nc, ${FIXevs}/masks/Bukovsky_G212_CONUS_Central.nc, ${FIXevs}/masks/Bukovsky_G212_CONUS_South.nc"
    elif [ "$NEST" = "ak" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="06"
        export GRID="G216"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.ak.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Alaska_G216.nc"
    elif [ "$NEST" = "hi" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="00"
        export GRID="G196"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.hi.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Hawaii_G196.nc"
    elif [ "$NEST" = "pr" ]; then
        export VERIF_TYPE="mrms"
        export COMINobs=${EVSINmrms}
        export OBSNAME="mrms"
        export MIN_IHOUR="06"
        export GRID="G194"
        export MODEL_INPUT_TEMPLATE="hiresw.{init?fmt=%Y%m%d}/hiresw.t{init?fmt=%2H}z.arw_5km.f{lead?fmt=%2H}.pr.grib2"
        export MASK_POLY_LIST="${FIXevs}/masks/Puerto_Rico_G194.nc"
    else
        err_exit "The provided NEST, $NEST, is not supported for $MODELNAME"
    fi
fi

echo "END: $(basename ${BASH_SOURCE[0]})"