[Met_help] [rt.rap.ucar.edu #62959] History for POINT_STAT Error APCP

John Halley Gotway via RT met_help at ucar.edu
Wed Sep 11 09:52:08 MDT 2013


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

Hello MET Helpdesk

I am trying to use point_stat to verify WRF against rain gauge data.
I am using METv4.0.
When I run 
point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc ./PointStatConfig_default
I get an error message:

DEBUG 1: Default Config File: /media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_default
ERROR  : 
ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field abbreviation 'APCP_surface' for table version 2
ERROR  :

I have post processed WRF with UPPv1.1
Which created a grib file called WRFPRS_d01.023.
The total precipitation variable in the grib file is called "APCP_surface".

If I try with fcst variable APCP, I get an error message:
ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24 found in file: ./WRFPRS_d01.023
If I try A1 instead of A24 I get.
ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1 found in file: ./WRFPRS_d01.023

I would lke to use MET to verify rainfall but it seams impossible.
Have I done something wrong with running UPPv1.1?
>From my wrf_cntrl.parm
" (ACM TOTAL PRECIP    ) SCAL=( 5.0)
 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000)"

Any help would be appreciated.

Thanks
Michael Weston

My config files looks like this:
// Forecast and observation fields to be verified
//
fcst = {
   wind_thresh  = [ NA ];
   message_type = [ "ADPSFC" ];

   field = [
      {
        name       = "APCP_surface";
        level      = "A24";
        cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
      }

   ];

};
obs = {
   wind_thresh  = [ NA ];
   message_type = [ "ADPSFC" ];

   field = [
      {
        name       = "APCP";
        level      = [ "A24" ];
        cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
      }

   ];

};



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

Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
From: John Halley Gotway
Time: Mon Sep 09 10:25:30 2013

Michael,

Can you please try the following first... just edit your config file
like this:

//
// Forecast and observation fields to be verified
//
fcst = {
    wind_thresh  = [ NA ];
    message_type = [ "ADPSFC" ];

    field = [
       {
         name       = "APCP";
         level      = "A24";
         cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
       }

    ];

};
obs = fcst;

Then try running Point-Stat.  If that doesn't fix the problem, can you
please send us the data you're passing to Point-Stat?  You can follow
the instructions listed here to post data to our anonymous
ftp site:
    http://www.dtcenter.org/met/users/support/met_help.php#ftp

Please send us:
(1) GRIB forecast file: WRFPRS_d01.023
(2) NetCDF point observation file: ascii2nc_lieb_2000.nc
(3) Point-Stat configuration file: PointStatConfig_default

I'll try running your data here and figure out what's going on.

Thanks,
John

On 09/09/2013 09:50 AM, Michael Weston via RT wrote:
>
> Mon Sep 09 09:50:52 2013: Request 62959 was acted upon.
> Transaction: Ticket created by michael at escience.co.za
>         Queue: met_help
>       Subject: POINT_STAT Error APCP
>         Owner: Nobody
>    Requestors: michael at escience.co.za
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>
>
> Hello MET Helpdesk
>
> I am trying to use point_stat to verify WRF against rain gauge data.
> I am using METv4.0.
> When I run
> point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc
./PointStatConfig_default
> I get an error message:
>
> DEBUG 1: Default Config File:
/media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_default
> ERROR  :
> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'APCP_surface' for table version 2
> ERROR  :
>
> I have post processed WRF with UPPv1.1
> Which created a grib file called WRFPRS_d01.023.
> The total precipitation variable in the grib file is called
"APCP_surface".
>
> If I try with fcst variable APCP, I get an error message:
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: ./WRFPRS_d01.023
> If I try A1 instead of A24 I get.
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1
found in file: ./WRFPRS_d01.023
>
> I would lke to use MET to verify rainfall but it seams impossible.
> Have I done something wrong with running UPPv1.1?
>>From my wrf_cntrl.parm
> " (ACM TOTAL PRECIP    ) SCAL=( 5.0)
>   L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000 00000 00000 00000)"
>
> Any help would be appreciated.
>
> Thanks
> Michael Weston
>
> My config files looks like this:
> // Forecast and observation fields to be verified
> //
> fcst = {
>     wind_thresh  = [ NA ];
>     message_type = [ "ADPSFC" ];
>
>     field = [
>        {
>          name       = "APCP_surface";
>          level      = "A24";
>          cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>        }
>
>     ];
>
> };
> obs = {
>     wind_thresh  = [ NA ];
>     message_type = [ "ADPSFC" ];
>
>     field = [
>        {
>          name       = "APCP";
>          level      = [ "A24" ];
>          cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>        }
>
>     ];
>
> };
>

