[gradsusr] UKMET

Fri Mar 9 14:51:15 EST 2012

Thanks.

Jennifer/Wesley,

I tried to break up a UKMET grib file, as per tricks using wgrib (ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib/tricks.wgrib) using #15 (separating grids by GDS).  So, I did the command:

wgrib -GDS fh.0006_tl.precp and I got the result below:

1:40:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff00490000000249f080015f9081d4c0ffff04e240000000000049004900490049004900490049004900480048004800470047004700460046004500450044004300430042004100410040003f003e003d003c003c003b003a003900380037003600340033003200310030002f002d002c002b002a0028002700260024002300210020001e001d001c001a0019001700160014001300110010000e000c000b000900080006000500030002
2:3759:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff0049815f9000ea60800000000249f0ffff04e24000000000000200030005000600080009000b000c000e0010001100130014001600170019001a001c001d001e0020002100230024002600270028002a002b002c002d002f003000310032003300340036003700380039003a003b003c003c003d003e003f0040004100410042004300430044004500450046004600470047004700480048004800490049004900490049004900490049
3:7478:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff0049815f900249f08000000081d4c0ffff04e24000000000000200030005000600080009000b000c000e0010001100130014001600170019001a001c001d001e0020002100230024002600270028002a002b002c002d002f003000310032003300340036003700380039003a003b003c003c003d003e003f0040004100410042004300430044004500450046004600470047004700480048004800490049004900490049004900490049
4:11629:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff0049815f9081d4c080000000807530ffff04e24000000000000200030005000600080009000b000c000e0010001100130014001600170019001a001c001d001e0020002100230024002600270028002a002b002c002d002f003000310032003300340036003700380039003a003b003c003c003d003e003f0040004100410042004300430044004500450046004600470047004700480048004800490049004900490049004900490049
5:15348:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff004900000000ea6080015f900249f0ffff04e240000000000049004900490049004900490049004900480048004800470047004700460046004500450044004300430042004100410040003f003e003d003c003c003b003a003900380037003600340033003200310030002f002d002c002b002a0028002700260024002300210020001e001d001c001a0019001700160014001300110010000e000c000b000900080006000500030002
6:19067:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff004900000080753080015f9000ea60ffff04e240000000000049004900490049004900490049004900480048004800470047004700460046004500450044004300430042004100410040003f003e003d003c003c003b003a003900380037003600340033003200310030002f002d002c002b002a0028002700260024002300210020001e001d001c001a0019001700160014001300110010000e000c000b000900080006000500030002
7:23218:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff004900000081d4c080015f90807530ffff04e240000000000049004900490049004900490049004900480048004800470047004700460046004500450044004300430042004100410040003f003e003d003c003c003b003a003900380037003600340033003200310030002f002d002c002b002a0028002700260024002300210020001e001d001c001a0019001700160014001300110010000e000c000b000900080006000500030002
8:27369:d=12030900:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-6hr acc:NAve=2:GDS=0000b2002100ffff0049815f908075308000000000ea60ffff04e24000000000000200030005000600080009000b000c000e0010001100130014001600170019001a001c001d001e0020002100230024002600270028002a002b002c002d002f003000310032003300340036003700380039003a003b003c003c003d003e003f0040004100410042004300430044004500450046004600470047004700480048004800490049004900490049004900490049

So, to parse everything but the first entry, I did the trick command:

wgrib fh.0006_tl.precp -GDS | grep "GDS=0000b2002100ffff00490000000249f080015f9081d4c0ffff04e240000000000049004900490049004900490049004900480048004800470047004700460046004500450044004300430042004100410040003f003e003d003c003c003b003a003900380037003600340033003200310030002f002d002c002b002a0028002700260024002300210020001e001d001c001a0019001700160014001300110010000e000c000b000900080006000500030002" | wgrib fh.0006_tl.precp -I -grib -o UKMET.grib

where UKMET.grib is the name of the new grib file

It completes the command, but says "bad input from stdin".

The resulting grib file does only contain this grid, but when I run grib2ctl on the file, it looks like this:

