[Met_help] [rt.rap.ucar.edu #66204] History for point_stat problem

[John Halley Gotway via RT]

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Greetings.
I'm busy testing MET software using the data from my region(i.e. South Africa). But when I ran the command below I get text files containing only the headers. I have already tested grid_stat and its working well with the same data. I have placed the data in your ftp server under the folder named Phaduli_data. Can you please test it on your side.

/home/verification/METv4.1/bin/point_stat 2014031400_nmm_d01_tmp.grb gdas1.t12z.prepbufr.nc PointStatConfig_default

Thank you in advance......


----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: Re: [rt.rap.ucar.edu #66204] point_stat problem
From: John Halley Gotway
Time: Mon Apr 14 09:57:38 2014

Elelwani,

First, I noticed the following setting in your config file:
    sid  = [ "SA.stations" ];

I see that you've defined a list of stations over which to compute
statistics.  I'd suggest including the full path to that file so that
Point-Stat can find it.

I ran the data you sent and see the several log messages like this
printed to the screen:
    DEBUG 2: Processing TMP/Z2 versus TMP/Z2, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using 0
pairs.

You're getting text files with only headers because Point-Stat is
finding zero matched pairs.  With zero matched pairs, there are no
statistics to compute and thus, no output.  This is a pretty common
question that comes into MET-Help - why am I getting zero matched
pairs from Point-Stat?  There are several possible reasons, but to
better diagnose the issue, try setting the Point-Stat verbosity
level to 4 (-v 4).

At verbosity level 4, Point-Stat dumps out counts of reasons why
observations were not used for each verification task:
    DEBUG 2: Processing TMP/Z2 versus TMP/Z2, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using 0
pairs.
    DEBUG 3: Number of matched pairs  = 0
    DEBUG 3: Observations processed   = 1166
    DEBUG 3: Rejected: GRIB code      = 1166
    DEBUG 3: Rejected: valid time     = 0
    DEBUG 3: Rejected: bad obs value  = 0
    DEBUG 3: Rejected: off the grid   = 0
    DEBUG 3: Rejected: level mismatch = 0
    DEBUG 3: Rejected: quality marker = 0
    DEBUG 3: Rejected: message type   = 0
    DEBUG 3: Rejected: masking region = 0
    DEBUG 3: Rejected: bad fcst value = 0
    DEBUG 3: Rejected: duplicates     = 0

The reason codes are listed in the order that Point-Stat applies the
logic.  In your case, Point-Stat processed 1166 observations from the
gdas1.t12z.prepbufr.nc file and skipped all of them because
of the GRIB code.  So that file doesn't contain any observations of
temperature.

Take a step back and look at how you generated that file by running
PB2NC.  Look in the PB2NC config file - did you include temperature in
the list of observations types to retain?

Thanks,
John


On 04/14/2014 12:37 AM, Elelwani Phaduli via RT wrote:
>
> Mon Apr 14 00:37:57 2014: Request 66204 was acted upon.
> Transaction: Ticket created by Elelwani.Phaduli at weathersa.co.za
>         Queue: met_help
>       Subject: point_stat problem
>         Owner: Nobody
>    Requestors: Elelwani.Phaduli at weathersa.co.za
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=66204 >
>
>
> Greetings.
> I'm busy testing MET software using the data from my region(i.e.
South Africa). But when I ran the command below I get text files
containing only the headers. I have already tested grid_stat and its
working well with the same data. I have placed the data in your ftp
server under the folder named Phaduli_data. Can you please test it on
your side.
>
> /home/verification/METv4.1/bin/point_stat 2014031400_nmm_d01_tmp.grb
gdas1.t12z.prepbufr.nc PointStatConfig_default
>
> Thank you in advance......
>

------------------------------------------------