------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
From: Michael Weston
Time: Tue Sep 10 02:48:55 2013

Hi John,

Thanks for your reply.
After editing the config file so that obs=fcst I receive the following
error message:

ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: WRFPRS_d01.023

I have uploaded my data to /incoming/irap/met_help/weston_data/
I think the units in my ascii file are incorrect for now. But I am
just trying to get the point_stat to run at this stage.

Regards
Michael

EScience Associates (Pty) Ltd
E-mail: michael at escience.co.za
Web:   www.escience.co.za
PO Box 2950, Saxonwold, 2132
9 Victoria Street, Oaklands, Johannesburg, 2192
Tel:  +27 (0)11 718 6380
Fax: 086 527 5241
VAT No: 473 025 4416
Reg No: 2009/014472/07
________________________________________
From: John Halley Gotway via RT [met_help at ucar.edu]
Sent: Monday, September 09, 2013 6:25 PM
To: Michael Weston
Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP

Michael,

Can you please try the following first... just edit your config file
like this:

//
// Forecast and observation fields to be verified
//
fcst = {
    wind_thresh  = [ NA ];
    message_type = [ "ADPSFC" ];

    field = [
       {
         name       = "APCP";
         level      = "A24";
         cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
       }

    ];

};
obs = fcst;

Then try running Point-Stat.  If that doesn't fix the problem, can you
please send us the data you're passing to Point-Stat?  You can follow
the instructions listed here to post data to our anonymous
ftp site:
    http://www.dtcenter.org/met/users/support/met_help.php#ftp

Please send us:
(1) GRIB forecast file: WRFPRS_d01.023
(2) NetCDF point observation file: ascii2nc_lieb_2000.nc
(3) Point-Stat configuration file: PointStatConfig_default

I'll try running your data here and figure out what's going on.

Thanks,
John

On 09/09/2013 09:50 AM, Michael Weston via RT wrote:
>
> Mon Sep 09 09:50:52 2013: Request 62959 was acted upon.
> Transaction: Ticket created by michael at escience.co.za
>         Queue: met_help
>       Subject: POINT_STAT Error APCP
>         Owner: Nobody
>    Requestors: michael at escience.co.za
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>
>
> Hello MET Helpdesk
>
> I am trying to use point_stat to verify WRF against rain gauge data.
> I am using METv4.0.
> When I run
> point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc
./PointStatConfig_default
> I get an error message:
>
> DEBUG 1: Default Config File:
/media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_default
> ERROR  :
> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'APCP_surface' for table version 2
> ERROR  :
>
> I have post processed WRF with UPPv1.1
> Which created a grib file called WRFPRS_d01.023.
> The total precipitation variable in the grib file is called
"APCP_surface".
>
> If I try with fcst variable APCP, I get an error message:
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: ./WRFPRS_d01.023
> If I try A1 instead of A24 I get.
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1
found in file: ./WRFPRS_d01.023
>
> I would lke to use MET to verify rainfall but it seams impossible.
> Have I done something wrong with running UPPv1.1?
>>From my wrf_cntrl.parm
> " (ACM TOTAL PRECIP    ) SCAL=( 5.0)
>   L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000 00000 00000 00000)"
>
> Any help would be appreciated.
>
> Thanks
> Michael Weston
>
> My config files looks like this:
> // Forecast and observation fields to be verified
> //
> fcst = {
>     wind_thresh  = [ NA ];
>     message_type = [ "ADPSFC" ];
>
>     field = [
>        {
>          name       = "APCP_surface";
>          level      = "A24";
>          cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>        }
>
>     ];
>
> };
> obs = {
>     wind_thresh  = [ NA ];
>     message_type = [ "ADPSFC" ];
>
>     field = [
>        {
>          name       = "APCP";
>          level      = [ "A24" ];
>          cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>        }
>
>     ];
>
> };
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
From: Julie Prestopnik
Time: Tue Sep 10 10:51:41 2013

Hi Michael.

Thank you for providing your data.

Using "wgrib" (wgrib WRFPRS_d01.023 | grep APCP) to view the fields in
forecast file, I see that the file contains the 23 hour forecast:

278:9741704:d=00092800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
acc:NAve=0

