[gradsusr] Cannot Contour Grid - all undefined values error

Wed Jun 7 14:37:12 EDT 2017

Wesley,

Can't say it enough, thanks again for helping.

When I start Grads here is what is says:

So I guess Version 2.1.a3.p4

Grid Analysis and Display System (GrADS) Version 2.1.a3.p4
Copyright (C) 1988-2015 by the Institute for Global Environment and Society
(IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v2.1.a3.p4 little-endian readline grib2 netcdf hdf4-sds hdf5
opendap-grids,stn geotiff shapefile cairo
Issue 'q config' command for more detailed configuration information
GX Package Initialization: Size = 11 8.5
Running in Batch mode

On Wed, Jun 7, 2017 at 10:42 AM, Wesley Ebisuzaki - NOAA Federal <
wesley.ebisuzaki at noaa.gov> wrote:

> Jennifer, Bill,
>
> This is very strange:
>
> > gribmap -v -i pctile.2017060700 <%28201%29%20706-0700>.000_to_024hr
> s.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_024hr
>> s.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
>>
>>
>
> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170607/f4205bc1/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/f4205bc1/attachment-0001.png 