[Met_help] [rt.rap.ucar.edu #77989] History for About pcp_combine
Julie Prestopnik via RT
met_help at ucar.edu
Tue Oct 18 14:07:17 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear all,
I received the following errors message while running the pcp_combine of met-5.0. I have put the test data on your ftp site at incoming/irap/met_help/liu_data directory.
DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
ERROR :
ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not found in grib record
ERROR :
ERROR :
ERROR : GribFile::seek_record(int) -> range check error
ERROR :
DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1
ERROR :
ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not found in grib record
ERROR :
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1".
WARNING:
ERROR :
ERROR : get_field() -> can't get data plane from file "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1"
ERROR :
Kindly advise.
Thanks in advance.
Yuanpu
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Tue Sep 20 01:57:24 2016
Dear all,
I received the following errors message while running the pcp_combine
of met-5.0. I have put the test data on your ftp site at
incoming/irap/met_help/liu_data directory.
DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
ERROR :
ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found in grib record
ERROR :
ERROR :
ERROR : GribFile::seek_record(int) -> range check error
ERROR :
DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1
ERROR :
ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found in grib record
ERROR :
WARNING:
WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
VarInfo "TP/Z0" in GRIB file
"/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1".
WARNING:
ERROR :
ERROR : get_field() -> can't get data plane from file
"/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_9_grid_simple.grib1"
ERROR :
Kindly advise.
Thanks in advance.
Yuanpu
------------------------------------------------
Subject: About pcp_combine
From: Julie Prestopnik
Time: Tue Sep 20 09:57:31 2016
Hi Yuanpu.
Could you please send me the command the you are running?
I think part of the problem with the GRIB data you sent to us might be
an
encoding issue. I tried to use the standard NCEP wgrib tool to list
out
the records present in the GRIB file, but it gives the following
error:
jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
Len of grib message is inconsistent.
missing end section
46 f 46 f
Are you able to list the records in that GRIB file using wgrib? If
so,
perhaps the file is being corrupted by the FTP transfer. Would you
mind
putting all the files you uploaded into a tarball and then uploading
that?
That should preserve the binary files and help avoid corruption.
Sorry for
the trouble.
Thanks,
Julie Prestopnik
On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> Tue Sep 20 01:27:52 2016: Request 77989 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=77989 >
>
>
>
> Dear all,
>
> I received the following errors message while running the
pcp_combine of
> met-5.0. I have put the test data on your ftp site at
> incoming/irap/met_help/liu_data directory.
>
> DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> 2016091900/sfc_20160919_6_grid_simple.grib1
> ERROR :
> ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found
> in grib record
> ERROR :
> ERROR :
> ERROR : GribFile::seek_record(int) -> range check error
> ERROR :
> DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> simple/2016091900/sfc_20160919_9_grid_simple.grib1
> ERROR :
> ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found
> in grib record
> ERROR :
> WARNING:
> WARNING: MetGrib1DataFile::data_plane() -> No exact match found for
> VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_
> simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> WARNING:
> ERROR :
> ERROR : get_field() -> can't get data plane from file
> "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> 20160919_9_grid_simple.grib1"
> 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: Tue Sep 20 20:03:48 2016
Hi Julie
Thank you for your reply.
I used the pcp_combine commands like follows, and when I used wgrib
tool to list the recoreds, I can get the right records as follows. So,
I have put the test tar file on your ftp site at
incoming/irap/met_help/liu_data directory.
COMMAND LINES:
../met-install/bin/pcp_combine -subtract \
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
'name="TP";level="Z0";GRIB1_ptv=128;' \
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
'name="TP";level="Z0";GRIB1_ptv=128;' \
${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_grid_simple.nc
[/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
0-6hr:NAve=0
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
3-6hr:NAve=0
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
.............................
Thanks,
yuanpu
> -----原始邮件-----
> 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> 发送时间: 2016年9月20日 星期二
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
>
> Hi Yuanpu.
>
> Could you please send me the command the you are running?
>
> I think part of the problem with the GRIB data you sent to us might
be an
> encoding issue. I tried to use the standard NCEP wgrib tool to list
out
> the records present in the GRIB file, but it gives the following
error:
>
> jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
>
> Len of grib message is inconsistent.
>
>
> missing end section
> 46 f 46 f
>
> Are you able to list the records in that GRIB file using wgrib? If
so,
> perhaps the file is being corrupted by the FTP transfer. Would you
mind
> putting all the files you uploaded into a tarball and then uploading
that?
> That should preserve the binary files and help avoid corruption.
Sorry for
> the trouble.
>
> Thanks,
> Julie Prestopnik
>
> On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > Tue Sep 20 01:27:52 2016: Request 77989 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=77989 >
> >
> >
> >
> > Dear all,
> >
> > I received the following errors message while running the
pcp_combine of
> > met-5.0. I have put the test data on your ftp site at
> > incoming/irap/met_help/liu_data directory.
> >
> > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> > 2016091900/sfc_20160919_6_grid_simple.grib1
> > ERROR :
> > ERROR : GribFile::read_record(GribRecord &) -> trailing "7777"
not found
> > in grib record
> > ERROR :
> > ERROR :
> > ERROR : GribFile::seek_record(int) -> range check error
> > ERROR :
> > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > ERROR :
> > ERROR : GribFile::read_record(GribRecord &) -> trailing "7777"
not found
> > in grib record
> > ERROR :
> > WARNING:
> > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_
> > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > WARNING:
> > ERROR :
> > ERROR : get_field() -> can't get data plane from file
> > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > 20160919_9_grid_simple.grib1"
> > 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: Wed Sep 21 10:10:49 2016
Hi Yuanpu.
Thank you for creating a tar file with the data.
Can you please tell me how the GRIB1 data that you are trying to use
with
MET was generated? Is it the output of the WRF Unified Post Processor
(UPP)?
It appear there may be an encoding problem that is causing this error.
Unfortunately, I am still unable to list the records for
sfc_20160919_6_grid_simple.grib1 using wgrib. I am able to use wgrib
to
view the records for sfc_20160919_9_grid_simple.grib1 but the same
issue
exists as when I use wgrib on the sfc_20160919_6_grid_simple.grib1
file:
Undefined parameter table (center 98-0 table 228), using NCEP-opn
Len of grib message is inconsistent.
missing end section
46 f 46 f
1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
In that output, I do not see the field "TP" that you wish to get the
values
for. Since wgrib doesn't know anything about table number 228 for
center
98, we can't trust the abbreviations wgrib spits out. It may be that
the
field exists, but due to an encoding issue, I am (and the MET code is)
unable to see it.
Thanks,
Julie
On Tue, Sep 20, 2016 at 8:03 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
>
> Hi Julie
>
> Thank you for your reply.
>
> I used the pcp_combine commands like follows, and when I used wgrib
tool
> to list the recoreds, I can get the right records as follows. So, I
have
> put the test tar file on your ftp site at
incoming/irap/met_help/liu_data
> directory.
>
>
> COMMAND LINES:
>
> ../met-install/bin/pcp_combine -subtract \
> ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1
\
> 'name="TP";level="Z0";GRIB1_ptv=128;' \
> ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
> 'name="TP";level="Z0";GRIB1_ptv=128;' \
> ${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_
> grid_simple.nc
>
> [/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
>
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> 0-6hr:NAve=0
>
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> 3-6hr:NAve=0
>
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
> .............................
>
> Thanks,
>
> yuanpu
> > -----原始邮件-----
> > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > 发送时间: 2016年9月20日 星期二
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> >
> > Hi Yuanpu.
> >
> > Could you please send me the command the you are running?
> >
> > I think part of the problem with the GRIB data you sent to us
might be an
> > encoding issue. I tried to use the standard NCEP wgrib tool to
list out
> > the records present in the GRIB file, but it gives the following
error:
> >
> > jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
> >
> > Len of grib message is inconsistent.
> >
> >
> > missing end section
> > 46 f 46 f
> >
> > Are you able to list the records in that GRIB file using wgrib?
If so,
> > perhaps the file is being corrupted by the FTP transfer. Would
you mind
> > putting all the files you uploaded into a tarball and then
uploading
> that?
> > That should preserve the binary files and help avoid corruption.
Sorry
> for
> > the trouble.
> >
> > Thanks,
> > Julie Prestopnik
> >
> > On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > Tue Sep 20 01:27:52 2016: Request 77989 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=77989
> >
> > >
> > >
> > >
> > > Dear all,
> > >
> > > I received the following errors message while running the
pcp_combine
> of
> > > met-5.0. I have put the test data on your ftp site at
> > > incoming/irap/met_help/liu_data directory.
> > >
> > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/
> > > 2016091900/sfc_20160919_6_grid_simple.grib1
> > > ERROR :
> > > ERROR : GribFile::read_record(GribRecord &) -> trailing "7777"
not
> found
> > > in grib record
> > > ERROR :
> > > ERROR :
> > > ERROR : GribFile::seek_record(int) -> range check error
> > > ERROR :
> > > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> > > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > > ERROR :
> > > ERROR : GribFile::read_record(GribRecord &) -> trailing "7777"
not
> found
> > > in grib record
> > > ERROR :
> > > WARNING:
> > > WARNING: MetGrib1DataFile::data_plane() -> No exact match found
for
> > > VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_
> > > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > > WARNING:
> > > ERROR :
> > > ERROR : get_field() -> can't get data plane from file
> > > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > > 20160919_9_grid_simple.grib1"
> > > 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: Wed Sep 21 19:17:44 2016
Hi Julie,
Thanks for your help!
The grib1 data that I used was generated with grib_api tool, and the
original data was ECMWF forcast results.The original ECMWF data
enclude grib1 and grib2 format data in one file, and the grib_api tool
can decompose that data to girb1 file and grib2 file. When I used the
wgrib tool to list the records(as following) of test file ,I can see
the field "TP".
[/]#wgrib sfc_20160919_6_grid_simple.grib1
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
0-6hr:NAve=0
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
3-6hr:NAve=0
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
6:35794866:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
7:42020070:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
8:48245274:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
9:54470478:d=16091900:var216:kpds5=216:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
10:62770722:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
11:71070966:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
12:77296170:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
13:85596414:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
14:93896658:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
15:102196902:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
16:108422106:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
3-6hr:NAve=0
17:114647310:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
18:118797474:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
19:127097718:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
3-6hr:NAve=0
20:135397962:d=16091900:MN2T6:kpds5=122:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
0-6hr:NAve=0
21:141623166:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
22:147848370:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
23:154073574:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
24:160298778:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
25:166523982:d=16091900:var247:kpds5=247:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
26:174824226:d=16091900:var246:kpds5=246:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
27:183124470:d=16091900:MX2T6:kpds5=121:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
0-6hr:NAve=0
28:189349674:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
29:197649918:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
30:205950162:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
31:214250406:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
32:222550650:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
fcst:NAve=0
[/]#wgrib sfc_20160919_9_grid_simple.grib1
1:0:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
2:8300244:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
3:16600488:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
4:22825692:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
5:29050896:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
14:91302852:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
15:100121862:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
16:108422106:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
17:112572270:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
18:120872514:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
19:127097718:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
20:133322922:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
21:141623166:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
22:147848370:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
23:154073574:d=16091900:var216:kpds5=216:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
24:162373818:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
25:166523982:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
26:172749186:d=16091900:var246:kpds5=246:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
27:181049430:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
28:189349674:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
29:197649918:d=16091900:var247:kpds5=247:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
Thanks
yuanpu
> -----原始邮件-----
> 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> 发送时间: 2016年9月22日 星期四
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
>
> Hi Yuanpu.
>
> Thank you for creating a tar file with the data.
>
> Can you please tell me how the GRIB1 data that you are trying to use
with
> MET was generated? Is it the output of the WRF Unified Post
Processor
> (UPP)?
>
> It appear there may be an encoding problem that is causing this
error.
>
> Unfortunately, I am still unable to list the records for
> sfc_20160919_6_grid_simple.grib1 using wgrib. I am able to use
wgrib to
> view the records for sfc_20160919_9_grid_simple.grib1 but the same
issue
> exists as when I use wgrib on the sfc_20160919_6_grid_simple.grib1
file:
>
> Undefined parameter table (center 98-0 table 228), using NCEP-opn
> Len of grib message is inconsistent.
>
> missing end section
> 46 f 46 f
>
1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
> In that output, I do not see the field "TP" that you wish to get the
values
> for. Since wgrib doesn't know anything about table number 228 for
center
> 98, we can't trust the abbreviations wgrib spits out. It may be
that the
> field exists, but due to an encoding issue, I am (and the MET code
is)
> unable to see it.
>
> Thanks,
> Julie
>
> On Tue, Sep 20, 2016 at 8:03 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
> >
> > Hi Julie
> >
> > Thank you for your reply.
> >
> > I used the pcp_combine commands like follows, and when I used
wgrib tool
> > to list the recoreds, I can get the right records as follows. So,
I have
> > put the test tar file on your ftp site at
incoming/irap/met_help/liu_data
> > directory.
> >
> >
> > COMMAND LINES:
> >
> > ../met-install/bin/pcp_combine -subtract \
> >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
> > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
> > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > ${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_
> > grid_simple.nc
> >
> > [/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
> >
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> > 0-6hr:NAve=0
> >
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> > 3-6hr:NAve=0
> >
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> > .............................
> >
> > Thanks,
> >
> > yuanpu
> > > -----原始邮件-----
> > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > 发送时间: 2016年9月20日 星期二
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > >
> > > Hi Yuanpu.
> > >
> > > Could you please send me the command the you are running?
> > >
> > > I think part of the problem with the GRIB data you sent to us
might be an
> > > encoding issue. I tried to use the standard NCEP wgrib tool to
list out
> > > the records present in the GRIB file, but it gives the following
error:
> > >
> > > jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
> > >
> > > Len of grib message is inconsistent.
> > >
> > >
> > > missing end section
> > > 46 f 46 f
> > >
> > > Are you able to list the records in that GRIB file using wgrib?
If so,
> > > perhaps the file is being corrupted by the FTP transfer. Would
you mind
> > > putting all the files you uploaded into a tarball and then
uploading
> > that?
> > > That should preserve the binary files and help avoid corruption.
Sorry
> > for
> > > the trouble.
> > >
> > > Thanks,
> > > Julie Prestopnik
> > >
> > > On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > Tue Sep 20 01:27:52 2016: Request 77989 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=77989
> > >
> > > >
> > > >
> > > >
> > > > Dear all,
> > > >
> > > > I received the following errors message while running the
pcp_combine
> > of
> > > > met-5.0. I have put the test data on your ftp site at
> > > > incoming/irap/met_help/liu_data directory.
> > > >
> > > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/
> > > > 2016091900/sfc_20160919_6_grid_simple.grib1
> > > > ERROR :
> > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777" not
> > found
> > > > in grib record
> > > > ERROR :
> > > > ERROR :
> > > > ERROR : GribFile::seek_record(int) -> range check error
> > > > ERROR :
> > > > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > > > ERROR :
> > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777" not
> > found
> > > > in grib record
> > > > ERROR :
> > > > WARNING:
> > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found for
> > > > VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_
> > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > > > WARNING:
> > > > ERROR :
> > > > ERROR : get_field() -> can't get data plane from file
> > > > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > > > 20160919_9_grid_simple.grib1"
> > > > 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: John Halley Gotway
Time: Thu Sep 22 09:41:39 2016
Yuanpu,
Great, if wgrib is able to read your GRIB1 files then hopefully the
MET
tools will be able to as well. However there are some issues to
consider...
(1) The default settings for pcp_combine expect accumulated
precipitation
to be stored in GRIB code 61 (for APCP). Your total precip (TP)
variable
likely is not stored in GRIB code 61. Based on the wgrib output you
sent
it's stored in GRIB code 228 (i.e. kpds5=228).
(2) In your original email, you mentioned MET version 5.0. The MET
version
5.2 release included some pretty major updates to GRIB tables... so
using
5.0 versus 5.2 will have a pretty big impact. The GRIB table updates
in
5.2 are described here:
http://www.dtcenter.org/met/users/support/release_notes/METv5.2_release_notes.php
The GRIB tables in MET 5.2 include information about center 98 and
table
number 228, whereas MET 5.0 is limited to NCEP's GRIB table usage.
Using 5.0, setting 'name="TP";' will likely not work while it may work
in
5.2.
Regardless, instead of specifying the name, let's try explicitly
specifying
the GRIB code and parameter table version number. Does the following
command work?
../met-install/bin/pcp_combine -subtract \
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_grid_simple.
nc
Note that it's just a coincidence that the GRIB code of 228 for total
precip and the parameter table version number of 228 just happen to
the be
the same value.
Please let me know if that works.
I also find it helpful to visualize data like this using the plot data
plane tool:
../met-install/bin/plot_data_plane \
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
plot_TP.ps \
'GRIB1_code=228; GRIB1_ptv=228; level="Z0";'
And then open the plot_TP.ps PostScript file.
Thanks,
John Halley Gotway
On Wed, Sep 21, 2016 at 7:17 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
>
> Hi Julie,
>
> Thanks for your help!
>
> The grib1 data that I used was generated with grib_api tool, and the
> original data was ECMWF forcast results.The original ECMWF data
enclude
> grib1 and grib2 format data in one file, and the grib_api tool can
> decompose that data to girb1 file and grib2 file. When I used the
wgrib
> tool to list the records(as following) of test file ,I can see the
field
> "TP".
>
>
> [/]#wgrib sfc_20160919_6_grid_simple.grib1
>
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> 0-6hr:NAve=0
>
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> 3-6hr:NAve=0
>
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
6:35794866:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
7:42020070:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
8:48245274:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
9:54470478:d=16091900:var216:kpds5=216:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
10:62770722:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
11:71070966:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
12:77296170:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
13:85596414:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
14:93896658:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
15:102196902:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
16:108422106:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> 3-6hr:NAve=0
>
17:114647310:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
18:118797474:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
19:127097718:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> 3-6hr:NAve=0
>
20:135397962:d=16091900:MN2T6:kpds5=122:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> 0-6hr:NAve=0
>
21:141623166:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
22:147848370:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
23:154073574:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
24:160298778:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
> 25:166523982:d=16091900:var247:kpds5=247:kpds6=1:
> kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
> 26:174824226:d=16091900:var246:kpds5=246:kpds6=1:
> kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
>
27:183124470:d=16091900:MX2T6:kpds5=121:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> 0-6hr:NAve=0
>
28:189349674:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
29:197649918:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
30:205950162:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
31:214250406:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
32:222550650:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
> [/]#wgrib sfc_20160919_9_grid_simple.grib1
>
1:0:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
2:8300244:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
3:16600488:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
4:22825692:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> 6-9hr:NAve=0
>
5:29050896:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
14:91302852:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
15:100121862:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
16:108422106:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
17:112572270:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
18:120872514:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
19:127097718:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
20:133322922:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
21:141623166:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
22:147848370:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
> 23:154073574:d=16091900:var216:kpds5=216:kpds6=1:
> kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
>
24:162373818:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
25:166523982:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
> 26:172749186:d=16091900:var246:kpds5=246:kpds6=1:
> kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
>
27:181049430:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
>
28:189349674:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> fcst:NAve=0
> 29:197649918:d=16091900:var247:kpds5=247:kpds6=1:
> kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
>
> Thanks
>
> yuanpu
>
> > -----原始邮件-----
> > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > 发送时间: 2016年9月22日 星期四
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> >
> > Hi Yuanpu.
> >
> > Thank you for creating a tar file with the data.
> >
> > Can you please tell me how the GRIB1 data that you are trying to
use with
> > MET was generated? Is it the output of the WRF Unified Post
Processor
> > (UPP)?
> >
> > It appear there may be an encoding problem that is causing this
error.
> >
> > Unfortunately, I am still unable to list the records for
> > sfc_20160919_6_grid_simple.grib1 using wgrib. I am able to use
wgrib to
> > view the records for sfc_20160919_9_grid_simple.grib1 but the same
issue
> > exists as when I use wgrib on the sfc_20160919_6_grid_simple.grib1
file:
> >
> > Undefined parameter table (center 98-0 table 228), using NCEP-opn
> > Len of grib message is inconsistent.
> >
> > missing end section
> > 46 f 46 f
> > 1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=
> 9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> > 2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:
> P1=6:P2=9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> > 4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:
> P1=6:P2=9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> > 5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
> > In that output, I do not see the field "TP" that you wish to get
the
> values
> > for. Since wgrib doesn't know anything about table number 228 for
center
> > 98, we can't trust the abbreviations wgrib spits out. It may be
that the
> > field exists, but due to an encoding issue, I am (and the MET code
is)
> > unable to see it.
> >
> > Thanks,
> > Julie
> >
> > On Tue, Sep 20, 2016 at 8:03 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
> > >
> > > Hi Julie
> > >
> > > Thank you for your reply.
> > >
> > > I used the pcp_combine commands like follows, and when I used
wgrib
> tool
> > > to list the recoreds, I can get the right records as follows.
So, I
> have
> > > put the test tar file on your ftp site at
incoming/irap/met_help/liu_
> data
> > > directory.
> > >
> > >
> > > COMMAND LINES:
> > >
> > > ../met-install/bin/pcp_combine -subtract \
> > >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
> > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
> \
> > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > > ${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_
> > > grid_simple.nc
> > >
> > > [/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
> > > 1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=
> 0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=
> 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > 0-6hr:NAve=0
> > > 3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:
> P1=3:P2=6:TimeU=1:sfc:valid
> > > 3-6hr:NAve=0
> > > 5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > .............................
> > >
> > > Thanks,
> > >
> > > yuanpu
> > > > -----原始邮件-----
> > > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > > 发送时间: 2016年9月20日 星期二
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > > >
> > > > Hi Yuanpu.
> > > >
> > > > Could you please send me the command the you are running?
> > > >
> > > > I think part of the problem with the GRIB data you sent to us
might
> be an
> > > > encoding issue. I tried to use the standard NCEP wgrib tool
to list
> out
> > > > the records present in the GRIB file, but it gives the
following
> error:
> > > >
> > > > jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
> > > >
> > > > Len of grib message is inconsistent.
> > > >
> > > >
> > > > missing end section
> > > > 46 f 46 f
> > > >
> > > > Are you able to list the records in that GRIB file using
wgrib? If
> so,
> > > > perhaps the file is being corrupted by the FTP transfer.
Would you
> mind
> > > > putting all the files you uploaded into a tarball and then
uploading
> > > that?
> > > > That should preserve the binary files and help avoid
corruption.
> Sorry
> > > for
> > > > the trouble.
> > > >
> > > > Thanks,
> > > > Julie Prestopnik
> > > >
> > > > On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > Tue Sep 20 01:27:52 2016: Request 77989 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=77989
> > > >
> > > > >
> > > > >
> > > > >
> > > > > Dear all,
> > > > >
> > > > > I received the following errors message while running the
> pcp_combine
> > > of
> > > > > met-5.0. I have put the test data on your ftp site at
> > > > > incoming/irap/met_help/liu_data directory.
> > > > >
> > > > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/
> > > > > 2016091900/sfc_20160919_6_grid_simple.grib1
> > > > > ERROR :
> > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777"
> not
> > > found
> > > > > in grib record
> > > > > ERROR :
> > > > > ERROR :
> > > > > ERROR : GribFile::seek_record(int) -> range check error
> > > > > ERROR :
> > > > > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > > > > ERROR :
> > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777"
> not
> > > found
> > > > > in grib record
> > > > > ERROR :
> > > > > WARNING:
> > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found for
> > > > > VarInfo "TP/Z0" in GRIB file "/data1/home/liuyp/DATA/ECMWF_
> > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > > > > WARNING:
> > > > > ERROR :
> > > > > ERROR : get_field() -> can't get data plane from file
> > > > > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > > > > 20160919_9_grid_simple.grib1"
> > > > > 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: Thu Sep 22 21:35:07 2016
Hi John,
I am sad that when I test the data according to your suggestion, and
used the new version met-5.2, the erros still exist like follows:
DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
ERROR :
ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found in grib record
ERROR :
ERROR :
ERROR : GribFile::seek_record(int) -> range check error
ERROR :
While the commands lines in met-5.2 like following:
../met-install/bin/pcp_combine -subtract \
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
06 \
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_3_grid_simple.grib1
03 \
./20160919.nc
And the commands lines in met-5.0 like this followsing:
../met-install/bin/pcp_combine -subtract \
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
\
'GRIB1_code=228;level="Z0";GRIB1_ptv=128;' \
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_3_grid_simple.grib1
\
'GRIB1_code=228;level="Z0";GRIB1_ptv=128;' \
./20160919.nc
Thanks,
yuanpu
> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2016年9月22日 星期四
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
>
> Yuanpu,
>
> Great, if wgrib is able to read your GRIB1 files then hopefully the
MET
> tools will be able to as well. However there are some issues to
consider...
>
> (1) The default settings for pcp_combine expect accumulated
precipitation
> to be stored in GRIB code 61 (for APCP). Your total precip (TP)
variable
> likely is not stored in GRIB code 61. Based on the wgrib output you
sent
> it's stored in GRIB code 228 (i.e. kpds5=228).
>
> (2) In your original email, you mentioned MET version 5.0. The MET
version
> 5.2 release included some pretty major updates to GRIB tables... so
using
> 5.0 versus 5.2 will have a pretty big impact. The GRIB table
updates in
> 5.2 are described here:
>
>
http://www.dtcenter.org/met/users/support/release_notes/METv5.2_release_notes.php
>
> The GRIB tables in MET 5.2 include information about center 98 and
table
> number 228, whereas MET 5.0 is limited to NCEP's GRIB table usage.
>
> Using 5.0, setting 'name="TP";' will likely not work while it may
work in
> 5.2.
>
> Regardless, instead of specifying the name, let's try explicitly
specifying
> the GRIB code and parameter table version number. Does the
following
> command work?
>
> ../met-install/bin/pcp_combine -subtract \
> ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1
\
> 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
> ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
> 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
>
${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_grid_simple.
> nc
>
> Note that it's just a coincidence that the GRIB code of 228 for
total
> precip and the parameter table version number of 228 just happen to
the be
> the same value.
>
> Please let me know if that works.
>
> I also find it helpful to visualize data like this using the plot
data
> plane tool:
>
> ../met-install/bin/plot_data_plane \
> ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1
\
> plot_TP.ps \
> 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";'
>
> And then open the plot_TP.ps PostScript file.
>
> Thanks,
> John Halley Gotway
>
>
> On Wed, Sep 21, 2016 at 7:17 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
> >
> > Hi Julie,
> >
> > Thanks for your help!
> >
> > The grib1 data that I used was generated with grib_api tool, and
the
> > original data was ECMWF forcast results.The original ECMWF data
enclude
> > grib1 and grib2 format data in one file, and the grib_api tool can
> > decompose that data to girb1 file and grib2 file. When I used the
wgrib
> > tool to list the records(as following) of test file ,I can see the
field
> > "TP".
> >
> >
> > [/]#wgrib sfc_20160919_6_grid_simple.grib1
> >
1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> > 0-6hr:NAve=0
> >
3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> > 3-6hr:NAve=0
> >
5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
6:35794866:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
7:42020070:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
8:48245274:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
9:54470478:d=16091900:var216:kpds5=216:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
10:62770722:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
11:71070966:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
12:77296170:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
13:85596414:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
14:93896658:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
15:102196902:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
16:108422106:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> > 3-6hr:NAve=0
> >
17:114647310:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
18:118797474:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
19:127097718:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=3:P2=6:TimeU=1:sfc:valid
> > 3-6hr:NAve=0
> >
20:135397962:d=16091900:MN2T6:kpds5=122:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> > 0-6hr:NAve=0
> >
21:141623166:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
22:147848370:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
23:154073574:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
24:160298778:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> > 25:166523982:d=16091900:var247:kpds5=247:kpds6=1:
> > kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
> > 26:174824226:d=16091900:var246:kpds5=246:kpds6=1:
> > kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
> >
27:183124470:d=16091900:MX2T6:kpds5=121:kpds6=1:kpds7=0:TR=2:P1=0:P2=6:TimeU=1:sfc:valid
> > 0-6hr:NAve=0
> >
28:189349674:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
29:197649918:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
30:205950162:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
31:214250406:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> >
32:222550650:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> > fcst:NAve=0
> > [/]#wgrib sfc_20160919_9_grid_simple.grib1
> >
1:0:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> >
2:8300244:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
3:16600488:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> >
4:22825692:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
> > 6-9hr:NAve=0
> >
5:29050896:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
14:91302852:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
15:100121862:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
16:108422106:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
17:112572270:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
18:120872514:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
19:127097718:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
20:133322922:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
21:141623166:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
22:147848370:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 23:154073574:d=16091900:var216:kpds5=216:kpds6=1:
> > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> >
24:162373818:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
25:166523982:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 26:172749186:d=16091900:var246:kpds5=246:kpds6=1:
> > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> >
27:181049430:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> >
28:189349674:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
> > fcst:NAve=0
> > 29:197649918:d=16091900:var247:kpds5=247:kpds6=1:
> > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> >
> > Thanks
> >
> > yuanpu
> >
> > > -----原始邮件-----
> > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > 发送时间: 2016年9月22日 星期四
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > >
> > > Hi Yuanpu.
> > >
> > > Thank you for creating a tar file with the data.
> > >
> > > Can you please tell me how the GRIB1 data that you are trying to
use with
> > > MET was generated? Is it the output of the WRF Unified Post
Processor
> > > (UPP)?
> > >
> > > It appear there may be an encoding problem that is causing this
error.
> > >
> > > Unfortunately, I am still unable to list the records for
> > > sfc_20160919_6_grid_simple.grib1 using wgrib. I am able to use
wgrib to
> > > view the records for sfc_20160919_9_grid_simple.grib1 but the
same issue
> > > exists as when I use wgrib on the
sfc_20160919_6_grid_simple.grib1 file:
> > >
> > > Undefined parameter table (center 98-0 table 228), using NCEP-
opn
> > > Len of grib message is inconsistent.
> > >
> > > missing end section
> > > 46 f 46 f
> > > 1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=
> > 9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:
> > P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:
> > P1=6:P2=9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:
> > P1=6:P2=9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:
> > P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=
> > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:
> > P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > >
> > > In that output, I do not see the field "TP" that you wish to get
the
> > values
> > > for. Since wgrib doesn't know anything about table number 228
for center
> > > 98, we can't trust the abbreviations wgrib spits out. It may be
that the
> > > field exists, but due to an encoding issue, I am (and the MET
code is)
> > > unable to see it.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Tue, Sep 20, 2016 at 8:03 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989
>
> > > >
> > > > Hi Julie
> > > >
> > > > Thank you for your reply.
> > > >
> > > > I used the pcp_combine commands like follows, and when I used
wgrib
> > tool
> > > > to list the recoreds, I can get the right records as follows.
So, I
> > have
> > > > put the test tar file on your ftp site at
incoming/irap/met_help/liu_
> > data
> > > > directory.
> > > >
> > > >
> > > > COMMAND LINES:
> > > >
> > > > ../met-install/bin/pcp_combine -subtract \
> > > >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
> > > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > > >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
> > \
> > > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > > >
${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_
> > > > grid_simple.nc
> > > >
> > > > [/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
> > > > 1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=
> > 0:TimeU=1:sfc:6hr
> > > > fcst:NAve=0
> > > > 2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=
> > 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > > 0-6hr:NAve=0
> > > > 3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> > 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > > fcst:NAve=0
> > > > 4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:
> > P1=3:P2=6:TimeU=1:sfc:valid
> > > > 3-6hr:NAve=0
> > > > 5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> > 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > > fcst:NAve=0
> > > > .............................
> > > >
> > > > Thanks,
> > > >
> > > > yuanpu
> > > > > -----原始邮件-----
> > > > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > > > 发送时间: 2016年9月20日 星期二
> > > > > 收件人: lyp at lzb.ac.cn
> > > > > 抄送:
> > > > > 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > > > >
> > > > > Hi Yuanpu.
> > > > >
> > > > > Could you please send me the command the you are running?
> > > > >
> > > > > I think part of the problem with the GRIB data you sent to
us might
> > be an
> > > > > encoding issue. I tried to use the standard NCEP wgrib tool
to list
> > out
> > > > > the records present in the GRIB file, but it gives the
following
> > error:
> > > > >
> > > > > jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
> > > > >
> > > > > Len of grib message is inconsistent.
> > > > >
> > > > >
> > > > > missing end section
> > > > > 46 f 46 f
> > > > >
> > > > > Are you able to list the records in that GRIB file using
wgrib? If
> > so,
> > > > > perhaps the file is being corrupted by the FTP transfer.
Would you
> > mind
> > > > > putting all the files you uploaded into a tarball and then
uploading
> > > > that?
> > > > > That should preserve the binary files and help avoid
corruption.
> > Sorry
> > > > for
> > > > > the trouble.
> > > > >
> > > > > Thanks,
> > > > > Julie Prestopnik
> > > > >
> > > > > On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT
<met_help at ucar.edu>
> > wrote:
> > > > >
> > > > > >
> > > > > > Tue Sep 20 01:27:52 2016: Request 77989 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=77989
> > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > Dear all,
> > > > > >
> > > > > > I received the following errors message while running the
> > pcp_combine
> > > > of
> > > > > > met-5.0. I have put the test data on your ftp site at
> > > > > > incoming/irap/met_help/liu_data directory.
> > > > > >
> > > > > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/
> > > > > > 2016091900/sfc_20160919_6_grid_simple.grib1
> > > > > > ERROR :
> > > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777"
> > not
> > > > found
> > > > > > in grib record
> > > > > > ERROR :
> > > > > > ERROR :
> > > > > > ERROR : GribFile::seek_record(int) -> range check error
> > > > > > ERROR :
> > > > > > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_
> > > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > > > > > ERROR :
> > > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
"7777"
> > not
> > > > found
> > > > > > in grib record
> > > > > > ERROR :
> > > > > > WARNING:
> > > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact match
found for
> > > > > > VarInfo "TP/Z0" in GRIB file
"/data1/home/liuyp/DATA/ECMWF_
> > > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > > > > > WARNING:
> > > > > > ERROR :
> > > > > > ERROR : get_field() -> can't get data plane from file
> > > > > > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > > > > > 20160919_9_grid_simple.grib1"
> > > > > > 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: John Halley Gotway
Time: Fri Sep 23 10:30:14 2016
Yuanpu,
That matches the behavior we saw when trying to run the sample GRIB
files
you sent us through MET. As Julie pointed out, the version of wgrib
we
have here wasn't able to list any of the records in the 6-hour file
and
only 13 from the 9-hour file, as shown at the bottom of this message.
However, in order to demonstrate that MET can in fact read some of
your
data, you could run the following command:
met-5.2/bin/plot_data_plane \
sfc_20160919_9_grid_simple.grib1 plot_gc44.ps \
'GRIB1_code = 44; GRIB1_ptv = 228; level="L0";'
I converted the resulting output PostScript file to png and have
attached
it to this message.
Unfortunately, I do think there's a problem with the way your GRIB
file is
encoded, and we are not well-equipped to assist you with that piece.
Sorry we can't be of more help.
John Halley Gotway
[johnhg at number5]% wgrib sfc_20160919_6_grid_simple.grib1
Len of grib message is inconsistent.
missing end section
46 f 46 f
[johnhg at number5]% wgrib sfc_20160919_9_grid_simple.grib1
Undefined parameter table (center 98-0 table 228), using NCEP-opn
1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:P1=6:P2=9:TimeU=1:sfc:valid
6-9hr:NAve=0
5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr
fcst:NAve=0
Len of grib message is inconsistent.
missing end section
46 f 46 f
On Thu, Sep 22, 2016 at 9:35 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
>
> Hi John,
>
> I am sad that when I test the data according to your suggestion, and
used
> the new version met-5.2, the erros still exist like follows:
> DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> 2016091900/sfc_20160919_6_grid_simple.grib1
> ERROR :
> ERROR : GribFile::read_record(GribRecord &) -> trailing "7777" not
found
> in grib record
> ERROR :
> ERROR :
> ERROR : GribFile::seek_record(int) -> range check error
> ERROR :
>
>
> While the commands lines in met-5.2 like following:
>
> ../met-install/bin/pcp_combine -subtract \
>
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
> 06 \
>
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_3_grid_simple.grib1
> 03 \
> ./20160919.nc
>
>
> And the commands lines in met-5.0 like this followsing:
>
> ../met-install/bin/pcp_combine -subtract \
>
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_6_grid_simple.grib1
> \
> 'GRIB1_code=228;level="Z0";GRIB1_ptv=128;' \
>
/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_20160919_3_grid_simple.grib1
> \
> 'GRIB1_code=228;level="Z0";GRIB1_ptv=128;' \
> ./20160919.nc
>
> Thanks,
>
> yuanpu
> > -----原始邮件-----
> > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > 发送时间: 2016年9月22日 星期四
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> >
> > Yuanpu,
> >
> > Great, if wgrib is able to read your GRIB1 files then hopefully
the MET
> > tools will be able to as well. However there are some issues to
> consider...
> >
> > (1) The default settings for pcp_combine expect accumulated
precipitation
> > to be stored in GRIB code 61 (for APCP). Your total precip (TP)
variable
> > likely is not stored in GRIB code 61. Based on the wgrib output
you sent
> > it's stored in GRIB code 228 (i.e. kpds5=228).
> >
> > (2) In your original email, you mentioned MET version 5.0. The
MET
> version
> > 5.2 release included some pretty major updates to GRIB tables...
so using
> > 5.0 versus 5.2 will have a pretty big impact. The GRIB table
updates in
> > 5.2 are described here:
> >
> > http://www.dtcenter.org/met/users/support/release_notes/
> METv5.2_release_notes.php
> >
> > The GRIB tables in MET 5.2 include information about center 98 and
table
> > number 228, whereas MET 5.0 is limited to NCEP's GRIB table usage.
> >
> > Using 5.0, setting 'name="TP";' will likely not work while it may
work in
> > 5.2.
> >
> > Regardless, instead of specifying the name, let's try explicitly
> specifying
> > the GRIB code and parameter table version number. Does the
following
> > command work?
> >
> > ../met-install/bin/pcp_combine -subtract \
> >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
> > 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
> >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_simple.grib1
\
> > 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";' \
> > ${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${
> datebefore}_${hh}_grid_simple.
> > nc
> >
> > Note that it's just a coincidence that the GRIB code of 228 for
total
> > precip and the parameter table version number of 228 just happen
to the
> be
> > the same value.
> >
> > Please let me know if that works.
> >
> > I also find it helpful to visualize data like this using the plot
data
> > plane tool:
> >
> > ../met-install/bin/plot_data_plane \
> >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1 \
> > plot_TP.ps \
> > 'GRIB1_code=228; GRIB1_ptv=228; level="Z0";'
> >
> > And then open the plot_TP.ps PostScript file.
> >
> > Thanks,
> > John Halley Gotway
> >
> >
> > On Wed, Sep 21, 2016 at 7:17 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
> > >
> > > Hi Julie,
> > >
> > > Thanks for your help!
> > >
> > > The grib1 data that I used was generated with grib_api tool, and
the
> > > original data was ECMWF forcast results.The original ECMWF data
enclude
> > > grib1 and grib2 format data in one file, and the grib_api tool
can
> > > decompose that data to girb1 file and grib2 file. When I used
the wgrib
> > > tool to list the records(as following) of test file ,I can see
the
> field
> > > "TP".
> > >
> > >
> > > [/]#wgrib sfc_20160919_6_grid_simple.grib1
> > > 1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=
> 0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=
> 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > 0-6hr:NAve=0
> > > 3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:
> P1=3:P2=6:TimeU=1:sfc:valid
> > > 3-6hr:NAve=0
> > > 5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 6:35794866:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 7:42020070:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 8:48245274:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 9:54470478:d=16091900:var216:kpds5=216:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 10:62770722:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 11:71070966:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 12:77296170:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 13:85596414:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:
> P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 14:93896658:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:
> P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 15:102196902:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 16:108422106:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:
> P1=3:P2=6:TimeU=1:sfc:valid
> > > 3-6hr:NAve=0
> > > 17:114647310:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 18:118797474:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:
> P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 19:127097718:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:
> P1=3:P2=6:TimeU=1:sfc:valid
> > > 3-6hr:NAve=0
> > > 20:135397962:d=16091900:MN2T6:kpds5=122:kpds6=1:kpds7=0:TR=
> 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > 0-6hr:NAve=0
> > > 21:141623166:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 22:147848370:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 23:154073574:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 24:160298778:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 25:166523982:d=16091900:var247:kpds5=247:kpds6=1:
> > > kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
> > > 26:174824226:d=16091900:var246:kpds5=246:kpds6=1:
> > > kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
> > > 27:183124470:d=16091900:MX2T6:kpds5=121:kpds6=1:kpds7=0:TR=
> 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > 0-6hr:NAve=0
> > > 28:189349674:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:
> P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 29:197649918:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:
> P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 30:205950162:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 31:214250406:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > 32:222550650:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=
> 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > fcst:NAve=0
> > > [/]#wgrib sfc_20160919_9_grid_simple.grib1
> > > 1:0:d=16091900:var28:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=
> 9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 2:8300244:d=16091900:var24:kpds5=24:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 3:16600488:d=16091900:var26:kpds5=26:kpds6=1:kpds7=0:TR=2:
> P1=6:P2=9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 4:22825692:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:
> P1=6:P2=9:TimeU=1:sfc:valid
> > > 6-9hr:NAve=0
> > > 5:29050896:d=16091900:var44:kpds5=44:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 14:91302852:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 15:100121862:d=16091900:SD:kpds5=141:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 16:108422106:d=16091900:LCC:kpds5=186:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 17:112572270:d=16091900:CAPE:kpds5=59:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 18:120872514:d=16091900:10V:kpds5=166:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 19:127097718:d=16091900:2T:kpds5=167:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 20:133322922:d=16091900:TP:kpds5=228:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 21:141623166:d=16091900:10U:kpds5=165:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 22:147848370:d=16091900:SKT:kpds5=235:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 23:154073574:d=16091900:var216:kpds5=216:kpds6=1:
> > > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> > > 24:162373818:d=16091900:FAL:kpds5=243:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 25:166523982:d=16091900:SP:kpds5=134:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 26:172749186:d=16091900:var246:kpds5=246:kpds6=1:
> > > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> > > 27:181049430:d=16091900:LSP:kpds5=142:kpds6=1:kpds7=0:TR=
> 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 28:189349674:d=16091900:RSN:kpds5=33:kpds6=1:kpds7=0:TR=0:
> P1=9:P2=0:TimeU=1:sfc:9hr
> > > fcst:NAve=0
> > > 29:197649918:d=16091900:var247:kpds5=247:kpds6=1:
> > > kpds7=0:TR=0:P1=9:P2=0:TimeU=1:sfc:9hr fcst:NAve=0
> > >
> > > Thanks
> > >
> > > yuanpu
> > >
> > > > -----原始邮件-----
> > > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > > 发送时间: 2016年9月22日 星期四
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > > >
> > > > Hi Yuanpu.
> > > >
> > > > Thank you for creating a tar file with the data.
> > > >
> > > > Can you please tell me how the GRIB1 data that you are trying
to use
> with
> > > > MET was generated? Is it the output of the WRF Unified Post
Processor
> > > > (UPP)?
> > > >
> > > > It appear there may be an encoding problem that is causing
this
> error.
> > > >
> > > > Unfortunately, I am still unable to list the records for
> > > > sfc_20160919_6_grid_simple.grib1 using wgrib. I am able to
use
> wgrib to
> > > > view the records for sfc_20160919_9_grid_simple.grib1 but the
same
> issue
> > > > exists as when I use wgrib on the
sfc_20160919_6_grid_simple.grib1
> file:
> > > >
> > > > Undefined parameter table (center 98-0 table 228), using NCEP-
opn
> > > > Len of grib message is inconsistent.
> > > >
> > > > missing end section
> > > > 46 f 46 f
> > > > 1:0:d=16091900:WVSP1:kpds5=28:kpds6=1:kpds7=0:TR=2:P1=6:P2=
> > > 9:TimeU=1:sfc:valid
> > > > 6-9hr:NAve=0
> > > > 2:8300244:d=16091900:PLI:kpds5=24:kpds6=1:kpds7=0:TR=0:
> > > P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 3:16600488:d=16091900:PRESA:kpds5=26:kpds6=1:kpds7=0:TR=2:
> > > P1=6:P2=9:TimeU=1:sfc:valid
> > > > 6-9hr:NAve=0
> > > > 4:22825692:d=16091900:GPA:kpds5=27:kpds6=1:kpds7=0:TR=2:
> > > P1=6:P2=9:TimeU=1:sfc:valid
> > > > 6-9hr:NAve=0
> > > > 5:29050896:d=16091900:RELD:kpds5=44:kpds6=1:kpds7=0:TR=0:
> > > P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 6:37351140:d=16091900:TCWV:kpds5=137:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 7:43576344:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 8:51876588:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 9:56026752:d=16091900:2D:kpds5=168:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 10:62251956:d=16091900:MSL:kpds5=151:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 11:68477160:d=16091900:SF:kpds5=144:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 12:76777404:d=16091900:TCW:kpds5=136:kpds6=1:kpds7=0:TR=
> > > 0:P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > > 13:83002608:d=16091900:VIS:kpds5=20:kpds6=1:kpds7=0:TR=0:
> > > P1=9:P2=0:TimeU=1:sfc:9hr
> > > > fcst:NAve=0
> > > >
> > > > In that output, I do not see the field "TP" that you wish to
get the
> > > values
> > > > for. Since wgrib doesn't know anything about table number 228
for
> center
> > > > 98, we can't trust the abbreviations wgrib spits out. It may
be
> that the
> > > > field exists, but due to an encoding issue, I am (and the MET
code
> is)
> > > > unable to see it.
> > > >
> > > > Thanks,
> > > > Julie
> > > >
> > > > On Tue, Sep 20, 2016 at 8:03 PM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77989 >
> > > > >
> > > > > Hi Julie
> > > > >
> > > > > Thank you for your reply.
> > > > >
> > > > > I used the pcp_combine commands like follows, and when I
used wgrib
> > > tool
> > > > > to list the recoreds, I can get the right records as
follows. So, I
> > > have
> > > > > put the test tar file on your ftp site at
> incoming/irap/met_help/liu_
> > > data
> > > > > directory.
> > > > >
> > > > >
> > > > > COMMAND LINES:
> > > > >
> > > > > ../met-install/bin/pcp_combine -subtract \
> > > > >
${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh}_grid_simple.grib1
> \
> > > > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > > > > ${TP_dir}/${datebefore}00/sfc_${datebefore}_${hh1}_grid_
> simple.grib1
> > > \
> > > > > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > > > >
${pcp_dir}/output_EC/${datebefore}00/3hr_TP_${datebefore}_${hh}_
> > > > > grid_simple.nc
> > > > >
> > > > > [/]#wgrib sfc_20160919_6_grid_simple.grib1 |more
> > > > > 1:0:d=16091900:SSTK:kpds5=34:kpds6=1:kpds7=0:TR=0:P1=6:P2=
> > > 0:TimeU=1:sfc:6hr
> > > > > fcst:NAve=0
> > > > > 2:8819010:d=16091900:10FG6:kpds5=123:kpds6=1:kpds7=0:TR=
> > > 2:P1=0:P2=6:TimeU=1:sfc:valid
> > > > > 0-6hr:NAve=0
> > > > > 3:17119254:d=16091900:CP:kpds5=143:kpds6=1:kpds7=0:TR=
> > > 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > > > fcst:NAve=0
> > > > > 4:25419498:d=16091900:var27:kpds5=27:kpds6=1:kpds7=0:TR=2:
> > > P1=3:P2=6:TimeU=1:sfc:valid
> > > > > 3-6hr:NAve=0
> > > > > 5:31644702:d=16091900:TCC:kpds5=164:kpds6=1:kpds7=0:TR=
> > > 0:P1=6:P2=0:TimeU=1:sfc:6hr
> > > > > fcst:NAve=0
> > > > > .............................
> > > > >
> > > > > Thanks,
> > > > >
> > > > > yuanpu
> > > > > > -----原始邮件-----
> > > > > > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > > > > > 发送时间: 2016年9月20日 星期二
> > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > 抄送:
> > > > > > 主题: Re: [rt.rap.ucar.edu #77989] About pcp_combine
> > > > > >
> > > > > > Hi Yuanpu.
> > > > > >
> > > > > > Could you please send me the command the you are running?
> > > > > >
> > > > > > I think part of the problem with the GRIB data you sent to
us
> might
> > > be an
> > > > > > encoding issue. I tried to use the standard NCEP wgrib
tool to
> list
> > > out
> > > > > > the records present in the GRIB file, but it gives the
following
> > > error:
> > > > > >
> > > > > > jpresto% wgrib sfc_20160919_6_grid_simple.grib1 | more
> > > > > >
> > > > > > Len of grib message is inconsistent.
> > > > > >
> > > > > >
> > > > > > missing end section
> > > > > > 46 f 46 f
> > > > > >
> > > > > > Are you able to list the records in that GRIB file using
wgrib?
> If
> > > so,
> > > > > > perhaps the file is being corrupted by the FTP transfer.
Would
> you
> > > mind
> > > > > > putting all the files you uploaded into a tarball and then
> uploading
> > > > > that?
> > > > > > That should preserve the binary files and help avoid
corruption.
> > > Sorry
> > > > > for
> > > > > > the trouble.
> > > > > >
> > > > > > Thanks,
> > > > > > Julie Prestopnik
> > > > > >
> > > > > > On Tue, Sep 20, 2016 at 1:27 AM, 柳媛普 via RT
<met_help at ucar.edu>
> > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Tue Sep 20 01:27:52 2016: Request 77989 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=77989
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Dear all,
> > > > > > >
> > > > > > > I received the following errors message while running
the
> > > pcp_combine
> > > > > of
> > > > > > > met-5.0. I have put the test data on your ftp site at
> > > > > > > incoming/irap/met_help/liu_data directory.
> > > > > > >
> > > > > > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_
> simple/
> > > > > > > 2016091900/sfc_20160919_6_grid_simple.grib1
> > > > > > > ERROR :
> > > > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
> "7777"
> > > not
> > > > > found
> > > > > > > in grib record
> > > > > > > ERROR :
> > > > > > > ERROR :
> > > > > > > ERROR : GribFile::seek_record(int) -> range check error
> > > > > > > ERROR :
> > > > > > > DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_
> > > > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1
> > > > > > > ERROR :
> > > > > > > ERROR : GribFile::read_record(GribRecord &) -> trailing
> "7777"
> > > not
> > > > > found
> > > > > > > in grib record
> > > > > > > ERROR :
> > > > > > > WARNING:
> > > > > > > WARNING: MetGrib1DataFile::data_plane() -> No exact
match
> found for
> > > > > > > VarInfo "TP/Z0" in GRIB file
"/data1/home/liuyp/DATA/ECMWF_
> > > > > > > simple/2016091900/sfc_20160919_9_grid_simple.grib1".
> > > > > > > WARNING:
> > > > > > > ERROR :
> > > > > > > ERROR : get_field() -> can't get data plane from file
> > > > > > > "/data1/home/liuyp/DATA/ECMWF_simple/2016091900/sfc_
> > > > > > > 20160919_9_grid_simple.grib1"
> > > > > > > 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
> > > >
> > >
> > >
> > >
> > >
> > >
> >
>
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list