By default, TC-Stat writes its job output to the screen but can easily be redirected to a file using the -out option. The -out option is specific to each job. When run on the command line, TC-Stat executes a single job. Therefore using the -out option once is sufficient. When run with a configuration file, the -out option should be specified once for each job defined in the jobs array.
The MET tarball also includes an Rscript which automates calls to TC-Stat and the generation of graphics. The scripts/Rscripts/plot_tcmpr.R Rscript performs two main steps: call TC-Stat to filter track data and create plots of the filtered track data. If Rscript or the R boot package are not available on your system, you may skip this step.
Plotting TCMPR lines
The plot_tcmpr.R Rscript requires that the MET_BUILD_BASE environment variable be set to the top-level MET source code directory and that TC-Stat be found in your path. Run the following commands based on whether you are using C-Shell or the Bourne Shell:
- Determine your shell:
- For C-Shell:
setenv PATH `pwd`/bin:$PATH
- For Bourne Shell:
export PATH=`pwd`/bin:$PATH
Listed below are several examples of running this Rscript. The output printed to the screen lists the name of the output files. You may display them using a web browser, the display command, or any other graphics program.
- Run the Rscript with no arguments to see the usage:
- Specify an input file and output directory:
-lookin $MET_TUTORIAL_DATA/output/tc_pairs/tc_pairs_sandy.tcst \
-outdir $MET_TUTORIAL_DATA/output/tc_pairs
- By default, the script creates boxplots of event equalized track errors (TK_ERR). Display the output image:
- Turn off event equalization, plot only 3 models (OFCL, OCD5, and HWRF), and create several plot types for track and intensity errors:
-lookin $MET_TUTORIAL_DATA/output/tc_pairs/tc_pairs_sandy.tcst \
-outdir $MET_TUTORIAL_DATA/output/tc_pairs \
-no_ee \
-filter '-amodel OFCL,OCD5,HWRF' \
-dep TK_ERR,ABS\(AMAX_WIND-BMAX_WIND\) \
-plot BOXPLOT,MEAN,MEDIAN
The -no_ee option disables event equalization. The -filter options are passed directly to the TC-Stat filtering job. The -dep options specifies the column(s) of data to be plotted. The -plot specifies the desired plot types. Display the output images.
- Plot HWRF track errors broken down by storm intensity level:
-lookin $MET_TUTORIAL_DATA/output/tc_pairs/tc_pairs_sandy.tcst \
-outdir $MET_TUTORIAL_DATA/output/tc_pairs \
-no_ee \
-filter '-amodel HWRF' \
-dep TK_ERR \
-plot MEAN \
-series LEVEL \
-title "HWRF Mean Track Error for Sandy by Intensity Level"
The -series option overrides the default plotting by the AMODEL column. Use -title to override the default plot title. Display the output image:
These examples are meant as a brief introduction to this plotting script. We encourage you to run it many more times and try plotting different columns of data or pass different filtering options.
