[Met_help] [rt.rap.ucar.edu #81470] History for problem with pcp_combine
John Halley Gotway via RT
met_help at ucar.edu
Wed Sep 6 12:11:44 MDT 2017
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
To whom it may concern,
I had the following problem pcp_combine. I supplied two input grib files
from WPC and pcp_combine could not read the second one (see output below),
both files came from the ncep site at:
http://ftp.wpc.ncep.noaa.gov/5km_qpf/
Any assistance would be appreciated,
Brad Diehl
NWS/MDL
(g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add p06m_2017080300f006.grb 6
p06m_2017080312f006.grb 6 p06_add_out.nc
DEBUG 1: Reading input file: p06m_2017080300f006.grb
WARNING:
WARNING: grd_file_type() -> the file type indicated by the suffix
"FileType_Gb1" does not match the file type indicated by the data
"FileType_Gb2".
WARNING:
ERROR :
ERROR : GribFile::read() -> file read error ... requested -8 bytes, got -1
ERROR :
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Fri Aug 04 11:29:46 2017
Upon using degrib, noticed that files were in grib2 and give the files
the
correct extensions and did the following pcp_combine:
(g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080200f006.grb2 6
p06m_2017080212f006.grb2 6 p06_add_out.nc
DEBUG 1: Reading input file: p06m_2017080200f006.grb2
WARNING:
WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not currently
support Lambert Conformal grids where dx (5.07862) != dy (5.08049) and
may
produce unexpected results!
WARNING:
DEBUG 1: Reading input file: p06m_2017080212f006.grb2
WARNING:
WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not currently
support Lambert Conformal grids where dx (5.07862) != dy (5.08049) and
may
produce unexpected results!
WARNING:
DEBUG 1: Writing output file: p06_add_out.nc
So, I guess we can close this ticket, sorry for the
trouble........Brad
On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT
<met_help at ucar.edu>
wrote:
> NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17 AND
5/30/17.
> RESPONSES MAY BE DELAYED.
>
> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
> "problem with pcp_combine",
> a summary of which appears below.
>
> There is no need to reply to this message right now. Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #81470].
>
> Please include the string:
>
> [rt.rap.ucar.edu #81470]
>
> in the subject line of all future correspondence about this issue.
To do
> so,
> you may reply to this message.
>
> Thank you,
> met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> To whom it may concern,
>
> I had the following problem pcp_combine. I supplied two input grib
files
> from WPC and pcp_combine could not read the second one (see output
below),
> both files came from the ncep site at:
> http://ftp.wpc.ncep.noaa.gov/5km_qpf/
>
> Any assistance would be appreciated,
>
> Brad Diehl
> NWS/MDL
>
> (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080300f006.grb 6
> p06m_2017080312f006.grb 6 p06_add_out.nc
> DEBUG 1: Reading input file: p06m_2017080300f006.grb
> WARNING:
> WARNING: grd_file_type() -> the file type indicated by the suffix
> "FileType_Gb1" does not match the file type indicated by the data
> "FileType_Gb2".
> WARNING:
> ERROR :
> ERROR : GribFile::read() -> file read error ... requested -8 bytes,
got -1
> ERROR :
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: John Halley Gotway
Time: Fri Aug 04 11:44:00 2017
Hi Brad,
I was writing up a response when I saw your email...
I see that you're having trouble getting pcp_combine to run on some
GRIB
data. As you may have guessed based on the error message, MET is
complaining that you're processing GRIB2 files that are named with a
".grb"
suffix, which it assumes to be GRIB1.
When you pass MET a gridded data file, it needs to figure out how to
process it. The first step is looking at the filename suffix.
Suffixes
like .grb, .grib, and .GRIB are assumed to be GRIB1. Suffixes like
.grb2,
grib2, and .GRIB2 are assumed to be GRIB2. With your data, it's
assumed
your file is GRIB1 because of the .grb suffix. And that's the reason
for
the error.
To get it to work, you could either...
(1) Rename the .grb files to .grib2 ... or ...
(2) Use command line arguments to tell MET that it's really GRIB2.
Here's an example of running MET's plot_data_plane tool, passing to it
a
GRIB2 file that named .grb:
met-6.0/bin/plot_data_plane \
p06m_2017080200f006.grb tmp_p500.ps \
'name="TMP"; level="P500"; file_type = GRIB2;'
The "file_type" option explicitly tells MET to interpret this .grb
file as
a GRIB2 file.
John
On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>
> Upon using degrib, noticed that files were in grib2 and give the
files the
> correct extensions and did the following pcp_combine:
>
> (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080200f006.grb2 6
> p06m_2017080212f006.grb2 6 p06_add_out.nc
> DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> DEBUG 1: Writing output file: p06_add_out.nc
>
> So, I guess we can close this ticket, sorry for the
trouble........Brad
>
>
> On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> met_help at ucar.edu>
> wrote:
>
> > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17 AND
5/30/17.
> > RESPONSES MAY BE DELAYED.
> >
> > Greetings,
> >
> > This message has been automatically generated in response to the
> > creation of a trouble ticket regarding:
> > "problem with pcp_combine",
> > a summary of which appears below.
> >
> > There is no need to reply to this message right now. Your ticket
has
> been
> > assigned an ID of [rt.rap.ucar.edu #81470].
> >
> > Please include the string:
> >
> > [rt.rap.ucar.edu #81470]
> >
> > in the subject line of all future correspondence about this issue.
To do
> > so,
> > you may reply to this message.
> >
> > Thank you,
> > met_help at ucar.edu
> >
> > ------------------------------------------------------------
> -------------
> > To whom it may concern,
> >
> > I had the following problem pcp_combine. I supplied two input grib
files
> > from WPC and pcp_combine could not read the second one (see output
> below),
> > both files came from the ncep site at:
> > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> >
> > Any assistance would be appreciated,
> >
> > Brad Diehl
> > NWS/MDL
> >
> > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080300f006.grb
> 6
> > p06m_2017080312f006.grb 6 p06_add_out.nc
> > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > WARNING:
> > WARNING: grd_file_type() -> the file type indicated by the suffix
> > "FileType_Gb1" does not match the file type indicated by the data
> > "FileType_Gb2".
> > WARNING:
> > ERROR :
> > ERROR : GribFile::read() -> file read error ... requested -8
bytes, got
> -1
> > ERROR :
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Fri Aug 04 11:52:21 2017
John, thanks for your message and help, I'll keep that around for
reference.
Quick question, have you used matplotlib to plot results of MET or are
you
basically using the MET Viewer?
If you have used matplotlib, if you don't mind I have a couple of
questions, basically about the netcdf file and variables.
Thanks, Brad
On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Hi Brad,
>
> I was writing up a response when I saw your email...
>
> I see that you're having trouble getting pcp_combine to run on some
GRIB
> data. As you may have guessed based on the error message, MET is
> complaining that you're processing GRIB2 files that are named with a
".grb"
> suffix, which it assumes to be GRIB1.
>
> When you pass MET a gridded data file, it needs to figure out how to
> process it. The first step is looking at the filename suffix.
Suffixes
> like .grb, .grib, and .GRIB are assumed to be GRIB1. Suffixes like
.grb2,
> grib2, and .GRIB2 are assumed to be GRIB2. With your data, it's
assumed
> your file is GRIB1 because of the .grb suffix. And that's the
reason for
> the error.
>
> To get it to work, you could either...
> (1) Rename the .grb files to .grib2 ... or ...
> (2) Use command line arguments to tell MET that it's really GRIB2.
>
> Here's an example of running MET's plot_data_plane tool, passing to
it a
> GRIB2 file that named .grb:
>
> met-6.0/bin/plot_data_plane \
> p06m_2017080200f006.grb tmp_p500.ps \
> 'name="TMP"; level="P500"; file_type = GRIB2;'
>
> The "file_type" option explicitly tells MET to interpret this .grb
file as
> a GRIB2 file.
>
> John
>
> On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> >
> > Upon using degrib, noticed that files were in grib2 and give the
files
> the
> > correct extensions and did the following pcp_combine:
> >
> > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> p06m_2017080200f006.grb2 6
> > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > WARNING:
> > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> > support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and
> may
> > produce unexpected results!
> > WARNING:
> > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > WARNING:
> > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> > support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and
> may
> > produce unexpected results!
> > WARNING:
> > DEBUG 1: Writing output file: p06_add_out.nc
> >
> > So, I guess we can close this ticket, sorry for the
trouble........Brad
> >
> >
> > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17 AND
> 5/30/17.
> > > RESPONSES MAY BE DELAYED.
> > >
> > > Greetings,
> > >
> > > This message has been automatically generated in response to the
> > > creation of a trouble ticket regarding:
> > > "problem with pcp_combine",
> > > a summary of which appears below.
> > >
> > > There is no need to reply to this message right now. Your
ticket has
> > been
> > > assigned an ID of [rt.rap.ucar.edu #81470].
> > >
> > > Please include the string:
> > >
> > > [rt.rap.ucar.edu #81470]
> > >
> > > in the subject line of all future correspondence about this
issue. To
> do
> > > so,
> > > you may reply to this message.
> > >
> > > Thank you,
> > > met_help at ucar.edu
> > >
> > > ------------------------------------------------------------
> > -------------
> > > To whom it may concern,
> > >
> > > I had the following problem pcp_combine. I supplied two input
grib
> files
> > > from WPC and pcp_combine could not read the second one (see
output
> > below),
> > > both files came from the ncep site at:
> > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > >
> > > Any assistance would be appreciated,
> > >
> > > Brad Diehl
> > > NWS/MDL
> > >
> > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> p06m_2017080300f006.grb
> > 6
> > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > WARNING:
> > > WARNING: grd_file_type() -> the file type indicated by the
suffix
> > > "FileType_Gb1" does not match the file type indicated by the
data
> > > "FileType_Gb2".
> > > WARNING:
> > > ERROR :
> > > ERROR : GribFile::read() -> file read error ... requested -8
bytes,
> got
> > -1
> > > ERROR :
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: John Halley Gotway
Time: Fri Aug 04 12:00:15 2017
Brad,
I have only used "matplotlib" in a class... still spinning up on
python.
But this email is cc'ed to other folks on the MET team. Minna has
done
much more Python plotting than I have and she may be able to answer
your
question. I may be able to answer questions about NetCDF variable
names.
So feel free to shoot.
John
On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>
> John, thanks for your message and help, I'll keep that around for
> reference.
>
> Quick question, have you used matplotlib to plot results of MET or
are you
> basically using the MET Viewer?
> If you have used matplotlib, if you don't mind I have a couple of
> questions, basically about the netcdf file and variables.
>
> Thanks, Brad
>
> On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Brad,
> >
> > I was writing up a response when I saw your email...
> >
> > I see that you're having trouble getting pcp_combine to run on
some GRIB
> > data. As you may have guessed based on the error message, MET is
> > complaining that you're processing GRIB2 files that are named with
a
> ".grb"
> > suffix, which it assumes to be GRIB1.
> >
> > When you pass MET a gridded data file, it needs to figure out how
to
> > process it. The first step is looking at the filename suffix.
Suffixes
> > like .grb, .grib, and .GRIB are assumed to be GRIB1. Suffixes
like
> .grb2,
> > grib2, and .GRIB2 are assumed to be GRIB2. With your data, it's
assumed
> > your file is GRIB1 because of the .grb suffix. And that's the
reason for
> > the error.
> >
> > To get it to work, you could either...
> > (1) Rename the .grb files to .grib2 ... or ...
> > (2) Use command line arguments to tell MET that it's really GRIB2.
> >
> > Here's an example of running MET's plot_data_plane tool, passing
to it a
> > GRIB2 file that named .grb:
> >
> > met-6.0/bin/plot_data_plane \
> > p06m_2017080200f006.grb tmp_p500.ps \
> > 'name="TMP"; level="P500"; file_type = GRIB2;'
> >
> > The "file_type" option explicitly tells MET to interpret this .grb
file
> as
> > a GRIB2 file.
> >
> > John
> >
> > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > >
> > > Upon using degrib, noticed that files were in grib2 and give the
files
> > the
> > > correct extensions and did the following pcp_combine:
> > >
> > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > p06m_2017080200f006.grb2 6
> > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > WARNING:
> > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049) and
> > may
> > > produce unexpected results!
> > > WARNING:
> > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > WARNING:
> > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049) and
> > may
> > > produce unexpected results!
> > > WARNING:
> > > DEBUG 1: Writing output file: p06_add_out.nc
> > >
> > > So, I guess we can close this ticket, sorry for the
trouble........Brad
> > >
> > >
> > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17
AND
> > 5/30/17.
> > > > RESPONSES MAY BE DELAYED.
> > > >
> > > > Greetings,
> > > >
> > > > This message has been automatically generated in response to
the
> > > > creation of a trouble ticket regarding:
> > > > "problem with pcp_combine",
> > > > a summary of which appears below.
> > > >
> > > > There is no need to reply to this message right now. Your
ticket has
> > > been
> > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > >
> > > > Please include the string:
> > > >
> > > > [rt.rap.ucar.edu #81470]
> > > >
> > > > in the subject line of all future correspondence about this
issue. To
> > do
> > > > so,
> > > > you may reply to this message.
> > > >
> > > > Thank you,
> > > > met_help at ucar.edu
> > > >
> > > > ------------------------------------------------------------
> > > -------------
> > > > To whom it may concern,
> > > >
> > > > I had the following problem pcp_combine. I supplied two input
grib
> > files
> > > > from WPC and pcp_combine could not read the second one (see
output
> > > below),
> > > > both files came from the ncep site at:
> > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > >
> > > > Any assistance would be appreciated,
> > > >
> > > > Brad Diehl
> > > > NWS/MDL
> > > >
> > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > p06m_2017080300f006.grb
> > > 6
> > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > WARNING:
> > > > WARNING: grd_file_type() -> the file type indicated by the
suffix
> > > > "FileType_Gb1" does not match the file type indicated by the
data
> > > > "FileType_Gb2".
> > > > WARNING:
> > > > ERROR :
> > > > ERROR : GribFile::read() -> file read error ... requested -8
bytes,
> > got
> > > -1
> > > > ERROR :
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Fri Aug 04 12:26:52 2017
John
I did'nt see a cc on that last message otherwise I would have sent it
straight to Minna, if you could pass it on it would be appreciated.
Anyway, the issue is this:
In the sample matplotlib code I was using, before calling basemap,
there
was the following:
lon_0 = -nc.variables['true_lon'].getValue()
lat_0 = nc.variables['true_lat'].getValue()
Yet when I did print nc.variables in the python code, there were no
netcdf
variables named true_lon or true_lat.
So a little unsure how to define lat_0, lon_0.
On Fri, Aug 4, 2017 at 6:00 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Brad,
>
> I have only used "matplotlib" in a class... still spinning up on
python.
> But this email is cc'ed to other folks on the MET team. Minna has
done
> much more Python plotting than I have and she may be able to answer
your
> question. I may be able to answer questions about NetCDF variable
names.
> So feel free to shoot.
>
> John
>
> On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> >
> > John, thanks for your message and help, I'll keep that around for
> > reference.
> >
> > Quick question, have you used matplotlib to plot results of MET or
are
> you
> > basically using the MET Viewer?
> > If you have used matplotlib, if you don't mind I have a couple of
> > questions, basically about the netcdf file and variables.
> >
> > Thanks, Brad
> >
> > On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Brad,
> > >
> > > I was writing up a response when I saw your email...
> > >
> > > I see that you're having trouble getting pcp_combine to run on
some
> GRIB
> > > data. As you may have guessed based on the error message, MET
is
> > > complaining that you're processing GRIB2 files that are named
with a
> > ".grb"
> > > suffix, which it assumes to be GRIB1.
> > >
> > > When you pass MET a gridded data file, it needs to figure out
how to
> > > process it. The first step is looking at the filename suffix.
> Suffixes
> > > like .grb, .grib, and .GRIB are assumed to be GRIB1. Suffixes
like
> > .grb2,
> > > grib2, and .GRIB2 are assumed to be GRIB2. With your data, it's
> assumed
> > > your file is GRIB1 because of the .grb suffix. And that's the
reason
> for
> > > the error.
> > >
> > > To get it to work, you could either...
> > > (1) Rename the .grb files to .grib2 ... or ...
> > > (2) Use command line arguments to tell MET that it's really
GRIB2.
> > >
> > > Here's an example of running MET's plot_data_plane tool, passing
to it
> a
> > > GRIB2 file that named .grb:
> > >
> > > met-6.0/bin/plot_data_plane \
> > > p06m_2017080200f006.grb tmp_p500.ps \
> > > 'name="TMP"; level="P500"; file_type = GRIB2;'
> > >
> > > The "file_type" option explicitly tells MET to interpret this
.grb file
> > as
> > > a GRIB2 file.
> > >
> > > John
> > >
> > > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470
>
> > > >
> > > > Upon using degrib, noticed that files were in grib2 and give
the
> files
> > > the
> > > > correct extensions and did the following pcp_combine:
> > > >
> > > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > p06m_2017080200f006.grb2 6
> > > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > > WARNING:
> > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> currently
> > > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049)
> and
> > > may
> > > > produce unexpected results!
> > > > WARNING:
> > > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > > WARNING:
> > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> currently
> > > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049)
> and
> > > may
> > > > produce unexpected results!
> > > > WARNING:
> > > > DEBUG 1: Writing output file: p06_add_out.nc
> > > >
> > > > So, I guess we can close this ticket, sorry for the
> trouble........Brad
> > > >
> > > >
> > > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> > > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17
AND
> > > 5/30/17.
> > > > > RESPONSES MAY BE DELAYED.
> > > > >
> > > > > Greetings,
> > > > >
> > > > > This message has been automatically generated in response to
the
> > > > > creation of a trouble ticket regarding:
> > > > > "problem with pcp_combine",
> > > > > a summary of which appears below.
> > > > >
> > > > > There is no need to reply to this message right now. Your
ticket
> has
> > > > been
> > > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > > >
> > > > > Please include the string:
> > > > >
> > > > > [rt.rap.ucar.edu #81470]
> > > > >
> > > > > in the subject line of all future correspondence about this
issue.
> To
> > > do
> > > > > so,
> > > > > you may reply to this message.
> > > > >
> > > > > Thank you,
> > > > > met_help at ucar.edu
> > > > >
> > > > > ------------------------------------------------------------
> > > > -------------
> > > > > To whom it may concern,
> > > > >
> > > > > I had the following problem pcp_combine. I supplied two
input grib
> > > files
> > > > > from WPC and pcp_combine could not read the second one (see
output
> > > > below),
> > > > > both files came from the ncep site at:
> > > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > > >
> > > > > Any assistance would be appreciated,
> > > > >
> > > > > Brad Diehl
> > > > > NWS/MDL
> > > > >
> > > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > p06m_2017080300f006.grb
> > > > 6
> > > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > > WARNING:
> > > > > WARNING: grd_file_type() -> the file type indicated by the
suffix
> > > > > "FileType_Gb1" does not match the file type indicated by the
data
> > > > > "FileType_Gb2".
> > > > > WARNING:
> > > > > ERROR :
> > > > > ERROR : GribFile::read() -> file read error ... requested
-8
> bytes,
> > > got
> > > > -1
> > > > > ERROR :
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Minna Win
Time: Fri Aug 04 13:17:18 2017
Hi Brad,
>From what I can see, it looks like you are trying to do the correct
thing-
assign a list of longitudes to lon_0 and a list of latitudes to lat_0.
Is
it possible for me to see the rest of your script?
Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Fri, Aug 4, 2017 at 6:26 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>
> John
>
> I did'nt see a cc on that last message otherwise I would have sent
it
> straight to Minna, if you could pass it on it would be appreciated.
>
> Anyway, the issue is this:
>
> In the sample matplotlib code I was using, before calling basemap,
there
> was the following:
>
> lon_0 = -nc.variables['true_lon'].getValue()
> lat_0 = nc.variables['true_lat'].getValue()
>
> Yet when I did print nc.variables in the python code, there were no
netcdf
> variables named true_lon or true_lat.
>
> So a little unsure how to define lat_0, lon_0.
>
> On Fri, Aug 4, 2017 at 6:00 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> > wrote:
>
> > Brad,
> >
> > I have only used "matplotlib" in a class... still spinning up on
python.
> > But this email is cc'ed to other folks on the MET team. Minna has
done
> > much more Python plotting than I have and she may be able to
answer your
> > question. I may be able to answer questions about NetCDF variable
names.
> > So feel free to shoot.
> >
> > John
> >
> > On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > >
> > > John, thanks for your message and help, I'll keep that around
for
> > > reference.
> > >
> > > Quick question, have you used matplotlib to plot results of MET
or are
> > you
> > > basically using the MET Viewer?
> > > If you have used matplotlib, if you don't mind I have a couple
of
> > > questions, basically about the netcdf file and variables.
> > >
> > > Thanks, Brad
> > >
> > > On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Brad,
> > > >
> > > > I was writing up a response when I saw your email...
> > > >
> > > > I see that you're having trouble getting pcp_combine to run on
some
> > GRIB
> > > > data. As you may have guessed based on the error message, MET
is
> > > > complaining that you're processing GRIB2 files that are named
with a
> > > ".grb"
> > > > suffix, which it assumes to be GRIB1.
> > > >
> > > > When you pass MET a gridded data file, it needs to figure out
how to
> > > > process it. The first step is looking at the filename suffix.
> > Suffixes
> > > > like .grb, .grib, and .GRIB are assumed to be GRIB1. Suffixes
like
> > > .grb2,
> > > > grib2, and .GRIB2 are assumed to be GRIB2. With your data,
it's
> > assumed
> > > > your file is GRIB1 because of the .grb suffix. And that's the
reason
> > for
> > > > the error.
> > > >
> > > > To get it to work, you could either...
> > > > (1) Rename the .grb files to .grib2 ... or ...
> > > > (2) Use command line arguments to tell MET that it's really
GRIB2.
> > > >
> > > > Here's an example of running MET's plot_data_plane tool,
passing to
> it
> > a
> > > > GRIB2 file that named .grb:
> > > >
> > > > met-6.0/bin/plot_data_plane \
> > > > p06m_2017080200f006.grb tmp_p500.ps \
> > > > 'name="TMP"; level="P500"; file_type = GRIB2;'
> > > >
> > > > The "file_type" option explicitly tells MET to interpret this
.grb
> file
> > > as
> > > > a GRIB2 file.
> > > >
> > > > John
> > > >
> > > > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > > > >
> > > > > Upon using degrib, noticed that files were in grib2 and give
the
> > files
> > > > the
> > > > > correct extensions and did the following pcp_combine:
> > > > >
> > > > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > p06m_2017080200f006.grb2 6
> > > > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > > > WARNING:
> > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> > currently
> > > > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049)
> > and
> > > > may
> > > > > produce unexpected results!
> > > > > WARNING:
> > > > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > > > WARNING:
> > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> > currently
> > > > > support Lambert Conformal grids where dx (5.07862) != dy
(5.08049)
> > and
> > > > may
> > > > > produce unexpected results!
> > > > > WARNING:
> > > > > DEBUG 1: Writing output file: p06_add_out.nc
> > > > >
> > > > > So, I guess we can close this ticket, sorry for the
> > trouble........Brad
> > > > >
> > > > >
> > > > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> > > > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN
5/23/17 AND
> > > > 5/30/17.
> > > > > > RESPONSES MAY BE DELAYED.
> > > > > >
> > > > > > Greetings,
> > > > > >
> > > > > > This message has been automatically generated in response
to the
> > > > > > creation of a trouble ticket regarding:
> > > > > > "problem with pcp_combine",
> > > > > > a summary of which appears below.
> > > > > >
> > > > > > There is no need to reply to this message right now. Your
ticket
> > has
> > > > > been
> > > > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > > > >
> > > > > > Please include the string:
> > > > > >
> > > > > > [rt.rap.ucar.edu #81470]
> > > > > >
> > > > > > in the subject line of all future correspondence about
this
> issue.
> > To
> > > > do
> > > > > > so,
> > > > > > you may reply to this message.
> > > > > >
> > > > > > Thank you,
> > > > > > met_help at ucar.edu
> > > > > >
> > > > > >
------------------------------------------------------------
> > > > > -------------
> > > > > > To whom it may concern,
> > > > > >
> > > > > > I had the following problem pcp_combine. I supplied two
input
> grib
> > > > files
> > > > > > from WPC and pcp_combine could not read the second one
(see
> output
> > > > > below),
> > > > > > both files came from the ncep site at:
> > > > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > > > >
> > > > > > Any assistance would be appreciated,
> > > > > >
> > > > > > Brad Diehl
> > > > > > NWS/MDL
> > > > > >
> > > > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > p06m_2017080300f006.grb
> > > > > 6
> > > > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > > > WARNING:
> > > > > > WARNING: grd_file_type() -> the file type indicated by the
suffix
> > > > > > "FileType_Gb1" does not match the file type indicated by
the data
> > > > > > "FileType_Gb2".
> > > > > > WARNING:
> > > > > > ERROR :
> > > > > > ERROR : GribFile::read() -> file read error ... requested
-8
> > bytes,
> > > > got
> > > > > -1
> > > > > > ERROR :
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Fri Aug 04 14:07:07 2017
Minna,
Thanks for your help. What I discovered was I cannot do the simplest
plot
in matplotlib, as soon as I call any basemap or even one simple plot
via
plt (see below), it all implodes. Should say I am pretty new to all
this.
With people's help advice, as background info, via minianaconda I
built
matplotlib, python and mpl_toolkits, including basemap. So I'm running
everything in python 2.7. I heard all the buzz about matplotlib and
was
looking forward to using it, but feeling pretty frustrated at this
point.
Any help would be appreciated even though this is not directly related
to
MET.
Thanks, Brad
>>> import matplotlib.pyplot as plt
>>> plt.plot([0, 1, 2, 3, 4], [0, 3, 5, 9, 11])
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File
"/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/python2.7/site-
packages/matplotlib/pyplot.py",
line 3306, in plot
ax = gca()
File
"/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/python2.7/site-
packages/matplotlib/pyplot.py",
line 950, in gca
return gcf().gca(**kwargs)
File
"/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/python2.7/site-
packages/matplotlib/pyplot.py",
line 586, in gcf
return figure()
File
"/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/python2.7/site-
packages/matplotlib/pyplot.py",
line 535, in figure
**kwargs)
File
"/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/python2.7/site-
packages/matplotlib/backends/backend_qt5agg.py",
line 44, in new_figure_manager
On Fri, Aug 4, 2017 at 7:17 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Brad,
>
> From what I can see, it looks like you are trying to do the correct
thing-
> assign a list of longitudes to lon_0 and a list of latitudes to
lat_0. Is
> it possible for me to see the rest of your script?
>
> Thanks,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Fri, Aug 4, 2017 at 6:26 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> >
> > John
> >
> > I did'nt see a cc on that last message otherwise I would have sent
it
> > straight to Minna, if you could pass it on it would be
appreciated.
> >
> > Anyway, the issue is this:
> >
> > In the sample matplotlib code I was using, before calling basemap,
there
> > was the following:
> >
> > lon_0 = -nc.variables['true_lon'].getValue()
> > lat_0 = nc.variables['true_lat'].getValue()
> >
> > Yet when I did print nc.variables in the python code, there were
no
> netcdf
> > variables named true_lon or true_lat.
> >
> > So a little unsure how to define lat_0, lon_0.
> >
> > On Fri, Aug 4, 2017 at 6:00 PM, John Halley Gotway via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Brad,
> > >
> > > I have only used "matplotlib" in a class... still spinning up on
> python.
> > > But this email is cc'ed to other folks on the MET team. Minna
has done
> > > much more Python plotting than I have and she may be able to
answer
> your
> > > question. I may be able to answer questions about NetCDF
variable
> names.
> > > So feel free to shoot.
> > >
> > > John
> > >
> > > On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470
>
> > > >
> > > > John, thanks for your message and help, I'll keep that around
for
> > > > reference.
> > > >
> > > > Quick question, have you used matplotlib to plot results of
MET or
> are
> > > you
> > > > basically using the MET Viewer?
> > > > If you have used matplotlib, if you don't mind I have a couple
of
> > > > questions, basically about the netcdf file and variables.
> > > >
> > > > Thanks, Brad
> > > >
> > > > On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > > Hi Brad,
> > > > >
> > > > > I was writing up a response when I saw your email...
> > > > >
> > > > > I see that you're having trouble getting pcp_combine to run
on some
> > > GRIB
> > > > > data. As you may have guessed based on the error message,
MET is
> > > > > complaining that you're processing GRIB2 files that are
named with
> a
> > > > ".grb"
> > > > > suffix, which it assumes to be GRIB1.
> > > > >
> > > > > When you pass MET a gridded data file, it needs to figure
out how
> to
> > > > > process it. The first step is looking at the filename
suffix.
> > > Suffixes
> > > > > like .grb, .grib, and .GRIB are assumed to be GRIB1.
Suffixes like
> > > > .grb2,
> > > > > grib2, and .GRIB2 are assumed to be GRIB2. With your data,
it's
> > > assumed
> > > > > your file is GRIB1 because of the .grb suffix. And that's
the
> reason
> > > for
> > > > > the error.
> > > > >
> > > > > To get it to work, you could either...
> > > > > (1) Rename the .grb files to .grib2 ... or ...
> > > > > (2) Use command line arguments to tell MET that it's really
GRIB2.
> > > > >
> > > > > Here's an example of running MET's plot_data_plane tool,
passing to
> > it
> > > a
> > > > > GRIB2 file that named .grb:
> > > > >
> > > > > met-6.0/bin/plot_data_plane \
> > > > > p06m_2017080200f006.grb tmp_p500.ps \
> > > > > 'name="TMP"; level="P500"; file_type = GRIB2;'
> > > > >
> > > > > The "file_type" option explicitly tells MET to interpret
this .grb
> > file
> > > > as
> > > > > a GRIB2 file.
> > > > >
> > > > > John
> > > > >
> > > > > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > > > > >
> > > > > > Upon using degrib, noticed that files were in grib2 and
give the
> > > files
> > > > > the
> > > > > > correct extensions and did the following pcp_combine:
> > > > > >
> > > > > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > p06m_2017080200f006.grb2 6
> > > > > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > > > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > > > > WARNING:
> > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > > currently
> > > > > > support Lambert Conformal grids where dx (5.07862) != dy
> (5.08049)
> > > and
> > > > > may
> > > > > > produce unexpected results!
> > > > > > WARNING:
> > > > > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > > > > WARNING:
> > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > > currently
> > > > > > support Lambert Conformal grids where dx (5.07862) != dy
> (5.08049)
> > > and
> > > > > may
> > > > > > produce unexpected results!
> > > > > > WARNING:
> > > > > > DEBUG 1: Writing output file: p06_add_out.nc
> > > > > >
> > > > > > So, I guess we can close this ticket, sorry for the
> > > trouble........Brad
> > > > > >
> > > > > >
> > > > > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT <
> > > > > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN
5/23/17 AND
> > > > > 5/30/17.
> > > > > > > RESPONSES MAY BE DELAYED.
> > > > > > >
> > > > > > > Greetings,
> > > > > > >
> > > > > > > This message has been automatically generated in
response to
> the
> > > > > > > creation of a trouble ticket regarding:
> > > > > > > "problem with pcp_combine",
> > > > > > > a summary of which appears below.
> > > > > > >
> > > > > > > There is no need to reply to this message right now.
Your
> ticket
> > > has
> > > > > > been
> > > > > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > > > > >
> > > > > > > Please include the string:
> > > > > > >
> > > > > > > [rt.rap.ucar.edu #81470]
> > > > > > >
> > > > > > > in the subject line of all future correspondence about
this
> > issue.
> > > To
> > > > > do
> > > > > > > so,
> > > > > > > you may reply to this message.
> > > > > > >
> > > > > > > Thank you,
> > > > > > > met_help at ucar.edu
> > > > > > >
> > > > > > >
------------------------------------------------------------
> > > > > > -------------
> > > > > > > To whom it may concern,
> > > > > > >
> > > > > > > I had the following problem pcp_combine. I supplied two
input
> > grib
> > > > > files
> > > > > > > from WPC and pcp_combine could not read the second one
(see
> > output
> > > > > > below),
> > > > > > > both files came from the ncep site at:
> > > > > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > > > > >
> > > > > > > Any assistance would be appreciated,
> > > > > > >
> > > > > > > Brad Diehl
> > > > > > > NWS/MDL
> > > > > > >
> > > > > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > p06m_2017080300f006.grb
> > > > > > 6
> > > > > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > > > > WARNING:
> > > > > > > WARNING: grd_file_type() -> the file type indicated by
the
> suffix
> > > > > > > "FileType_Gb1" does not match the file type indicated by
the
> data
> > > > > > > "FileType_Gb2".
> > > > > > > WARNING:
> > > > > > > ERROR :
> > > > > > > ERROR : GribFile::read() -> file read error ...
requested -8
> > > bytes,
> > > > > got
> > > > > > -1
> > > > > > > ERROR :
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Minna Win
Time: Fri Aug 04 15:39:09 2017
Hi Brad,
It looks like something went awry with you package installation. You
might
want to try installing again. There is a conda troubleshooting page
you
might find useful: https://conda.io/docs/troubleshooting.html
One more thing, do you have proj (http://proj4.org/) and GEOS (
https://trac.osgeo.org/geos/) installed-you will need these for
Basemap? A
'conda install' won't be able to install those dependencies for you so
you'll have to install (or have sys admin) install that for you.
Unfortunately, we can't really be of much additional help since this
is not
directly related to MET. Please don't hesitate to contact us if you
have
any further MET related questions.
Best of luck,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Fri, Aug 4, 2017 at 8:07 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>
> Minna,
>
> Thanks for your help. What I discovered was I cannot do the simplest
plot
> in matplotlib, as soon as I call any basemap or even one simple plot
via
> plt (see below), it all implodes. Should say I am pretty new to all
this.
> With people's help advice, as background info, via minianaconda I
built
> matplotlib, python and mpl_toolkits, including basemap. So I'm
running
> everything in python 2.7. I heard all the buzz about matplotlib and
was
> looking forward to using it, but feeling pretty frustrated at this
point.
> Any help would be appreciated even though this is not directly
related to
> MET.
>
> Thanks, Brad
>
> >>> import matplotlib.pyplot as plt
> >>> plt.plot([0, 1, 2, 3, 4], [0, 3, 5, 9, 11])
> Traceback (most recent call last):
> File "<stdin>", line 1, in <module>
> File
> "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> python2.7/site-packages/matplotlib/pyplot.py",
> line 3306, in plot
> ax = gca()
> File
> "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> python2.7/site-packages/matplotlib/pyplot.py",
> line 950, in gca
> return gcf().gca(**kwargs)
> File
> "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> python2.7/site-packages/matplotlib/pyplot.py",
> line 586, in gcf
> return figure()
> File
> "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> python2.7/site-packages/matplotlib/pyplot.py",
> line 535, in figure
> **kwargs)
> File
> "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> python2.7/site-packages/matplotlib/backends/backend_qt5agg.py",
> line 44, in new_figure_manager
>
>
> On Fri, Aug 4, 2017 at 7:17 PM, Minna Win via RT <met_help at ucar.edu>
> wrote:
>
> > Hi Brad,
> >
> > From what I can see, it looks like you are trying to do the
correct
> thing-
> > assign a list of longitudes to lon_0 and a list of latitudes to
lat_0.
> Is
> > it possible for me to see the rest of your script?
> >
> > Thanks,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Fri, Aug 4, 2017 at 6:26 PM, Brad Diehl - NOAA Federal via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > >
> > > John
> > >
> > > I did'nt see a cc on that last message otherwise I would have
sent it
> > > straight to Minna, if you could pass it on it would be
appreciated.
> > >
> > > Anyway, the issue is this:
> > >
> > > In the sample matplotlib code I was using, before calling
basemap,
> there
> > > was the following:
> > >
> > > lon_0 = -nc.variables['true_lon'].getValue()
> > > lat_0 = nc.variables['true_lat'].getValue()
> > >
> > > Yet when I did print nc.variables in the python code, there were
no
> > netcdf
> > > variables named true_lon or true_lat.
> > >
> > > So a little unsure how to define lat_0, lon_0.
> > >
> > > On Fri, Aug 4, 2017 at 6:00 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Brad,
> > > >
> > > > I have only used "matplotlib" in a class... still spinning up
on
> > python.
> > > > But this email is cc'ed to other folks on the MET team. Minna
has
> done
> > > > much more Python plotting than I have and she may be able to
answer
> > your
> > > > question. I may be able to answer questions about NetCDF
variable
> > names.
> > > > So feel free to shoot.
> > > >
> > > > John
> > > >
> > > > On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > > > >
> > > > > John, thanks for your message and help, I'll keep that
around for
> > > > > reference.
> > > > >
> > > > > Quick question, have you used matplotlib to plot results of
MET or
> > are
> > > > you
> > > > > basically using the MET Viewer?
> > > > > If you have used matplotlib, if you don't mind I have a
couple of
> > > > > questions, basically about the netcdf file and variables.
> > > > >
> > > > > Thanks, Brad
> > > > >
> > > > > On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > > Hi Brad,
> > > > > >
> > > > > > I was writing up a response when I saw your email...
> > > > > >
> > > > > > I see that you're having trouble getting pcp_combine to
run on
> some
> > > > GRIB
> > > > > > data. As you may have guessed based on the error message,
MET is
> > > > > > complaining that you're processing GRIB2 files that are
named
> with
> > a
> > > > > ".grb"
> > > > > > suffix, which it assumes to be GRIB1.
> > > > > >
> > > > > > When you pass MET a gridded data file, it needs to figure
out how
> > to
> > > > > > process it. The first step is looking at the filename
suffix.
> > > > Suffixes
> > > > > > like .grb, .grib, and .GRIB are assumed to be GRIB1.
Suffixes
> like
> > > > > .grb2,
> > > > > > grib2, and .GRIB2 are assumed to be GRIB2. With your
data, it's
> > > > assumed
> > > > > > your file is GRIB1 because of the .grb suffix. And that's
the
> > reason
> > > > for
> > > > > > the error.
> > > > > >
> > > > > > To get it to work, you could either...
> > > > > > (1) Rename the .grb files to .grib2 ... or ...
> > > > > > (2) Use command line arguments to tell MET that it's
really
> GRIB2.
> > > > > >
> > > > > > Here's an example of running MET's plot_data_plane tool,
passing
> to
> > > it
> > > > a
> > > > > > GRIB2 file that named .grb:
> > > > > >
> > > > > > met-6.0/bin/plot_data_plane \
> > > > > > p06m_2017080200f006.grb tmp_p500.ps \
> > > > > > 'name="TMP"; level="P500"; file_type = GRIB2;'
> > > > > >
> > > > > > The "file_type" option explicitly tells MET to interpret
this
> .grb
> > > file
> > > > > as
> > > > > > a GRIB2 file.
> > > > > >
> > > > > > John
> > > > > >
> > > > > > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA Federal
via
> RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470
> >
> > > > > > >
> > > > > > > Upon using degrib, noticed that files were in grib2 and
give
> the
> > > > files
> > > > > > the
> > > > > > > correct extensions and did the following pcp_combine:
> > > > > > >
> > > > > > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > > p06m_2017080200f006.grb2 6
> > > > > > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > > > > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > > > > > WARNING:
> > > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > > > currently
> > > > > > > support Lambert Conformal grids where dx (5.07862) != dy
> > (5.08049)
> > > > and
> > > > > > may
> > > > > > > produce unexpected results!
> > > > > > > WARNING:
> > > > > > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > > > > > WARNING:
> > > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > > > currently
> > > > > > > support Lambert Conformal grids where dx (5.07862) != dy
> > (5.08049)
> > > > and
> > > > > > may
> > > > > > > produce unexpected results!
> > > > > > > WARNING:
> > > > > > > DEBUG 1: Writing output file: p06_add_out.nc
> > > > > > >
> > > > > > > So, I guess we can close this ticket, sorry for the
> > > > trouble........Brad
> > > > > > >
> > > > > > >
> > > > > > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via RT
<
> > > > > > > met_help at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN
5/23/17
> AND
> > > > > > 5/30/17.
> > > > > > > > RESPONSES MAY BE DELAYED.
> > > > > > > >
> > > > > > > > Greetings,
> > > > > > > >
> > > > > > > > This message has been automatically generated in
response to
> > the
> > > > > > > > creation of a trouble ticket regarding:
> > > > > > > > "problem with pcp_combine",
> > > > > > > > a summary of which appears below.
> > > > > > > >
> > > > > > > > There is no need to reply to this message right now.
Your
> > ticket
> > > > has
> > > > > > > been
> > > > > > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > > > > > >
> > > > > > > > Please include the string:
> > > > > > > >
> > > > > > > > [rt.rap.ucar.edu #81470]
> > > > > > > >
> > > > > > > > in the subject line of all future correspondence about
this
> > > issue.
> > > > To
> > > > > > do
> > > > > > > > so,
> > > > > > > > you may reply to this message.
> > > > > > > >
> > > > > > > > Thank you,
> > > > > > > > met_help at ucar.edu
> > > > > > > >
> > > > > > > >
------------------------------------------------------------
> > > > > > > -------------
> > > > > > > > To whom it may concern,
> > > > > > > >
> > > > > > > > I had the following problem pcp_combine. I supplied
two input
> > > grib
> > > > > > files
> > > > > > > > from WPC and pcp_combine could not read the second one
(see
> > > output
> > > > > > > below),
> > > > > > > > both files came from the ncep site at:
> > > > > > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > > > > > >
> > > > > > > > Any assistance would be appreciated,
> > > > > > > >
> > > > > > > > Brad Diehl
> > > > > > > > NWS/MDL
> > > > > > > >
> > > > > > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > > p06m_2017080300f006.grb
> > > > > > > 6
> > > > > > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > > > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > > > > > WARNING:
> > > > > > > > WARNING: grd_file_type() -> the file type indicated by
the
> > suffix
> > > > > > > > "FileType_Gb1" does not match the file type indicated
by the
> > data
> > > > > > > > "FileType_Gb2".
> > > > > > > > WARNING:
> > > > > > > > ERROR :
> > > > > > > > ERROR : GribFile::read() -> file read error ...
requested -8
> > > > bytes,
> > > > > > got
> > > > > > > -1
> > > > > > > > ERROR :
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Julie Prestopnik
Time: Mon Aug 07 13:13:42 2017
I see that you are having some trouble running some files through
pcp_combine. I believe the problem is due to the suffix of thsoe
files.
MET interprets .gb, .grb, and .grib to be GRIB1 files and .gb2, .grb2,
and
.grib2 to be GRIB2 files. So it is trying to read a GRIB2 file as a
GRIB1
files.
When I change the .grb extension to .grb2, I am able to get a
successful
run, however with some WARNING notices:
/usr/local/met-6.0/bin/pcp_combine -add p06m_2017080300f006.grb2 6
p06m_2017080312f006.grb2 6 p06_add_out.nc
DEBUG 1: Reading input file: p06m_2017080300f006.grb2
WARNING:
WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not currently
support Lambert Conformal grids where dx (5.07862) != dy (5.08049) and
may
produce unexpected results!
WARNING:
WARNING:
WARNING: Multiple GRIB2 table entries match lookup criteria (parm_name
=
APCP):
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr =
0,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr =
57,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr =
57,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
= 0, grib2_ltab = 0, index_b = 1, index_c = 8
WARNING:
DEBUG 1: Reading input file: p06m_2017080312f006.grb2
WARNING:
WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not currently
support Lambert Conformal grids where dx (5.07862) != dy (5.08049) and
may
produce unexpected results!
WARNING:
WARNING:
WARNING: Multiple GRIB2 table entries match lookup criteria (parm_name
=
APCP):
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr =
0,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr =
57,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr =
57,
grib2_ltab = 0, index_b = 1, index_c = 8
WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
= 0, grib2_ltab = 0, index_b = 1, index_c = 8
WARNING:
DEBUG 1: Writing output file: p06_add_out.nc
MET expects the dx and dy values to be the same, however, it looks
like
there are different values for dx and dy in these files. MET also
gives a
warning about multiple GRIB2 table entries matching the lookup
criteria.
Please take a look and verify MET is using the correct one.
I hope that helps!
Thanks,
Julie
On Fri, Aug 4, 2017 at 8:04 AM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> Fri Aug 04 08:04:41 2017: Request 81470 was acted upon.
> Transaction: Ticket created by brad.diehl at noaa.gov
> Queue: met_help
> Subject: problem with pcp_combine
> Owner: Nobody
> Requestors: brad.diehl at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>
>
> To whom it may concern,
>
> I had the following problem pcp_combine. I supplied two input grib
files
> from WPC and pcp_combine could not read the second one (see output
below),
> both files came from the ncep site at:
> http://ftp.wpc.ncep.noaa.gov/5km_qpf/
>
> Any assistance would be appreciated,
>
> Brad Diehl
> NWS/MDL
>
> (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080300f006.grb 6
> p06m_2017080312f006.grb 6 p06_add_out.nc
> DEBUG 1: Reading input file: p06m_2017080300f006.grb
> WARNING:
> WARNING: grd_file_type() -> the file type indicated by the suffix
> "FileType_Gb1" does not match the file type indicated by the data
> "FileType_Gb2".
> WARNING:
> ERROR :
> ERROR : GribFile::read() -> file read error ... requested -8 bytes,
got -1
> ERROR :
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Julie Prestopnik
Time: Mon Aug 07 13:14:41 2017
My apologies, I didn't see the other correspondence before responding
to
this ticket.
Julie
On Mon, Aug 7, 2017 at 1:13 PM, Julie Prestopnik <jpresto at ucar.edu>
wrote:
> I see that you are having some trouble running some files through
> pcp_combine. I believe the problem is due to the suffix of thsoe
files.
> MET interprets .gb, .grb, and .grib to be GRIB1 files and .gb2,
.grb2, and
> .grib2 to be GRIB2 files. So it is trying to read a GRIB2 file as a
GRIB1
> files.
>
> When I change the .grb extension to .grb2, I am able to get a
successful
> run, however with some WARNING notices:
>
> /usr/local/met-6.0/bin/pcp_combine -add p06m_2017080300f006.grb2 6
> p06m_2017080312f006.grb2 6 p06_add_out.nc
> DEBUG 1: Reading input file: p06m_2017080300f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> WARNING:
> WARNING: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> APCP):
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
> = 0, grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING:
> DEBUG 1: Reading input file: p06m_2017080312f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> WARNING:
> WARNING: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> APCP):
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
> = 0, grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING:
> DEBUG 1: Writing output file: p06_add_out.nc
>
> MET expects the dx and dy values to be the same, however, it looks
like
> there are different values for dx and dy in these files. MET also
gives a
> warning about multiple GRIB2 table entries matching the lookup
criteria.
> Please take a look and verify MET is using the correct one.
>
> I hope that helps!
>
> Thanks,
> Julie
>
>
>
>
>
> On Fri, Aug 4, 2017 at 8:04 AM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
>>
>> Fri Aug 04 08:04:41 2017: Request 81470 was acted upon.
>> Transaction: Ticket created by brad.diehl at noaa.gov
>> Queue: met_help
>> Subject: problem with pcp_combine
>> Owner: Nobody
>> Requestors: brad.diehl at noaa.gov
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
>>
>>
>> To whom it may concern,
>>
>> I had the following problem pcp_combine. I supplied two input grib
files
>> from WPC and pcp_combine could not read the second one (see output
below),
>> both files came from the ncep site at:
>> http://ftp.wpc.ncep.noaa.gov/5km_qpf/
>>
>> Any assistance would be appreciated,
>>
>> Brad Diehl
>> NWS/MDL
>>
>> (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080300f006.grb 6
>> p06m_2017080312f006.grb 6 p06_add_out.nc
>> DEBUG 1: Reading input file: p06m_2017080300f006.grb
>> WARNING:
>> WARNING: grd_file_type() -> the file type indicated by the suffix
>> "FileType_Gb1" does not match the file type indicated by the data
>> "FileType_Gb2".
>> WARNING:
>> ERROR :
>> ERROR : GribFile::read() -> file read error ... requested -8
bytes, got
>> -1
>> ERROR :
>>
>>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Mon Aug 07 14:36:49 2017
Julie,
I did eventually see the issue with grib2, then everything ran fine. I
thought I closed the trouble ticket.
But I appreciate you trying to help.
Thanks, Brad
On Mon, Aug 7, 2017 at 7:13 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> I see that you are having some trouble running some files through
> pcp_combine. I believe the problem is due to the suffix of thsoe
files.
> MET interprets .gb, .grb, and .grib to be GRIB1 files and .gb2,
.grb2, and
> .grib2 to be GRIB2 files. So it is trying to read a GRIB2 file as a
GRIB1
> files.
>
> When I change the .grb extension to .grb2, I am able to get a
successful
> run, however with some WARNING notices:
>
> /usr/local/met-6.0/bin/pcp_combine -add p06m_2017080300f006.grb2 6
> p06m_2017080312f006.grb2 6 p06_add_out.nc
> DEBUG 1: Reading input file: p06m_2017080300f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> WARNING:
> WARNING: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> APCP):
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
> = 0, grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING:
> DEBUG 1: Reading input file: p06m_2017080312f006.grb2
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (5.07862) != dy (5.08049)
and may
> produce unexpected results!
> WARNING:
> WARNING:
> WARNING: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> APCP):
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: parm_name: APCP, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING: Using: parm_name: APCP, index_a = 0, grib2_mtab = 1,
grib2_cntr
> = 0, grib2_ltab = 0, index_b = 1, index_c = 8
> WARNING:
> DEBUG 1: Writing output file: p06_add_out.nc
>
> MET expects the dx and dy values to be the same, however, it looks
like
> there are different values for dx and dy in these files. MET also
gives a
> warning about multiple GRIB2 table entries matching the lookup
criteria.
> Please take a look and verify MET is using the correct one.
>
> I hope that helps!
>
> Thanks,
> Julie
>
>
>
>
>
> On Fri, Aug 4, 2017 at 8:04 AM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Fri Aug 04 08:04:41 2017: Request 81470 was acted upon.
> > Transaction: Ticket created by brad.diehl at noaa.gov
> > Queue: met_help
> > Subject: problem with pcp_combine
> > Owner: Nobody
> > Requestors: brad.diehl at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> >
> >
> > To whom it may concern,
> >
> > I had the following problem pcp_combine. I supplied two input grib
files
> > from WPC and pcp_combine could not read the second one (see output
> below),
> > both files came from the ncep site at:
> > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> >
> > Any assistance would be appreciated,
> >
> > Brad Diehl
> > NWS/MDL
> >
> > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
p06m_2017080300f006.grb
> 6
> > p06m_2017080312f006.grb 6 p06_add_out.nc
> > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > WARNING:
> > WARNING: grd_file_type() -> the file type indicated by the suffix
> > "FileType_Gb1" does not match the file type indicated by the data
> > "FileType_Gb2".
> > WARNING:
> > ERROR :
> > ERROR : GribFile::read() -> file read error ... requested -8
bytes, got
> -1
> > ERROR :
> >
> >
>
>
------------------------------------------------
Subject: problem with pcp_combine
From: Brad Diehl - NOAA Federal
Time: Tue Aug 08 15:28:09 2017
Minna,
I solved the problems I had building matplotlib, netcdf, etc., via
conda
and can plot things via matplotlib, read netcdf4 files, etc. Basically
good
to go. But I'm hung up on plotting a contour to the basemap. This is
indirectly related to MET-6.0 in that I used MET to add two NCEP QPF
files
and just simply want to plot the results of the MET ncdf file. I know
what
I need to know, but here's where I'm hung up: How can I get lat_1,
lon_1
and lon_0 and lat_0 for the matplotlib basemap call? I look at ncdump
on
the file and there are no obvious variables of that name (the ncdump
results are below). When I run my matplotlib code (attached) there
are no
errors, yet I know that the map projection is not quite right. At
first I
thought maybe lat_pin and lon_pin might be lon_0 and lat_0 and used
those
in my basemap call.
If you can offer any help, great, if not I understand and appreciate
your
original help.
Thanks, Brad
netcdf p06_add_out {
dimensions:
lat = 689 ;
lon = 1073 ;
variables:
float lat(lat, lon) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
float lon(lat, lon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:standard_name = "longitude" ;
float APCP_12(lat, lon) ;
APCP_12:name = "APCP_12" ;
APCP_12:long_name = "Total Precipitation" ;
APCP_12:level = "A12" ;
APCP_12:units = "kg/m^2" ;
APCP_12:_FillValue = -9999.f ;
APCP_12:init_time = "20170802_180000" ;
APCP_12:init_time_ut = "1501696800" ;
APCP_12:valid_time = "20170802_180000" ;
APCP_12:valid_time_ut = "1501696800" ;
APCP_12:accum_time = "120000" ;
APCP_12:accum_time_sec = 43200 ;
// global attributes:
:FileOrigins = "File p06_add_out.nc generated
20170804_172418 UTC on host g14a1 by the MET pcp_combine tool" ;
:MET_version = "V6.0" ;
:MET_tool = "pcp_combine" ;
:RunCommand = "Addition: 2 files." ;
:Projection = "Lambert Conformal" ;
:scale_lat_1 = "25.000000" ;
:scale_lat_2 = "25.000000" ;
:lat_pin = "20.190001" ;
:lon_pin = "-121.550004" ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = "-95.000000" ;
:d_km = "5.078616" ;
:r_km = "6371.200000" ;
:nx = "1073" ;
:ny = "689 grid_points" ;
On Fri, Aug 4, 2017 at 9:39 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Brad,
>
> It looks like something went awry with you package installation.
You might
> want to try installing again. There is a conda troubleshooting page
you
> might find useful: https://conda.io/docs/troubleshooting.html
>
> One more thing, do you have proj (http://proj4.org/) and GEOS (
> https://trac.osgeo.org/geos/) installed-you will need these for
Basemap? A
> 'conda install' won't be able to install those dependencies for you
so
> you'll have to install (or have sys admin) install that for you.
> Unfortunately, we can't really be of much additional help since this
is not
> directly related to MET. Please don't hesitate to contact us if you
have
> any further MET related questions.
>
> Best of luck,
> Minna
>
>
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Fri, Aug 4, 2017 at 8:07 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> >
> > Minna,
> >
> > Thanks for your help. What I discovered was I cannot do the
simplest plot
> > in matplotlib, as soon as I call any basemap or even one simple
plot via
> > plt (see below), it all implodes. Should say I am pretty new to
all this.
> > With people's help advice, as background info, via minianaconda I
built
> > matplotlib, python and mpl_toolkits, including basemap. So I'm
running
> > everything in python 2.7. I heard all the buzz about matplotlib
and was
> > looking forward to using it, but feeling pretty frustrated at this
point.
> > Any help would be appreciated even though this is not directly
related to
> > MET.
> >
> > Thanks, Brad
> >
> > >>> import matplotlib.pyplot as plt
> > >>> plt.plot([0, 1, 2, 3, 4], [0, 3, 5, 9, 11])
> > Traceback (most recent call last):
> > File "<stdin>", line 1, in <module>
> > File
> > "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> > python2.7/site-packages/matplotlib/pyplot.py",
> > line 3306, in plot
> > ax = gca()
> > File
> > "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> > python2.7/site-packages/matplotlib/pyplot.py",
> > line 950, in gca
> > return gcf().gca(**kwargs)
> > File
> > "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> > python2.7/site-packages/matplotlib/pyplot.py",
> > line 586, in gcf
> > return figure()
> > File
> > "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> > python2.7/site-packages/matplotlib/pyplot.py",
> > line 535, in figure
> > **kwargs)
> > File
> > "/mdlstat/save/usr/Brad.Diehl/miniconda2/envs/QPF/lib/
> > python2.7/site-packages/matplotlib/backends/backend_qt5agg.py",
> > line 44, in new_figure_manager
> >
> >
> > On Fri, Aug 4, 2017 at 7:17 PM, Minna Win via RT
<met_help at ucar.edu>
> > wrote:
> >
> > > Hi Brad,
> > >
> > > From what I can see, it looks like you are trying to do the
correct
> > thing-
> > > assign a list of longitudes to lon_0 and a list of latitudes to
lat_0.
> > Is
> > > it possible for me to see the rest of your script?
> > >
> > > Thanks,
> > > Minna
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423
> > > Fax: 302-497-8401
> > >
> > >
> > > On Fri, Aug 4, 2017 at 6:26 PM, Brad Diehl - NOAA Federal via RT
<
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470
>
> > > >
> > > > John
> > > >
> > > > I did'nt see a cc on that last message otherwise I would have
sent it
> > > > straight to Minna, if you could pass it on it would be
appreciated.
> > > >
> > > > Anyway, the issue is this:
> > > >
> > > > In the sample matplotlib code I was using, before calling
basemap,
> > there
> > > > was the following:
> > > >
> > > > lon_0 = -nc.variables['true_lon'].getValue()
> > > > lat_0 = nc.variables['true_lat'].getValue()
> > > >
> > > > Yet when I did print nc.variables in the python code, there
were no
> > > netcdf
> > > > variables named true_lon or true_lat.
> > > >
> > > > So a little unsure how to define lat_0, lon_0.
> > > >
> > > > On Fri, Aug 4, 2017 at 6:00 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > > Brad,
> > > > >
> > > > > I have only used "matplotlib" in a class... still spinning
up on
> > > python.
> > > > > But this email is cc'ed to other folks on the MET team.
Minna has
> > done
> > > > > much more Python plotting than I have and she may be able to
answer
> > > your
> > > > > question. I may be able to answer questions about NetCDF
variable
> > > names.
> > > > > So feel free to shoot.
> > > > >
> > > > > John
> > > > >
> > > > > On Fri, Aug 4, 2017 at 11:52 AM, Brad Diehl - NOAA Federal
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81470 >
> > > > > >
> > > > > > John, thanks for your message and help, I'll keep that
around
> for
> > > > > > reference.
> > > > > >
> > > > > > Quick question, have you used matplotlib to plot results
of MET
> or
> > > are
> > > > > you
> > > > > > basically using the MET Viewer?
> > > > > > If you have used matplotlib, if you don't mind I have a
couple of
> > > > > > questions, basically about the netcdf file and variables.
> > > > > >
> > > > > > Thanks, Brad
> > > > > >
> > > > > > On Fri, Aug 4, 2017 at 5:44 PM, John Halley Gotway via RT
<
> > > > > > met_help at ucar.edu
> > > > > > > wrote:
> > > > > >
> > > > > > > Hi Brad,
> > > > > > >
> > > > > > > I was writing up a response when I saw your email...
> > > > > > >
> > > > > > > I see that you're having trouble getting pcp_combine to
run on
> > some
> > > > > GRIB
> > > > > > > data. As you may have guessed based on the error
message, MET
> is
> > > > > > > complaining that you're processing GRIB2 files that are
named
> > with
> > > a
> > > > > > ".grb"
> > > > > > > suffix, which it assumes to be GRIB1.
> > > > > > >
> > > > > > > When you pass MET a gridded data file, it needs to
figure out
> how
> > > to
> > > > > > > process it. The first step is looking at the filename
suffix.
> > > > > Suffixes
> > > > > > > like .grb, .grib, and .GRIB are assumed to be GRIB1.
Suffixes
> > like
> > > > > > .grb2,
> > > > > > > grib2, and .GRIB2 are assumed to be GRIB2. With your
data,
> it's
> > > > > assumed
> > > > > > > your file is GRIB1 because of the .grb suffix. And
that's the
> > > reason
> > > > > for
> > > > > > > the error.
> > > > > > >
> > > > > > > To get it to work, you could either...
> > > > > > > (1) Rename the .grb files to .grib2 ... or ...
> > > > > > > (2) Use command line arguments to tell MET that it's
really
> > GRIB2.
> > > > > > >
> > > > > > > Here's an example of running MET's plot_data_plane tool,
> passing
> > to
> > > > it
> > > > > a
> > > > > > > GRIB2 file that named .grb:
> > > > > > >
> > > > > > > met-6.0/bin/plot_data_plane \
> > > > > > > p06m_2017080200f006.grb tmp_p500.ps \
> > > > > > > 'name="TMP"; level="P500"; file_type = GRIB2;'
> > > > > > >
> > > > > > > The "file_type" option explicitly tells MET to interpret
this
> > .grb
> > > > file
> > > > > > as
> > > > > > > a GRIB2 file.
> > > > > > >
> > > > > > > John
> > > > > > >
> > > > > > > On Fri, Aug 4, 2017 at 11:29 AM, Brad Diehl - NOAA
Federal via
> > RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81470
> > >
> > > > > > > >
> > > > > > > > Upon using degrib, noticed that files were in grib2
and give
> > the
> > > > > files
> > > > > > > the
> > > > > > > > correct extensions and did the following pcp_combine:
> > > > > > > >
> > > > > > > > (g14a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > > > p06m_2017080200f006.grb2 6
> > > > > > > > p06m_2017080212f006.grb2 6 p06_add_out.nc
> > > > > > > > DEBUG 1: Reading input file: p06m_2017080200f006.grb2
> > > > > > > > WARNING:
> > > > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET
does not
> > > > > currently
> > > > > > > > support Lambert Conformal grids where dx (5.07862) !=
dy
> > > (5.08049)
> > > > > and
> > > > > > > may
> > > > > > > > produce unexpected results!
> > > > > > > > WARNING:
> > > > > > > > DEBUG 1: Reading input file: p06m_2017080212f006.grb2
> > > > > > > > WARNING:
> > > > > > > > WARNING: MetGrib2DataFile::read_grib2_grid() -> MET
does not
> > > > > currently
> > > > > > > > support Lambert Conformal grids where dx (5.07862) !=
dy
> > > (5.08049)
> > > > > and
> > > > > > > may
> > > > > > > > produce unexpected results!
> > > > > > > > WARNING:
> > > > > > > > DEBUG 1: Writing output file: p06_add_out.nc
> > > > > > > >
> > > > > > > > So, I guess we can close this ticket, sorry for the
> > > > > trouble........Brad
> > > > > > > >
> > > > > > > >
> > > > > > > > On Fri, Aug 4, 2017 at 2:04 PM, met_help at ucar.edu via
RT <
> > > > > > > > met_help at ucar.edu>
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > > NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN
5/23/17
> > AND
> > > > > > > 5/30/17.
> > > > > > > > > RESPONSES MAY BE DELAYED.
> > > > > > > > >
> > > > > > > > > Greetings,
> > > > > > > > >
> > > > > > > > > This message has been automatically generated in
response
> to
> > > the
> > > > > > > > > creation of a trouble ticket regarding:
> > > > > > > > > "problem with pcp_combine",
> > > > > > > > > a summary of which appears below.
> > > > > > > > >
> > > > > > > > > There is no need to reply to this message right now.
Your
> > > ticket
> > > > > has
> > > > > > > > been
> > > > > > > > > assigned an ID of [rt.rap.ucar.edu #81470].
> > > > > > > > >
> > > > > > > > > Please include the string:
> > > > > > > > >
> > > > > > > > > [rt.rap.ucar.edu #81470]
> > > > > > > > >
> > > > > > > > > in the subject line of all future correspondence
about this
> > > > issue.
> > > > > To
> > > > > > > do
> > > > > > > > > so,
> > > > > > > > > you may reply to this message.
> > > > > > > > >
> > > > > > > > > Thank you,
> > > > > > > > > met_help at ucar.edu
> > > > > > > > >
> > > > > > > > > ------------------------------
> ------------------------------
> > > > > > > > -------------
> > > > > > > > > To whom it may concern,
> > > > > > > > >
> > > > > > > > > I had the following problem pcp_combine. I supplied
two
> input
> > > > grib
> > > > > > > files
> > > > > > > > > from WPC and pcp_combine could not read the second
one (see
> > > > output
> > > > > > > > below),
> > > > > > > > > both files came from the ncep site at:
> > > > > > > > > http://ftp.wpc.ncep.noaa.gov/5km_qpf/
> > > > > > > > >
> > > > > > > > > Any assistance would be appreciated,
> > > > > > > > >
> > > > > > > > > Brad Diehl
> > > > > > > > > NWS/MDL
> > > > > > > > >
> > > > > > > > > (g10a1) /u/Brad.Diehl/map_plot$ pcp_combine -add
> > > > > > > p06m_2017080300f006.grb
> > > > > > > > 6
> > > > > > > > > p06m_2017080312f006.grb 6 p06_add_out.nc
> > > > > > > > > DEBUG 1: Reading input file: p06m_2017080300f006.grb
> > > > > > > > > WARNING:
> > > > > > > > > WARNING: grd_file_type() -> the file type indicated
by the
> > > suffix
> > > > > > > > > "FileType_Gb1" does not match the file type
indicated by
> the
> > > data
> > > > > > > > > "FileType_Gb2".
> > > > > > > > > WARNING:
> > > > > > > > > ERROR :
> > > > > > > > > ERROR : GribFile::read() -> file read error ...
requested
> -8
> > > > > bytes,
> > > > > > > got
> > > > > > > > -1
> > > > > > > > > ERROR :
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list