[gradsusr] Decoding ECMWF Data

Sat Aug 10 19:14:06 EDT 2013

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130810/517184e7/attachment-0003.html 