[Met_help] [rt.rap.ucar.edu #74609] History for Problems running grid_stat precipitation analyses
John Halley Gotway via RT
met_help at ucar.edu
Fri Jan 8 09:45:03 MST 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Met help,
I'm trying to run MET precipitation analyses for NMMB, and get the following error:
DEBUG 1: Forecast File: nmbprs_d01.12
DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
WARNING:
WARNING:
WARNING: process_scores() -> APCP/A03 not found in file: nmbprs_d01.12
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
The namelist has the following:
field = [
{
name = "APCP";
level = [ "A03" ];
cat_thresh = [ >=0,>=12, >=24 ];
GRIB1_ptv = 129;
I also tried "L0", "Z0" , "A12" and "A06"
Any help will be greatly appreciated.
Jose M. Garcia-Rivera
This e-mail, including attachments, may include confidential and/or proprietary information, and may be used only by the person or entity to which it is addressed. If the reader of this e-mail is not the intended recipient or his or her authorized agent, the reader is hereby notified that any dissemination, distribution or copying of this e-mail is prohibited. If you have received this e-mail in error, please notify the sender by replying to this message and delete this e-mail immediately.
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 09:49:39 2016
Jose,
I see that you're having a hard time extracting 3-hourly precip from
an
NMMB GRIB file. The question is if and how the accumulated
precipitation
is stored in that GRIB file. Try running the following wgrib command
to
see what APCP records it contains:
wgrib nmbprs_d01.12 | grep APCP
Please send me the wgrib output, and that should tell us more about
what
your GRIB file contains.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> Transaction: Ticket created by GarciaJ at imsg.com
> Queue: met_help
> Subject: Problems running grid_stat precipitation analyses
> Owner: Nobody
> Requestors: GarciaJ at imsg.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
>
> Dear Met help,
>
>
> I'm trying to run MET precipitation analyses for NMMB, and get the
> following error:
>
> DEBUG 1: Forecast File: nmbprs_d01.12
> DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> WARNING:
> WARNING:
> WARNING: process_scores() -> APCP/A03 not found in file:
nmbprs_d01.12
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
>
>
> The namelist has the following:
>
> field = [
> {
> name = "APCP";
> level = [ "A03" ];
> cat_thresh = [ >=0,>=12, >=24 ];
> GRIB1_ptv = 129;
>
> I also tried "L0", "Z0" , "A12" and "A06"
>
> Any help will be greatly appreciated.
>
>
> Jose M. Garcia-Rivera
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Wed Jan 06 09:52:17 2016
Dear John,
Thanks for the reply. I've been playing with the var info using wgrib,
but can't make grid_stat work. Here's the info you want:
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:TimeU=1:sfc:9-
12hr acc:NAve=0
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 11:50 AM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
Jose,
I see that you're having a hard time extracting 3-hourly precip from
an NMMB GRIB file. The question is if and how the accumulated
precipitation is stored in that GRIB file. Try running the following
wgrib command to see what APCP records it contains:
wgrib nmbprs_d01.12 | grep APCP
Please send me the wgrib output, and that should tell us more about
what your GRIB file contains.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> Transaction: Ticket created by GarciaJ at imsg.com
> Queue: met_help
> Subject: Problems running grid_stat precipitation analyses
> Owner: Nobody
> Requestors: GarciaJ at imsg.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> >
>
>
> Dear Met help,
>
>
> I'm trying to run MET precipitation analyses for NMMB, and get the
> following error:
>
> DEBUG 1: Forecast File: nmbprs_d01.12
> DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> WARNING:
> WARNING:
> WARNING: process_scores() -> APCP/A03 not found in file:
nmbprs_d01.12
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
>
>
> The namelist has the following:
>
> field = [
> {
> name = "APCP";
> level = [ "A03" ];
> cat_thresh = [ >=0,>=12, >=24 ];
> GRIB1_ptv = 129;
>
> I also tried "L0", "Z0" , "A12" and "A06"
>
> Any help will be greatly appreciated.
>
>
> Jose M. Garcia-Rivera
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 10:50:44 2016
Jose,
Please just try removing the "GRIB1_ptv = 129;" setting or change it
to the
default PTV value of 2. Then it should work.
The reason why it's not working now is that your GRIB record is most
likely
encoded using parameter table version number 2. For example, try
running
this command:
wgrib -d 304 -V -PDS10 nmbprs_d01.12
You'll see in the output something like this:
PDS10(1..28)=
1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8 7
0
0 1 0 3
21: 4 0 0 0 21 0 0 3
The PDS is the product description section and the 4th value listed is
probably the number 2, as in the example above. Looking at the PDS
definition table
(http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
the 4th value is the parameter table version number.
When you set "GRIB1_ptv = 129;", MET will skip over that record
because the
PTV value in the data (2) doesn't match the value you specified (129).
Hope that helps.
Thanks,
Joihn
On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> Thanks for the reply. I've been playing with the var info using
wgrib, but
> can't make grid_stat work. Here's the info you want:
>
>
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:TimeU=1:sfc:9-
12hr
> acc:NAve=0
>
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 11:50 AM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> I see that you're having a hard time extracting 3-hourly precip from
an
> NMMB GRIB file. The question is if and how the accumulated
precipitation
> is stored in that GRIB file. Try running the following wgrib
command to
> see what APCP records it contains:
>
> wgrib nmbprs_d01.12 | grep APCP
>
> Please send me the wgrib output, and that should tell us more about
what
> your GRIB file contains.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > Transaction: Ticket created by GarciaJ at imsg.com
> > Queue: met_help
> > Subject: Problems running grid_stat precipitation analyses
> > Owner: Nobody
> > Requestors: GarciaJ at imsg.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > >
> >
> >
> > Dear Met help,
> >
> >
> > I'm trying to run MET precipitation analyses for NMMB, and get the
> > following error:
> >
> > DEBUG 1: Forecast File: nmbprs_d01.12
> > DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > WARNING:
> > WARNING:
> > WARNING: process_scores() -> APCP/A03 not found in file:
nmbprs_d01.12
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> >
> >
> > The namelist has the following:
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "A03" ];
> > cat_thresh = [ >=0,>=12, >=24 ];
> > GRIB1_ptv = 129;
> >
> > I also tried "L0", "Z0" , "A12" and "A06"
> >
> > Any help will be greatly appreciated.
> >
> >
> > Jose M. Garcia-Rivera
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Wed Jan 06 11:06:43 2016
Dear John,
Thanks for the help, now it's working for the NMMB file. However, the
error now shows up with the NARR data file:
DEBUG 1: Forecast File: nmbprs_d01.12
DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
VarInfo "APCP/A03" in GRIB file
"narr.2015102812.RS.sfc_regridded.grb".
WARNING:
WARNING:
WARNING: process_scores() -> APCP/A03 not found in file:
narr.2015102812.RS.sfc_regridded.grb
Running wgrib for NARR:
wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr acc:NAve=0
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr acc:NAve=0
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 12:51 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
Jose,
Please just try removing the "GRIB1_ptv = 129;" setting or change it
to the default PTV value of 2. Then it should work.
The reason why it's not working now is that your GRIB record is most
likely encoded using parameter table version number 2. For example,
try running this command:
wgrib -d 304 -V -PDS10 nmbprs_d01.12
You'll see in the output something like this:
PDS10(1..28)=
1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8 7
0
0 1 0 3
21: 4 0 0 0 21 0 0 3
The PDS is the product description section and the 4th value listed is
probably the number 2, as in the example above. Looking at the PDS
definition table
(http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
the 4th value is the parameter table version number.
When you set "GRIB1_ptv = 129;", MET will skip over that record
because the PTV value in the data (2) doesn't match the value you
specified (129).
Hope that helps.
Thanks,
Joihn
On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> Thanks for the reply. I've been playing with the var info using
wgrib,
> but can't make grid_stat work. Here's the info you want:
>
>
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> TimeU=1:sfc:9-12hr
> acc:NAve=0
>
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 11:50 AM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> I see that you're having a hard time extracting 3-hourly precip from
> an NMMB GRIB file. The question is if and how the accumulated
> precipitation is stored in that GRIB file. Try running the
following
> wgrib command to see what APCP records it contains:
>
> wgrib nmbprs_d01.12 | grep APCP
>
> Please send me the wgrib output, and that should tell us more about
> what your GRIB file contains.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > Transaction: Ticket created by GarciaJ at imsg.com
> > Queue: met_help
> > Subject: Problems running grid_stat precipitation analyses
> > Owner: Nobody
> > Requestors: GarciaJ at imsg.com
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > >
> >
> >
> > Dear Met help,
> >
> >
> > I'm trying to run MET precipitation analyses for NMMB, and get the
> > following error:
> >
> > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > narr.2015102812.RS.sfc_regridded.grb
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > WARNING:
> > WARNING:
> > WARNING: process_scores() -> APCP/A03 not found in file:
> > nmbprs_d01.12
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------
> > --
> > ----------
> > DEBUG 2:
> >
> >
> > The namelist has the following:
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "A03" ];
> > cat_thresh = [ >=0,>=12, >=24 ];
> > GRIB1_ptv = 129;
> >
> > I also tried "L0", "Z0" , "A12" and "A06"
> >
> > Any help will be greatly appreciated.
> >
> >
> > Jose M. Garcia-Rivera
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 11:42:01 2016
Jose,
It's probably be easiest to just have you send me sample files.
Please
post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb, and your
Grid-Stat configuration file to our anonymous ftp server following
these
instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Once I hear back from you that the data is posted, I'll go grab it and
check for issues.
Thanks,
John
On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> Thanks for the help, now it's working for the NMMB file. However,
the
> error now shows up with the NARR data file:
>
> DEBUG 1: Forecast File: nmbprs_d01.12
> DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "APCP/A03" in GRIB file
"narr.2015102812.RS.sfc_regridded.grb".
> WARNING:
> WARNING:
> WARNING: process_scores() -> APCP/A03 not found in file:
narr.2015102812
> .RS.sfc_regridded.grb
>
>
> Running wgrib for NARR:
> wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
>
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
> acc:NAve=0
>
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
> acc:NAve=0
>
> Best Regards,
> Jose
>
>
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 12:51 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> Please just try removing the "GRIB1_ptv = 129;" setting or change it
to
> the default PTV value of 2. Then it should work.
>
> The reason why it's not working now is that your GRIB record is most
> likely encoded using parameter table version number 2. For example,
try
> running this command:
>
> wgrib -d 304 -V -PDS10 nmbprs_d01.12
>
> You'll see in the output something like this:
> PDS10(1..28)=
> 1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8
7 0
> 0 1 0 3
> 21: 4 0 0 0 21 0 0 3
>
> The PDS is the product description section and the 4th value listed
is
> probably the number 2, as in the example above. Looking at the PDS
> definition table (
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> the 4th value is the parameter table version number.
>
> When you set "GRIB1_ptv = 129;", MET will skip over that record
because
> the PTV value in the data (2) doesn't match the value you specified
(129).
>
> Hope that helps.
>
> Thanks,
> Joihn
>
> On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > Thanks for the reply. I've been playing with the var info using
wgrib,
> > but can't make grid_stat work. Here's the info you want:
> >
> >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > TimeU=1:sfc:9-12hr
> > acc:NAve=0
> >
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 11:50 AM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > I see that you're having a hard time extracting 3-hourly precip
from
> > an NMMB GRIB file. The question is if and how the accumulated
> > precipitation is stored in that GRIB file. Try running the
following
> > wgrib command to see what APCP records it contains:
> >
> > wgrib nmbprs_d01.12 | grep APCP
> >
> > Please send me the wgrib output, and that should tell us more
about
> > what your GRIB file contains.
> >
> > Thanks,
> > John Halley Gotway
> > met_help at ucar.edu
> >
> > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > Transaction: Ticket created by GarciaJ at imsg.com
> > > Queue: met_help
> > > Subject: Problems running grid_stat precipitation analyses
> > > Owner: Nobody
> > > Requestors: GarciaJ at imsg.com
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > >
> > >
> > >
> > > Dear Met help,
> > >
> > >
> > > I'm trying to run MET precipitation analyses for NMMB, and get
the
> > > following error:
> > >
> > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > > narr.2015102812.RS.sfc_regridded.grb
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > > VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > WARNING:
> > > WARNING:
> > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > nmbprs_d01.12
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
--------------------------------------------------------------------
> > > --
> > > ----------
> > > DEBUG 2:
> > >
> > >
> > > The namelist has the following:
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "A03" ];
> > > cat_thresh = [ >=0,>=12, >=24 ];
> > > GRIB1_ptv = 129;
> > >
> > > I also tried "L0", "Z0" , "A12" and "A06"
> > >
> > > Any help will be greatly appreciated.
> > >
> > >
> > > Jose M. Garcia-Rivera
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Wed Jan 06 12:18:23 2016
Dear John,
The files are now at 'garcia_data_new'
Thanks for the help,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 1:42 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
Jose,
It's probably be easiest to just have you send me sample files.
Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb, and
your Grid-Stat configuration file to our anonymous ftp server
following these
instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Once I hear back from you that the data is posted, I'll go grab it and
check for issues.
Thanks,
John
On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> Thanks for the help, now it's working for the NMMB file. However,
the
> error now shows up with the NARR data file:
>
> DEBUG 1: Forecast File: nmbprs_d01.12
> DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "APCP/A03" in GRIB file
"narr.2015102812.RS.sfc_regridded.grb".
> WARNING:
> WARNING:
> WARNING: process_scores() -> APCP/A03 not found in file:
> narr.2015102812 .RS.sfc_regridded.grb
>
>
> Running wgrib for NARR:
> wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
>
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:Tim
> eU=1:sfc:0-3hr
> acc:NAve=0
>
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> imeU=1:sfc:0-3hr
> acc:NAve=0
>
> Best Regards,
> Jose
>
>
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 12:51 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> Please just try removing the "GRIB1_ptv = 129;" setting or change it
> to the default PTV value of 2. Then it should work.
>
> The reason why it's not working now is that your GRIB record is most
> likely encoded using parameter table version number 2. For example,
> try running this command:
>
> wgrib -d 304 -V -PDS10 nmbprs_d01.12
>
> You'll see in the output something like this:
> PDS10(1..28)=
> 1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8
7 0
> 0 1 0 3
> 21: 4 0 0 0 21 0 0 3
>
> The PDS is the product description section and the 4th value listed
is
> probably the number 2, as in the example above. Looking at the PDS
> definition table (
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> the 4th value is the parameter table version number.
>
> When you set "GRIB1_ptv = 129;", MET will skip over that record
> because the PTV value in the data (2) doesn't match the value you
specified (129).
>
> Hope that helps.
>
> Thanks,
> Joihn
>
> On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > Thanks for the reply. I've been playing with the var info using
> > wgrib, but can't make grid_stat work. Here's the info you want:
> >
> >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > TimeU=1:sfc:9-12hr
> > acc:NAve=0
> >
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 11:50 AM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > I see that you're having a hard time extracting 3-hourly precip
from
> > an NMMB GRIB file. The question is if and how the accumulated
> > precipitation is stored in that GRIB file. Try running the
> > following wgrib command to see what APCP records it contains:
> >
> > wgrib nmbprs_d01.12 | grep APCP
> >
> > Please send me the wgrib output, and that should tell us more
about
> > what your GRIB file contains.
> >
> > Thanks,
> > John Halley Gotway
> > met_help at ucar.edu
> >
> > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > Transaction: Ticket created by GarciaJ at imsg.com
> > > Queue: met_help
> > > Subject: Problems running grid_stat precipitation analyses
> > > Owner: Nobody
> > > Requestors: GarciaJ at imsg.com
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > >
> > >
> > >
> > > Dear Met help,
> > >
> > >
> > > I'm trying to run MET precipitation analyses for NMMB, and get
the
> > > following error:
> > >
> > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > > narr.2015102812.RS.sfc_regridded.grb
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
> > > for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > WARNING:
> > > WARNING:
> > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > nmbprs_d01.12
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
------------------------------------------------------------------
> > > --
> > > --
> > > ----------
> > > DEBUG 2:
> > >
> > >
> > > The namelist has the following:
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "A03" ];
> > > cat_thresh = [ >=0,>=12, >=24 ];
> > > GRIB1_ptv = 129;
> > >
> > > I also tried "L0", "Z0" , "A12" and "A06"
> > >
> > > Any help will be greatly appreciated.
> > >
> > >
> > > Jose M. Garcia-Rivera
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 13:56:53 2016
Jose,
For whatever reason that "narr.2015102812.RS.sfc_regridded.grb" file
encodes APCP using GRIB parameter table version number 131, as shown
in the
wgrib command below:
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1 0 0 15 10 28 12 0 1 0
3 4 0
0 0 21 15 0 0
Notice that the 4th value in "PDS10=" list is 131.
I often use the plot_data_plane utility to test out this sort of
thing.
Try running the following commands:
met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
'name="APCP";
level="A03";
met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
GRIB1_ptv =
131;'
The resulting images are attached.
In your Grid-Stat config file, rather than saying "obs = fcst" just
set
"obs" to the same thing as the "fcst" but add "GRIB1_ptv = 131;".
Hopefully that does the trick.
John
On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> The files are now at 'garcia_data_new'
>
>
> Thanks for the help,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 1:42 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> It's probably be easiest to just have you send me sample files.
Please
> post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb, and your
> Grid-Stat configuration file to our anonymous ftp server following
these
> instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Once I hear back from you that the data is posted, I'll go grab it
and
> check for issues.
>
> Thanks,
> John
>
> On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > Thanks for the help, now it's working for the NMMB file. However,
the
> > error now shows up with the NARR data file:
> >
> > DEBUG 1: Forecast File: nmbprs_d01.12
> > DEBUG 1: Observation File: narr.2015102812.RS.sfc_regridded.grb
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "APCP/A03" in GRIB file
"narr.2015102812.RS.sfc_regridded.grb".
> > WARNING:
> > WARNING:
> > WARNING: process_scores() -> APCP/A03 not found in file:
> > narr.2015102812 .RS.sfc_regridded.grb
> >
> >
> > Running wgrib for NARR:
> > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:Tim
> > eU=1:sfc:0-3hr
> > acc:NAve=0
> >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> > imeU=1:sfc:0-3hr
> > acc:NAve=0
> >
> > Best Regards,
> > Jose
> >
> >
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 12:51 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > Please just try removing the "GRIB1_ptv = 129;" setting or change
it
> > to the default PTV value of 2. Then it should work.
> >
> > The reason why it's not working now is that your GRIB record is
most
> > likely encoded using parameter table version number 2. For
example,
> > try running this command:
> >
> > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> >
> > You'll see in the output something like this:
> > PDS10(1..28)=
> > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8
7 0
> > 0 1 0 3
> > 21: 4 0 0 0 21 0 0 3
> >
> > The PDS is the product description section and the 4th value
listed is
> > probably the number 2, as in the example above. Looking at the
PDS
> > definition table (
> > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > the 4th value is the parameter table version number.
> >
> > When you set "GRIB1_ptv = 129;", MET will skip over that record
> > because the PTV value in the data (2) doesn't match the value you
> specified (129).
> >
> > Hope that helps.
> >
> > Thanks,
> > Joihn
> >
> > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > Thanks for the reply. I've been playing with the var info using
> > > wgrib, but can't make grid_stat work. Here's the info you want:
> > >
> > >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > TimeU=1:sfc:9-12hr
> > > acc:NAve=0
> > >
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > I see that you're having a hard time extracting 3-hourly precip
from
> > > an NMMB GRIB file. The question is if and how the accumulated
> > > precipitation is stored in that GRIB file. Try running the
> > > following wgrib command to see what APCP records it contains:
> > >
> > > wgrib nmbprs_d01.12 | grep APCP
> > >
> > > Please send me the wgrib output, and that should tell us more
about
> > > what your GRIB file contains.
> > >
> > > Thanks,
> > > John Halley Gotway
> > > met_help at ucar.edu
> > >
> > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > Queue: met_help
> > > > Subject: Problems running grid_stat precipitation
analyses
> > > > Owner: Nobody
> > > > Requestors: GarciaJ at imsg.com
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > >
> > > >
> > > >
> > > > Dear Met help,
> > > >
> > > >
> > > > I'm trying to run MET precipitation analyses for NMMB, and get
the
> > > > following error:
> > > >
> > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > narr.2015102812.RS.sfc_regridded.grb
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > WARNING:
> > > > WARNING:
> > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > nmbprs_d01.12
> > > > WARNING:
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
------------------------------------------------------------------
> > > > --
> > > > --
> > > > ----------
> > > > DEBUG 2:
> > > >
> > > >
> > > > The namelist has the following:
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "A03" ];
> > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > GRIB1_ptv = 129;
> > > >
> > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > >
> > > > Any help will be greatly appreciated.
> > > >
> > > >
> > > > Jose M. Garcia-Rivera
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Wed Jan 06 15:03:21 2016
Dear John,
It worked! Your help is greatly appreciated.
I have a question regarding neighborhood methods:
We've been changing the parameters, and we keep getting the same HSS
values. How can we activate neighborhood methods in the namelist?
For example, we changed for Reflectivity to:
// Neighborhood methods
//
nbrhd = {
vld_thresh = 0.1;
width = [ 7 ];
cov_thresh = [ >=0.1 ];
}
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 3:57 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
Jose,
For whatever reason that "narr.2015102812.RS.sfc_regridded.grb" file
encodes APCP using GRIB parameter table version number 131, as shown
in the wgrib command below:
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1 0 0 15 10 28 12 0 1 0
3 4 0
0 0 21 15 0 0
Notice that the 4th value in "PDS10=" list is 131.
I often use the plot_data_plane utility to test out this sort of
thing.
Try running the following commands:
met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
'name="APCP"; level="A03";
met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
GRIB1_ptv = 131;'
The resulting images are attached.
In your Grid-Stat config file, rather than saying "obs = fcst" just
set "obs" to the same thing as the "fcst" but add "GRIB1_ptv = 131;".
Hopefully that does the trick.
John
On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> The files are now at 'garcia_data_new'
>
>
> Thanks for the help,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 1:42 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> It's probably be easiest to just have you send me sample files.
> Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb, and
> your Grid-Stat configuration file to our anonymous ftp server
> following these
> instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Once I hear back from you that the data is posted, I'll go grab it
and
> check for issues.
>
> Thanks,
> John
>
> On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > Thanks for the help, now it's working for the NMMB file. However,
> > the error now shows up with the NARR data file:
> >
> > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > narr.2015102812.RS.sfc_regridded.grb
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "APCP/A03" in GRIB file
"narr.2015102812.RS.sfc_regridded.grb".
> > WARNING:
> > WARNING:
> > WARNING: process_scores() -> APCP/A03 not found in file:
> > narr.2015102812 .RS.sfc_regridded.grb
> >
> >
> > Running wgrib for NARR:
> > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> > im
> > eU=1:sfc:0-3hr
> > acc:NAve=0
> >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3
> > :T
> > imeU=1:sfc:0-3hr
> > acc:NAve=0
> >
> > Best Regards,
> > Jose
> >
> >
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 12:51 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > Please just try removing the "GRIB1_ptv = 129;" setting or change
it
> > to the default PTV value of 2. Then it should work.
> >
> > The reason why it's not working now is that your GRIB record is
most
> > likely encoded using parameter table version number 2. For
example,
> > try running this command:
> >
> > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> >
> > You'll see in the output something like this:
> > PDS10(1..28)=
> > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5 8
7 0
> > 0 1 0 3
> > 21: 4 0 0 0 21 0 0 3
> >
> > The PDS is the product description section and the 4th value
listed
> > is probably the number 2, as in the example above. Looking at the
> > PDS definition table (
> > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > the 4th value is the parameter table version number.
> >
> > When you set "GRIB1_ptv = 129;", MET will skip over that record
> > because the PTV value in the data (2) doesn't match the value you
> specified (129).
> >
> > Hope that helps.
> >
> > Thanks,
> > Joihn
> >
> > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > Thanks for the reply. I've been playing with the var info using
> > > wgrib, but can't make grid_stat work. Here's the info you want:
> > >
> > >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > TimeU=1:sfc:9-12hr
> > > acc:NAve=0
> > >
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > I see that you're having a hard time extracting 3-hourly precip
> > > from an NMMB GRIB file. The question is if and how the
> > > accumulated precipitation is stored in that GRIB file. Try
> > > running the following wgrib command to see what APCP records it
contains:
> > >
> > > wgrib nmbprs_d01.12 | grep APCP
> > >
> > > Please send me the wgrib output, and that should tell us more
> > > about what your GRIB file contains.
> > >
> > > Thanks,
> > > John Halley Gotway
> > > met_help at ucar.edu
> > >
> > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > Queue: met_help
> > > > Subject: Problems running grid_stat precipitation
analyses
> > > > Owner: Nobody
> > > > Requestors: GarciaJ at imsg.com
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > >
> > > >
> > > >
> > > > Dear Met help,
> > > >
> > > >
> > > > I'm trying to run MET precipitation analyses for NMMB, and get
> > > > the following error:
> > > >
> > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > narr.2015102812.RS.sfc_regridded.grb
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > WARNING:
> > > > WARNING:
> > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > nmbprs_d01.12
> > > > WARNING:
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
----------------------------------------------------------------
> > > > --
> > > > --
> > > > --
> > > > ----------
> > > > DEBUG 2:
> > > >
> > > >
> > > > The namelist has the following:
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "A03" ];
> > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > GRIB1_ptv = 129;
> > > >
> > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > >
> > > > Any help will be greatly appreciated.
> > > >
> > > >
> > > > Jose M. Garcia-Rivera
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
> > > > person or entity to which it is addressed. If the reader of
this
> > > > e-mail is not the intended recipient or his or her authorized
> > > > agent, the reader is hereby notified that any dissemination,
> > > > distribution or copying of this e-mail is prohibited. If you
> > > > have received this e-mail in error, please notify the sender
by
> > > > replying to this message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 15:33:50 2016
Jose,
The neighborhood methods are used primarily to compute Fractions Skill
Score (FSS) which is found in the NBRCNT output line type.
The Heidke Skill Score can be found in the CTS (contingency table
statistics) output line type. It is derived from a 2x2 contingency
table.
If you'd like to look at FSS, I'd suggest:
- Turn on the NBRCNT output line type in the Grid-Stat config file:
nbrcnt
= STAT or BOTH;
- Specify several "width" settings in the "nbrhd" section:
nbrhd = {
vld_thresh = 0.1;
width = [ 3, 5, 7, 9, 11, 13 ];
cov_thresh = [ >=0.1 ];
}
People often see how FSS changes as you increase the neighborhood
size.
Hope that helps.
Thanks,
John
On Wed, Jan 6, 2016 at 3:03 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> It worked! Your help is greatly appreciated.
>
> I have a question regarding neighborhood methods:
> We've been changing the parameters, and we keep getting the same HSS
> values. How can we activate neighborhood methods in the namelist?
>
> For example, we changed for Reflectivity to:
>
> // Neighborhood methods
> //
> nbrhd = {
> vld_thresh = 0.1;
> width = [ 7 ];
> cov_thresh = [ >=0.1 ];
> }
>
>
> Best Regards,
> Jose
>
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 3:57 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> For whatever reason that "narr.2015102812.RS.sfc_regridded.grb" file
> encodes APCP using GRIB parameter table version number 131, as shown
in the
> wgrib command below:
>
>
>
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
> acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1 0 0 15 10 28 12 0 1
0 3 4 0
> 0 0 21 15 0 0
>
> Notice that the 4th value in "PDS10=" list is 131.
>
> I often use the plot_data_plane utility to test out this sort of
thing.
> Try running the following commands:
>
> met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
'name="APCP";
> level="A03";
>
> met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
> narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
GRIB1_ptv
> = 131;'
>
> The resulting images are attached.
>
> In your Grid-Stat config file, rather than saying "obs = fcst" just
set
> "obs" to the same thing as the "fcst" but add "GRIB1_ptv = 131;".
>
> Hopefully that does the trick.
>
> John
>
>
>
> On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > The files are now at 'garcia_data_new'
> >
> >
> > Thanks for the help,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 1:42 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > It's probably be easiest to just have you send me sample files.
> > Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb,
and
> > your Grid-Stat configuration file to our anonymous ftp server
> > following these
> > instructions:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Once I hear back from you that the data is posted, I'll go grab it
and
> > check for issues.
> >
> > Thanks,
> > John
> >
> > On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > Thanks for the help, now it's working for the NMMB file.
However,
> > > the error now shows up with the NARR data file:
> > >
> > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > > narr.2015102812.RS.sfc_regridded.grb
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > > VarInfo "APCP/A03" in GRIB file "narr.2015102812
> .RS.sfc_regridded.grb".
> > > WARNING:
> > > WARNING:
> > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > narr.2015102812 .RS.sfc_regridded.grb
> > >
> > >
> > > Running wgrib for NARR:
> > > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> > > im
> > > eU=1:sfc:0-3hr
> > > acc:NAve=0
> > >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2=3
> > > :T
> > > imeU=1:sfc:0-3hr
> > > acc:NAve=0
> > >
> > > Best Regards,
> > > Jose
> > >
> > >
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 12:51 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > Please just try removing the "GRIB1_ptv = 129;" setting or
change it
> > > to the default PTV value of 2. Then it should work.
> > >
> > > The reason why it's not working now is that your GRIB record is
most
> > > likely encoded using parameter table version number 2. For
example,
> > > try running this command:
> > >
> > > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> > >
> > > You'll see in the output something like this:
> > > PDS10(1..28)=
> > > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5
8 7
> 0
> > > 0 1 0 3
> > > 21: 4 0 0 0 21 0 0 3
> > >
> > > The PDS is the product description section and the 4th value
listed
> > > is probably the number 2, as in the example above. Looking at
the
> > > PDS definition table (
> > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > > the 4th value is the parameter table version number.
> > >
> > > When you set "GRIB1_ptv = 129;", MET will skip over that record
> > > because the PTV value in the data (2) doesn't match the value
you
> > specified (129).
> > >
> > > Hope that helps.
> > >
> > > Thanks,
> > > Joihn
> > >
> > > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
>
> > > >
> > > > Dear John,
> > > >
> > > > Thanks for the reply. I've been playing with the var info
using
> > > > wgrib, but can't make grid_stat work. Here's the info you
want:
> > > >
> > > >
>
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > > TimeU=1:sfc:9-12hr
> > > > acc:NAve=0
> > > >
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > precipitation analyses
> > > >
> > > > Jose,
> > > >
> > > > I see that you're having a hard time extracting 3-hourly
precip
> > > > from an NMMB GRIB file. The question is if and how the
> > > > accumulated precipitation is stored in that GRIB file. Try
> > > > running the following wgrib command to see what APCP records
it
> contains:
> > > >
> > > > wgrib nmbprs_d01.12 | grep APCP
> > > >
> > > > Please send me the wgrib output, and that should tell us more
> > > > about what your GRIB file contains.
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > > met_help at ucar.edu
> > > >
> > > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > Queue: met_help
> > > > > Subject: Problems running grid_stat precipitation
analyses
> > > > > Owner: Nobody
> > > > > Requestors: GarciaJ at imsg.com
> > > > > Status: new
> > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > >
> > > > >
> > > > >
> > > > > Dear Met help,
> > > > >
> > > > >
> > > > > I'm trying to run MET precipitation analyses for NMMB, and
get
> > > > > the following error:
> > > > >
> > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > WARNING:
> > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > > for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > > WARNING:
> > > > > WARNING:
> > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > nmbprs_d01.12
> > > > > WARNING:
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
----------------------------------------------------------------
> > > > > --
> > > > > --
> > > > > --
> > > > > ----------
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > > > The namelist has the following:
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A03" ];
> > > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > > GRIB1_ptv = 129;
> > > > >
> > > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > > >
> > > > > Any help will be greatly appreciated.
> > > > >
> > > > >
> > > > > Jose M. Garcia-Rivera
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
this
> > > > > e-mail is not the intended recipient or his or her
authorized
> > > > > agent, the reader is hereby notified that any dissemination,
> > > > > distribution or copying of this e-mail is prohibited. If you
> > > > > have received this e-mail in error, please notify the sender
by
> > > > > replying to this message and delete this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Wed Jan 06 16:10:39 2016
Dear John,
That cleared a lot of doubts. Does the width in the namelist is
measured in gridpoints? i.e. 3 = one gridpoint to the left and one to
the right of the center,,,, or is it in km, or m?
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 5:34 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
Jose,
The neighborhood methods are used primarily to compute Fractions Skill
Score (FSS) which is found in the NBRCNT output line type.
The Heidke Skill Score can be found in the CTS (contingency table
statistics) output line type. It is derived from a 2x2 contingency
table.
If you'd like to look at FSS, I'd suggest:
- Turn on the NBRCNT output line type in the Grid-Stat config file:
nbrcnt = STAT or BOTH;
- Specify several "width" settings in the "nbrhd" section:
nbrhd = {
vld_thresh = 0.1;
width = [ 3, 5, 7, 9, 11, 13 ];
cov_thresh = [ >=0.1 ];
}
People often see how FSS changes as you increase the neighborhood
size.
Hope that helps.
Thanks,
John
On Wed, Jan 6, 2016 at 3:03 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> It worked! Your help is greatly appreciated.
>
> I have a question regarding neighborhood methods:
> We've been changing the parameters, and we keep getting the same HSS
> values. How can we activate neighborhood methods in the namelist?
>
> For example, we changed for Reflectivity to:
>
> // Neighborhood methods
> //
> nbrhd = {
> vld_thresh = 0.1;
> width = [ 7 ];
> cov_thresh = [ >=0.1 ];
> }
>
>
> Best Regards,
> Jose
>
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 3:57 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> For whatever reason that "narr.2015102812.RS.sfc_regridded.grb" file
> encodes APCP using GRIB parameter table version number 131, as shown
> in the wgrib command below:
>
>
>
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:Tim
> eU=1:sfc:0-3hr acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1 0 0
15
> 10 28 12 0 1 0 3 4 0
> 0 0 21 15 0 0
>
> Notice that the 4th value in "PDS10=" list is 131.
>
> I often use the plot_data_plane utility to test out this sort of
thing.
> Try running the following commands:
>
> met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
> 'name="APCP"; level="A03";
>
> met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
> narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
> GRIB1_ptv = 131;'
>
> The resulting images are attached.
>
> In your Grid-Stat config file, rather than saying "obs = fcst" just
> set "obs" to the same thing as the "fcst" but add "GRIB1_ptv =
131;".
>
> Hopefully that does the trick.
>
> John
>
>
>
> On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > The files are now at 'garcia_data_new'
> >
> >
> > Thanks for the help,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 1:42 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > It's probably be easiest to just have you send me sample files.
> > Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb,
and
> > your Grid-Stat configuration file to our anonymous ftp server
> > following these
> > instructions:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Once I hear back from you that the data is posted, I'll go grab it
> > and check for issues.
> >
> > Thanks,
> > John
> >
> > On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > Thanks for the help, now it's working for the NMMB file.
However,
> > > the error now shows up with the NARR data file:
> > >
> > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation File:
> > > narr.2015102812.RS.sfc_regridded.grb
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
> > > for VarInfo "APCP/A03" in GRIB file "narr.2015102812
> .RS.sfc_regridded.grb".
> > > WARNING:
> > > WARNING:
> > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > narr.2015102812 .RS.sfc_regridded.grb
> > >
> > >
> > > Running wgrib for NARR:
> > > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3
> > > :T
> > > im
> > > eU=1:sfc:0-3hr
> > > acc:NAve=0
> > >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2
> > > =3
> > > :T
> > > imeU=1:sfc:0-3hr
> > > acc:NAve=0
> > >
> > > Best Regards,
> > > Jose
> > >
> > >
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 12:51 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > Please just try removing the "GRIB1_ptv = 129;" setting or
change
> > > it to the default PTV value of 2. Then it should work.
> > >
> > > The reason why it's not working now is that your GRIB record is
> > > most likely encoded using parameter table version number 2. For
> > > example, try running this command:
> > >
> > > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> > >
> > > You'll see in the output something like this:
> > > PDS10(1..28)=
> > > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5
8 7
> 0
> > > 0 1 0 3
> > > 21: 4 0 0 0 21 0 0 3
> > >
> > > The PDS is the product description section and the 4th value
> > > listed is probably the number 2, as in the example above.
Looking
> > > at the PDS definition table (
> > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > > the 4th value is the parameter table version number.
> > >
> > > When you set "GRIB1_ptv = 129;", MET will skip over that record
> > > because the PTV value in the data (2) doesn't match the value
you
> > specified (129).
> > >
> > > Hope that helps.
> > >
> > > Thanks,
> > > Joihn
> > >
> > > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
>
> > > >
> > > > Dear John,
> > > >
> > > > Thanks for the reply. I've been playing with the var info
using
> > > > wgrib, but can't make grid_stat work. Here's the info you
want:
> > > >
> > > >
>
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > > TimeU=1:sfc:9-12hr
> > > > acc:NAve=0
> > > >
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > precipitation analyses
> > > >
> > > > Jose,
> > > >
> > > > I see that you're having a hard time extracting 3-hourly
precip
> > > > from an NMMB GRIB file. The question is if and how the
> > > > accumulated precipitation is stored in that GRIB file. Try
> > > > running the following wgrib command to see what APCP records
it
> contains:
> > > >
> > > > wgrib nmbprs_d01.12 | grep APCP
> > > >
> > > > Please send me the wgrib output, and that should tell us more
> > > > about what your GRIB file contains.
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > > met_help at ucar.edu
> > > >
> > > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > Queue: met_help
> > > > > Subject: Problems running grid_stat precipitation
analyses
> > > > > Owner: Nobody
> > > > > Requestors: GarciaJ at imsg.com
> > > > > Status: new
> > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > >
> > > > >
> > > > >
> > > > > Dear Met help,
> > > > >
> > > > >
> > > > > I'm trying to run MET precipitation analyses for NMMB, and
get
> > > > > the following error:
> > > > >
> > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > WARNING:
> > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
> > > > > found for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > > WARNING:
> > > > > WARNING:
> > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > nmbprs_d01.12
> > > > > WARNING:
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
--------------------------------------------------------------
> > > > > --
> > > > > --
> > > > > --
> > > > > --
> > > > > ----------
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > > > The namelist has the following:
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A03" ];
> > > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > > GRIB1_ptv = 129;
> > > > >
> > > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > > >
> > > > > Any help will be greatly appreciated.
> > > > >
> > > > >
> > > > > Jose M. Garcia-Rivera
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
> > > > > this e-mail is not the intended recipient or his or her
> > > > > authorized agent, the reader is hereby notified that any
> > > > > dissemination, distribution or copying of this e-mail is
> > > > > prohibited. If you have received this e-mail in error,
please
> > > > > notify the sender by replying to this message and delete
this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
> > > > person or entity to which it is addressed. If the reader of
this
> > > > e-mail is not the intended recipient or his or her authorized
> > > > agent, the reader is hereby notified that any dissemination,
> > > > distribution or copying of this e-mail is prohibited. If you
> > > > have received this e-mail in error, please notify the sender
by
> > > > replying to this message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Wed Jan 06 16:59:08 2016
The width is in grid squares. A width of 3 means a 3x3 box centered
around
each point. 5 means a 5x5 box, and so on.
In the output you'll actually see that the INTERP_PNTS columns says 9
for a
3x3 box and 25 for a 5x5 box.
John
On Wed, Jan 6, 2016 at 4:10 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> That cleared a lot of doubts. Does the width in the namelist is
measured
> in gridpoints? i.e. 3 = one gridpoint to the left and one to the
right of
> the center,,,, or is it in km, or m?
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 5:34 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> The neighborhood methods are used primarily to compute Fractions
Skill
> Score (FSS) which is found in the NBRCNT output line type.
>
> The Heidke Skill Score can be found in the CTS (contingency table
> statistics) output line type. It is derived from a 2x2 contingency
table.
>
> If you'd like to look at FSS, I'd suggest:
> - Turn on the NBRCNT output line type in the Grid-Stat config file:
> nbrcnt = STAT or BOTH;
> - Specify several "width" settings in the "nbrhd" section:
>
> nbrhd = {
> vld_thresh = 0.1;
> width = [ 3, 5, 7, 9, 11, 13 ];
> cov_thresh = [ >=0.1 ];
> }
>
> People often see how FSS changes as you increase the neighborhood
size.
>
> Hope that helps.
>
> Thanks,
> John
>
>
> On Wed, Jan 6, 2016 at 3:03 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > It worked! Your help is greatly appreciated.
> >
> > I have a question regarding neighborhood methods:
> > We've been changing the parameters, and we keep getting the same
HSS
> > values. How can we activate neighborhood methods in the namelist?
> >
> > For example, we changed for Reflectivity to:
> >
> > // Neighborhood methods
> > //
> > nbrhd = {
> > vld_thresh = 0.1;
> > width = [ 7 ];
> > cov_thresh = [ >=0.1 ];
> > }
> >
> >
> > Best Regards,
> > Jose
> >
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 3:57 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > For whatever reason that "narr.2015102812.RS.sfc_regridded.grb"
file
> > encodes APCP using GRIB parameter table version number 131, as
shown
> > in the wgrib command below:
> >
> >
> >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:Tim
> > eU=1:sfc:0-3hr acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1 0 0
15
> > 10 28 12 0 1 0 3 4 0
> > 0 0 21 15 0 0
> >
> > Notice that the 4th value in "PDS10=" list is 131.
> >
> > I often use the plot_data_plane utility to test out this sort of
thing.
> > Try running the following commands:
> >
> > met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
> > 'name="APCP"; level="A03";
> >
> > met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
> > narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
> > GRIB1_ptv = 131;'
> >
> > The resulting images are attached.
> >
> > In your Grid-Stat config file, rather than saying "obs = fcst"
just
> > set "obs" to the same thing as the "fcst" but add "GRIB1_ptv =
131;".
> >
> > Hopefully that does the trick.
> >
> > John
> >
> >
> >
> > On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > The files are now at 'garcia_data_new'
> > >
> > >
> > > Thanks for the help,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 1:42 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > It's probably be easiest to just have you send me sample files.
> > > Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb,
and
> > > your Grid-Stat configuration file to our anonymous ftp server
> > > following these
> > > instructions:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Once I hear back from you that the data is posted, I'll go grab
it
> > > and check for issues.
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
>
> > > >
> > > > Dear John,
> > > >
> > > > Thanks for the help, now it's working for the NMMB file.
However,
> > > > the error now shows up with the NARR data file:
> > > >
> > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > narr.2015102812.RS.sfc_regridded.grb
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > for VarInfo "APCP/A03" in GRIB file "narr.2015102812
> > .RS.sfc_regridded.grb".
> > > > WARNING:
> > > > WARNING:
> > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > narr.2015102812 .RS.sfc_regridded.grb
> > > >
> > > >
> > > > Running wgrib for NARR:
> > > > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> > > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3
> > > > :T
> > > > im
> > > > eU=1:sfc:0-3hr
> > > > acc:NAve=0
> > > >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:P2
> > > > =3
> > > > :T
> > > > imeU=1:sfc:0-3hr
> > > > acc:NAve=0
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > >
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Wednesday, January 06, 2016 12:51 PM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > precipitation analyses
> > > >
> > > > Jose,
> > > >
> > > > Please just try removing the "GRIB1_ptv = 129;" setting or
change
> > > > it to the default PTV value of 2. Then it should work.
> > > >
> > > > The reason why it's not working now is that your GRIB record
is
> > > > most likely encoded using parameter table version number 2.
For
> > > > example, try running this command:
> > > >
> > > > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> > > >
> > > > You'll see in the output something like this:
> > > > PDS10(1..28)=
> > > > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5
8 7
> > 0
> > > > 0 1 0 3
> > > > 21: 4 0 0 0 21 0 0 3
> > > >
> > > > The PDS is the product description section and the 4th value
> > > > listed is probably the number 2, as in the example above.
Looking
> > > > at the PDS definition table (
> > > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > > > the 4th value is the parameter table version number.
> > > >
> > > > When you set "GRIB1_ptv = 129;", MET will skip over that
record
> > > > because the PTV value in the data (2) doesn't match the value
you
> > > specified (129).
> > > >
> > > > Hope that helps.
> > > >
> > > > Thanks,
> > > > Joihn
> > > >
> > > > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > > > >
> > > > > Dear John,
> > > > >
> > > > > Thanks for the reply. I've been playing with the var info
using
> > > > > wgrib, but can't make grid_stat work. Here's the info you
want:
> > > > >
> > > > >
> >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > > > TimeU=1:sfc:9-12hr
> > > > > acc:NAve=0
> > > > >
> > > > >
> > > > > Best Regards,
> > > > > Jose
> > > > >
> > > > > -----Original Message-----
> > > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > > > To: Jose Garcia-Rivera
> > > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > > precipitation analyses
> > > > >
> > > > > Jose,
> > > > >
> > > > > I see that you're having a hard time extracting 3-hourly
precip
> > > > > from an NMMB GRIB file. The question is if and how the
> > > > > accumulated precipitation is stored in that GRIB file. Try
> > > > > running the following wgrib command to see what APCP records
it
> > contains:
> > > > >
> > > > > wgrib nmbprs_d01.12 | grep APCP
> > > > >
> > > > > Please send me the wgrib output, and that should tell us
more
> > > > > about what your GRIB file contains.
> > > > >
> > > > > Thanks,
> > > > > John Halley Gotway
> > > > > met_help at ucar.edu
> > > > >
> > > > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > >
> > > > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > > Queue: met_help
> > > > > > Subject: Problems running grid_stat precipitation
analyses
> > > > > > Owner: Nobody
> > > > > > Requestors: GarciaJ at imsg.com
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > > >
> > > > > >
> > > > > >
> > > > > > Dear Met help,
> > > > > >
> > > > > >
> > > > > > I'm trying to run MET precipitation analyses for NMMB, and
get
> > > > > > the following error:
> > > > > >
> > > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > > WARNING:
> > > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
> > > > > > found for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > > > WARNING:
> > > > > > WARNING:
> > > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > > nmbprs_d01.12
> > > > > > WARNING:
> > > > > > DEBUG 2:
> > > > > > DEBUG 2:
> > > > > >
--------------------------------------------------------------
> > > > > > --
> > > > > > --
> > > > > > --
> > > > > > --
> > > > > > ----------
> > > > > > DEBUG 2:
> > > > > >
> > > > > >
> > > > > > The namelist has the following:
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "A03" ];
> > > > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > > > GRIB1_ptv = 129;
> > > > > >
> > > > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > > > >
> > > > > > Any help will be greatly appreciated.
> > > > > >
> > > > > >
> > > > > > Jose M. Garcia-Rivera
> > > > > > This e-mail, including attachments, may include
confidential
> > > > > > and/or proprietary information, and may be used only by
the
> > > > > > person or entity to which it is addressed. If the reader
of
> > > > > > this e-mail is not the intended recipient or his or her
> > > > > > authorized agent, the reader is hereby notified that any
> > > > > > dissemination, distribution or copying of this e-mail is
> > > > > > prohibited. If you have received this e-mail in error,
please
> > > > > > notify the sender by replying to this message and delete
this
> > > > > e-mail immediately.
> > > > > >
> > > > > >
> > > > >
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
this
> > > > > e-mail is not the intended recipient or his or her
authorized
> > > > > agent, the reader is hereby notified that any dissemination,
> > > > > distribution or copying of this e-mail is prohibited. If you
> > > > > have received this e-mail in error, please notify the sender
by
> > > > > replying to this message and delete this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74609] Problems running grid_stat precipitation analyses
From: Jose Garcia-Rivera
Time: Thu Jan 07 13:45:00 2016
Dear John,
Thank you so much for all of the information. It was very helpful.
Best Regards
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, January 06, 2016 6:59 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
precipitation analyses
The width is in grid squares. A width of 3 means a 3x3 box centered
around each point. 5 means a 5x5 box, and so on.
In the output you'll actually see that the INTERP_PNTS columns says 9
for a
3x3 box and 25 for a 5x5 box.
John
On Wed, Jan 6, 2016 at 4:10 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> That cleared a lot of doubts. Does the width in the namelist is
> measured in gridpoints? i.e. 3 = one gridpoint to the left and one
to
> the right of the center,,,, or is it in km, or m?
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 5:34 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> Jose,
>
> The neighborhood methods are used primarily to compute Fractions
Skill
> Score (FSS) which is found in the NBRCNT output line type.
>
> The Heidke Skill Score can be found in the CTS (contingency table
> statistics) output line type. It is derived from a 2x2 contingency
table.
>
> If you'd like to look at FSS, I'd suggest:
> - Turn on the NBRCNT output line type in the Grid-Stat config file:
> nbrcnt = STAT or BOTH;
> - Specify several "width" settings in the "nbrhd" section:
>
> nbrhd = {
> vld_thresh = 0.1;
> width = [ 3, 5, 7, 9, 11, 13 ];
> cov_thresh = [ >=0.1 ];
> }
>
> People often see how FSS changes as you increase the neighborhood
size.
>
> Hope that helps.
>
> Thanks,
> John
>
>
> On Wed, Jan 6, 2016 at 3:03 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > It worked! Your help is greatly appreciated.
> >
> > I have a question regarding neighborhood methods:
> > We've been changing the parameters, and we keep getting the same
HSS
> > values. How can we activate neighborhood methods in the namelist?
> >
> > For example, we changed for Reflectivity to:
> >
> > // Neighborhood methods
> > //
> > nbrhd = {
> > vld_thresh = 0.1;
> > width = [ 7 ];
> > cov_thresh = [ >=0.1 ];
> > }
> >
> >
> > Best Regards,
> > Jose
> >
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 3:57 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > For whatever reason that "narr.2015102812.RS.sfc_regridded.grb"
file
> > encodes APCP using GRIB parameter table version number 131, as
shown
> > in the wgrib command below:
> >
> >
> >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> > im eU=1:sfc:0-3hr acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61 1
0
> > 0 15
> > 10 28 12 0 1 0 3 4 0
> > 0 0 21 15 0 0
> >
> > Notice that the 4th value in "PDS10=" list is 131.
> >
> > I often use the plot_data_plane utility to test out this sort of
thing.
> > Try running the following commands:
> >
> > met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
> > 'name="APCP"; level="A03";
> >
> > met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
> > narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
> > GRIB1_ptv = 131;'
> >
> > The resulting images are attached.
> >
> > In your Grid-Stat config file, rather than saying "obs = fcst"
just
> > set "obs" to the same thing as the "fcst" but add "GRIB1_ptv =
131;".
> >
> > Hopefully that does the trick.
> >
> > John
> >
> >
> >
> > On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > The files are now at 'garcia_data_new'
> > >
> > >
> > > Thanks for the help,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 1:42 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > It's probably be easiest to just have you send me sample files.
> > > Please post nmbprs_d01.12, narr.2015102812.RS.sfc_regridded.grb,
> > > and your Grid-Stat configuration file to our anonymous ftp
server
> > > following these
> > > instructions:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Once I hear back from you that the data is posted, I'll go grab
it
> > > and check for issues.
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
>
> > > >
> > > > Dear John,
> > > >
> > > > Thanks for the help, now it's working for the NMMB file.
> > > > However, the error now shows up with the NARR data file:
> > > >
> > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > narr.2015102812.RS.sfc_regridded.grb
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > for VarInfo "APCP/A03" in GRIB file "narr.2015102812
> > .RS.sfc_regridded.grb".
> > > > WARNING:
> > > > WARNING:
> > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > narr.2015102812 .RS.sfc_regridded.grb
> > > >
> > > >
> > > > Running wgrib for NARR:
> > > > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> > > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2
> > > > =3
> > > > :T
> > > > im
> > > > eU=1:sfc:0-3hr
> > > > acc:NAve=0
> > > >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:
> > > > P2
> > > > =3
> > > > :T
> > > > imeU=1:sfc:0-3hr
> > > > acc:NAve=0
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > >
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Wednesday, January 06, 2016 12:51 PM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > precipitation analyses
> > > >
> > > > Jose,
> > > >
> > > > Please just try removing the "GRIB1_ptv = 129;" setting or
> > > > change it to the default PTV value of 2. Then it should work.
> > > >
> > > > The reason why it's not working now is that your GRIB record
is
> > > > most likely encoded using parameter table version number 2.
For
> > > > example, try running this command:
> > > >
> > > > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> > > >
> > > > You'll see in the output something like this:
> > > > PDS10(1..28)=
> > > > 1: 0 0 28 2 59 89 212 192 61 1 0 0 5
8 7
> > 0
> > > > 0 1 0 3
> > > > 21: 4 0 0 0 21 0 0 3
> > > >
> > > > The PDS is the product description section and the 4th value
> > > > listed is probably the number 2, as in the example above.
> > > > Looking at the PDS definition table (
> > > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > > > the 4th value is the parameter table version number.
> > > >
> > > > When you set "GRIB1_ptv = 129;", MET will skip over that
record
> > > > because the PTV value in the data (2) doesn't match the value
> > > > you
> > > specified (129).
> > > >
> > > > Hope that helps.
> > > >
> > > > Thanks,
> > > > Joihn
> > > >
> > > > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > >
> > > > >
> > > > > Dear John,
> > > > >
> > > > > Thanks for the reply. I've been playing with the var info
> > > > > using wgrib, but can't make grid_stat work. Here's the info
you want:
> > > > >
> > > > >
> >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > > > TimeU=1:sfc:9-12hr
> > > > > acc:NAve=0
> > > > >
> > > > >
> > > > > Best Regards,
> > > > > Jose
> > > > >
> > > > > -----Original Message-----
> > > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > > > To: Jose Garcia-Rivera
> > > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
> > > > > grid_stat precipitation analyses
> > > > >
> > > > > Jose,
> > > > >
> > > > > I see that you're having a hard time extracting 3-hourly
> > > > > precip from an NMMB GRIB file. The question is if and how
the
> > > > > accumulated precipitation is stored in that GRIB file. Try
> > > > > running the following wgrib command to see what APCP records
> > > > > it
> > contains:
> > > > >
> > > > > wgrib nmbprs_d01.12 | grep APCP
> > > > >
> > > > > Please send me the wgrib output, and that should tell us
more
> > > > > about what your GRIB file contains.
> > > > >
> > > > > Thanks,
> > > > > John Halley Gotway
> > > > > met_help at ucar.edu
> > > > >
> > > > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > >
> > > > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > > Queue: met_help
> > > > > > Subject: Problems running grid_stat precipitation
analyses
> > > > > > Owner: Nobody
> > > > > > Requestors: GarciaJ at imsg.com
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > > >
> > > > > >
> > > > > >
> > > > > > Dear Met help,
> > > > > >
> > > > > >
> > > > > > I'm trying to run MET precipitation analyses for NMMB, and
> > > > > > get the following error:
> > > > > >
> > > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > > WARNING:
> > > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
> > > > > > found for VarInfo "APCP/A03" in GRIB file "nmbprs_d01.12".
> > > > > > WARNING:
> > > > > > WARNING:
> > > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > > nmbprs_d01.12
> > > > > > WARNING:
> > > > > > DEBUG 2:
> > > > > > DEBUG 2:
> > > > > >
------------------------------------------------------------
> > > > > > --
> > > > > > --
> > > > > > --
> > > > > > --
> > > > > > --
> > > > > > ----------
> > > > > > DEBUG 2:
> > > > > >
> > > > > >
> > > > > > The namelist has the following:
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "A03" ];
> > > > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > > > GRIB1_ptv = 129;
> > > > > >
> > > > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > > > >
> > > > > > Any help will be greatly appreciated.
> > > > > >
> > > > > >
> > > > > > Jose M. Garcia-Rivera
> > > > > > This e-mail, including attachments, may include
confidential
> > > > > > and/or proprietary information, and may be used only by
the
> > > > > > person or entity to which it is addressed. If the reader
of
> > > > > > this e-mail is not the intended recipient or his or her
> > > > > > authorized agent, the reader is hereby notified that any
> > > > > > dissemination, distribution or copying of this e-mail is
> > > > > > prohibited. If you have received this e-mail in error,
> > > > > > please notify the sender by replying to this message and
> > > > > > delete this
> > > > > e-mail immediately.
> > > > > >
> > > > > >
> > > > >
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
> > > > > this e-mail is not the intended recipient or his or her
> > > > > authorized agent, the reader is hereby notified that any
> > > > > dissemination, distribution or copying of this e-mail is
> > > > > prohibited. If you have received this e-mail in error,
please
> > > > > notify the sender by replying to this message and delete
this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
> > > > person or entity to which it is addressed. If the reader of
this
> > > > e-mail is not the intended recipient or his or her authorized
> > > > agent, the reader is hereby notified that any dissemination,
> > > > distribution or copying of this e-mail is prohibited. If you
> > > > have received this e-mail in error, please notify the sender
by
> > > > replying to this message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Problems running grid_stat precipitation analyses
From: John Halley Gotway
Time: Fri Jan 08 09:42:16 2016
Happy to help. I'll go ahead and resolve this ticket now.
Just let us know if more issues or questions arise in your use of MET.
Thanks,
John
On Thu, Jan 7, 2016 at 1:45 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
>
> Dear John,
>
> Thank you so much for all of the information. It was very helpful.
>
> Best Regards
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, January 06, 2016 6:59 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> precipitation analyses
>
> The width is in grid squares. A width of 3 means a 3x3 box centered
> around each point. 5 means a 5x5 box, and so on.
>
> In the output you'll actually see that the INTERP_PNTS columns says
9 for a
> 3x3 box and 25 for a 5x5 box.
>
> John
>
> On Wed, Jan 6, 2016 at 4:10 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> >
> > Dear John,
> >
> > That cleared a lot of doubts. Does the width in the namelist is
> > measured in gridpoints? i.e. 3 = one gridpoint to the left and one
to
> > the right of the center,,,, or is it in km, or m?
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, January 06, 2016 5:34 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > precipitation analyses
> >
> > Jose,
> >
> > The neighborhood methods are used primarily to compute Fractions
Skill
> > Score (FSS) which is found in the NBRCNT output line type.
> >
> > The Heidke Skill Score can be found in the CTS (contingency table
> > statistics) output line type. It is derived from a 2x2
contingency
> table.
> >
> > If you'd like to look at FSS, I'd suggest:
> > - Turn on the NBRCNT output line type in the Grid-Stat config
file:
> > nbrcnt = STAT or BOTH;
> > - Specify several "width" settings in the "nbrhd" section:
> >
> > nbrhd = {
> > vld_thresh = 0.1;
> > width = [ 3, 5, 7, 9, 11, 13 ];
> > cov_thresh = [ >=0.1 ];
> > }
> >
> > People often see how FSS changes as you increase the neighborhood
size.
> >
> > Hope that helps.
> >
> > Thanks,
> > John
> >
> >
> > On Wed, Jan 6, 2016 at 3:03 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > >
> > > Dear John,
> > >
> > > It worked! Your help is greatly appreciated.
> > >
> > > I have a question regarding neighborhood methods:
> > > We've been changing the parameters, and we keep getting the same
HSS
> > > values. How can we activate neighborhood methods in the
namelist?
> > >
> > > For example, we changed for Reflectivity to:
> > >
> > > // Neighborhood methods
> > > //
> > > nbrhd = {
> > > vld_thresh = 0.1;
> > > width = [ 7 ];
> > > cov_thresh = [ >=0.1 ];
> > > }
> > >
> > >
> > > Best Regards,
> > > Jose
> > >
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Wednesday, January 06, 2016 3:57 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running grid_stat
> > > precipitation analyses
> > >
> > > Jose,
> > >
> > > For whatever reason that "narr.2015102812.RS.sfc_regridded.grb"
file
> > > encodes APCP using GRIB parameter table version number 131, as
shown
> > > in the wgrib command below:
> > >
> > >
> > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:T
> > > im eU=1:sfc:0-3hr acc:NAve=0:PDS10= 0 0 28 131 7 140 255 128 61
1 0
> > > 0 15
> > > 10 28 12 0 1 0 3 4 0
> > > 0 0 21 15 0 0
> > >
> > > Notice that the 4th value in "PDS10=" list is 131.
> > >
> > > I often use the plot_data_plane utility to test out this sort of
thing.
> > > Try running the following commands:
> > >
> > > met-5.1/bin/plot_data_plane nmbprs_d01.12 nmbprs_d01.12.ps
> > > 'name="APCP"; level="A03";
> > >
> > > met-5.1/bin/plot_data_plane narr.2015102812.RS.sfc_regridded.grb
> > > narr.2015102812.RS.sfc_regridded.ps 'name="APCP"; level="A03";
> > > GRIB1_ptv = 131;'
> > >
> > > The resulting images are attached.
> > >
> > > In your Grid-Stat config file, rather than saying "obs = fcst"
just
> > > set "obs" to the same thing as the "fcst" but add "GRIB1_ptv =
131;".
> > >
> > > Hopefully that does the trick.
> > >
> > > John
> > >
> > >
> > >
> > > On Wed, Jan 6, 2016 at 12:18 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
>
> > > >
> > > > Dear John,
> > > >
> > > > The files are now at 'garcia_data_new'
> > > >
> > > >
> > > > Thanks for the help,
> > > > Jose
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Wednesday, January 06, 2016 1:42 PM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > precipitation analyses
> > > >
> > > > Jose,
> > > >
> > > > It's probably be easiest to just have you send me sample
files.
> > > > Please post nmbprs_d01.12,
narr.2015102812.RS.sfc_regridded.grb,
> > > > and your Grid-Stat configuration file to our anonymous ftp
server
> > > > following these
> > > > instructions:
> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Once I hear back from you that the data is posted, I'll go
grab it
> > > > and check for issues.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Wed, Jan 6, 2016 at 11:06 AM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609 >
> > > > >
> > > > > Dear John,
> > > > >
> > > > > Thanks for the help, now it's working for the NMMB file.
> > > > > However, the error now shows up with the NARR data file:
> > > > >
> > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1: Observation
File:
> > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > WARNING:
> > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found
> > > > > for VarInfo "APCP/A03" in GRIB file "narr.2015102812
> > > .RS.sfc_regridded.grb".
> > > > > WARNING:
> > > > > WARNING:
> > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > narr.2015102812 .RS.sfc_regridded.grb
> > > > >
> > > > >
> > > > > Running wgrib for NARR:
> > > > > wgrib narr.2015102812.RS.sfc_regridded.grb | grep APCP
> > > > >
24:3248052:d=15102812:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2
> > > > > =3
> > > > > :T
> > > > > im
> > > > > eU=1:sfc:0-3hr
> > > > > acc:NAve=0
> > > > >
25:3348146:d=15102812:APCPN:kpds5=202:kpds6=1:kpds7=0:TR=4:P1=0:
> > > > > P2
> > > > > =3
> > > > > :T
> > > > > imeU=1:sfc:0-3hr
> > > > > acc:NAve=0
> > > > >
> > > > > Best Regards,
> > > > > Jose
> > > > >
> > > > >
> > > > >
> > > > > -----Original Message-----
> > > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > > Sent: Wednesday, January 06, 2016 12:51 PM
> > > > > To: Jose Garcia-Rivera
> > > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
grid_stat
> > > > > precipitation analyses
> > > > >
> > > > > Jose,
> > > > >
> > > > > Please just try removing the "GRIB1_ptv = 129;" setting or
> > > > > change it to the default PTV value of 2. Then it should
work.
> > > > >
> > > > > The reason why it's not working now is that your GRIB record
is
> > > > > most likely encoded using parameter table version number 2.
For
> > > > > example, try running this command:
> > > > >
> > > > > wgrib -d 304 -V -PDS10 nmbprs_d01.12
> > > > >
> > > > > You'll see in the output something like this:
> > > > > PDS10(1..28)=
> > > > > 1: 0 0 28 2 59 89 212 192 61 1 0 0
5 8
> 7
> > > 0
> > > > > 0 1 0 3
> > > > > 21: 4 0 0 0 21 0 0 3
> > > > >
> > > > > The PDS is the product description section and the 4th value
> > > > > listed is probably the number 2, as in the example above.
> > > > > Looking at the PDS definition table (
> > > > > http://www.nco.ncep.noaa.gov/pmb/docs/on388/section1.html)
> > > > > the 4th value is the parameter table version number.
> > > > >
> > > > > When you set "GRIB1_ptv = 129;", MET will skip over that
record
> > > > > because the PTV value in the data (2) doesn't match the
value
> > > > > you
> > > > specified (129).
> > > > >
> > > > > Hope that helps.
> > > > >
> > > > > Thanks,
> > > > > Joihn
> > > > >
> > > > > On Wed, Jan 6, 2016 at 9:52 AM, Jose Garcia-Rivera via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > > >
> > > > > >
> > > > > > Dear John,
> > > > > >
> > > > > > Thanks for the reply. I've been playing with the var info
> > > > > > using wgrib, but can't make grid_stat work. Here's the
info you
> want:
> > > > > >
> > > > > >
> > >
304:77038000:d=15102800:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:
> > > > > > TimeU=1:sfc:9-12hr
> > > > > > acc:NAve=0
> > > > > >
> > > > > >
> > > > > > Best Regards,
> > > > > > Jose
> > > > > >
> > > > > > -----Original Message-----
> > > > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > > > Sent: Wednesday, January 06, 2016 11:50 AM
> > > > > > To: Jose Garcia-Rivera
> > > > > > Subject: Re: [rt.rap.ucar.edu #74609] Problems running
> > > > > > grid_stat precipitation analyses
> > > > > >
> > > > > > Jose,
> > > > > >
> > > > > > I see that you're having a hard time extracting 3-hourly
> > > > > > precip from an NMMB GRIB file. The question is if and how
the
> > > > > > accumulated precipitation is stored in that GRIB file.
Try
> > > > > > running the following wgrib command to see what APCP
records
> > > > > > it
> > > contains:
> > > > > >
> > > > > > wgrib nmbprs_d01.12 | grep APCP
> > > > > >
> > > > > > Please send me the wgrib output, and that should tell us
more
> > > > > > about what your GRIB file contains.
> > > > > >
> > > > > > Thanks,
> > > > > > John Halley Gotway
> > > > > > met_help at ucar.edu
> > > > > >
> > > > > > On Tue, Jan 5, 2016 at 4:28 PM, Jose Garcia-Rivera via RT
<
> > > > > > met_help at ucar.edu
> > > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Tue Jan 05 16:28:24 2016: Request 74609 was acted upon.
> > > > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > > > Queue: met_help
> > > > > > > Subject: Problems running grid_stat precipitation
analyses
> > > > > > > Owner: Nobody
> > > > > > > Requestors: GarciaJ at imsg.com
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74609
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Dear Met help,
> > > > > > >
> > > > > > >
> > > > > > > I'm trying to run MET precipitation analyses for NMMB,
and
> > > > > > > get the following error:
> > > > > > >
> > > > > > > DEBUG 1: Forecast File: nmbprs_d01.12 DEBUG 1:
Observation
> File:
> > > > > > > narr.2015102812.RS.sfc_regridded.grb
> > > > > > > WARNING:
> > > > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact
match
> > > > > > > found for VarInfo "APCP/A03" in GRIB file
"nmbprs_d01.12".
> > > > > > > WARNING:
> > > > > > > WARNING:
> > > > > > > WARNING: process_scores() -> APCP/A03 not found in file:
> > > > > > > nmbprs_d01.12
> > > > > > > WARNING:
> > > > > > > DEBUG 2:
> > > > > > > DEBUG 2:
> > > > > > >
------------------------------------------------------------
> > > > > > > --
> > > > > > > --
> > > > > > > --
> > > > > > > --
> > > > > > > --
> > > > > > > ----------
> > > > > > > DEBUG 2:
> > > > > > >
> > > > > > >
> > > > > > > The namelist has the following:
> > > > > > >
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "APCP";
> > > > > > > level = [ "A03" ];
> > > > > > > cat_thresh = [ >=0,>=12, >=24 ];
> > > > > > > GRIB1_ptv = 129;
> > > > > > >
> > > > > > > I also tried "L0", "Z0" , "A12" and "A06"
> > > > > > >
> > > > > > > Any help will be greatly appreciated.
> > > > > > >
> > > > > > >
> > > > > > > Jose M. Garcia-Rivera
> > > > > > > This e-mail, including attachments, may include
confidential
> > > > > > > and/or proprietary information, and may be used only by
the
> > > > > > > person or entity to which it is addressed. If the reader
of
> > > > > > > this e-mail is not the intended recipient or his or her
> > > > > > > authorized agent, the reader is hereby notified that any
> > > > > > > dissemination, distribution or copying of this e-mail is
> > > > > > > prohibited. If you have received this e-mail in error,
> > > > > > > please notify the sender by replying to this message and
> > > > > > > delete this
> > > > > > e-mail immediately.
> > > > > > >
> > > > > > >
> > > > > >
> > > > > > This e-mail, including attachments, may include
confidential
> > > > > > and/or proprietary information, and may be used only by
the
> > > > > > person or entity to which it is addressed. If the reader
of
> > > > > > this e-mail is not the intended recipient or his or her
> > > > > > authorized agent, the reader is hereby notified that any
> > > > > > dissemination, distribution or copying of this e-mail is
> > > > > > prohibited. If you have received this e-mail in error,
please
> > > > > > notify the sender by replying to this message and delete
this
> > > > > e-mail immediately.
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
this
> > > > > e-mail is not the intended recipient or his or her
authorized
> > > > > agent, the reader is hereby notified that any dissemination,
> > > > > distribution or copying of this e-mail is prohibited. If you
> > > > > have received this e-mail in error, please notify the sender
by
> > > > > replying to this message and delete this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
More information about the Met_help
mailing list