Given that information, you would want change the name from
"APCP_surface" to "APCP" in the fcst section of the config file and
provide the 24 hour forecast file instead or update the name and the
level if you wish to compare with the 23 hour forecast

name       = "APCP_surface";
level      = "A24";

to

name       = "APCP";
level      = "A23";

Using "ncdump" (ncdump -v obs_arr ascii2nc_lieb_2000.nc | more) to
view
the obs_arr data, I see that the file contains a 1 hour (3600s)
accumulation interval:

   0, 61, 3600, -9999, 0,

The level in your config file in the obs section was

level      = [ "A24" ];

but needs to be

level      = [ "A1" ];

While you could compare the 24 hour forecast with the one hour
observations, it probably doesn't make much sense to do that.  I would
suggest verifying against the one hour forecast file or summing up
your
point observations into 24 hour accumulations.

I hope this helps.

Julie


On 09/10/2013 02:48 AM, Michael Weston via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>
> Hi John,
>
> Thanks for your reply.
> After editing the config file so that obs=fcst I receive the
following error message:
>
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: WRFPRS_d01.023
>
> I have uploaded my data to /incoming/irap/met_help/weston_data/
> I think the units in my ascii file are incorrect for now. But I am
just trying to get the point_stat to run at this stage.
>
> Regards
> Michael
>
> EScience Associates (Pty) Ltd
> E-mail: michael at escience.co.za
> Web:   www.escience.co.za
> PO Box 2950, Saxonwold, 2132
> 9 Victoria Street, Oaklands, Johannesburg, 2192
> Tel:  +27 (0)11 718 6380
> Fax: 086 527 5241
> VAT No: 473 025 4416
> Reg No: 2009/014472/07
> ________________________________________
> From: John Halley Gotway via RT [met_help at ucar.edu]
> Sent: Monday, September 09, 2013 6:25 PM
> To: Michael Weston
> Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
>
> Michael,
>
> Can you please try the following first... just edit your config file
like this:
>
> //
> // Forecast and observation fields to be verified
> //
> fcst = {
>      wind_thresh  = [ NA ];
>      message_type = [ "ADPSFC" ];
>
>      field = [
>         {
>           name       = "APCP";
>           level      = "A24";
>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>         }
>
>      ];
>
> };
> obs = fcst;
>
> Then try running Point-Stat.  If that doesn't fix the problem, can
you please send us the data you're passing to Point-Stat?  You can
follow the instructions listed here to post data to our anonymous
> ftp site:
>      http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please send us:
> (1) GRIB forecast file: WRFPRS_d01.023
> (2) NetCDF point observation file: ascii2nc_lieb_2000.nc
> (3) Point-Stat configuration file: PointStatConfig_default
>
> I'll try running your data here and figure out what's going on.
>
> Thanks,
> John
>
> On 09/09/2013 09:50 AM, Michael Weston via RT wrote:
>>
>> Mon Sep 09 09:50:52 2013: Request 62959 was acted upon.
>> Transaction: Ticket created by michael at escience.co.za
>>          Queue: met_help
>>        Subject: POINT_STAT Error APCP
>>          Owner: Nobody
>>     Requestors: michael at escience.co.za
>>         Status: new
>>    Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>>
>>
>> Hello MET Helpdesk
>>
>> I am trying to use point_stat to verify WRF against rain gauge
data.
>> I am using METv4.0.
>> When I run
>> point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc
./PointStatConfig_default
>> I get an error message:
>>
>> DEBUG 1: Default Config File:
/media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
>> DEBUG 1: User Config File: ./PointStatConfig_default
>> ERROR  :
>> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'APCP_surface' for table version 2
>> ERROR  :
>>
>> I have post processed WRF with UPPv1.1
>> Which created a grib file called WRFPRS_d01.023.
>> The total precipitation variable in the grib file is called
"APCP_surface".
>>
>> If I try with fcst variable APCP, I get an error message:
>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: ./WRFPRS_d01.023
>> If I try A1 instead of A24 I get.
>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1
found in file: ./WRFPRS_d01.023
>>
>> I would lke to use MET to verify rainfall but it seams impossible.
>> Have I done something wrong with running UPPv1.1?
>> >From my wrf_cntrl.parm
>> " (ACM TOTAL PRECIP    ) SCAL=( 5.0)
>>    L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000 00000 00000 00000)"
>>
>> Any help would be appreciated.
>>
>> Thanks
>> Michael Weston
>>
>> My config files looks like this:
>> // Forecast and observation fields to be verified
>> //
>> fcst = {
>>      wind_thresh  = [ NA ];
>>      message_type = [ "ADPSFC" ];
>>
>>      field = [
>>         {
>>           name       = "APCP_surface";
>>           level      = "A24";
>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>         }
>>
>>      ];
>>
>> };
>> obs = {
>>      wind_thresh  = [ NA ];
>>      message_type = [ "ADPSFC" ];
>>
>>      field = [
>>         {
>>           name       = "APCP";
>>           level      = [ "A24" ];
>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>         }
>>
>>      ];
>>
>> };
>>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
From: Michael Weston
Time: Wed Sep 11 07:54:06 2013

