Run

Run cindyhg Thu, 04/25/2019 - 10:45

Next, run Wavelet-Stat on the command line using the following command:

wavelet_stat \
$MET_TUTORIAL_DATA/output/pcp_combine/sample_fcst_24L_2005080800V_12A.nc \
$MET_TUTORIAL_DATA/output/pcp_combine/sample_obs_2005080800V_12A.nc \
$MET_TUTORIAL_DATA/config/WaveletStatConfig_tutorial \
-outdir $MET_TUTORIAL_DATA/output/wavelet_stat \
-v 2

Wavelet-Stat is now performing the verification task we requested in the configuration file. It should take several seconds to run. Generally, the Wavelet-Stat tool runs pretty quickly.

When Wavelet-Stat is finished, it will have created four files: two ASCII statistics files, a NetCDF scale decomposition file, and a PostScript summary plot. Open up the PostScript summary plot using the PostScript viewer of your choice, gv, or Ghostview, for example:

gv $MET_TUTORIAL_DATA/output/wavelet_stat/wavelet_stat_240000L_20050808_000000V.ps

This PostScript summary plot contains five pages. The first page summarizes the definition of the tile(s) in the domain. The remaining pages show the difference field (f-o) for each decomposed scale and the statistics for each scale.

Now, let's modify the configuration file and rerun this case. Again, open up the $MET_TUTORIAL_DATA/config/WaveletStatConfig_tutorial file and edit it as follows:

  • Set grid_decomp_flag = AUTO;
    To let the Wavelet-Stat tool automatically define the largest 2^n by 2^n tile that fits in the center of the domain.

Now, rerun the Wavelet-Stat command listed above, and when it is finished, reload the PostScript plot. On the first page of the PostScript plot, note the following:

  • A tile of dimension 128 by 128 was chosen in the center of the domain. Since the dimension increased from 64 (= 2^6) to 128 (= 2^7) the number of scales has increased by 1.
  • The tile chosen includes a large amount of missing data in the observation field. You should try to avoid including missing data when running the Wavelet-Stat tool as it will cause misleading results. Missing data is replaced with a value of 0 for preciptation fields or the mean value of the valid data for non-preciptation fields.

Close that PostScript file.