[Met_help] [rt.rap.ucar.edu #83485] History for point_stat problem of probability forecast

Thu Jan 4 15:35:54 MST 2018

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

Hi All,

     I try to run point_stat mode by using ensemble_stat results ,

the NC output information of ensemble is like:

calculating min and maxes for APCP_A1_ENS_FREQ_gt5...

calculating min and maxes for APCP_A1_ENS_FREQ_ge15.0...

calculating min and maxes for APCP_A1_ENS_FREQ_ge25...

calculating min and maxes for APCP_A1_ENS_FREQ_ge50...

******************************************

But when I run the point_stat, the error likes:

*** Running Point-Stat ***

DEBUG 1: Default Config File: /home-gs/users/SZQXG152/METOOLS/met/share/met/config/PointStatConfig_default

DEBUG 1: User Config File: /home-gs/users/SZQXG152/metenm_test/program/src/PointEnsmConfig_acc01

ERROR  :

ERROR  : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field abbreviation 'APCP_A1_ENS_FREQ_ge15.0' for table version 2

*****************************

Do you know how to fix this problem?

Thank you very much!

Mandy

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

Subject: point_stat problem of probability forecast
From: John Halley Gotway
Time: Thu Jan 04 15:12:42 2018

Hi Mandy,

I would need to see your point_stat config file to be sure, but I'm
guessing you have this set in there:

    obs = fcst;

When processing NetCDF input files (i.e. the NetCDF output of
ensemble-stat), that won't work.

Instead, you'll need to define a full "obs" dictionary to define which
point observations should be used for verification.  By convention,
you
specify the point observations as if they were GRIB1 files.  So for
example...

fcst = {
  field = [ { name="APCP_A1_ENS_FREQ_gt5"; level="(*,*)"; prob = TRUE;
cat_thresh = ==0.1; } ];
}
obs = {
   field = [ { name="APCP"; level="A1"; cat_thresh = >5; } ];
}

Here we're comparing a forecast of the probability of 1-hour precip
greater
than 5 to the observation of 1-hour precip greater than 5.

Hope this helps.

Thanks,
John Halley Gotway

On Thu, Jan 4, 2018 at 8:15 AM, Zhang, Mandy via RT
<met_help at ucar.edu>
wrote:

>
> Thu Jan 04 08:15:02 2018: Request 83485 was acted upon.
> Transaction: Ticket created by Man.Zhang-1 at ou.edu
>        Queue: met_help
>      Subject: point_stat problem of probability forecast
>        Owner: Nobody
>   Requestors: Man.Zhang-1 at ou.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83485 >
>
>
> Hi All,
>
>      I try to run point_stat mode by using ensemble_stat results ,
>
> the NC output information of ensemble is like:
>
> calculating min and maxes for APCP_A1_ENS_FREQ_gt5...
>
> calculating min and maxes for APCP_A1_ENS_FREQ_ge15.0...
>
> calculating min and maxes for APCP_A1_ENS_FREQ_ge25...
>
> calculating min and maxes for APCP_A1_ENS_FREQ_ge50...
>
> ******************************************
>
> But when I run the point_stat, the error likes:
>
>
> *** Running Point-Stat ***
>
> DEBUG 1: Default Config File: /home-gs/users/SZQXG152/
> METOOLS/met/share/met/config/PointStatConfig_default
>
> DEBUG 1: User Config File: /home-gs/users/SZQXG152/
> metenm_test/program/src/PointEnsmConfig_acc01
>
> ERROR  :
>
> ERROR  : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
> abbreviation 'APCP_A1_ENS_FREQ_ge15.0' for table version 2
>
> *****************************
>
> Do you know how to fix this problem?
>
> Thank you very much!
>
> Mandy
>
>

------------------------------------------------
Subject: point_stat problem of probability forecast
From: Zhang, Mandy
Time: Thu Jan 04 15:32:16 2018

Hi John,