Skip to content

Commit

Permalink
Mods to METplus vx to address multiple issues (ufs-community#826)
Browse files Browse the repository at this point in the history
* Mods to METplus conf files: TCDC specifications, correction to level specificatons in point-stat mean and prob files, and added functionality to make METplus output dirs in ex scripts.

* Updated comments in MET ex-scripts for creating output directories.

* Fixed minor formatting issue in exregional_run_gridstatvx.sh
  • Loading branch information
michelleharrold authored Aug 25, 2022
1 parent 83472a2 commit f870a53
Show file tree
Hide file tree
Showing 13 changed files with 152 additions and 57 deletions.
20 changes: 20 additions & 0 deletions scripts/exregional_run_ensgridvx.sh
Original file line number Diff line number Diff line change
Expand Up @@ -117,6 +117,26 @@ else
LOG_SUFFIX=ensgrid_${CDATE}_${VAR}
fi

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat" # Output directory for ensemble_stat tool.

#
# If the variable is accumulated precipitation for a time interval
# (bucket) other than 1 hour, the MET/METplus tools called below will
# include pcp_combine. In that case, create (if necessary) directories
# needed by pcp_combine.
#
if [ "${VAR}" = "APCP" ] && [ "${ACCUM: -1}" != "1" ]; then
mkdir_vrfy -p "${EXPTDIR}/metprd/pcp_combine" # For observations
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/pcp_combine" # For forecast
fi

#
#-----------------------------------------------------------------------
#
Expand Down
9 changes: 9 additions & 0 deletions scripts/exregional_run_ensgridvx_mean.sh
Original file line number Diff line number Diff line change
Expand Up @@ -115,6 +115,15 @@ else
LOG_SUFFIX=ensgrid_mean_${CDATE}_${VAR}
fi

#
#-----------------------------------------------------------------------
#
# Make sure directory in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat_mean" # Output directory for grid_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
9 changes: 9 additions & 0 deletions scripts/exregional_run_ensgridvx_prob.sh
Original file line number Diff line number Diff line change
Expand Up @@ -115,6 +115,15 @@ else
LOG_SUFFIX=ensgrid_prob_${CDATE}_${VAR}
fi

#
#-----------------------------------------------------------------------
#
# Make sure directory in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat_prob" # Output directory for grid_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
10 changes: 10 additions & 0 deletions scripts/exregional_run_enspointvx.sh
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,16 @@ export fhr_list
#
LOG_SUFFIX=enspoint_${CDATE}

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/metprd/pb2nc" # Output directory for pb2nc tool.
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat" # Output directory for ensemble_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
10 changes: 10 additions & 0 deletions scripts/exregional_run_enspointvx_mean.sh
Original file line number Diff line number Diff line change
Expand Up @@ -108,6 +108,16 @@ export fhr_list
INPUT_BASE=${EXPTDIR}/${CDATE}/metprd/ensemble_stat
LOG_SUFFIX=enspoint_mean_${CDATE}

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/metprd/pb2nc" # Output directory for pb2nc tool.
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat_mean" # Output directory for point_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
10 changes: 10 additions & 0 deletions scripts/exregional_run_enspointvx_prob.sh
Original file line number Diff line number Diff line change
Expand Up @@ -108,6 +108,16 @@ export fhr_list
INPUT_BASE=${EXPTDIR}/${CDATE}/metprd/ensemble_stat
LOG_SUFFIX=enspoint_prob_${CDATE}

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/metprd/pb2nc" # Output directory for pb2nc tool.
mkdir_vrfy -p "${EXPTDIR}/${CDATE}/metprd/ensemble_stat_prob" # Output directory for point_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
19 changes: 19 additions & 0 deletions scripts/exregional_run_gridstatvx.sh
Original file line number Diff line number Diff line change
Expand Up @@ -118,6 +118,25 @@ elif [[ ${DO_ENSEMBLE} == "TRUE" ]]; then
fi
fi

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${OUTPUT_BASE}/metprd/grid_stat"
#
# If the variable is accumulated precipitation for a time interval
# (bucket) other than 1 hour, the MET/METplus tools called below will
# include pcp_combine. In that case, create (if necessary) directories
# needed by pcp_combine.
#
if [ "${VAR}" = "APCP" ] && [ "${ACCUM: -1}" != "1" ]; then
mkdir_vrfy -p "${EXPTDIR}/metprd/pcp_combine" # For observations
mkdir_vrfy -p "${OUTPUT_BASE}/metprd/pcp_combine" # For forecast
fi

#
#-----------------------------------------------------------------------
#
Expand Down
10 changes: 10 additions & 0 deletions scripts/exregional_run_pointstatvx.sh
Original file line number Diff line number Diff line change
Expand Up @@ -110,6 +110,16 @@ elif [[ ${DO_ENSEMBLE} == "TRUE" ]]; then
LOG_SUFFIX=pointstat_${CDATE}_${ENSMEM}
fi

#
#-----------------------------------------------------------------------
#
# Make sure directories in which output files will be placed exist.
#
#-----------------------------------------------------------------------
#
mkdir_vrfy -p "${EXPTDIR}/metprd/pb2nc" # Output directory for pb2nc tool.
mkdir_vrfy -p "${OUTPUT_BASE}/metprd/point_stat" # Output directory for point_stat tool.

#
#-----------------------------------------------------------------------
#
Expand Down
5 changes: 2 additions & 3 deletions ush/templates/parm/metplus/EnsembleStat_conus_sfc.conf
Original file line number Diff line number Diff line change
Expand Up @@ -215,7 +215,7 @@ ENS_VAR3_OPTIONS = GRIB2_pdt = 0; ;; derive instantaneous 10-m wind from U/V com
ENS_VAR4_NAME = TCDC
ENS_VAR4_LEVELS = L0
ENS_VAR4_THRESH = <25, >75
ENS_VAR4_OPTIONS = GRIB_lvl_typ = 200; GRIB2_ipdtmpl_index = [ 14 ]; GRIB2_ipdtmpl_val = [ 0 ]; interp = { type = [ { method = NEAREST; width = 1; } ]; }
ENS_VAR4_OPTIONS = GRIB_lvl_typ = 200; GRIB2_ipdtmpl_index=[27]; GRIB2_ipdtmpl_val=[255]; interp = { type = [ { method = NEAREST; width = 1; } ]; }

ENS_VAR5_NAME = VIS
ENS_VAR5_LEVELS = L0
Expand Down Expand Up @@ -261,11 +261,10 @@ OBS_VAR3_OPTIONS = obs_error = { flag = TRUE; dist_type = NONE; dist_parm = [];
FCST_VAR4_NAME = TCDC
FCST_VAR4_LEVELS = L0
FCST_VAR4_THRESH = <25, >75
FCST_VAR4_OPTIONS = GRIB_lvl_typ = 200; GRIB2_ipdtmpl_index = 8; GRIB2_ipdtmpl_val = {lead?fmt=%H};
FCST_VAR4_OPTIONS = GRIB_lvl_typ = 200; GRIB2_ipdtmpl_index=[27]; GRIB2_ipdtmpl_val=[255]; interp = { type = [ { method = NEAREST; width = 1; } ]; }
OBS_VAR4_NAME = TCDC
OBS_VAR4_LEVELS = L0
OBS_VAR4_THRESH = <25, >75
OBS_VAR4_OPTIONS = interp = { type = [ { method = NEAREST; width = 1; } ]; }

BOTH_VAR5_NAME = VIS
BOTH_VAR5_LEVELS = L0
Expand Down
3 changes: 1 addition & 2 deletions ush/templates/parm/metplus/PointStat_conus_sfc.conf
Original file line number Diff line number Diff line change
Expand Up @@ -212,10 +212,9 @@ BOTH_VAR7_LEVELS = Z0

FCST_VAR8_NAME = TCDC
FCST_VAR8_LEVELS = L0
FCST_VAR8_OPTIONS = GRIB_lvl_typ = 200; interp = { type = [ { method = NEAREST; width = 1; } ]; }
FCST_VAR8_OPTIONS = GRIB_lvl_typ = 200; GRIB2_ipdtmpl_index=[27]; GRIB2_ipdtmpl_val=[255]; interp = { type = [ { method = NEAREST; width = 1; } ]; }
OBS_VAR8_NAME = TCDC
OBS_VAR8_LEVELS = L0
OBS_VAR8_OPTIONS = interp = { type = [ { method = NEAREST; width = 1; } ]; }

BOTH_VAR9_NAME = VIS
BOTH_VAR9_LEVELS = L0
Expand Down
24 changes: 12 additions & 12 deletions ush/templates/parm/metplus/PointStat_conus_sfc_prob.conf
Original file line number Diff line number Diff line change
Expand Up @@ -185,95 +185,95 @@ POINT_STAT_ONCE_PER_FIELD = False
# Note: If FCST_VAR<n>_* is set, then a corresponding OBS_VAR<n>_* variable must be set
# To use one variables for both forecast and observation data, set BOTH_VAR<n>_* instead
FCST_VAR1_NAME = TMP_Z2_ENS_FREQ_ge268
FCST_VAR1_LEVELS = (*,*)
FCST_VAR1_LEVELS = Z2
FCST_VAR1_THRESH = ==0.1

OBS_VAR1_NAME = TMP
OBS_VAR1_LEVELS = Z2
OBS_VAR1_THRESH = >=268

FCST_VAR2_NAME = TMP_Z2_ENS_FREQ_ge273
FCST_VAR2_LEVELS = (*,*)
FCST_VAR2_LEVELS = Z2
FCST_VAR2_THRESH = ==0.1

OBS_VAR2_NAME = TMP
OBS_VAR2_LEVELS = Z2
OBS_VAR2_THRESH = >=273

FCST_VAR3_NAME = TMP_Z2_ENS_FREQ_ge278
FCST_VAR3_LEVELS = (*,*)
FCST_VAR3_LEVELS = Z2
FCST_VAR3_THRESH = ==0.1

OBS_VAR3_NAME = TMP
OBS_VAR3_LEVELS = Z2
OBS_VAR3_THRESH = >=278

FCST_VAR4_NAME = TMP_Z2_ENS_FREQ_ge293
FCST_VAR4_LEVELS = (*,*)
FCST_VAR4_LEVELS = Z2
FCST_VAR4_THRESH = ==0.1

OBS_VAR4_NAME = TMP
OBS_VAR4_LEVELS = Z2
OBS_VAR4_THRESH = >=293

FCST_VAR5_NAME = TMP_Z2_ENS_FREQ_ge298
FCST_VAR5_LEVELS = (*,*)
FCST_VAR5_LEVELS = Z2
FCST_VAR5_THRESH = ==0.1

OBS_VAR5_NAME = TMP
OBS_VAR5_LEVELS = Z2
OBS_VAR5_THRESH = >=298

FCST_VAR6_NAME = TMP_Z2_ENS_FREQ_ge303
FCST_VAR6_LEVELS = (*,*)
FCST_VAR6_LEVELS = Z2
FCST_VAR6_THRESH = ==0.1

OBS_VAR6_NAME = TMP
OBS_VAR6_LEVELS = Z2
OBS_VAR6_THRESH = >=303

FCST_VAR7_NAME = DPT_Z2_ENS_FREQ_ge263
FCST_VAR7_LEVELS = (*,*)
FCST_VAR7_LEVELS = Z2
FCST_VAR7_THRESH = ==0.1

OBS_VAR7_NAME = DPT
OBS_VAR7_LEVELS = Z2
OBS_VAR7_THRESH = >=263

FCST_VAR8_NAME = DPT_Z2_ENS_FREQ_ge268
FCST_VAR8_LEVELS = (*,*)
FCST_VAR8_LEVELS = Z2
FCST_VAR8_THRESH = ==0.1

OBS_VAR8_NAME = DPT
OBS_VAR8_LEVELS = Z2
OBS_VAR8_THRESH = >=268

FCST_VAR9_NAME = DPT_Z2_ENS_FREQ_ge273
FCST_VAR9_LEVELS = (*,*)
FCST_VAR9_LEVELS = Z2
FCST_VAR9_THRESH = ==0.1

OBS_VAR9_NAME = DPT
OBS_VAR9_LEVELS = Z2
OBS_VAR9_THRESH = >=273

FCST_VAR10_NAME = DPT_Z2_ENS_FREQ_ge288
FCST_VAR10_LEVELS = (*,*)
FCST_VAR10_LEVELS = Z2
FCST_VAR10_THRESH = ==0.1

OBS_VAR10_NAME = DPT
OBS_VAR10_LEVELS = Z2
OBS_VAR10_THRESH = >=288

FCST_VAR11_NAME = DPT_Z2_ENS_FREQ_ge293
FCST_VAR11_LEVELS = (*,*)
FCST_VAR11_LEVELS = Z2
FCST_VAR11_THRESH = ==0.1

OBS_VAR11_NAME = DPT
OBS_VAR11_LEVELS = Z2
OBS_VAR11_THRESH = >=293

FCST_VAR12_NAME = DPT_Z2_ENS_FREQ_ge298
FCST_VAR12_LEVELS = (*,*)
FCST_VAR12_LEVELS = Z2
FCST_VAR12_THRESH = ==0.1

OBS_VAR12_NAME = DPT
Expand Down
2 changes: 1 addition & 1 deletion ush/templates/parm/metplus/PointStat_upper_air_mean.conf
Original file line number Diff line number Diff line change
Expand Up @@ -216,7 +216,7 @@ OBS_VAR4_NAME = DPT
OBS_VAR4_LEVELS = P850
OBS_VAR4_THRESH = >=273, >=278, >=283

FCST_VAR5_NAME = DPT_P850_ENS_MEAN
FCST_VAR5_NAME = DPT_P700_ENS_MEAN
FCST_VAR5_LEVELS = P700
FCST_VAR5_THRESH = >=263, >=286, >=273

Expand Down
Loading

0 comments on commit f870a53

Please sign in to comment.