As mentioned on the previous page, the output of Wavelet-Stat typically consists of four files: two ASCII statistics files, one NetCDF scale decomposition file, and one PostScript summary plot. The output of these files may be disabled in the configuration file or using the appropriate command line argument. In this example, the output is written to the $MET_TUTORIAL_DATA/output/wavelet_stat directory as we requested on the command line.
The Wavelet-Stat output file naming convention is similar to that of the Point-Stat and Grid-Stat tools. The four Wavelet-Stat output files are described briefly below:
- The PostScript file was described on the previous page.
- The NetCDF scale decomposition file contains the raw, thresholded, and decomposed fields for each variable, tile, and threshold used. Note that while the PostScript plot only shows the difference (f-o) fields, the NetCDF file contains the actual forecast, observation, and difference fields decomposed for each scale.
- The ASCII ISC file contains the ISC line type with a header row for the column names.
- The ASCII STAT file contains only the ISC line type. Currently, the Wavelet-Stat tool only creates one output line type. So the ISC file and the STAT file are almost identical. In future versions of MET, the Wavelet-Stat tool may be enhanced to produce additional line types.
Open up the $MET_TUTORIAL_DATA/output/wavelet_stat/wavelet_stat_240000L_20050808_000000V_isc.txt ISC file using the text editor of your choice, and note the following:
- The header columns are identical to the other ASCII output files from Point-Stat and Grid-Stat.
- The LINE_TYPE column is set to ISC, indicating that the columns to follow contain information about the intensity-scale method.
- This file contains eights rows of data. The ISCALE column indicates the scale for that row. The row with ISCALE equal to 0 contains scores for the thresholded binary fields. The rows with ISCALE greater than 0 contain scores for the thresholded binary fields decomposed into separate scales.
- Looking carefully you'll see that the columns for MSE, FENERGY2, and OENERGY2 are additive across the scales. The sum of these values in the lines where ISCALE is greater than 0 equals the values in the line where ISCALE equals 0.
Close this file and use the ncdump and ncview utilities (if available on your machine) to view the NetCDF output of Wavelet-Stat:
When clicking through and displaying each variable, note that some have a dimension for scale. Click through the different scales to see the decompositions of those fields.