dset ^UKMET.grib
index ^UKMET.grib.idx
undef 9.999E+20
title UKMET.grib
*  produced by grib2ctl v0.9.12.5p32amec
dtype grib 40
tdef 1 linear 06Z09mar2012 1mo
zdef 1 linear 1 1
vars 1
APCPsfc  0 61,1, 0  ** surface Total precipitation [kg/m^2]
ENDVARS

So, it doesn't have any X/Y coordinates.  Any ideas?

Thanks

Brian

From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Jennifer Adams
Sent: Thursday, January 19, 2012 7:18 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] UKMET

I grabbed a sample file:
ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/RD.20120119/PT.grid_DF.bb/fh.0006_tl.press

There are 3 different variables (grib codes 6, 7, and 102). Each variable has 8 records, each with a different grid, except for variable 102, which has 9. The ninth grid for variable 102 is a latlon grid with 1813 points, all the others are thinned grids with nxny=3447. Grib2ctl creates a pdef file for the thinned grid, presumably based on the grid specs it finds in record #1. I ran wgrib -V on this file and then grepped on 'grid=' and got this:

rec 1:00040:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=44 tropopause 6hr fcst:
rec 2:05915:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=41 tropopause 6hr fcst:
rec 3:11790:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=43 tropopause 6hr fcst:
rec 4:17233:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=37 tropopause 6hr fcst:
rec 5:23108:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=39 tropopause 6hr fcst:
rec 6:28983:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=42 tropopause 6hr fcst:
rec 7:34858:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=38 tropopause 6hr fcst:
rec 8:40733:date 2012011912 PRES kpds5=1 kpds6=7 kpds7=0 levels=(0,0) grid=40 tropopause 6hr fcst:

When checking for matches for PDEF'd grids, gribmap will check the nxny value to make sure it matches what the PDEF entry says, but doesn't go beyond that. The one latlon grid in the file does spawn an error from gribmap, "GRIB grid size does not match descriptor". Each of the 8 thinned grids will appear identical to gribmap, so they will all match, and the one that the index file will point to will be the last record in the file to match, which will not likely correspond to the information used to create the pdef file. Wesley's suggestion to break up the file and separate out the nine different grids is the only way this data can be read with GrADS.
--Jennifer

On Dec 16, 2011, at 5:06 PM, Wesley Ebisuzaki wrote:

Brian,

Looks like a problem with gribmap in GrADS v2.

The files have several different grids in them. GrADS should check that
the grid
matches the grid in the ctl file (it does with grib1). It doesn't so all
the grids match
and you can plot the wrong data. :(

Options: grib1, delete all but one of the grids.

18.1: lev1=101,0 var=0,3,1 MATCH
19.1: lev1=101,0 var=0,3,1 MATCH
20.1: lev1=101,0 var=0,3,1 MATCH
21.1: lev1=101,0 var=0,3,1 MATCH
22.1: lev1=101,0 var=0,3,1 MATCH
23.1: lev1=101,0 var=0,3,1 MATCH
24.1: lev1=101,0 var=0,3,1 MATCH
25.1: lev1=101,0 var=0,3,1 MATCH
grib2map: Writing out the index file
Wesley

Walsh, Brian J wrote:

Hi,

I'm wondering if anyone has been able to plot the UKMET grib data
available on the TGFTP ftp site
(ftp://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/MT.bmoav_CY.12/) using
GrADS. It is a thinned lat-long grid and I am able to convert it to a
GRIB2 file using cnvgrib and run wgrib2 on the files, but it's not
plotting. Just wanted to see if anyone had success plotting these
files and how it was accomplished.

Thanks

Brian

The information contained in this e-mail is intended only for the
individual or entity to whom it is addressed.
Its contents (including any attachments) may contain confidential
and/or privileged information.
If you are not an intended recipient you must not use, disclose,
disseminate, copy or print its contents.
If you receive this e-mail in error, please notify the sender by reply
e-mail and delete and destroy the message.
------------------------------------------------------------------------

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org<mailto:jma at cola.iges.org>

________________________________
The information contained in this e-mail is intended only for the individual or entity to whom it is addressed.
Its contents (including any attachments) may contain confidential and/or privileged information.
If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents.
If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120309/cc329c49/attachment-0003.html 