[gradsusr] Decoding ECMWF Data

Sun Aug 11 14:32:01 EDT 2013

Hi,

   grib2ctl and g2ctl assume a vertical structure like NCEP
pgb files.  That is that there are no missing levels from the
lowest level to an upper level.  If you have missing levels,
then you should use the -no_profile option to grib2ctl/g2ctl.
That gives you a variable for each field and level. That
will keep you from getting confused like "Z is defined on 500
while U is defined on 200 and 850 and RH is defined on
750".

   Wesley

On Sat, Aug 10, 2013 at 7:14 PM, Andrew Revering <andy at f5data.com> wrote:

> Mohsen,
>
> I appreciate your help. However, I'm not convinced this is working
> properly. The other variable that goes with this problem is vprs... and if
> I do as instructed, I get undefined values when I try to display vprs. I
> feel like by changing the levels from 2 to 3 is just shifting the vprs for
> 850 to be extracted when I set the level to 500. Know what I mean?
>
> I should be able to get both 850 uprs & vprs as well as 500 uprs & vprs,
> but can't seem to get them all. yet again, if I were to open this grib1
> file in tkdegrib, I can see all of that data.
>
> Very strange. Its something in how the file is packed and how grib2ctl.plis reading the data when creating the ctl file. I greatly appreciate your
> assistance, and welcome any other thoughts you might have.
>
> Regards,
>
> Andrew Revering - Follow @andrewrevering<https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>on Twitter
> Convective Development, Inc.
> PO Box 511, Cambridge, MN 55008
> http://www.convectivedevelopment.com/
>
>
> On Sat, Aug 10, 2013 at 4:03 AM, Mohsen Soltani <soltani.clima at gmail.com>wrote:
>
>> Hello Andrew,
>>
>> I examined your ctl. and found that the only thing you need to do is to
>> change a number from 2 into 3 in the ctl file!
>> *****************
>> that is:
>> Uprs 2 131,100,0 ** (profile) U velocity [m s**-1]
>> should be changed into:
>> Uprs 3 131,100,0 ** (profile) U velocity [m s**-1]
>> ***********************************************
>> And then you will be able to display "uprs" in both 850-hPa and 500-hPa
>> levels. Therefore, do as follows:
>> *******************
>> open .....ctl
>> set lev 850
>> d uprs
>> c
>> set lev 500
>> d uprs
>> ************************
>> The attached are the outputs of each level.
>>
>>
>> Mohsen
>>
>> --
>> April showers bring May flowers!
>> --
>> *Mohsen Soltani, M.Sc.*
>> Climatology Grad Student,
>> Department of Climatology, Faculty of Geography,
>> University of Tehran, Tehran, Iran
>> voice: (+98) 9119772934
>> e-mail: soltani.clima at gmail.com
>>
>>
>> On Fri, Aug 9, 2013 at 11:28 PM, Andrew Revering <andy at f5data.com> wrote:
>>
>>> I appreciate the help, but I just don't get it.
>>> I played around with the CTL file and changing some values but couldn't
>>> get anything other than what appears to be 850mb data to display.
>>>
>>> Any other help is appreciated.
>>>
>>> Andrew Revering - Follow @andrewrevering<https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>on Twitter
>>> Convective Development, Inc.
>>> PO Box 511, Cambridge, MN 55008
>>> http://www.convectivedevelopment.com/
>>>
>>>
>>> On Wed, Aug 7, 2013 at 3:13 AM, Becker, Bernd <
>>> bernd.becker at metoffice.gov.uk> wrote:
>>>
>>>> **
>>>> try the following nuonces in the way you describe Z in the .ctl files:
>>>>
>>>>
>>>> Multilevel data:
>>>>
>>>> zdef 6 levels
>>>> 1000 850 700 500 200 10
>>>> tdef 33 linear 00Z18Jul2013 1dy
>>>> vars 1
>>>> Z         6  129,100,    0 ** geopotential             gsmtercboundary2
>>>> [gm]
>>>>
>>>> and single level
>>>>
>>>> z     0 129, 99 ,0   ** 500 hPa Geopotential
>>>> gsmobstercboundary2 [m**2 s**-2]
>>>>
>>>> and
>>>>
>>>> z     0 129, 100,500   ** 500 hPa Geopotential
>>>>
>>>> good luck!
>>>> Bernd
>>>>
>>>>
>>>> --
>>>> Bernd Becker   Climate Impacts Product Development
>>>> Met Office Hadley Centre    FitzRoy Road  Exeter   Devon EX1 3PB
>>>> United Kingdom
>>>> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050
>>>> E-mail:bernd.becker at metoffice.gov.uk -
>>>> http://www.metoffice.gov.uk/research/people/bernd-becker
>>>>
>>>>
>>>>  ------------------------------
>>>> *From:* gradsusr-bounces at gradsusr.org [mailto:
>>>> gradsusr-bounces at gradsusr.org] *On Behalf Of *Andrew Revering
>>>> *Sent:* 07 August 2013 06:21
>>>> *To:* GrADS User List
>>>> *Subject:* [gradsusr] Decoding ECMWF Data
>>>>
>>>>  I'm having a hard time with *some* ECMWF data that I'm subscribed to.
>>>> I am using grads (opengrads), gribmap & grib2ctl.pl to generate the
>>>> CTL & IDX files. The file format is Grib1.
>>>> - Using grib2ctl version 0.9.12.6
>>>> - Using gribmap "for grads version 2.0.1.oga.1"
>>>> - Using wgrib 1.8.1.0b
>>>>
>>>>  The problem exists with 'level' data, but not all level data. In the
>>>> example below, I can map the 850mb winds, but not the 500mb.
>>>> Try
>>>>
>>>>  set lev 500
>>>> d uprs
>>>>
>>>>  as an example. I get an 'entire grid undefined' message.
>>>>
>>>>  I was able to open the data using the tkdegrib package and output
>>>> uprs @ 500mb level in CSV format, and see that indeed the data is in there
>>>> and looks good. Just can't visualize it in Grads for some reason.
>>>>
>>>>  Here's a temporary link to example files:
>>>> http://www.f5data.net/temp/
>>>>
>>>>  Any thoughts?
>>>>
>>>> Andrew Revering - Follow @andrewrevering<https://twitter.com/intent/follow?screen_name=andrewrevering&source=email&variant=2.0&xd_token=ceae5c40173e5>on Twitter
>>>> Convective Development, Inc.
>>>> PO Box 511, Cambridge, MN 55008
>>>> http://www.convectivedevelopment.com/
>>>>
>>>> _______________________________________________
>>>> 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
>>>
>>>
>>
>> _______________________________________________
>> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130811/8d69b7a6/attachment-0003.html 