[Met_help] [rt.rap.ucar.edu #60895] History for Point Stat ANNN level

John Halley Gotway via RT met_help at ucar.edu
Thu May 23 14:36:12 MDT 2013


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

I was wondering what units of measure the point-stat level for accumulation
interval is?  Not sure if it was the same as ascii2nc.
Thanks,

-- 
Ronald David Leeper
Research Associate, CICS-NC
Quality Assurance Specialist, USCRN
National Climatic Data Center (NCDC)
151 Patton Ave.
Asheville, NC 28801-5001
ronald.leeper at noaa.gov
leeperd at gmail.com
Office:(828)257-3185
Fax:(828)271-4022


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

Subject: Re: [rt.rap.ucar.edu #60895] Point Stat ANNN level
From: John Halley Gotway
Time: Fri Apr 05 20:10:46 2013

Ronald,

Throughout MET, accumulation intervals should be specified in HH[MMSS]
format.  If that's not the case somewhere, please let me know.

We do specify time-offsets in seconds though.  For example, in the
"obs_window" setting "beg" and "end" are the offset in seconds
relative to
the valid time of the forecast:

//
// Point observation time window
//
obs_window = {
   beg = -5400;
   end =  5400;
}

Thanks,
John

>
> Fri Apr 05 09:38:01 2013: Request 60895 was acted upon.
> Transaction: Ticket created by ronald.leeper at noaa.gov
>        Queue: met_help
>      Subject: Point Stat ANNN level
>        Owner: Nobody
>   Requestors: ronald.leeper at noaa.gov
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>
>
> I was wondering what units of measure the point-stat level for
> accumulation
> interval is?  Not sure if it was the same as ascii2nc.
> Thanks,
>
> --
> Ronald David Leeper
> Research Associate, CICS-NC
> Quality Assurance Specialist, USCRN
> National Climatic Data Center (NCDC)
> 151 Patton Ave.
> Asheville, NC 28801-5001
> ronald.leeper at noaa.gov
> leeperd at gmail.com
> Office:(828)257-3185
> Fax:(828)271-4022
>



------------------------------------------------
Subject: Point Stat ANNN level
From: Ronald Leeper - NOAA Affiliate
Time: Sat Apr 06 12:42:50 2013

When I specify the 30 minute accumulation interval for level in
pointStat-config file, point-stat processes the first leadtime, but
fails
on the second, indicating no matches within the gribFile. I will load
the
Aug2008_allGribFiles shortly.

*DEBUG 1: Forecast File: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
*DEBUG 1: Climatology File: none*
*DEBUG 1: Observation File: ./KYmesoAug08_BWG.nc*
*DEBUG 2: *
*DEBUG 2:
--------------------------------------------------------------------------------
*
*DEBUG 2: *
*DEBUG 2: Reading data for APCP/A3000.*
*DEBUG 2: For APCP/A3000 found 1 forecast levels and 0 climatology
levels.*
*DEBUG 2: *
*DEBUG 2:
--------------------------------------------------------------------------------
*
*DEBUG 2: *
*DEBUG 2: Reading data for APCP/A3000.*
*ERROR  : *
*ERROR  : process_fcst_climo_files() -> no fields matching APCP/A3000
found
in file: ./Aug2008NARR6_4km/Aug2008_allGribFiles*


On Fri, Apr 5, 2013 at 10:10 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Ronald,
>
> Throughout MET, accumulation intervals should be specified in
HH[MMSS]
> format.  If that's not the case somewhere, please let me know.
>
> We do specify time-offsets in seconds though.  For example, in the
> "obs_window" setting "beg" and "end" are the offset in seconds
relative to
> the valid time of the forecast:
>
> //
> // Point observation time window
> //
> obs_window = {
>    beg = -5400;
>    end =  5400;
> }
>
> Thanks,
> John
>
> >
> > Fri Apr 05 09:38:01 2013: Request 60895 was acted upon.
> > Transaction: Ticket created by ronald.leeper at noaa.gov
> >        Queue: met_help
> >      Subject: Point Stat ANNN level
> >        Owner: Nobody
> >   Requestors: ronald.leeper at noaa.gov
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
> >
> >
> > I was wondering what units of measure the point-stat level for
> > accumulation
> > interval is?  Not sure if it was the same as ascii2nc.
> > Thanks,
> >
> > --
> > Ronald David Leeper
> > Research Associate, CICS-NC
> > Quality Assurance Specialist, USCRN
> > National Climatic Data Center (NCDC)
> > 151 Patton Ave.
> > Asheville, NC 28801-5001
> > ronald.leeper at noaa.gov
> > leeperd at gmail.com
> > Office:(828)257-3185
> > Fax:(828)271-4022
> >
>
>
>
>


--
Ronald David Leeper
Research Associate, CICS-NC
Quality Assurance Specialist, USCRN
National Climatic Data Center (NCDC)
151 Patton Ave.
Asheville, NC 28801-5001
ronald.leeper at noaa.gov
leeperd at gmail.com
Office:(828)257-3185
Fax:(828)271-4022

------------------------------------------------
Subject: Point Stat ANNN level
From: John Halley Gotway
Time: Mon Apr 08 10:11:06 2013

Ronald,

Thanks for sending the data.  It makes it a lot easier to debug.
Point-Stat can't find any data for the second verification task:
    {name="APCP"; level=[ "A003000" ]; lead_time="010000";}

You've requested an accumulation interval of 30 minutes for a lead
time of 1 hour.  So you're requesting the APCP accumulation from 30 -
60 minutes.  Point-Stat can't find any data in the file because
that data does not exist in the GRIB file you sent me.  I ran the
following wgrib command to see what APCP records the GRIB file you
sent contains:
    wgrib Aug2008_allGribFiles | grep APCP > APCP_wgrib.txt

Looking through the inventory, you can see there are no records for
the 30 - 60 minute accumulation.  There is a 0-1hr (record #13829) and
a 0-30min (record #16952), but no 30-60min.

If you need that 30-60min accumulation interval, you could use the
pcp_combine tool to compute it:
    METv4.0/bin/pcp_combine -subtract \
       Aug2008_allGribFiles 'name="APCP"; level="A010000";
lead_time="010000";' \
       Aug2008_allGribFiles 'name="APCP"; level="A003000";
lead_time="003000";' \
       apcp_30_60.nc -v 4

That results in a NetCDF file that you could then verify using Point-
Stat.  But it'd require a separate call to Point-Stat since it only
handle one input forecast file at a time.

All of these combinations of lead times and accumulation intervals get
tricky to keep straight!

Hope that helps.

Thanks,
John

On 04/06/2013 12:42 PM, Ronald Leeper - NOAA Affiliate via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>
> When I specify the 30 minute accumulation interval for level in
> pointStat-config file, point-stat processes the first leadtime, but
fails
> on the second, indicating no matches within the gribFile. I will
load the
> Aug2008_allGribFiles shortly.
>
> *DEBUG 1: Forecast File: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
> *DEBUG 1: Climatology File: none*
> *DEBUG 1: Observation File: ./KYmesoAug08_BWG.nc*
> *DEBUG 2: *
> *DEBUG 2:
>
--------------------------------------------------------------------------------
> *
> *DEBUG 2: *
> *DEBUG 2: Reading data for APCP/A3000.*
> *DEBUG 2: For APCP/A3000 found 1 forecast levels and 0 climatology
levels.*
> *DEBUG 2: *52
> *DEBUG 2:
>
--------------------------------------------------------------------------------
> *
> *DEBUG 2: *
> *DEBUG 2: Reading data for APCP/A3000.*
> *ERROR  : *
> *ERROR  : process_fcst_climo_files() -> no fields matching
APCP/A3000 found
> in file: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
>
>
> On Fri, Apr 5, 2013 at 10:10 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Ronald,
>>
>> Throughout MET, accumulation intervals should be specified in
HH[MMSS]
>> format.  If that's not the case somewhere, please let me know.
>>R
>> We do specify time-offsets in seconds though.  For example, in the
>> "obs_window" setting "beg" and "end" are the offset in seconds
relative to
>> the valid time of the forecast:
>>
>> //
>> // Point observation time window
>> //
>> obs_window = {
>>     beg = -5400;
>>     end =  5400;
>> }
>>
>> Thanks,
>> John
>>
>>>
>>> Fri Apr 05 09:38:01 2013: Request 60895 was acted upon.
>>> Transaction: Ticket created by ronald.leeper at noaa.gov
>>>         Queue: met_help
>>>       Subject: Point Stat ANNN level
>>>         Owner: Nobody
>>>    Requestors: ronald.leeper at noaa.gov
>>>        Status: new
>>>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>>>
>>>
>>> I was wondering what units of measure the point-stat level for
>>> accumulation
>>> interval is?  Not sure if it was the same as ascii2nc.
>>> Thanks,
>>>
>>> --
>>> Ronald David Leeper
>>> Research Associate, CICS-NC
>>> Quality Assurance Specialist, USCRN
>>> National Climatic Data Center (NCDC)
>>> 151 Patton Ave.
>>> Asheville, NC 28801-5001
>>> ronald.leeper at noaa.gov
>>> leeperd at gmail.com
>>> Office:(828)257-3185
>>> Fax:(828)271-4022
>>>
>>
>>
>>
>>
>
>

------------------------------------------------
Subject: Point Stat ANNN level
From: John Halley Gotway
Time: Mon Apr 08 10:11:06 2013

278:50110968:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=18:P2=19:TimeU=14:sfc:540-
570min acc:NAve=0
625:109702116:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=36:P2=37:TimeU=14:sfc:1080-
1110min acc:NAve=0
972:169042722:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr acc:NAve=0
1319:227941402:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=40:P2=41:TimeU=14:sfc:1200-
1230min acc:NAve=0
1684:287065214:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr acc:NAve=0
2031:346450380:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr acc:NAve=0
2378:405143250:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=46:P2=47:TimeU=14:sfc:1380-
1410min acc:NAve=0
2725:464218430:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr acc:NAve=0
3072:523456392:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr acc:NAve=0
3419:582195176:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr acc:NAve=0
3766:641683334:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr acc:NAve=0
4113:700611476:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=44:P2=45:TimeU=14:sfc:1320-
1350min acc:NAve=0
4460:759819290:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr acc:NAve=0
4807:819044736:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=16:P2=17:TimeU=14:sfc:480-
510min acc:NAve=0
5154:878503248:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=38:P2=39:TimeU=14:sfc:1140-
1170min acc:NAve=0
5501:937181038:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr acc:NAve=0
5848:996079362:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=42:P2=43:TimeU=14:sfc:1260-
1290min acc:NAve=0
6195:1054848648:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=20:P2=21:TimeU=14:sfc:600-
630min acc:NAve=0
6542:1114206724:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr acc:NAve=0
6889:1173211638:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr acc:NAve=0
7236:1232184166:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=2:P2=3:TimeU=14:sfc:60-
90min acc:NAve=0
7583:1291115056:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=10:P2=11:TimeU=14:sfc:300-
330min acc:NAve=0
7930:1350517006:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr acc:NAve=0
8277:1409827634:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr acc:NAve=0
8624:1469433330:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=28:P2=29:TimeU=14:sfc:840-
870min acc:NAve=0
8971:1528731600:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=26:P2=27:TimeU=14:sfc:780-
810min acc:NAve=0
9318:1588060712:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr acc:NAve=0
9665:1647074030:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr acc:NAve=0
10012:1701436486:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=4:P2=5:TimeU=14:sfc:120-
150min acc:NAve=0
10359:1760426152:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=12:P2=13:TimeU=14:sfc:360-
390min acc:NAve=0
10706:1819621452:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr acc:NAve=0
11053:1878813842:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr acc:NAve=0
11400:1938419382:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=34:P2=35:TimeU=14:sfc:1020-
1050min acc:NAve=0
11747:1997865886:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=14:P2=15:TimeU=14:sfc:420-
450min acc:NAve=0
12094:2057193482:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr acc:NAve=0
12441:2116505810:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr acc:NAve=0
12788:2176243976:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=32:P2=33:TimeU=14:sfc:960-
990min acc:NAve=0
13135:2235675398:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=24:P2=25:TimeU=14:sfc:720-
750min acc:NAve=0
13482:2294974866:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr acc:NAve=0
13829:2354652886:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr acc:NAve=0
14176:2413494666:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=6:P2=7:TimeU=14:sfc:180-
210min acc:NAve=0
14523:2472190114:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=8:P2=9:TimeU=14:sfc:240-
270min acc:NAve=0
14870:2531120656:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr acc:NAve=0
15217:2590711462:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr acc:NAve=0
15564:2650346122:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=30:P2=31:TimeU=14:sfc:900-
930min acc:NAve=0
15911:2709777878:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=22:P2=23:TimeU=14:sfc:660-
690min acc:NAve=0
16258:2769164746:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr acc:NAve=0
16605:2828565666:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr acc:NAve=0
16952:2887859328:d=08082406:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=14:sfc:0-
30min acc:NAve=0

------------------------------------------------
Subject: Point Stat ANNN level
From: Ronald Leeper - NOAA Affiliate
Time: Mon Apr 08 15:06:34 2013

I am checking into this currently.  We don't have wgrib on our systems
yet.
 Regardless, this makes sense.  So I would need to have the level set
to
"A010000" and then so on?

Thanks,

On Mon, Apr 8, 2013 at 12:11 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Ronald,
>
> Thanks for sending the data.  It makes it a lot easier to debug.
>  Point-Stat can't find any data for the second verification task:
>     {name="APCP"; level=[ "A003000" ]; lead_time="010000";}
>
> You've requested an accumulation interval of 30 minutes for a lead
time of
> 1 hour.  So you're requesting the APCP accumulation from 30 - 60
minutes.
>  Point-Stat can't find any data in the file because
> that data does not exist in the GRIB file you sent me.  I ran the
> following wgrib command to see what APCP records the GRIB file you
sent
> contains:
>     wgrib Aug2008_allGribFiles | grep APCP > APCP_wgrib.txt
>
> Looking through the inventory, you can see there are no records for
the 30
> - 60 minute accumulation.  There is a 0-1hr (record #13829) and a 0-
30min
> (record #16952), but no 30-60min.
>
> If you need that 30-60min accumulation interval, you could use the
> pcp_combine tool to compute it:
>     METv4.0/bin/pcp_combine -subtract \
>        Aug2008_allGribFiles 'name="APCP"; level="A010000";
> lead_time="010000";' \
>        Aug2008_allGribFiles 'name="APCP"; level="A003000";
> lead_time="003000";' \
>        apcp_30_60.nc -v 4
>
> That results in a NetCDF file that you could then verify using
Point-Stat.
>  But it'd require a separate call to Point-Stat since it only handle
one
> input forecast file at a time.
>
> All of these combinations of lead times and accumulation intervals
get
> tricky to keep straight!
>
> Hope that helps.
>
> Thanks,
> John
>
> On 04/06/2013 12:42 PM, Ronald Leeper - NOAA Affiliate via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
> >
> > When I specify the 30 minute accumulation interval for level in
> > pointStat-config file, point-stat processes the first leadtime,
but fails
> > on the second, indicating no matches within the gribFile. I will
load the
> > Aug2008_allGribFiles shortly.
> >
> > *DEBUG 1: Forecast File: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
> > *DEBUG 1: Climatology File: none*
> > *DEBUG 1: Observation File: ./KYmesoAug08_BWG.nc*
> > *DEBUG 2: *
> > *DEBUG 2:
> >
>
--------------------------------------------------------------------------------
> > *
> > *DEBUG 2: *
> > *DEBUG 2: Reading data for APCP/A3000.*
> > *DEBUG 2: For APCP/A3000 found 1 forecast levels and 0 climatology
> levels.*
> > *DEBUG 2: *52
> > *DEBUG 2:
> >
>
--------------------------------------------------------------------------------
> > *
> > *DEBUG 2: *
> > *DEBUG 2: Reading data for APCP/A3000.*
> > *ERROR  : *
> > *ERROR  : process_fcst_climo_files() -> no fields matching
APCP/A3000
> found
> > in file: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
> >
> >
> > On Fri, Apr 5, 2013 at 10:10 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Ronald,
> >>
> >> Throughout MET, accumulation intervals should be specified in
HH[MMSS]
> >> format.  If that's not the case somewhere, please let me know.
> >>R
> >> We do specify time-offsets in seconds though.  For example, in
the
> >> "obs_window" setting "beg" and "end" are the offset in seconds
relative
> to
> >> the valid time of the forecast:
> >>
> >> //
> >> // Point observation time window
> >> //
> >> obs_window = {
> >>     beg = -5400;
> >>     end =  5400;
> >> }
> >>
> >> Thanks,
> >> John
> >>
> >>>
> >>> Fri Apr 05 09:38:01 2013: Request 60895 was acted upon.
> >>> Transaction: Ticket created by ronald.leeper at noaa.gov
> >>>         Queue: met_help
> >>>       Subject: Point Stat ANNN level
> >>>         Owner: Nobody
> >>>    Requestors: ronald.leeper at noaa.gov
> >>>        Status: new
> >>>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895>
> >>>
> >>>
> >>> I was wondering what units of measure the point-stat level for
> >>> accumulation
> >>> interval is?  Not sure if it was the same as ascii2nc.
> >>> Thanks,
> >>>
> >>> --
> >>> Ronald David Leeper
> >>> Research Associate, CICS-NC
> >>> Quality Assurance Specialist, USCRN
> >>> National Climatic Data Center (NCDC)
> >>> 151 Patton Ave.
> >>> Asheville, NC 28801-5001
> >>> ronald.leeper at noaa.gov
> >>> leeperd at gmail.com
> >>> Office:(828)257-3185
> >>> Fax:(828)271-4022
> >>>
> >>
> >>
> >>
> >>
> >
> >
>
>


--
Ronald David Leeper
Research Associate, CICS-NC
Quality Assurance Specialist, USCRN
National Climatic Data Center (NCDC)
151 Patton Ave.
Asheville, NC 28801-5001
ronald.leeper at noaa.gov
leeperd at gmail.com
Office:(828)257-3185
Fax:(828)271-4022

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #60895] Point Stat ANNN level
From: John Halley Gotway
Time: Mon Apr 08 15:30:37 2013

Ronald,

wgrib is a very useful utility for dumping information about GRIB
files.  I'd highly recommend getting/using it.

Yes, your GRIB data does include a record for accumulated precip
(APCP) from 0 - 1 hours.  So you could use Point-Stat to verify that
field.  Keep in mind though that if you'd like to verify a
forecast of 1-hour accumulation, you'll need observations with that
same accumulation interval.  Otherwise, Point-Stat will find 0 matched
pairs for that task.

In my experience, the observation data is usually the limiting factor.
We often adjust the forecast accumulation intervals to match the
observation data that's available.

John

On 04/08/2013 03:06 PM, Ronald Leeper - NOAA Affiliate via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>
> I am checking into this currently.  We don't have wgrib on our
systems yet.
>   Regardless, this makes sense.  So I would need to have the level
set to
> "A010000" and then so on?
>
> Thanks,
>
> On Mon, Apr 8, 2013 at 12:11 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Ronald,
>>
>> Thanks for sending the data.  It makes it a lot easier to debug.
>>   Point-Stat can't find any data for the second verification task:
>>      {name="APCP"; level=[ "A003000" ]; lead_time="010000";}
>>
>> You've requested an accumulation interval of 30 minutes for a lead
time of
>> 1 hour.  So you're requesting the APCP accumulation from 30 - 60
minutes.
>>   Point-Stat can't find any data in the file because
>> that data does not exist in the GRIB file you sent me.  I ran the
>> following wgrib command to see what APCP records the GRIB file you
sent
>> contains:
>>      wgrib Aug2008_allGribFiles | grep APCP > APCP_wgrib.txt
>>
>> Looking through the inventory, you can see there are no records for
the 30
>> - 60 minute accumulation.  There is a 0-1hr (record #13829) and a
0-30min
>> (record #16952), but no 30-60min.
>>
>> If you need that 30-60min accumulation interval, you could use the
>> pcp_combine tool to compute it:
>>      METv4.0/bin/pcp_combine -subtract \
>>         Aug2008_allGribFiles 'name="APCP"; level="A010000";
>> lead_time="010000";' \
>>         Aug2008_allGribFiles 'name="APCP"; level="A003000";
>> lead_time="003000";' \
>>         apcp_30_60.nc -v 4
>>
>> That results in a NetCDF file that you could then verify using
Point-Stat.
>>   But it'd require a separate call to Point-Stat since it only
handle one
>> input forecast file at a time.
>>
>> All of these combinations of lead times and accumulation intervals
get
>> tricky to keep straight!
>>
>> Hope that helps.
>>
>> Thanks,
>> John
>>
>> On 04/06/2013 12:42 PM, Ronald Leeper - NOAA Affiliate via RT
wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>>>
>>> When I specify the 30 minute accumulation interval for level in
>>> pointStat-config file, point-stat processes the first leadtime,
but fails
>>> on the second, indicating no matches within the gribFile. I will
load the
>>> Aug2008_allGribFiles shortly.
>>>
>>> *DEBUG 1: Forecast File: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
>>> *DEBUG 1: Climatology File: none*
>>> *DEBUG 1: Observation File: ./KYmesoAug08_BWG.nc*
>>> *DEBUG 2: *
>>> *DEBUG 2:
>>>
>>
--------------------------------------------------------------------------------
>>> *
>>> *DEBUG 2: *
>>> *DEBUG 2: Reading data for APCP/A3000.*
>>> *DEBUG 2: For APCP/A3000 found 1 forecast levels and 0 climatology
>> levels.*
>>> *DEBUG 2: *52
>>> *DEBUG 2:
>>>
>>
--------------------------------------------------------------------------------
>>> *
>>> *DEBUG 2: *
>>> *DEBUG 2: Reading data for APCP/A3000.*
>>> *ERROR  : *
>>> *ERROR  : process_fcst_climo_files() -> no fields matching
APCP/A3000
>> found
>>> in file: ./Aug2008NARR6_4km/Aug2008_allGribFiles*
>>>
>>>
>>> On Fri, Apr 5, 2013 at 10:10 PM, John Halley Gotway via RT <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Ronald,
>>>>
>>>> Throughout MET, accumulation intervals should be specified in
HH[MMSS]
>>>> format.  If that's not the case somewhere, please let me know.
>>>> R
>>>> We do specify time-offsets in seconds though.  For example, in
the
>>>> "obs_window" setting "beg" and "end" are the offset in seconds
relative
>> to
>>>> the valid time of the forecast:
>>>>
>>>> //
>>>> // Point observation time window
>>>> //
>>>> obs_window = {
>>>>      beg = -5400;
>>>>      end =  5400;
>>>> }
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>>>
>>>>> Fri Apr 05 09:38:01 2013: Request 60895 was acted upon.
>>>>> Transaction: Ticket created by ronald.leeper at noaa.gov
>>>>>          Queue: met_help
>>>>>        Subject: Point Stat ANNN level
>>>>>          Owner: Nobody
>>>>>     Requestors: ronald.leeper at noaa.gov
>>>>>         Status: new
>>>>>    Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895>
>>>>>
>>>>>
>>>>> I was wondering what units of measure the point-stat level for
>>>>> accumulation
>>>>> interval is?  Not sure if it was the same as ascii2nc.
>>>>> Thanks,
>>>>>
>>>>> --
>>>>> Ronald David Leeper
>>>>> Research Associate, CICS-NC
>>>>> Quality Assurance Specialist, USCRN
>>>>> National Climatic Data Center (NCDC)
>>>>> 151 Patton Ave.
>>>>> Asheville, NC 28801-5001
>>>>> ronald.leeper at noaa.gov
>>>>> leeperd at gmail.com
>>>>> Office:(828)257-3185
>>>>> Fax:(828)271-4022
>>>>>
>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
>

------------------------------------------------
Subject: Point Stat ANNN level
From: Ronald Leeper - NOAA Affiliate
Time: Mon Apr 08 15:44:52 2013

Hopefully system Admins will find the installation of wgrib fairly
straightforward.  When I attempted to run point_stat with the
accumulation
level set to "A013000" it crased.  Was that period also missing?  I
guess I
should really just wait until wgrib is installed.

Thanks again for your help,


On Mon, Apr 8, 2013 at 5:30 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:

> n my experience, the observation data is usually the limiting
factor.  We
> often adjust the forecast accumulation intervals to match the
observation
> data that's available.





--
Ronald David Leeper
Research Associate, CICS-NC
Quality Assurance Specialist, USCRN
National Climatic Data Center (NCDC)
151 Patton Ave.
Asheville, NC 28801-5001
ronald.leeper at noaa.gov
leeperd at gmail.com
Office:(828)257-3185
Fax:(828)271-4022

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #60895] Point Stat ANNN level
From: John Halley Gotway
Time: Mon Apr 08 16:02:29 2013

Ronald,

No, I do not see a 90 minute accumulation in your file.  Yes, I think
wgrib will help you understand the data.

Thanks
John

On 04/08/2013 03:44 PM, Ronald Leeper - NOAA Affiliate via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=60895 >
>
> Hopefully system Admins will find the installation of wgrib fairly
> straightforward.  When I attempted to run point_stat with the
accumulation
> level set to "A013000" it crased.  Was that period also missing?  I
guess I
> should really just wait until wgrib is installed.
>
> Thanks again for your help,
>
>
> On Mon, Apr 8, 2013 at 5:30 PM, John Halley Gotway via RT
<met_help at ucar.edu
>> wrote:
>
>> n my experience, the observation data is usually the limiting
factor.  We
>> often adjust the forecast accumulation intervals to match the
observation
>> data that's available.
>
>
>
>
>

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


More information about the Met_help mailing list