Run

Series-Analysis Tool: Run

First, we need to prepare our observations by putting 1-hourly StageII precipitation forecasts into 3-hourly buckets. Create an output directory:

mkdir -p sample_obs/ST2ml_3h

Run the following PCP-Combine commands to prepare the observations:

pcp_combine -sum 00000000_000000 01 20050807_030000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080703V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_060000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080706V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_090000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080709V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_120000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080712V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_150000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080715V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_180000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080718V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050807_210000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080721V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

pcp_combine -sum 00000000_000000 01 20050808_000000 03 \

sample_obs/ST2ml_3h/sample_obs_2005080800V_03A.nc \

-pcpdir ${METPLUS_DATA}/met_test/data/sample_obs/ST2ml

Note that the previous set of PCP-Combine commands could easily be run by looping through times in METplus Wrappers! The MET tools are often run using METplus Wrappers rather than typing individual commands by hand. You'll learn more about automation using the METplus Wrappers throughout the tutorial.

Next, we'll run Series-Analysis using the following command:

series_analysis \

-fcst ${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_03.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_06.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_09.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_15.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_18.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_21.tm00_G212 \

${METPLUS_DATA}/met_test/data/sample_fcst/2005080700/wrfprs_ruc13_24.tm00_G212 \

-obs sample_obs/ST2ml_3h/sample_obs_2005080703V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080706V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080709V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080712V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080715V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080718V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080721V_03A.nc \

sample_obs/ST2ml_3h/sample_obs_2005080800V_03A.nc \

-out series_analysis_2005080700_2005080800_3A.nc \

-config SeriesAnalysisConfig_tutorial \

-v 2

The statistics we requested in the configuration file will be computed separately for each grid location and accumulated over a time series of eight three-hour accumulations over a 24-hour period. Each grid point will have up to 8 matched pair values.

Note how long this command line is. Imagine how long it would be for a series of 100 files! Instead of listing all of the input files on the command line, you can list them in an ASCII file and pass that to Series-Analysis using the -fcst and -obs options.

johnhg Fri, 07/26/2019 - 11:44