[gradsusr] Help scripts for WRF output in grads

Zilore Mumba zmumba at gmail.com
Wed Jun 13 08:22:51 EDT 2012 

Hello,
I have WRF model output converted to binary for use in grads.
Apparently WRF does not give humidity and temperature at different
levels, except T2, temperature at 2m.
I have to calculate rh at each level from
es = 6.112*EXP(17.67*(t-273.15)/(t-29.65))
e = qv * p/100./(qv + 0.622)
rh = e/es
But before that I have to get the temperature from potential temperature using
temperature = (T2 + 300) [(P/P0) ^ (R/Cp)]
I was wondering if there are users who visualise WRF output in grads,
and if they have scripts for the above conversions.
Thanks
Zilore
On Wed, Jun 13, 2012 at 1:27 AM,  <gradsusr-request at gradsusr.org> wrote:
> Send gradsusr mailing list submissions to
>        gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>        http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
>        gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
>        gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
>   1. Re: Problem posting nps grib2 data in grads (Wesley Ebisuzaki)
>   2. Re: How I can enlarge the font size of lat,       lon and legend
>      text? (Jayakrishnan P.R)
>   3. Re: How I can enlarge the font size of lat,       lon and legend
>      text? (Mekonnen Adnew)
>   4. Re: How I can enlarge the font size of lat,       lon and legend
>      text? (Jayakrishnan P.R)
>   5. Re: Fwd: dear sir.... (Kishore Babu)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 11 Jun 2012 11:26:48 -0400
> From: Wesley Ebisuzaki <wesley.ebisuzaki at noaa.gov>
> Subject: Re: [gradsusr] Problem posting nps grib2 data in grads
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
>        <CAMHOGenefM6pUGhwALCqSV-a1=BTQwJW0Op4uK79yYhQA+Wk2Q at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Ron,
>
> g2ctl can only handle grib files with one grid.  So you will have to
> separate the grids into their own files.
>
>  wgrib2  reg_00_024.grib2 -grid -one_line | grep  "493 x 399" | wgrib2 -i
> reg_00_024.grib2 -grib grid1
>  wgrib2  reg_00_024.grib2 -grid -one_line | grep  -v "493 x 399" | wgrib2
> -i reg_00_024.grib2 -grib grid2
>
> However, once you made your ctl files, you can modify the "dset" to point
> to the original file.
> gribmap is smart enough to handle grib files with different grids.
>
>           Wesley
>
>
>
> On Fri, Jun 8, 2012 at 4:14 PM, Goodson,Ron [Edm] <Ron.Goodson at ec.gc.ca>wrote:
>
>> **
>> Hi
>>
>> I already mentioned to Lewis (we work for the same organization) that his
>> problem is that his grib-file contains 2 different grids (depending upon
>> variable).
>>
>> I answer him to the entire list, though, in case somebody knows how to
>> handle this situation (other than breaking the file into two files .. each
>> with their own grid).
>>
>> tks
>>
>> ron
>>
>>  ------------------------------
>> *From:* gradsusr-bounces at gradsusr.org [mailto:
>> gradsusr-bounces at gradsusr.org] *On Behalf Of *Poulin,Lewis [CMC]
>> *Sent:* June 8, 2012 9:40 AM
>> *To:* gradsusr at gradsusr.org
>> *Subject:* [gradsusr] Problem posting nps grib2 data in grads
>>
>>  Hello
>>
>> I am having difficulty in posting in GrADS my grib2 data on a north polar
>> stereographic grid
>>
>> I placed all relevent files and information in the directory at:
>> *http://collaboration.cmc.ec.gc.ca/cmc/CMOI/grads/regional/*<http://collaboration.cmc.ec.gc.ca/cmc/CMOI/grads/regional/>
>>
>> From this site you can see my grib2 file: called reg_00_024.grib2
>>
>> I am using grads-2.0.a5
>>
>> I made the following ctl and idx files (available that above site)
>>   reg_00_024.grib2.ctl
>>   reg_00_024.grib2.idx
>>
>> When I run the following simple script:
>>
>> *start sample script
>>
>> 'open reg_00_024.grib2.ctl'
>>
>> 'clear'
>> 'set lon 230 310'
>> 'set lat 30 70'
>> 'set mpvals 235 305 35 65'
>> 'set mproj nps'
>>
>> 'set lev 850'
>> 'd HGTprs'
>>
>> 'set lev 500'
>> 'd TMPprs'
>>
>> *end sample script
>>
>>
>> I get the following messages:
>>
>> Notice: Implied interpolation for file reg_00_024.grib2.ctl
>>  Interpolation will be performed on any data displayed from this file
>> Cannot contour grid - all undefined values
>> Cannot contour grid - all undefined values
>>
>>
>> I cannot get my data to display.
>>
>> Thanks in advance for your suggestions and assistance.
>>
>> Lewis
>>
>> _______________________________________________
>> 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/20120611/54fa22f0/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Mon, 11 Jun 2012 21:05:59 +0530
> From: "Jayakrishnan P.R" <prjayakrishnan at gmail.com>
> Subject: Re: [gradsusr] How I can enlarge the font size of lat, lon
>        and legend text?
> To: Mekonnen Adnew <mekonnenadnew at yahoo.com>,   GrADS Users Forum
>        <gradsusr at gradsusr.org>
> Message-ID:
>        <CA+P5ZEEpAWj3Gz3MoO5J-oeWGON+rC4r2Ai7bo1UjScaK_37wQ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Adnew,
>             I think you can set the size options for x and y axis by these
> commands...
>
> 'set xlopts 1 4 0.18'
> 'set ylopts 1 4 0.18'
>
> Try this...Best of Luck....
>
> --
> Sincerely
> ***********************************************
> Jayakrishnan.P.R
> CSIR-Senior Research Fellow
> Department of Atmospheric Sciences
> Cochin University of Science and Technology (CUSAT), Cochin-682 016
> Kerala, India.
> Mob: 09895417565
>
> ***********************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120611/c3c4666e/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Tue, 12 Jun 2012 06:07:49 -0700 (PDT)
> From: Mekonnen Adnew <mekonnenadnew at yahoo.com>
> Subject: Re: [gradsusr] How I can enlarge the font size of lat, lon
>        and legend text?
> To: "Jayakrishnan P.R" <prjayakrishnan at gmail.com>,
>        "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
> Message-ID:
>        <1339506469.94710.YahooMailNeo at web121205.mail.ne1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Jayakrishnan:
> Thank you for your rapid and successful response to my question. I have successfully able to enlarge the font size of lat and lon in Grads. But I need also your assistant to enlarge the font size of the lagend text. Thank you very much.
> ?
> Best regards,
> Mekonnen Adnew Degefu
>
>
> ________________________________
>  From: Jayakrishnan P.R <prjayakrishnan at gmail.com>
> To: Mekonnen Adnew <mekonnenadnew at yahoo.com>; GrADS Users Forum <gradsusr at gradsusr.org>
> Sent: Monday, June 11, 2012 6:35 PM
> Subject: Re: [gradsusr] How I can enlarge the font size of lat, lon and legend text?
>
>
> Dear Adnew,
> ? ? ? ? ? ? ?I think you can set the size options for x and y axis by these commands...
>
> 'set xlopts 1 4 0.18'
> 'set ylopts 1 4 0.18'
>
> Try this...Best of Luck....
> --
> Sincerely
> ***********************************************
> Jayakrishnan.P.R
> CSIR-Senior Research Fellow
> Department of Atmospheric Sciences
> Cochin University of Science and Technology (CUSAT), Cochin-682 016
> Kerala, India.
> Mob: 09895417565
>
> ***********************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120612/9152857e/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Tue, 12 Jun 2012 21:19:26 +0530
> From: "Jayakrishnan P.R" <prjayakrishnan at gmail.com>
> Subject: Re: [gradsusr] How I can enlarge the font size of lat, lon
>        and legend text?
> To: Mekonnen Adnew <mekonnenadnew at yahoo.com>,   GrADS Users Forum
>        <gradsusr at gradsusr.org>
> Message-ID:
>        <CA+P5ZEFjEG4a+p+hk0WHO9ufoC3YSXr97XZDKU2SvOU7ftwUfw at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> *You can refer this.....**set digsiz*
>
> set digsiz *size*
>
> Size (in inches, or plotter units) of the numbers. 0.1 to 0.15 is usually a
> good range to use. Both of these options stay the same until changed.
>
> For further reference kindly use the following link....Best of luck...
>>
>
>
> http://www.iges.org/grads/gadoc/gadocindex.html
>
>
> --
> Sincerely
> ***********************************************
> Jayakrishnan.P.R
> CSIR-Senior Research Fellow
> Department of Atmospheric Sciences
> Cochin University of Science and Technology (CUSAT), Cochin-682 016
> Kerala, India.
> Mob: 09895417565
>
> ***********************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120612/b316f42b/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Wed, 13 Jun 2012 10:57:32 +0530
> From: Kishore Babu <kishoreragi at gmail.com>
> Subject: Re: [gradsusr] Fwd: dear sir....
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
>        <CAN4rXbRVysq0VCMLR-rNFaWTQN-wP_Se7_ky2N4CY4XxpvgKWQ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear sir,
>
> How can I use scale factor and add offset to netcdf file? Could you please
> give some reference to get the info.
>
> Thank you sir,
>
> Sincerely,
>
> Kishore
>
> On Sun, Jun 10, 2012 at 1:51 PM, Jayakrishnan P.R
> <prjayakrishnan at gmail.com>wrote:
>
>> Dear Kishore,
>>                  I will be able to give you some hint.....You can try
>> something like this....
>> 'reinit'
>> 'open trmm02.ctl'
>> 'set lat -30 30'
>> 'set lon 30 110'
>> 'define a=ave(rr,t=1,t=31,8)'
>> prompt 'key-in file name for rain :'
>> pull fname
>> 'set fwrite 'fname
>> 'set gxout fwrite'
>> 'd a'
>> 'disable fwrite'
>>
>>
>> This I tried with some errors......This is in the case of TRMM...we have 3
>> hourly rainfall....So what you have to do is to set the time once in 8
>> intervals...ie for getting 00 hour average you will have to get 00 hour
>> data only...so t=1,t=8 is given for time averaging....This is the rough
>> idea.....Hope you will be able to develop from this code....Best of luck....
>>
>>
>>
>> --
>> Sincerely
>> ***********************************************
>> Jayakrishnan.P.R
>> CSIR-Senior Research Fellow
>> Department of Atmospheric Sciences
>> Cochin University of Science and Technology (CUSAT), Cochin-682 016
>> Kerala, India.
>> Mob: 09895417565
>>
>> ***********************************************
>>
>> _______________________________________________
>> 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/20120613/1e631991/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 28, Issue 12
> ****************************************

More information about the gradsusr mailing list