[Met_help] [rt.rap.ucar.edu #73244] History for About pcp_combine
John Halley Gotway via RT
met_help at ucar.edu
Wed Sep 2 09:53:51 MDT 2015
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear all,
I saw the following errors when I run pcp_combine tool of met-5.0 to processed precipitation. The output nc file likes TP_Z0 =
_, _, _, _, _, . While the commands likes follwed. The name of "TP" indicat total pricipitation, and it's unit is m. I also has post the grib data on your ftp site in incoming/irap/met_help/liu_data directory.
Kindly advise.
Thanks in advance.
Yuanpu
COMMANDS:
../met-install/bin/pcp_combine -subtract \
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1 \
'name="TP";level="Z0";GRIB1_ptv=128;' \
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_48_grid_simple.grib1 \
'name="TP";level="Z0";GRIB1_ptv=128;' \
./72.nc
ERRORS:
DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1
DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_69_grid_simple.grib1
DEBUG 1: Writing output file: ./72.nc
NetCDF: Start+count exceeds dimension bound
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Mon Aug 31 21:51:36 2015
Dear all,
I saw the following errors when I run pcp_combine tool of met-5.0 to
processed precipitation. The output nc file likes
TP_Z0 =
_, _, _, _, _, .
While the commands likes follwed. The name of "TP" indicat total
pricipitation, and it's unit is m. I also has post the grib data on
your ftp site in incoming/irap/met_help/liu_data directory.
Kindly advise.
Thanks in advance.
Yuanpu
COMMANDS:
../met-install/bin/pcp_combine -subtract \
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1
\
'name="TP";level="Z0";GRIB1_ptv=128;' \
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_48_grid_simple.grib1
\
'name="TP";level="Z0";GRIB1_ptv=128;' \
./72.nc
ERRORS:
DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1
DEBUG 1: Reading input file:
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_69_grid_simple.grib1
DEBUG 1: Writing output file: ./72.nc
NetCDF: Start+count exceeds dimension bound
------------------------------------------------
Subject: About pcp_combine
From: John Halley Gotway
Time: Tue Sep 01 10:36:58 2015
MET assumes that all records in a GRIB file are defined on the same
grid.
This assumption is violated in the data you sent us, and that is the
source
of the problem.
MET extracts the grid information from the first record of that GRIB
file,
which is defined on a 361 x 281 grid:
wgrib -V -d 1 20150825_48_grid_simple.grib1
However, the TP data in record number 165 is defined on a 2880 x 1441
grid:
wgrib -V -d 165 20150825_48_grid_simple.grib1
pcp_combine has only allocated space to store data for that 361 x 281
grid,
and we get that error when we trying to write 2880 x 1441 values into
that
array.
This is the first example I've seen of people mixing grids in a GRIB
file.
Is this common practice at your institution? If so, we could look at
enhancing MET to better accommodate that in the future.
To make sure there aren't other problems, I ran the following commands
to
strip out those single GRIB records and run them through pcp_combine:
# Extract TP records
wgrib 20150825_48_grid_simple.grib1 | grep TP | wgrib -i -grib -o
20150825_48_grid_simple_TP.grib1 20150825_48_grid_simple.grib1
wgrib 20150825_72_grid_simple.grib1 | grep TP | wgrib -i -grib -o
20150825_72_grid_simple_TP.grib1 20150825_72_grid_simple.grib1
# Run pcp_combine subraction
met-5.0/bin/pcp_combine -subtract \
20150825_72_grid_simple_TP.grib1 \
'name="TP"; level="L0"; GRIB1_ptv=128;' \
20150825_48_grid_simple_TP.grib1 \
'name="TP"; level="L0"; GRIB1_ptv=128;' \
diff_72_48.nc
# Plot the result
met-5.0/bin/plot_data_plane diff_72_48.nc diff_72_48.ps 'name="TP_L0";
level="(*,*)";'
The resulting image is attached.
Hopefully that helps.
Thanks,
John Halley Gotway
On Mon, Aug 31, 2015 at 9:51 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73244 >
>
>
>
>
> Dear all,
>
> I saw the following errors when I run pcp_combine tool of met-5.0 to
> processed precipitation. The output nc file likes
>
> TP_Z0 =
> _, _, _, _, _, .
>
> While the commands likes follwed. The name of "TP" indicat total
> pricipitation, and it's unit is m. I also has post the grib data on
your
> ftp site in incoming/irap/met_help/liu_data directory.
>
> Kindly advise.
>
> Thanks in advance.
>
> Yuanpu
>
>
> COMMANDS:
>
> ../met-install/bin/pcp_combine -subtract \
>
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1
> \
> 'name="TP";level="Z0";GRIB1_ptv=128;' \
>
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_48_grid_simple.grib1
> \
> 'name="TP";level="Z0";GRIB1_ptv=128;' \
> ./72.nc
>
> ERRORS:
>
> DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> 2015082500/20150825_72_grid_simple.grib1
> DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> 2015082500/20150825_69_grid_simple.grib1
> DEBUG 1: Writing output file: ./72.nc
> NetCDF: Start+count exceeds dimension bound
>
>
>
>
------------------------------------------------
Subject: About pcp_combine
From: 柳媛普
Time: Tue Sep 01 21:04:41 2015
Hello John,
Thanks for your help. I appreciate it. I will extract the TP record
with wgrib first, and then run the pcp_combine tool. There will be no
errors.
The grib1 data were decoded from ECMWF complex packing data with
grib_api tool. There would be seperate to two files of different kinds
of grids or mixing grids file.
Yuanpu
> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年9月2日 星期三
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #73244] About pcp_combine
>
> MET assumes that all records in a GRIB file are defined on the same
grid.
> This assumption is violated in the data you sent us, and that is the
source
> of the problem.
>
> MET extracts the grid information from the first record of that GRIB
file,
> which is defined on a 361 x 281 grid:
> wgrib -V -d 1 20150825_48_grid_simple.grib1
>
> However, the TP data in record number 165 is defined on a 2880 x
1441 grid:
> wgrib -V -d 165 20150825_48_grid_simple.grib1
>
> pcp_combine has only allocated space to store data for that 361 x
281 grid,
> and we get that error when we trying to write 2880 x 1441 values
into that
> array.
>
> This is the first example I've seen of people mixing grids in a GRIB
file.
> Is this common practice at your institution? If so, we could look
at
> enhancing MET to better accommodate that in the future.
>
> To make sure there aren't other problems, I ran the following
commands to
> strip out those single GRIB records and run them through
pcp_combine:
>
> # Extract TP records
> wgrib 20150825_48_grid_simple.grib1 | grep TP | wgrib -i -grib -o
> 20150825_48_grid_simple_TP.grib1 20150825_48_grid_simple.grib1
>
> wgrib 20150825_72_grid_simple.grib1 | grep TP | wgrib -i -grib -o
> 20150825_72_grid_simple_TP.grib1 20150825_72_grid_simple.grib1
>
> # Run pcp_combine subraction
> met-5.0/bin/pcp_combine -subtract \
> 20150825_72_grid_simple_TP.grib1 \
> 'name="TP"; level="L0"; GRIB1_ptv=128;' \
> 20150825_48_grid_simple_TP.grib1 \
> 'name="TP"; level="L0"; GRIB1_ptv=128;' \
> diff_72_48.nc
>
> # Plot the result
> met-5.0/bin/plot_data_plane diff_72_48.nc diff_72_48.ps
'name="TP_L0";
> level="(*,*)";'
>
> The resulting image is attached.
>
> Hopefully that helps.
>
> Thanks,
> John Halley Gotway
>
>
>
>
> On Mon, Aug 31, 2015 at 9:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73244 >
> >
> >
> >
> >
> > Dear all,
> >
> > I saw the following errors when I run pcp_combine tool of met-5.0
to
> > processed precipitation. The output nc file likes
> >
> > TP_Z0 =
> > _, _, _, _, _, .
> >
> > While the commands likes follwed. The name of "TP" indicat total
> > pricipitation, and it's unit is m. I also has post the grib data
on your
> > ftp site in incoming/irap/met_help/liu_data directory.
> >
> > Kindly advise.
> >
> > Thanks in advance.
> >
> > Yuanpu
> >
> >
> > COMMANDS:
> >
> > ../met-install/bin/pcp_combine -subtract \
> >
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_72_grid_simple.grib1
> > \
> > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> >
/data1/home/liuyp/DATA/ECMWF_simple/2015082500/20150825_48_grid_simple.grib1
> > \
> > 'name="TP";level="Z0";GRIB1_ptv=128;' \
> > ./72.nc
> >
> > ERRORS:
> >
> > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> > 2015082500/20150825_72_grid_simple.grib1
> > DEBUG 1: Reading input file: /data1/home/liuyp/DATA/ECMWF_simple/
> > 2015082500/20150825_69_grid_simple.grib1
> > DEBUG 1: Writing output file: ./72.nc
> > NetCDF: Start+count exceeds dimension bound
> >
> >
> >
> >
>
------------------------------------------------
More information about the Met_help
mailing list