[Met_help] [rt.rap.ucar.edu #87647] History for ECMWF grib file question
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:07:17 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Good Day,
I am seeking to run MET on an ECMWF grib file I have for 100 meter winds. However, the ECMWF grib file I have list the UWND at 100m as var246 at level sfc and VWND as var247 at level sfc.
Is there a way for MET to read this and properly perform the calculations. Perhaps changing a grid table or configuration file in MET somewhere? Thanks.
Allan Huffman
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: ECMWF grib file question
From: John Halley Gotway
Time: Mon Nov 05 11:33:39 2018
Allan,
Please go ahead and try running your ECMWF data through MET since it
should
already include support for reading/plotting ECMWF GRIB data. When
trying
to figure out if MET reads a particular data source well, I suggest
running
the *plot_data_plane* tool. You request the data that you'd like on
the
command line, it reads in a single 2d field of data, and then creates
a
postscript plot of it.
Running plot_data_plane in this ways checks...
(1) That MET understands the GRIB code abbreviations you're using.
(2) That you're setting the level values correctly.
(3) Running at verbosity level 3 (-v 3), it tells you which GRIB
record
number its reading.
(4) Looking at the resulting image tells you that MET is placing the
data
in the correct spot on the earth.
The MET tarball includes some sample GRIB files... but none of them
are
ECMWF. So I'll use an NCEP one as an example instead, in which U and
V
wind are stored as UGRD and VGRD:
- First, run wgrib or (wgrib2 for GRIB2) to inventory your file...
* wgrib met-
8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 |
grep UGRD | grep "m above ground"*
44:555724:d=05080700:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
m above gnd:12hr fcst:NAve=0
84:1151750:d=05080700:RH:kpds5=52:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
m above gnd:12hr fcst:NAve=0
163:2230984:d=05080700:UGRD:kpds5=33:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
m above gnd:12hr fcst:NAve=0
164:2242920:d=05080700:VGRD:kpds5=34:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
m above gnd:12hr fcst:NAve=0
I have 2-meter temperature and relative humidity and 10-meter winds.
- Next run plot_data_plane to plot the 10-m UGRD field:
* /usr/local/met-8.0/bin/plot_data_plane \*
* met-8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \*
* UGRD_Z10.ps 'name="UGRD"; level="Z10";' -v 3*
Looking in the status messages, I see:
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for
VarInfo "UGRD/Z10" in GRIB record 163 of GRIB file
And indeed wgrib confirms that record number 163 contains 10 meter U
wind.
And use the "gv" or "ghostview" tool to view the resulting image:
* gv UGRD_Z10.ps*
Now, do the same steps with your ECMWF data, but you want to UWND at
100
meters, so use:
* 'name="UWND"; level="Z100";'*
Hopefully that helps get you going.
Thanks,
John Halley Gotway
On Sat, Nov 3, 2018 at 7:19 AM Allan Huffman via RT
<met_help at ucar.edu>
wrote:
>
> Sat Nov 03 07:19:20 2018: Request 87647 was acted upon.
> Transaction: Ticket created by awhuffma at yahoo.com
> Queue: met_help
> Subject: ECMWF grib file question
> Owner: Nobody
> Requestors: awhuffma at yahoo.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87647 >
>
>
> Good Day,
> I am seeking to run MET on an ECMWF grib file I have for 100 meter
winds.
> However, the ECMWF grib file I have list the UWND at 100m as var246
at
> level sfc and VWND as var247 at level sfc.
> Is there a way for MET to read this and properly perform the
calculations.
> Perhaps changing a grid table or configuration file in MET
somewhere?
> Thanks.
> Allan Huffman
>
------------------------------------------------
Subject: ECMWF grib file question
From: Allan Huffman
Time: Wed Nov 07 05:12:51 2018
Thanks John!
On Monday, November 5, 2018, 1:33:42 PM EST, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Allan,
Please go ahead and try running your ECMWF data through MET since it
should
already include support for reading/plotting ECMWF GRIB data. When
trying
to figure out if MET reads a particular data source well, I suggest
running
the *plot_data_plane* tool. You request the data that you'd like on
the
command line, it reads in a single 2d field of data, and then creates
a
postscript plot of it.
Running plot_data_plane in this ways checks...
(1) That MET understands the GRIB code abbreviations you're using.
(2) That you're setting the level values correctly.
(3) Running at verbosity level 3 (-v 3), it tells you which GRIB
record
number its reading.
(4) Looking at the resulting image tells you that MET is placing the
data
in the correct spot on the earth.
The MET tarball includes some sample GRIB files... but none of them
are
ECMWF. So I'll use an NCEP one as an example instead, in which U and
V
wind are stored as UGRD and VGRD:
- First, run wgrib or (wgrib2 for GRIB2) to inventory your file...
* wgrib met-8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
|
grep UGRD | grep "m above ground"*
44:555724:d=05080700:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
m above gnd:12hr fcst:NAve=0
84:1151750:d=05080700:RH:kpds5=52:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
m above gnd:12hr fcst:NAve=0
163:2230984:d=05080700:UGRD:kpds5=33:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
m above gnd:12hr fcst:NAve=0
164:2242920:d=05080700:VGRD:kpds5=34:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
m above gnd:12hr fcst:NAve=0
I have 2-meter temperature and relative humidity and 10-meter winds.
- Next run plot_data_plane to plot the 10-m UGRD field:
* /usr/local/met-8.0/bin/plot_data_plane \*
* met-8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \*
* UGRD_Z10.ps 'name="UGRD"; level="Z10";' -v 3*
Looking in the status messages, I see:
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for
VarInfo "UGRD/Z10" in GRIB record 163 of GRIB file
And indeed wgrib confirms that record number 163 contains 10 meter U
wind.
And use the "gv" or "ghostview" tool to view the resulting image:
* gv UGRD_Z10.ps*
Now, do the same steps with your ECMWF data, but you want to UWND at
100
meters, so use:
* 'name="UWND"; level="Z100";'*
Hopefully that helps get you going.
Thanks,
John Halley Gotway
On Sat, Nov 3, 2018 at 7:19 AM Allan Huffman via RT
<met_help at ucar.edu>
wrote:
>
> Sat Nov 03 07:19:20 2018: Request 87647 was acted upon.
> Transaction: Ticket created by awhuffma at yahoo.com
> Queue: met_help
> Subject: ECMWF grib file question
> Owner: Nobody
> Requestors: awhuffma at yahoo.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87647 >
>
>
> Good Day,
> I am seeking to run MET on an ECMWF grib file I have for 100 meter
winds.
> However, the ECMWF grib file I have list the UWND at 100m as var246
at
> level sfc and VWND as var247 at level sfc.
> Is there a way for MET to read this and properly perform the
calculations.
> Perhaps changing a grid table or configuration file in MET
somewhere?
> Thanks.
> Allan Huffman
>
------------------------------------------------
Subject: ECMWF grib file question
From: John Halley Gotway
Time: Wed Nov 07 08:13:27 2018
Allan,
Please let me know if/when you've been successful, at which point I'll
resolve this support ticket.
Thanks,
John
On Wed, Nov 7, 2018 at 5:13 AM Allan Huffman via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87647 >
>
> Thanks John!
> On Monday, November 5, 2018, 1:33:42 PM EST, John Halley Gotway
via RT
> <met_help at ucar.edu> wrote:
>
> Allan,
>
> Please go ahead and try running your ECMWF data through MET since it
should
> already include support for reading/plotting ECMWF GRIB data. When
trying
> to figure out if MET reads a particular data source well, I suggest
running
> the *plot_data_plane* tool. You request the data that you'd like on
the
> command line, it reads in a single 2d field of data, and then
creates a
> postscript plot of it.
>
> Running plot_data_plane in this ways checks...
> (1) That MET understands the GRIB code abbreviations you're using.
> (2) That you're setting the level values correctly.
> (3) Running at verbosity level 3 (-v 3), it tells you which GRIB
record
> number its reading.
> (4) Looking at the resulting image tells you that MET is placing the
data
> in the correct spot on the earth.
>
> The MET tarball includes some sample GRIB files... but none of them
are
> ECMWF. So I'll use an NCEP one as an example instead, in which U
and V
> wind are stored as UGRD and VGRD:
>
> - First, run wgrib or (wgrib2 for GRIB2) to inventory your file...
> * wgrib met-
8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 |
> grep UGRD | grep "m above ground"*
>
>
>
44:555724:d=05080700:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
> m above gnd:12hr fcst:NAve=0
>
>
84:1151750:d=05080700:RH:kpds5=52:kpds6=105:kpds7=2:TR=0:P1=12:P2=0:TimeU=1:2
> m above gnd:12hr fcst:NAve=0
>
>
163:2230984:d=05080700:UGRD:kpds5=33:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
> m above gnd:12hr fcst:NAve=0
>
>
164:2242920:d=05080700:VGRD:kpds5=34:kpds6=105:kpds7=10:TR=0:P1=12:P2=0:TimeU=1:10
> m above gnd:12hr fcst:NAve=0
> I have 2-meter temperature and relative humidity and 10-meter winds.
>
> - Next run plot_data_plane to plot the 10-m UGRD field:
> * /usr/local/met-8.0/bin/plot_data_plane \*
> * met-8.0/data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \*
> * UGRD_Z10.ps 'name="UGRD"; level="Z10";' -v 3*
>
> Looking in the status messages, I see:
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "UGRD/Z10" in GRIB record 163 of GRIB file
>
> And indeed wgrib confirms that record number 163 contains 10 meter U
wind.
>
> And use the "gv" or "ghostview" tool to view the resulting image:
> * gv UGRD_Z10.ps*
>
> Now, do the same steps with your ECMWF data, but you want to UWND at
100
> meters, so use:
> * 'name="UWND"; level="Z100";'*
>
> Hopefully that helps get you going.
>
> Thanks,
> John Halley Gotway
>
>
>
> On Sat, Nov 3, 2018 at 7:19 AM Allan Huffman via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Sat Nov 03 07:19:20 2018: Request 87647 was acted upon.
> > Transaction: Ticket created by awhuffma at yahoo.com
> > Queue: met_help
> > Subject: ECMWF grib file question
> > Owner: Nobody
> > Requestors: awhuffma at yahoo.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87647 >
> >
> >
> > Good Day,
> > I am seeking to run MET on an ECMWF grib file I have for 100 meter
winds.
> > However, the ECMWF grib file I have list the UWND at 100m as
var246 at
> > level sfc and VWND as var247 at level sfc.
> > Is there a way for MET to read this and properly perform the
> calculations.
> > Perhaps changing a grid table or configuration file in MET
somewhere?
> > Thanks.
> > Allan Huffman
> >
>
>
>
------------------------------------------------
More information about the Met_help
mailing list