[Met_help] [rt.rap.ucar.edu #75994] History for About pcp_combine
Julie Prestopnik via RT
met_help at ucar.edu
Mon May 2 15:00:59 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear all,
I received the following errors message while running the pcp_combine of met-5.0. But when I used plot_data_plane tool to plot the APCP variable, and I can get the plot correctly.I have put the test data on your ftp site at incoming/irap/met_help/liu_data directory.
DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
WARNING:
ERROR :
ERROR : get_field() -> can't get data plane from file "./rmf_gra_2016040600.24"
ERROR :
Kindly advise.
Thanks in advance.
Yuanpu
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Wed Apr 20 04:17:43 2016
Dear all,
I received the following errors message while running the pcp_combine
of met-5.0. But when I used plot_data_plane tool to plot the APCP
variable, and I can get the plot correctly.I have put the test data on
your ftp site at incoming/irap/met_help/liu_data directory.
DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
WARNING:
ERROR :
ERROR : get_field() -> can't get data plane from file
"./rmf_gra_2016040600.24"
ERROR :
Kindly advise.
Thanks in advance.
Yuanpu
------------------------------------------------
Subject: About pcp_combine
From: Julie Prestopnik
Time: Wed Apr 20 10:11:23 2016
Hi Yuanpu.
I see that you are having problems running pcp_combine from met-5.0.
Typically, accumulated precipitation is stored in GRIB files using an
accumulation interval with a "time range" indicator value of 4. Here
is 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 APCP in the GRIB files included in the
MET
tar ball:
wgrib met-5.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 5, not 4.
wgrib rmf_gra_2016040600.24 | grep APCP
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
diff:NAve=0
pcp_combine is looking in "rmf_gra_2016040600.24" for a 24 hour
*accumulation*, but since the time range indicator is no 4, it doesn't
find
a match.
If you have control over the generation of these GRIB files, you could
switch the time range indicator to 4. If not, there is 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 -add rmf_gra_2016040600.24 'name="APCP"; level="L0-24";'
rmf_gra_2016040600_APCP_00_24.nc
I hope that helps!
Julie
On Wed, Apr 20, 2016 at 3:59 AM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> Wed Apr 20 03:59:14 2016: Request 75994 was acted upon.
> Transaction: Ticket created by lyp at lzb.ac.cn
> Queue: met_help
> Subject: About pcp_combine
> Owner: Nobody
> Requestors: lyp at lzb.ac.cn
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
>
>
> Dear all,
>
> I received the following errors message while running the
pcp_combine of
> met-5.0. But when I used plot_data_plane tool to plot the APCP
variable,
> and I can get the plot correctly.I have put the test data on your
ftp site
> at incoming/irap/met_help/liu_data directory.
>
>
> DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
> WARNING:
> ERROR :
> ERROR : get_field() -> can't get data plane from file "./rmf_gra_
> 2016040600.24"
> ERROR :
>
>
> Kindly advise.
>
> Thanks in advance.
>
> Yuanpu
>
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Wed Apr 20 20:42:32 2016
Hello Julie
Thanks for your reply.
I used pcp_combine tool to processed the precipitation, and I want to
get 24 hours precipitation. When I used the following command, it
doen't work.
../met-install/bin/pcp_combine -subtract \
./rmf_gra_2016040600.24 \
'name="APCP";level="L0-24";'\
./rmf_gra_2016040600.00 \
'name="APCP";level="L0-0";'\
./rmf_gra_APCP24.024.nc
../met-install/bin/pcp_combine -subtract \
./rmf_gra_2016040600.48 \
'name="APCP";level="L0-48";'\
./rmf_gra_2016040600.24 \
'name="APCP";level="L0-24";'\
./rmf_gra_APCP24.048.nc
I also used ../met-install/bin/pcp_combine -add rmf_gra_2016040600.48
'name="APCP"; level="L0-48";' rmf_gra_2016040600_APCP_00_48.nc ,and
then used
../met-install/bin/pcp_combine -subtract \
./rmf_gra_2016040600_APCP_00_48.nc 48 \
./rmf_gra_2016040600_APCP_00_24.nc 24 \
./rmf_gra.48.nc
it appears follows errors:
DEBUG 1: Reading input file: ./rmf_gra_2016040600_APCP_00_48.nc
ERROR :
ERROR : get_field() -> can't get data plane from file
"./rmf_gra_2016040600_APCP_00_48.nc"
ERROR :
Kindly advise.
Thanks in advance.
Yuanpu
> -----原始邮件-----
> 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> 发送时间: 2016年4月21日 星期四
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #75994] About pcp_combine
>
> Hi Yuanpu.
>
> I see that you are having problems running pcp_combine from met-5.0.
>
> Typically, accumulated precipitation is stored in GRIB files using
an
> accumulation interval with a "time range" indicator value of 4.
Here is 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 APCP in the GRIB files included in
the MET
> tar ball:
>
> wgrib met-5.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 5, not 4.
>
> wgrib rmf_gra_2016040600.24 | grep APCP
>
>
>
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> diff:NAve=0
>
> pcp_combine is looking in "rmf_gra_2016040600.24" for a 24 hour
> *accumulation*, but since the time range indicator is no 4, it
doesn't find
> a match.
>
> If you have control over the generation of these GRIB files, you
could
> switch the time range indicator to 4. If not, there is 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 -add rmf_gra_2016040600.24 'name="APCP"; level="L0-24";'
> rmf_gra_2016040600_APCP_00_24.nc
>
> I hope that helps!
>
> Julie
>
> On Wed, Apr 20, 2016 at 3:59 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > Wed Apr 20 03:59:14 2016: Request 75994 was acted upon.
> > Transaction: Ticket created by lyp at lzb.ac.cn
> > Queue: met_help
> > Subject: About pcp_combine
> > Owner: Nobody
> > Requestors: lyp at lzb.ac.cn
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
> >
> >
> > Dear all,
> >
> > I received the following errors message while running the
pcp_combine of
> > met-5.0. But when I used plot_data_plane tool to plot the APCP
variable,
> > and I can get the plot correctly.I have put the test data on your
ftp site
> > at incoming/irap/met_help/liu_data directory.
> >
> >
> > DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
> > WARNING:
> > ERROR :
> > ERROR : get_field() -> can't get data plane from file "./rmf_gra_
> > 2016040600.24"
> > ERROR :
> >
> >
> > Kindly advise.
> >
> > Thanks in advance.
> >
> > Yuanpu
> >
> >
> >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
------------------------------------------------
Subject: About pcp_combine
From: Julie Prestopnik
Time: Thu Apr 21 09:45:00 2016
Hi Yuanpu.
I understand that you would like to get 24 hours of precipitation.
The
file, rmf_gra_2016040600.24, already contains the total for 24 hours
of
precipitation. By running:
wgrib rmf_gra_2016040600.24 | grep APCP
we can see in the "sfc" field in the output that the value is '0-
24hr':
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
diff:NAve=0
So that file already contains the 24 hours of precipitation. For
simplicity in scripting your calls to the other MET tools, you may
find it
easier to get the 24 hours of precipitation into NetCDF format as
well. To
do that, you could run the following as I suggested previously:
pcp_combine -add rmf_gra_2016040600.24 'name="APCP"; level="L0-24";'
rmf_gra_2016040600_APCP_00_24.nc
You mentioned that when you ran the following command it did not work:
../met-install/bin/pcp_combine -subtract \
./rmf_gra_2016040600.24 \
'name="APCP";level="L0-24";'\
./rmf_gra_2016040600.00 \
'name="APCP";level="L0-0";'\
./rmf_gra_APCP24.024.nc
Can you please provide more information about how it did not work?
What
error did you received? I was able to run the command successfully
here:
pcp_combine -subtract rmf_gra_2016040600.24 'name="APCP";level="L0-
24";'
rmf_gra_2016040600.00 'name="APCP";level="L0-0";'
rmf_gra_APCP24.024.nc
To be clear, the output of running that command produces the same file
as
the command I suggested running since the file, rmf_gra_2016040600.24,
already contains the total for 24 hours of precipitation.
In order to get the 24 hour precipitation by subtracting the 24 hour
from
the 48 hour file that you have, your command that you sent to generate
that
output looks correct to me:
../met-install/bin/pcp_combine -subtract \
./rmf_gra_2016040600.48 \
'name="APCP";level="L0-48";'\
./rmf_gra_2016040600.24 \
'name="APCP";level="L0-24";'\
./rmf_gra_APCP24.048.nc
If it is not working for you, please let us know what the error is and
please also upload the 48 hour file to the ftp server, so we can take
a
look and try to reproduce the results that you are getting.
I hope this helps.
Thanks,
Julie
On Wed, Apr 20, 2016 at 8:42 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
>
> Hello Julie
>
> Thanks for your reply.
>
> I used pcp_combine tool to processed the precipitation, and I want
to get
> 24 hours precipitation. When I used the following command, it doen't
work.
>
> ../met-install/bin/pcp_combine -subtract \
> ./rmf_gra_2016040600.24 \
> 'name="APCP";level="L0-24";'\
> ./rmf_gra_2016040600.00 \
> 'name="APCP";level="L0-0";'\
> ./rmf_gra_APCP24.024.nc
>
> ../met-install/bin/pcp_combine -subtract \
> ./rmf_gra_2016040600.48 \
> 'name="APCP";level="L0-48";'\
> ./rmf_gra_2016040600.24 \
> 'name="APCP";level="L0-24";'\
> ./rmf_gra_APCP24.048.nc
>
>
> I also used ../met-install/bin/pcp_combine -add
rmf_gra_2016040600.48
> 'name="APCP"; level="L0-48";' rmf_gra_2016040600_APCP_00_48.nc ,and
then
> used
>
> ../met-install/bin/pcp_combine -subtract \
> ./rmf_gra_2016040600_APCP_00_48.nc 48 \
> ./rmf_gra_2016040600_APCP_00_24.nc 24 \
> ./rmf_gra.48.nc
>
>
> it appears follows errors:
>
> DEBUG 1: Reading input file: ./rmf_gra_2016040600_APCP_00_48.nc
> ERROR :
> ERROR : get_field() -> can't get data plane from file "./rmf_gra_
> 2016040600_APCP_00_48.nc"
> ERROR :
>
> Kindly advise.
>
> Thanks in advance.
> Yuanpu
>
>
> > -----原始邮件-----
> > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > 发送时间: 2016年4月21日 星期四
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: [rt.rap.ucar.edu #75994] About pcp_combine
> >
> > Hi Yuanpu.
> >
> > I see that you are having problems running pcp_combine from met-
5.0.
> >
> > Typically, accumulated precipitation is stored in GRIB files using
an
> > accumulation interval with a "time range" indicator value of 4.
Here is
> 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 APCP in the GRIB files included in
the
> MET
> > tar ball:
> >
> > wgrib met-
5.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 5, not 4.
> >
> > wgrib rmf_gra_2016040600.24 | grep APCP
> >
> >
> >
>
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> > diff:NAve=0
> >
> > pcp_combine is looking in "rmf_gra_2016040600.24" for a 24 hour
> > *accumulation*, but since the time range indicator is no 4, it
doesn't
> find
> > a match.
> >
> > If you have control over the generation of these GRIB files, you
could
> > switch the time range indicator to 4. If not, there is 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 -add rmf_gra_2016040600.24 'name="APCP"; level="L0-
24";'
> > rmf_gra_2016040600_APCP_00_24.nc
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Wed, Apr 20, 2016 at 3:59 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > Wed Apr 20 03:59:14 2016: Request 75994 was acted upon.
> > > Transaction: Ticket created by lyp at lzb.ac.cn
> > > Queue: met_help
> > > Subject: About pcp_combine
> > > Owner: Nobody
> > > Requestors: lyp at lzb.ac.cn
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994
> >
> > >
> > >
> > > Dear all,
> > >
> > > I received the following errors message while running the
pcp_combine
> of
> > > met-5.0. But when I used plot_data_plane tool to plot the APCP
> variable,
> > > and I can get the plot correctly.I have put the test data on
your ftp
> site
> > > at incoming/irap/met_help/liu_data directory.
> > >
> > >
> > > DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > > VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
> > > WARNING:
> > > ERROR :
> > > ERROR : get_field() -> can't get data plane from file
"./rmf_gra_
> > > 2016040600.24"
> > > ERROR :
> > >
> > >
> > > Kindly advise.
> > >
> > > Thanks in advance.
> > >
> > > Yuanpu
> > >
> > >
> > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
>
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Fri Apr 22 01:14:20 2016
Hi Julie.
Thank you for your help.
I have resolved the problem. The reason is that I lost a spase on the
fllowing line of the command which I have sent to you.
'name="APCP";level="L0-24";'\
Thanks again,
yuanpu
> -----原始邮件-----
> 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> 发送时间: 2016年4月21日 星期四
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: [rt.rap.ucar.edu #75994] About pcp_combine
>
> Hi Yuanpu.
>
> I understand that you would like to get 24 hours of precipitation.
The
> file, rmf_gra_2016040600.24, already contains the total for 24 hours
of
> precipitation. By running:
>
> wgrib rmf_gra_2016040600.24 | grep APCP
>
> we can see in the "sfc" field in the output that the value is '0-
24hr':
>
>
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> diff:NAve=0
>
> So that file already contains the 24 hours of precipitation. For
> simplicity in scripting your calls to the other MET tools, you may
find it
> easier to get the 24 hours of precipitation into NetCDF format as
well. To
> do that, you could run the following as I suggested previously:
>
> pcp_combine -add rmf_gra_2016040600.24 'name="APCP"; level="L0-24";'
> rmf_gra_2016040600_APCP_00_24.nc
>
> You mentioned that when you ran the following command it did not
work:
>
> ../met-install/bin/pcp_combine -subtract \
> ./rmf_gra_2016040600.24 \
> 'name="APCP";level="L0-24";'\
> ./rmf_gra_2016040600.00 \
> 'name="APCP";level="L0-0";'\
> ./rmf_gra_APCP24.024.nc
>
> Can you please provide more information about how it did not work?
What
> error did you received? I was able to run the command successfully
here:
>
> pcp_combine -subtract rmf_gra_2016040600.24 'name="APCP";level="L0-
24";'
> rmf_gra_2016040600.00 'name="APCP";level="L0-0";'
rmf_gra_APCP24.024.nc
>
> To be clear, the output of running that command produces the same
file as
> the command I suggested running since the file,
rmf_gra_2016040600.24,
> already contains the total for 24 hours of precipitation.
>
> In order to get the 24 hour precipitation by subtracting the 24 hour
from
> the 48 hour file that you have, your command that you sent to
generate that
> output looks correct to me:
>
> ../met-install/bin/pcp_combine -subtract \
> ./rmf_gra_2016040600.48 \
> 'name="APCP";level="L0-48";'\
> ./rmf_gra_2016040600.24 \
> 'name="APCP";level="L0-24";'\
> ./rmf_gra_APCP24.048.nc
>
> If it is not working for you, please let us know what the error is
and
> please also upload the 48 hour file to the ftp server, so we can
take a
> look and try to reproduce the results that you are getting.
>
> I hope this helps.
>
> Thanks,
> Julie
>
> On Wed, Apr 20, 2016 at 8:42 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
> >
> > Hello Julie
> >
> > Thanks for your reply.
> >
> > I used pcp_combine tool to processed the precipitation, and I want
to get
> > 24 hours precipitation. When I used the following command, it
doen't work.
> >
> > ../met-install/bin/pcp_combine -subtract \
> > ./rmf_gra_2016040600.24 \
> > 'name="APCP";level="L0-24";'\
> > ./rmf_gra_2016040600.00 \
> > 'name="APCP";level="L0-0";'\
> > ./rmf_gra_APCP24.024.nc
> >
> > ../met-install/bin/pcp_combine -subtract \
> > ./rmf_gra_2016040600.48 \
> > 'name="APCP";level="L0-48";'\
> > ./rmf_gra_2016040600.24 \
> > 'name="APCP";level="L0-24";'\
> > ./rmf_gra_APCP24.048.nc
> >
> >
> > I also used ../met-install/bin/pcp_combine -add
rmf_gra_2016040600.48
> > 'name="APCP"; level="L0-48";' rmf_gra_2016040600_APCP_00_48.nc
,and then
> > used
> >
> > ../met-install/bin/pcp_combine -subtract \
> > ./rmf_gra_2016040600_APCP_00_48.nc 48 \
> > ./rmf_gra_2016040600_APCP_00_24.nc 24 \
> > ./rmf_gra.48.nc
> >
> >
> > it appears follows errors:
> >
> > DEBUG 1: Reading input file: ./rmf_gra_2016040600_APCP_00_48.nc
> > ERROR :
> > ERROR : get_field() -> can't get data plane from file "./rmf_gra_
> > 2016040600_APCP_00_48.nc"
> > ERROR :
> >
> > Kindly advise.
> >
> > Thanks in advance.
> > Yuanpu
> >
> >
> > > -----原始邮件-----
> > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > 发送时间: 2016年4月21日 星期四
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: [rt.rap.ucar.edu #75994] About pcp_combine
> > >
> > > Hi Yuanpu.
> > >
> > > I see that you are having problems running pcp_combine from met-
5.0.
> > >
> > > Typically, accumulated precipitation is stored in GRIB files
using an
> > > accumulation interval with a "time range" indicator value of 4.
Here is
> > 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 APCP in the GRIB files included
in the
> > MET
> > > tar ball:
> > >
> > > wgrib met-
5.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 5, not 4.
> > >
> > > wgrib rmf_gra_2016040600.24 | grep APCP
> > >
> > >
> > >
> >
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> > > diff:NAve=0
> > >
> > > pcp_combine is looking in "rmf_gra_2016040600.24" for a 24 hour
> > > *accumulation*, but since the time range indicator is no 4, it
doesn't
> > find
> > > a match.
> > >
> > > If you have control over the generation of these GRIB files, you
could
> > > switch the time range indicator to 4. If not, there is 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 -add rmf_gra_2016040600.24 'name="APCP"; level="L0-
24";'
> > > rmf_gra_2016040600_APCP_00_24.nc
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > > On Wed, Apr 20, 2016 at 3:59 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > Wed Apr 20 03:59:14 2016: Request 75994 was acted upon.
> > > > Transaction: Ticket created by lyp at lzb.ac.cn
> > > > Queue: met_help
> > > > Subject: About pcp_combine
> > > > Owner: Nobody
> > > > Requestors: lyp at lzb.ac.cn
> > > > Status: new
> > > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994
> > >
> > > >
> > > >
> > > > Dear all,
> > > >
> > > > I received the following errors message while running the
pcp_combine
> > of
> > > > met-5.0. But when I used plot_data_plane tool to plot the APCP
> > variable,
> > > > and I can get the plot correctly.I have put the test data on
your ftp
> > site
> > > > at incoming/irap/met_help/liu_data directory.
> > > >
> > > >
> > > > DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found for
> > > > VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
> > > > WARNING:
> > > > ERROR :
> > > > ERROR : get_field() -> can't get data plane from file
"./rmf_gra_
> > > > 2016040600.24"
> > > > ERROR :
> > > >
> > > >
> > > > Kindly advise.
> > > >
> > > > Thanks in advance.
> > > >
> > > > Yuanpu
> > > >
> > > >
> > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> >
> >
> >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
------------------------------------------------
Subject: About pcp_combine
From: Julie Prestopnik
Time: Fri Apr 22 08:39:04 2016
That's great, Yuanpu. Thank you for letting us know.
Julie
On Fri, Apr 22, 2016 at 1:14 AM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
>
> Hi Julie.
>
> Thank you for your help.
>
> I have resolved the problem. The reason is that I lost a spase on
the
> fllowing line of the command which I have sent to you.
>
> 'name="APCP";level="L0-24";'\
>
> Thanks again,
>
> yuanpu
> > -----原始邮件-----
> > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > 发送时间: 2016年4月21日 星期四
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: [rt.rap.ucar.edu #75994] About pcp_combine
> >
> > Hi Yuanpu.
> >
> > I understand that you would like to get 24 hours of precipitation.
The
> > file, rmf_gra_2016040600.24, already contains the total for 24
hours of
> > precipitation. By running:
> >
> > wgrib rmf_gra_2016040600.24 | grep APCP
> >
> > we can see in the "sfc" field in the output that the value is '0-
24hr':
> >
> >
>
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> > diff:NAve=0
> >
> > So that file already contains the 24 hours of precipitation. For
> > simplicity in scripting your calls to the other MET tools, you may
find
> it
> > easier to get the 24 hours of precipitation into NetCDF format as
well.
> To
> > do that, you could run the following as I suggested previously:
> >
> > pcp_combine -add rmf_gra_2016040600.24 'name="APCP"; level="L0-
24";'
> > rmf_gra_2016040600_APCP_00_24.nc
> >
> > You mentioned that when you ran the following command it did not
work:
> >
> > ../met-install/bin/pcp_combine -subtract \
> > ./rmf_gra_2016040600.24 \
> > 'name="APCP";level="L0-24";'\
> > ./rmf_gra_2016040600.00 \
> > 'name="APCP";level="L0-0";'\
> > ./rmf_gra_APCP24.024.nc
> >
> > Can you please provide more information about how it did not work?
What
> > error did you received? I was able to run the command
successfully here:
> >
> > pcp_combine -subtract rmf_gra_2016040600.24
'name="APCP";level="L0-24";'
> > rmf_gra_2016040600.00 'name="APCP";level="L0-0";'
rmf_gra_APCP24.024.nc
> >
> > To be clear, the output of running that command produces the same
file as
> > the command I suggested running since the file,
rmf_gra_2016040600.24,
> > already contains the total for 24 hours of precipitation.
> >
> > In order to get the 24 hour precipitation by subtracting the 24
hour from
> > the 48 hour file that you have, your command that you sent to
generate
> that
> > output looks correct to me:
> >
> > ../met-install/bin/pcp_combine -subtract \
> > ./rmf_gra_2016040600.48 \
> > 'name="APCP";level="L0-48";'\
> > ./rmf_gra_2016040600.24 \
> > 'name="APCP";level="L0-24";'\
> > ./rmf_gra_APCP24.048.nc
> >
> > If it is not working for you, please let us know what the error is
and
> > please also upload the 48 hour file to the ftp server, so we can
take a
> > look and try to reproduce the results that you are getting.
> >
> > I hope this helps.
> >
> > Thanks,
> > Julie
> >
> > On Wed, Apr 20, 2016 at 8:42 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994 >
> > >
> > > Hello Julie
> > >
> > > Thanks for your reply.
> > >
> > > I used pcp_combine tool to processed the precipitation, and I
want to
> get
> > > 24 hours precipitation. When I used the following command, it
doen't
> work.
> > >
> > > ../met-install/bin/pcp_combine -subtract \
> > > ./rmf_gra_2016040600.24 \
> > > 'name="APCP";level="L0-24";'\
> > > ./rmf_gra_2016040600.00 \
> > > 'name="APCP";level="L0-0";'\
> > > ./rmf_gra_APCP24.024.nc
> > >
> > > ../met-install/bin/pcp_combine -subtract \
> > > ./rmf_gra_2016040600.48 \
> > > 'name="APCP";level="L0-48";'\
> > > ./rmf_gra_2016040600.24 \
> > > 'name="APCP";level="L0-24";'\
> > > ./rmf_gra_APCP24.048.nc
> > >
> > >
> > > I also used ../met-install/bin/pcp_combine -add
rmf_gra_2016040600.48
> > > 'name="APCP"; level="L0-48";' rmf_gra_2016040600_APCP_00_48.nc
,and
> then
> > > used
> > >
> > > ../met-install/bin/pcp_combine -subtract \
> > > ./rmf_gra_2016040600_APCP_00_48.nc 48 \
> > > ./rmf_gra_2016040600_APCP_00_24.nc 24 \
> > > ./rmf_gra.48.nc
> > >
> > >
> > > it appears follows errors:
> > >
> > > DEBUG 1: Reading input file: ./rmf_gra_2016040600_APCP_00_48.nc
> > > ERROR :
> > > ERROR : get_field() -> can't get data plane from file
"./rmf_gra_
> > > 2016040600_APCP_00_48.nc"
> > > ERROR :
> > >
> > > Kindly advise.
> > >
> > > Thanks in advance.
> > > Yuanpu
> > >
> > >
> > > > -----原始邮件-----
> > > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > > 发送时间: 2016年4月21日 星期四
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: [rt.rap.ucar.edu #75994] About pcp_combine
> > > >
> > > > Hi Yuanpu.
> > > >
> > > > I see that you are having problems running pcp_combine from
met-5.0.
> > > >
> > > > Typically, accumulated precipitation is stored in GRIB files
using an
> > > > accumulation interval with a "time range" indicator value of
4.
> Here is
> > > 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 APCP in the GRIB files
included in
> the
> > > MET
> > > > tar ball:
> > > >
> > > > wgrib met-
5.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 5, not
4.
> > > >
> > > > wgrib rmf_gra_2016040600.24 | grep APCP
> > > >
> > > >
> > > >
> > >
>
291:28360360:d=16040600:APCP:kpds5=61:kpds6=1:kpds7=0:TR=5:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> > > > diff:NAve=0
> > > >
> > > > pcp_combine is looking in "rmf_gra_2016040600.24" for a 24
hour
> > > > *accumulation*, but since the time range indicator is no 4, it
> doesn't
> > > find
> > > > a match.
> > > >
> > > > If you have control over the generation of these GRIB files,
you
> could
> > > > switch the time range indicator to 4. If not, there is
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 -add rmf_gra_2016040600.24 'name="APCP";
level="L0-24";'
> > > > rmf_gra_2016040600_APCP_00_24.nc
> > > >
> > > > I hope that helps!
> > > >
> > > > Julie
> > > >
> > > > On Wed, Apr 20, 2016 at 3:59 AM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > Wed Apr 20 03:59:14 2016: Request 75994 was acted upon.
> > > > > Transaction: Ticket created by lyp at lzb.ac.cn
> > > > > Queue: met_help
> > > > > Subject: About pcp_combine
> > > > > Owner: Nobody
> > > > > Requestors: lyp at lzb.ac.cn
> > > > > Status: new
> > > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75994
> > > >
> > > > >
> > > > >
> > > > > Dear all,
> > > > >
> > > > > I received the following errors message while running the
> pcp_combine
> > > of
> > > > > met-5.0. But when I used plot_data_plane tool to plot the
APCP
> > > variable,
> > > > > and I can get the plot correctly.I have put the test data on
your
> ftp
> > > site
> > > > > at incoming/irap/met_help/liu_data directory.
> > > > >
> > > > >
> > > > > DEBUG 1: Reading input file: ./rmf_gra_2016040600.24
> > > > > WARNING:
> > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found for
> > > > > VarInfo "APCP/A24" in GRIB file "./rmf_gra_2016040600.24".
> > > > > WARNING:
> > > > > ERROR :
> > > > > ERROR : get_field() -> can't get data plane from file
"./rmf_gra_
> > > > > 2016040600.24"
> > > > > ERROR :
> > > > >
> > > > >
> > > > > Kindly advise.
> > > > >
> > > > > Thanks in advance.
> > > > >
> > > > > Yuanpu
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > >
> > >
> > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
>
>
>
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
More information about the Met_help
mailing list