[Met_help] [rt.rap.ucar.edu #64170] History for a question for pcp_combin

John Halley Gotway via RT met_help at ucar.edu
Mon Nov 11 14:26:04 MST 2013


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

hi!

I am runing met4.0 and use pcp_combine to subtract 192hr accum precipation and 168  accum precipation  to get 24hr precipitation.

but when I run it as : 

pcp_combine -subtract pgbf192.grapes.2011070612 192 pgbf168.grapes.2011070612 168 zb.nc 

there are always say that 

WARNING: 
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for VarInfo "APCP/A192" in GRIB file "pgbf192.grapes.2011070612".
WARNING: 
ERROR  : 
ERROR  : get_field() -> can't get data plane from file "pgbf192.grapes.2011070612"
ERROR  : 

the grib index for 192hr is 

352:2522142:d=11070612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=192:TimeU=1:sfc:192hr fcst:NAve=0

what wrong have I met?

thanks alot 

zhaobin

2013-11-06 



zhaob 


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

Subject: Re: [rt.rap.ucar.edu #64170] a question for pcp_combin
From: John Halley Gotway
Time: Thu Nov 07 11:14:03 2013

Zhaob,

Sorry for the delay in getting back to you.  I think I know what the
problem is.  Typically, accumulated precipitation is stored in GRIB
files using an accumulation interval with a "time range"
indicator value of 4.  Here's a description of the different time
range indicator values and meanings:
    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table5.html

For example, take a look at the the APCP in the GRIB files included in
the MET tar ball:
   wgrib METv4.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
| grep APCP
     1:0:d=05080700:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr acc:NAve=0
     2:31408:d=05080700:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:TimeU=1:sfc:9-
12hr acc:NAve=0

The "TR=4" indicates that these records contain an accumulation
between times P1 and P2.  In the first record, the precip is
accumulated between 0 and 12 hours.  In the second record, the precip
is
accumulated between 9 and 12 hours.

However, your GRIB data uses a time range indicator of 10, not 4.
pcp_combine is looking in "pgbf192.grapes.2011070612" for a 192 hour
*accumulation*.  But since the time range indicator is not 4, it
doesn't find a match.  Perhaps you know that the this data really is a
runtime accumulation from 0 to 192 hours, but according to the GRIB
flags, it's a 192 hour instantaneous forecast, not an
accumulation.

If you have control over the generation of these GRIB files, you could
switch the time range indicator to 4.  If not, there's another
workaround.  Instead of telling pcp_combine to look for a
particular accumulation interval, you can give it a more complete
description of the field you'd like to use from each file.  Please try
running the following:

   pcp_combine -subtract \
     pgbf192.grapes.2011070612 'name="APCP"; level="Z0";
lead_time="192";' \
     pgbf168.grapes.2011070612 'name="APCP"; level="Z0";
lead_time="168";' \
     zb.nc

Hope that helps.

Thanks,
John Halley Gotway
met_help at ucar.edu

On 11/06/2013 01:22 AM, zhaob via RT wrote:
>
> Wed Nov 06 01:22:23 2013: Request 64170 was acted upon.
> Transaction: Ticket created by zhaob at cma.gov.cn
>         Queue: met_help
>       Subject: a question for pcp_combin
>         Owner: Nobody
>    Requestors: zhaob at cma.gov.cn
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=64170 >
>
>
> hi!
>
> I am runing met4.0 and use pcp_combine to subtract 192hr accum
precipation and 168  accum precipation  to get 24hr precipitation.
>
> but when I run it as :
>
> pcp_combine -subtract pgbf192.grapes.2011070612 192
pgbf168.grapes.2011070612 168 zb.nc
>
> there are always say that
>
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
VarInfo "APCP/A192" in GRIB file "pgbf192.grapes.2011070612".
> WARNING:
> ERROR  :
> ERROR  : get_field() -> can't get data plane from file
"pgbf192.grapes.2011070612"
> ERROR  :
>
> the grib index for 192hr is
>
>
352:2522142:d=11070612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=192:TimeU=1:sfc:192hr
fcst:NAve=0
>
> what wrong have I met?
>
> thanks alot
>
> zhaobin
>
> 2013-11-06
>
>
>
> zhaob
>

------------------------------------------------
Subject: Re: Re: [rt.rap.ucar.edu #64170] a question for pcp_combin
From: zhaob
Time: Sun Nov 10 16:58:24 2013

thanks alot for your help.

I can I can make it by your method.

yours

zhaobin


> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2013-11-08 02:14:03 (星期五)
> 收件人: zhaob at cma.gov.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #64170] a question for pcp_combin
>
> Zhaob,
>
> Sorry for the delay in getting back to you.  I think I know what the
problem is.  Typically, accumulated precipitation is stored in GRIB
files using an accumulation interval with a "time range"
> indicator value of 4.  Here's a description of the different time
range indicator values and meanings:
>     http://www.nco.ncep.noaa.gov/pmb/docs/on388/table5.html
>
> For example, take a look at the the APCP in the GRIB files included
in the MET tar ball:
>    wgrib
METv4.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 | grep
APCP
>
1:0:d=05080700:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr acc:NAve=0
>
2:31408:d=05080700:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=9:P2=12:TimeU=1:sfc:9-
12hr acc:NAve=0
>
> The "TR=4" indicates that these records contain an accumulation
between times P1 and P2.  In the first record, the precip is
accumulated between 0 and 12 hours.  In the second record, the precip
is
> accumulated between 9 and 12 hours.
>
> However, your GRIB data uses a time range indicator of 10, not 4.
pcp_combine is looking in "pgbf192.grapes.2011070612" for a 192 hour
*accumulation*.  But since the time range indicator is not 4, it
> doesn't find a match.  Perhaps you know that the this data really is
a runtime accumulation from 0 to 192 hours, but according to the GRIB
flags, it's a 192 hour instantaneous forecast, not an
> accumulation.
>
> If you have control over the generation of these GRIB files, you
could switch the time range indicator to 4.  If not, there's another
workaround.  Instead of telling pcp_combine to look for a
> particular accumulation interval, you can give it a more complete
description of the field you'd like to use from each file.  Please try
running the following:
>
>    pcp_combine -subtract \
>      pgbf192.grapes.2011070612 'name="APCP"; level="Z0";
lead_time="192";' \
>      pgbf168.grapes.2011070612 'name="APCP"; level="Z0";
lead_time="168";' \
>      zb.nc
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On 11/06/2013 01:22 AM, zhaob via RT wrote:
> >
> > Wed Nov 06 01:22:23 2013: Request 64170 was acted upon.
> > Transaction: Ticket created by zhaob at cma.gov.cn
> >         Queue: met_help
> >       Subject: a question for pcp_combin
> >         Owner: Nobody
> >    Requestors: zhaob at cma.gov.cn
> >        Status: new
> >   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=64170 >
> >
> >
> > hi!
> >
> > I am runing met4.0 and use pcp_combine to subtract 192hr accum
precipation and 168  accum precipation  to get 24hr precipitation.
> >
> > but when I run it as :
> >
> > pcp_combine -subtract pgbf192.grapes.2011070612 192
pgbf168.grapes.2011070612 168 zb.nc
> >
> > there are always say that
> >
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for VarInfo "APCP/A192" in GRIB file "pgbf192.grapes.2011070612".
> > WARNING:
> > ERROR  :
> > ERROR  : get_field() -> can't get data plane from file
"pgbf192.grapes.2011070612"
> > ERROR  :
> >
> > the grib index for 192hr is
> >
> >
352:2522142:d=11070612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=192:TimeU=1:sfc:192hr
fcst:NAve=0
> >
> > what wrong have I met?
> >
> > thanks alot
> >
> > zhaobin
> >
> > 2013-11-06
> >
> >
> >
> > zhaob
> >
>


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


More information about the Met_help mailing list