Thanks for your help guys!
On the right track now.

-Michael
________________________________________
From: Julie Prestopnik via RT [met_help at ucar.edu]
Sent: Tuesday, September 10, 2013 6:51 PM
To: Michael Weston
Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP

Hi Michael.

Thank you for providing your data.

Using "wgrib" (wgrib WRFPRS_d01.023 | grep APCP) to view the fields in
forecast file, I see that the file contains the 23 hour forecast:

278:9741704:d=00092800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
acc:NAve=0

Given that information, you would want change the name from
"APCP_surface" to "APCP" in the fcst section of the config file and
provide the 24 hour forecast file instead or update the name and the
level if you wish to compare with the 23 hour forecast

name       = "APCP_surface";
level      = "A24";

to

name       = "APCP";
level      = "A23";

Using "ncdump" (ncdump -v obs_arr ascii2nc_lieb_2000.nc | more) to
view
the obs_arr data, I see that the file contains a 1 hour (3600s)
accumulation interval:

   0, 61, 3600, -9999, 0,

The level in your config file in the obs section was

level      = [ "A24" ];

but needs to be

level      = [ "A1" ];

While you could compare the 24 hour forecast with the one hour
observations, it probably doesn't make much sense to do that.  I would
suggest verifying against the one hour forecast file or summing up
your
point observations into 24 hour accumulations.

I hope this helps.

Julie


On 09/10/2013 02:48 AM, Michael Weston via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>
> Hi John,
>
> Thanks for your reply.
> After editing the config file so that obs=fcst I receive the
following error message:
>
> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: WRFPRS_d01.023
>
> I have uploaded my data to /incoming/irap/met_help/weston_data/
> I think the units in my ascii file are incorrect for now. But I am
just trying to get the point_stat to run at this stage.
>
> Regards
> Michael
>
> EScience Associates (Pty) Ltd
> E-mail: michael at escience.co.za
> Web:   www.escience.co.za
> PO Box 2950, Saxonwold, 2132
> 9 Victoria Street, Oaklands, Johannesburg, 2192
> Tel:  +27 (0)11 718 6380
> Fax: 086 527 5241
> VAT No: 473 025 4416
> Reg No: 2009/014472/07
> ________________________________________
> From: John Halley Gotway via RT [met_help at ucar.edu]
> Sent: Monday, September 09, 2013 6:25 PM
> To: Michael Weston
> Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
>
> Michael,
>
> Can you please try the following first... just edit your config file
like this:
>
> //
> // Forecast and observation fields to be verified
> //
> fcst = {
>      wind_thresh  = [ NA ];
>      message_type = [ "ADPSFC" ];
>
>      field = [
>         {
>           name       = "APCP";
>           level      = "A24";
>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>         }
>
>      ];
>
> };
> obs = fcst;
>
> Then try running Point-Stat.  If that doesn't fix the problem, can
you please send us the data you're passing to Point-Stat?  You can
follow the instructions listed here to post data to our anonymous
> ftp site:
>      http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please send us:
> (1) GRIB forecast file: WRFPRS_d01.023
> (2) NetCDF point observation file: ascii2nc_lieb_2000.nc
> (3) Point-Stat configuration file: PointStatConfig_default
>
> I'll try running your data here and figure out what's going on.
>
> Thanks,
> John
>
> On 09/09/2013 09:50 AM, Michael Weston via RT wrote:
>>
>> Mon Sep 09 09:50:52 2013: Request 62959 was acted upon.
>> Transaction: Ticket created by michael at escience.co.za
>>          Queue: met_help
>>        Subject: POINT_STAT Error APCP
>>          Owner: Nobody
>>     Requestors: michael at escience.co.za
>>         Status: new
>>    Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>>
>>
>> Hello MET Helpdesk
>>
>> I am trying to use point_stat to verify WRF against rain gauge
data.
>> I am using METv4.0.
>> When I run
>> point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc
./PointStatConfig_default
>> I get an error message:
>>
>> DEBUG 1: Default Config File:
/media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
>> DEBUG 1: User Config File: ./PointStatConfig_default
>> ERROR  :
>> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'APCP_surface' for table version 2
>> ERROR  :
>>
>> I have post processed WRF with UPPv1.1
>> Which created a grib file called WRFPRS_d01.023.
>> The total precipitation variable in the grib file is called
"APCP_surface".
>>
>> If I try with fcst variable APCP, I get an error message:
>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
found in file: ./WRFPRS_d01.023
>> If I try A1 instead of A24 I get.
>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1
found in file: ./WRFPRS_d01.023
>>
>> I would lke to use MET to verify rainfall but it seams impossible.
>> Have I done something wrong with running UPPv1.1?
>> >From my wrf_cntrl.parm
>> " (ACM TOTAL PRECIP    ) SCAL=( 5.0)
>>    L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000 00000 00000 00000)"
>>
>> Any help would be appreciated.
>>
>> Thanks
>> Michael Weston
>>
>> My config files looks like this:
>> // Forecast and observation fields to be verified
>> //
>> fcst = {
>>      wind_thresh  = [ NA ];
>>      message_type = [ "ADPSFC" ];
>>
>>      field = [
>>         {
>>           name       = "APCP_surface";
>>           level      = "A24";
>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>         }
>>
>>      ];
>>
>> };
>> obs = {
>>      wind_thresh  = [ NA ];
>>      message_type = [ "ADPSFC" ];
>>
>>      field = [
>>         {
>>           name       = "APCP";
>>           level      = [ "A24" ];
>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>         }
>>
>>      ];
>>
>> };
>>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu



