[gradsusr] Cannot Contour Grid - all undefined values error

Wed Jun 7 13:42:32 EDT 2017

Jennifer, Bill,

This is very strange:

> gribmap -v -i pctile.2017060700 <%28201%29%20706-0700>.000_to_
024hrs.grib2.ctl
gribmap: scanning GRIB2 file: pctile.2017060700 <%28201%29%20706-0700>
.000_to_024hrs.grib2
1.1: lev1=1,0 var=0,1,8,4 ens=197
2.1: lev1=103,2 var=0,0,0,0 ens=197
3.1: lev1=103,10 var=0,2,1,0 ens=197
gribmap: Writing out the GRIB2 index file (version 3)

Why does gribmap say the file has "ens=197"?  This is not
an ensemble file.  I don't get "ens=197" and I get a match.

Bill what version of GrADS are you using?

Wesley

On Wed, Jun 7, 2017 at 12:01 PM, William Rasch - NOAA Federal <
william.rasch at noaa.gov> wrote:

> Jennifer/Wesley,
>
> The level of support you both are providing is crazy great, so thanks.
> Ugh!  But I'm still having issues!  Really strange.
>
> So here is my latest,
>
> I do the g2ctl on the grib2 file that I downloaded.  I then did a cat
> command on that control file...below...there are some lines different, not
> sure they are significant?  In the control file there is ENDVARS, Jennifer,
> yours has ENDEDEF for one.
>
> > g2ctl pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2 >
> pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2.ctl
> > cat pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2.ctl
> dset ^pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2
> index ^pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2.idx
> undef 9.999E+20
> title pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2
> * produced by g2ctl v0.1.1
> * command line options: pctile.2017060700 <%28201%29%20706-0700>.000_to_
> 024hrs.grib2
> * griddef=1:0:(131 x 66):grid_template=0:winds(N/S): lat-lon grid:(131 x
> 66) units 1e-06 input WE:SN output WE:SN res 48 lat 10.000000 to 75.000000
> by 1.000000 lon 180.000000 to 310.000000 by 1.000000 #points=8646:winds(N/S)
>
> dtype grib2
> ydef 66 linear 10.000000 1
> xdef 131 linear 180.000000 1.000000
> tdef 1 linear 00Z08jun2017 1mo
> zdef 1 linear 1 1
> vars 3
> APCPsfc   0,1,0   0,1,8,4 ** surface Total Precipitation [kg/m^2]
> TMP2m   0,103,2   0,0,0,0 ** 2 m above ground Temperature [K]
> WIND10m   0,103,10   0,2,1,0 ** 10 m above ground Wind Speed [m/s]
> ENDVARS
>
> So now I do the index command and don't get any matches...
>
> > gribmap -v -i pctile.2017060700 <%28201%29%20706-0700>.000_to_
> 024hrs.grib2.ctl
> gribmap: scanning GRIB2 file: pctile.2017060700 <%28201%29%20706-0700>
> .000_to_024hrs.grib2
> 1.1: lev1=1,0 var=0,1,8,4 ens=197
> 2.1: lev1=103,2 var=0,0,0,0 ens=197
> 3.1: lev1=103,10 var=0,2,1,0 ens=197
> gribmap: Writing out the GRIB2 index file (version 3)
>
> When I then do a:
> > grads -l
> > open pctile.2017060700 <%28201%29%20706-0700>.000_to_024hrs.grib2.ctl
> > d tmp2m
>
> I get the UNDEFINED ERROR!!!
>
> Ugh!
>
> Thanks so much if you have ANY other ideas!
>
> Bill
>
>
>
>
>
> On Tue, Jun 6, 2017 at 2:03 PM, Jennifer M Adams <jadams21 at gmu.edu> wrote:
>
>>
>> On Jun 6, 2017, at 4:07 PM, William Rasch - NOAA Federal <
>> william.rasch at noaa.gov> wrote:
>>
>> Hi Jennifer,
>>
>> I don't know what my descriptor file is.
>>
>>
>> It is the output from g2ctl. It is a text file that organizes metadata in
>> a format that GrADS (and humans) can easily understand. I copied it into a
>> file (I named it pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2.ctl)
>> like this:
>>
>> > g2ctl pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2 >
>> pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2.ctl
>>
>> And here is what is in it:
>>
>> dset ^pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2
>> index ^pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2.idx
>> undef 9.999E+20
>> title pctile.2017060500 <%28201%29%20706-0500>.000_to_024hrs.grib2
>> * produced by g2ctl v0.1.1
>> * command line options: pctile.2017060500 <%28201%29%20706-0500>
>> .000_to_024hrs.grib2
>> * griddef=1:0:(131 x 66):grid_template=0:winds(N/S): lat-lon grid:(131 x
>> 66) units 1e-06 input WE:SN output WE:SN res 48 lat 10.000000 to 75.000000
>> by 1.000000 lon 180.000000 to 310.000000 by 1.000000 #points=8646:winds(N/S)
>> dtype grib2
>> ydef 66 linear 10.000000 1
>> xdef 131 linear 180.000000 1.000000
>> tdef 1 linear 00Z06jun2017 1mo
>> zdef 1 linear 1 1
>> vars 3
>> APCPsfc   0,1,0   0,1,8,4 ** surface Total Precipitation [kg/m^2]
>> TMP2m   0,103,2   0,0,0,0 ** 2 m above ground Temperature [K]
>> WIND10m   0,103,10   0,2,1,0 ** 10 m above ground Wind Speed [m/s]
>> ENDVARS
>> EDEF 1
>> E197 1  00Z06jun2017 197
>> ENDEDEF
>>
>>
>> Then I ran gribmap on it, like this:
>>
>> >  gribmap -v -i pctile.2017060500 <%28201%29%20706-0500>.000_to_024hr
>> s.grib2.ctl
>> gribmap: scanning GRIB2 file: pctile.2017060500 <%28201%29%20706-0500>
>> .000_to_024hrs.grib2
>> 1.1: lev1=1,0 var=0,1,8,4 ens=197 valid at 2017-06-06-00:00 (t=1)  MATCH
>> 2.1: lev1=103,2 var=0,0,0,0 ens=197 valid at 2017-06-06-00:00 (t=1)
>>  MATCH
>> 3.1: lev1=103,10 var=0,2,1,0 ens=197 valid at 2017-06-06-00:00 (t=1)
>>  MATCH
>> gribmap: Writing out the GRIB2 index file (version 3)
>>
>>
>> I could open this file with GrADS and display tmp2m:
>>
>>
>>
>>
>>
>> Bill
>>
>>
>> On Tue, Jun 6, 2017 at 8:50 AM, Jennifer M Adams <jadams21 at gmu.edu> wrote
>> :
>> There is no need to apologize for being a novice. You may find this doc
>> page on how GrADS deals with GRIB data useful:
>> http://cola.gmu.edu/grads/gadoc/grib.html
>> Please see more comments inline.
>>
>>
>> > On Jun 6, 2017, at 11:34 AM, William Rasch - NOAA Federal <
>> william.rasch at noaa.gov> wrote:
>> > I grabbed one of the files, only one, then cat'ed that to the efi.master
>>
>> I would haved run g2ctl on just the one file, before cat’ing it to
>> efi.master.
>>
>> >
>> > Ran:  g2ctl  efi.master.ctl
>> >
>> > Then, did a gribmap -v efi.master.ctl
>> >
>> > It then prompted me the "descriptor file", and I entered:
>>  efi.master.ctl
>> >
>> > The feedback from that was this:
>>
>> In general, it would be better to send the actual text instead of an
>> image of it.
>>
>> >
>> > <image.png>
>> It doesn’t look you got a MATCH for any of the records in your file.
>> What is in your descriptor file?
>>
>> —Jennifer
>>
>> >
>> > I then tried to display tmp2m and got the same error.
>> >
>> > Thanks for any additional help and again, sorry for my incompetence
>> with this.
>> >
>> > Bill
>> >
>> > On Mon, Jun 5, 2017 at 6:20 PM, Jennifer M Adams <jadams21 at gmu.edu>
>> wrote:
>> > Hi, Bill -
>> > I find that the g2ctl output often needs a few tweaks to get evertyhing
>> right and some variables need a fourth grib code in the units field in
>> order to match properly. You also appear to be concatenating files with
>> different time stamps and this may be confusing the issue. Try setting up
>> your descriptor with only one file, and run gribmap in verbose mode (with
>> the -v) and look at what is matching/not matching. If records don’t match
>> what’s in your descriptor, gribmap will show you what the entry should
>> contain. Once you are matching all the records in one file, then add the
>> templating for multiple time steps. If it’s still not working please send
>> your descriptor and the gribmap verbose output.
>> > —Jennifer
>> >
>> >
>> >
>> > On Jun 5, 2017, at 11:33 AM, William Rasch - NOAA Federal <
>> william.rasch at noaa.gov> wrote:
>> >
>> > > Hi,
>> > >
>> > > I have what would seem to be a very simple data-set and "variable"
>> I'm trying to plot wiht GrADS.
>> > >
>> > > The grib files are located here...
>> > >
>> > > ftp://ftp.cdc.noaa.gov/Public/gbates/reforecast2/efi_grib2
>> /pctile.2017060500.000_to_024hrs.grib2
>> > >
>> > > I am simply downloading several of these files, then do this:
>> > >
>> > > cat * > efi.master
>> > > g2ctl efi.master > efi.master.ctl
>> > > gribmap -i efi.master.ctl
>> > >
>> > > I'm simply trying to plot either:
>> > >
>> > > temp2m
>> > > wind10m
>> > > apcpsfc
>> > >
>> > > When I try to display the data, I get the error:
>> > >
>> > > Cannot Contour Grid - all undefined values
>> > >
>> > > Thanks for any help....
>> > >
>> > > Bill
>> > > _______________________________________________
>> > > gradsusr mailing list
>> > > gradsusr at gradsusr.org
>> > > http://gradsusr.org/mailman/listinfo/gradsusr
>> >
>> > --
>> > Jennifer Miletta Adams
>> > Center for Ocean-Land-Atmosphere Studies (COLA)
>> > George Mason University
>> >
>> >
>> >
>> >
>> > _______________________________________________
>> > gradsusr mailing list
>> > gradsusr at gradsusr.org
>> > http://gradsusr.org/mailman/listinfo/gradsusr
>> >
>> >
>> >
>> > --
>> >
>> > William Rasch
>> > Science and Operations Officer
>> > Weather Forecast Office
>> > Sacramento, CA
>> > _______________________________________________
>> > gradsusr mailing list
>> > gradsusr at gradsusr.org
>> > http://gradsusr.org/mailman/listinfo/gradsusr
>>
>> --
>> Jennifer Miletta Adams
>> Center for Ocean-Land-Atmosphere Studies (COLA)
>> George Mason University
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>>
>> --
>>
>> William Rasch
>> Science and Operations Officer
>> Weather Forecast Office
>> Sacramento, CA
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>> --
>> Jennifer Miletta Adams
>> Center for Ocean-Land-Atmosphere Studies (COLA)
>> George Mason University
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> <https://twitter.com/NWSSacramento>
> <https://www.facebook.com/NWS.Sacramento>
> <http://www.youtube.com/nwssacramento>
> William Rasch
> Science and Operations Officer
> Weather Forecast Office
> Sacramento, CA
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170607/38bb3894/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pctile.2017060500.000_to_024hrs.grib2.png
Type: image/png
Size: 166929 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20170607/38bb3894/attachment-0001.png 