------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
From: Julie Prestopnik
Time: Wed Sep 11 09:28:09 2013

Great, Michael!  Thanks for the update.

Julie

On 2013-09-11 07:54, Michael Weston via RT wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>
> Thanks for your help guys!
> On the right track now.
>
> -Michael
> ________________________________________
> From: Julie Prestopnik via RT [met_help at ucar.edu]
> Sent: Tuesday, September 10, 2013 6:51 PM
> To: Michael Weston
> Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
>
> Hi Michael.
>
> Thank you for providing your data.
>
> Using "wgrib" (wgrib WRFPRS_d01.023 | grep APCP) to view the fields
> in
> forecast file, I see that the file contains the 23 hour forecast:
>
>
>
278:9741704:d=00092800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
> acc:NAve=0
>
> Given that information, you would want change the name from
> "APCP_surface" to "APCP" in the fcst section of the config file and
> provide the 24 hour forecast file instead or update the name and the
> level if you wish to compare with the 23 hour forecast
>
> name       = "APCP_surface";
> level      = "A24";
>
> to
>
> name       = "APCP";
> level      = "A23";
>
> Using "ncdump" (ncdump -v obs_arr ascii2nc_lieb_2000.nc | more) to
> view
> the obs_arr data, I see that the file contains a 1 hour (3600s)
> accumulation interval:
>
>    0, 61, 3600, -9999, 0,
>
> The level in your config file in the obs section was
>
> level      = [ "A24" ];
>
> but needs to be
>
> level      = [ "A1" ];
>
> While you could compare the 24 hour forecast with the one hour
> observations, it probably doesn't make much sense to do that.  I
> would
> suggest verifying against the one hour forecast file or summing up
> your
> point observations into 24 hour accumulations.
>
> I hope this helps.
>
> Julie
>
>
> On 09/10/2013 02:48 AM, Michael Weston via RT wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>>
>> Hi John,
>>
>> Thanks for your reply.
>> After editing the config file so that obs=fcst I receive the
>> following error message:
>>
>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
>> found in file: WRFPRS_d01.023
>>
>> I have uploaded my data to /incoming/irap/met_help/weston_data/
>> I think the units in my ascii file are incorrect for now. But I am
>> just trying to get the point_stat to run at this stage.
>>
>> Regards
>> Michael
>>
>> EScience Associates (Pty) Ltd
>> E-mail: michael at escience.co.za
>> Web:   www.escience.co.za
>> PO Box 2950, Saxonwold, 2132
>> 9 Victoria Street, Oaklands, Johannesburg, 2192
>> Tel:  +27 (0)11 718 6380
>> Fax: 086 527 5241
>> VAT No: 473 025 4416
>> Reg No: 2009/014472/07
>> ________________________________________
>> From: John Halley Gotway via RT [met_help at ucar.edu]
>> Sent: Monday, September 09, 2013 6:25 PM
>> To: Michael Weston
>> Subject: Re: [rt.rap.ucar.edu #62959] POINT_STAT Error APCP
>>
>> Michael,
>>
>> Can you please try the following first... just edit your config
file
>> like this:
>>
>> //
>> // Forecast and observation fields to be verified
>> //
>> fcst = {
>>      wind_thresh  = [ NA ];
>>      message_type = [ "ADPSFC" ];
>>
>>      field = [
>>         {
>>           name       = "APCP";
>>           level      = "A24";
>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>         }
>>
>>      ];
>>
>> };
>> obs = fcst;
>>
>> Then try running Point-Stat.  If that doesn't fix the problem, can
>> you please send us the data you're passing to Point-Stat?  You can
>> follow the instructions listed here to post data to our anonymous
>> ftp site:
>>      http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> Please send us:
>> (1) GRIB forecast file: WRFPRS_d01.023
>> (2) NetCDF point observation file: ascii2nc_lieb_2000.nc
>> (3) Point-Stat configuration file: PointStatConfig_default
>>
>> I'll try running your data here and figure out what's going on.
>>
>> Thanks,
>> John
>>
>> On 09/09/2013 09:50 AM, Michael Weston via RT wrote:
>>>
>>> Mon Sep 09 09:50:52 2013: Request 62959 was acted upon.
>>> Transaction: Ticket created by michael at escience.co.za
>>>          Queue: met_help
>>>        Subject: POINT_STAT Error APCP
>>>          Owner: Nobody
>>>     Requestors: michael at escience.co.za
>>>         Status: new
>>>    Ticket <URL:
>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62959 >
>>>
>>>
>>> Hello MET Helpdesk
>>>
>>> I am trying to use point_stat to verify WRF against rain gauge
>>> data.
>>> I am using METv4.0.
>>> When I run
>>> point_stat ./WRFPRS_d01.023 ./ascii2nc_lieb_2000.nc
>>> ./PointStatConfig_default
>>> I get an error message:
>>>
>>> DEBUG 1: Default Config File:
>>>
/media/prod/chris/WRFv3.4.1/METv4.0/data/config/PointStatConfig_default
>>> DEBUG 1: User Config File: ./PointStatConfig_default
>>> ERROR  :
>>> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>>> abbreviation 'APCP_surface' for table version 2
>>> ERROR  :
>>>
>>> I have post processed WRF with UPPv1.1
>>> Which created a grib file called WRFPRS_d01.023.
>>> The total precipitation variable in the grib file is called
>>> "APCP_surface".
>>>
>>> If I try with fcst variable APCP, I get an error message:
>>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A24
>>> found in file: ./WRFPRS_d01.023
>>> If I try A1 instead of A24 I get.
>>> ERROR  : process_fcst_climo_files() -> no fields matching APCP/A1
>>> found in file: ./WRFPRS_d01.023
>>>
>>> I would lke to use MET to verify rainfall but it seams impossible.
>>> Have I done something wrong with running UPPv1.1?
>>> >From my wrf_cntrl.parm
>>> " (ACM TOTAL PRECIP    ) SCAL=( 5.0)
>>>    L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
>>> 00000 00000 00000 00000)"
>>>
>>> Any help would be appreciated.
>>>
>>> Thanks
>>> Michael Weston
>>>
>>> My config files looks like this:
>>> // Forecast and observation fields to be verified
>>> //
>>> fcst = {
>>>      wind_thresh  = [ NA ];
>>>      message_type = [ "ADPSFC" ];
>>>
>>>      field = [
>>>         {
>>>           name       = "APCP_surface";
>>>           level      = "A24";
>>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>>         }
>>>
>>>      ];
>>>
>>> };
>>> obs = {
>>>      wind_thresh  = [ NA ];
>>>      message_type = [ "ADPSFC" ];
>>>
>>>      field = [
>>>         {
>>>           name       = "APCP";
>>>           level      = [ "A24" ];
>>>           cat_thresh = [ >0, >1, >5, >10, >25, >25 ];
>>>         }
>>>
>>>      ];
>>>
>>> };
>>>
>>
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu


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


More information about the Met_help mailing list