[gradsusr] Help invoking user function
Zilore Mumba
zmumba at gmail.com
Tue May 22 01:06:09 EDT 2012
Help will be appreciated how to invoke a user function. I have to put
time of validity on model output. In the script attached my BaseDate
is 19/05/2012 which the script takes from "set t 1". Now I want to put
validity time of a forecast valid at 06Z on 20th which I want to take
from "set t 6".
My script is not able to take the result of t 6, instead it print
valid to 06Z on / /.
Please assist.
Zilore Mumba
On Mon, May 21, 2012 at 3:52 PM, <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. wgrib2 (Jeff Lake)
> 2. Re: Reg: re.gs file (Arlindo da Silva)
> 3. Re: Building 2.0.1 on Ubuntu 12.04 (hdf4, netcdf, shapefile
> not detected) (Arlindo da Silva)
> 4. Re: g2grb subregion (Mike Manning)
> 5. Re: Stippling/cross hatching for south polar (KODAMA Chihiro)
> 6. Area average only for California (rjmera)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 18 May 2012 14:16:23 -0400
> From: Jeff Lake <admin at michiganwxsystem.com>
> Subject: [gradsusr] wgrib2
> To: Grads <gradsusr at gradsusr.org>
> Message-ID: <4FB691F7.4010105 at michiganwxsystem.com>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> yes this is a GrADS list, but those of you that have experience with wgrib2
> I have a problem with a simple snippet ..
> Attempting to get a 12hr fcst_ave from the NWS NDFD grib files
> following the directions on the wgrib2 pages, this looks like it should work
>
> wgrib2 -vt -var -lev -misc ds.td.bin -fcst_ave 12hr newdew.bin
> but its not doing anything but making an exact copy of the ds.td.bin
>
> anybody used wgrib2 like this so I can use GrADS to make 12hr fcsts and
> not 3hr
>
> --
> -Jeff Lake
> MichiganWxSystem.com
> WeatherMichigan.net
> TheWeatherCenter.net
> GRLevelXStuff.com
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 19 May 2012 17:33:27 -0400
> From: Arlindo da Silva <dasilva at alum.mit.edu>
> Subject: Re: [gradsusr] Reg: re.gs file
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CALz7xZcRi-u7keX1hFqrSGo4UuwdPGF5mFo3t9_zZEAjw7kJvw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> On Mon, May 14, 2012 at 7:34 AM, ?? ??????? ?????? <nelli.narendra at gmail.com
>> wrote:
>
>> Dear friends
>>
>> Can any one send me re.gs file for regrid the rawdata in Opengrads.
>>
>>
>
> I am not sure what you mean by "raw data". Now, re is not a script but an
> opengrads extension (written in C). If you have the opengrads bundle it
> should be included. Try this:
>
> ga-> q udf
>
> Do you see re() list? Documentation is here:
>
> http://opengrads.org/doc/udxt/re/
>
> Arlindo
>
>
>
>
>
>
>
>> Thanks in advance
>>
>> Narendra Reddy
>>
>> --
>> Regards,
>>
>> N.Narendra Reddy
>> Research Scholar,
>> C/O Dr.Kusuma G Rao
>> Atmospheric Science Programme
>> Antariksh Bhavan, NEW BEL Road
>> ISRO HQ, B'lore, 560 094.
>> Office: +91 80 2217 2091, 2077
>> Home: +91 80 42035437
>> Mobile: 09986589807 (*18:30 - 21:30 IST*)
>> nelli.narendra at gmail.com
>> ************************************************************************
>> ???? ???? ???? ????? ???? ?????? ???? ??????? ????? ....
>>
>> *************************************************************************
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> Arlindo da Silva
> *dasilva at alum.mit.edu*
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120519/c1c7e26b/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Sat, 19 May 2012 17:42:30 -0400
> From: Arlindo da Silva <dasilva at alum.mit.edu>
> Subject: Re: [gradsusr] Building 2.0.1 on Ubuntu 12.04 (hdf4, netcdf,
> shapefile not detected)
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CALz7xZcj0dui5PknWX=1BQNudV9LbPzj8yHQPk2h1yq9gUOebA at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> On Sun, May 13, 2012 at 10:38 PM, klo uo <klonuo at gmail.com> wrote:
>
>> OK, sorry for previous message, but I'll take smarter route and
>> install what's available in repository
>> Building whole grads from source on Debian/Ubuntu is not trivial task (for
>> me)
>>
>>
> Building GrADS is not trivial even for us who have been doing this for
> years. GrADS itself builds very easily, the problem are those pesky
> dependencies...
>
> Arlindo
>
> --
> Arlindo da Silva
> *dasilva at alum.mit.edu*
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120519/af8b0571/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Mon, 21 May 2012 09:03:52 +1000
> From: Mike Manning <michael at bsch.au.com>
> Subject: Re: [gradsusr] g2grb subregion
> To: gradsusr at gradsusr.org
> Message-ID: <4FB97858.2030806 at bsch.au.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Thanks Wesley,
>
> I'll give that a try this week :)
>
> Cheers, Mike
>
>
> On 19/05/2012 1:07 AM, Wesley Ebisuzaki wrote:
>> Mike,
>>
>> The error message suggests that you haven't installed all the files
>> needed
>> for g2grb.gs <http://g2grb.gs> to run. You need
>>
>> template.grb2 in the $GADDIR directory
>> you need wgrib2 (v1.8.4+) on your path
>>
>>
>> BTW The simpler way of making regional grib2 subset is to use wgrib2.
>>
>> For lat-lon -> subset lat-lon: use the -small_grib option
>>
>> wgrib2 IN.grb2 -set_grib_type j -small_grib 100:170 -55:0 OUT.grb2
>>
>> You can also interpolate to a lat-lon grid by the -new_grid option.
>>
>> Yours,
>> Wesley
>>
>>
>>
>> On Fri, May 18, 2012 at 3:19 AM, Mike Manning <michael at bsch.au.com
>> <mailto:michael at bsch.au.com>> wrote:
>>
>> Hi all,
>>
>> I've got data downloaded from an archive now however I'd like to
>> extract
>> a subregion from the grib2 files using g2grb. I've created a test.ctl
>> files as follows:
>>
>> 'open gfs.00z.ctl'
>> 'set lat -55 0'
>> 'set lon 100 170'
>> 'g2grb LFTXsfc extrap.grb UGRD:1000 mb:anl:packing=j:center=7'
>> #'g2grb h(lev=925) +extrap.grb HGT:925 mb:anl:packing=j:center=7'
>> #'g2grb h(lev=850) +extrap.grb HGT:850 mb:anl:packing=j:center=7'
>> #'g2grb h(lev=700) +extrap.grb HGT:700 mb:anl:packing=j:center=7'
>> #'g2grb h(lev=600) +extrap.grb HGT:600 mb:anl:packing=j:center=7'
>> #'g2grb h(lev=500) +extrap.grb HGT:500 mb:anl:packing=j:center=7'
>> quit
>>
>>
>> I run it using grads -blc however it's throwing up some errors:
>> Starting "/grads2/Linux/Versions/2.0.a9.oga.1/i686/grads -blc
>> test.ctl
>> " ...
>>
>>
>> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1
>> Copyright (c) 1988-2010 by Brian Doty and the
>> Institute for Global Environment and Society (IGES)
>> GrADS comes with ABSOLUTELY NO WARRANTY
>> See file COPYRIGHT for more information
>>
>> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf
>> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile
>> Issue 'q config' command for more detailed configuration information
>> Loading User Defined Extensions table
>> </grads2/Linux/Versions/2.0.a9.oga.1/i686/gex/udxt> ... ok.
>> GX Package Initialization: Size = 11 8.5
>> Running in Batch mode
>> dlon=1 nx = 71 lon0 = 100
>> Y is varying Lat = -55 to 0 Y = 36 to 91
>> dlat=1 ny = 56 lat0 = -55
>> need to generate date
>> new inv=d=2011121000:UGRD:1000 mb:anl:packing=j:center=7
>> new grads found undef=-999000000.000000
>> making template file
>> ! /grads2/wgrib2 $GADDIR/template.grb2 -d 1 -lola 100:71:1 -55:56:1
>> extrap.grb.template grib
>> argument: -lola ????
>> argument: 100:71:1 ????
>> argument: -55:56:1 ????
>> argument: extrap.grb.template ????
>> argument: grib ????
>> missing GRIB record(s)
>> opening fwrite extrap.grb.tmp
>> going to display LFTXsfc
>> result is Wrote 3976 of 3976 elements to extrap.grb.tmp as Sequential
>> Big_Endian
>>
>> metadata = 1:0:d=2011121000:UGRD:1000 mb:anl:packing=j:center=7:
>> ! /grads2/wgrib2 extrap.grb.template -import_ieee extrap.grb.tmp
>> -set_metadata extrap.grb.inv -undefine_val -999000000.000000 -grib_out
>> extrap.grb
>> argument: -import_ieee ????
>> argument: extrap.grb.tmp ????
>> argument: -set_grib_type ????
>> argument: c3 ????
>> argument: -set_metadata ????
>> argument: extrap.grb.inv ????
>> argument: -undefine_val ????
>> argument: -999000000.000000 ????
>> argument: -grib_out ????
>> argument: extrap.grb ????
>> could not open file: extrap.grb.template
>>
>> !rm extrap.grb.tmp extrap.grb.inv extrap.grb.template
>> rm: cannot remove `extrap.grb.template': No such file or directory
>>
>>
>> I'm not sure why it's not working :( All I'd like to do is extract the
>> subregion from lat -55 to 0 and lon 100 to 170 for all variables
>> in the
>> grib data files I've downloaded.
>>
>> Cheers, Mike
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org <mailto: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/20120521/21f56adc/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Mon, 21 May 2012 18:37:30 +0900
> From: KODAMA Chihiro <kodamac at jamstec.go.jp>
> Subject: Re: [gradsusr] Stippling/cross hatching for south polar
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <4FBA0CDA.3050303 at jamstec.go.jp>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Dear Katherine,
>
> Just use maskout function, like:
>
> ga-> hatch maskout(var,-lat-20) min max
>
>
> Good luck!
>
> Chihiro Kodama
>
>
> (2012/05/17 13:27), Katherine Sadler wrote:
>> Thanks Chihiro, it now works.
>>
>> However it also draws cross-hatching outside the boundary i've chosen (i
>> want between the south pole and the equator, but it seems to draw
>> cross-hatching from the pole to about 20 degrees north), is there a way
>> to limit it?
>>
>> Thank you also for sharing hatch.gs <http://hatch.gs>, it has been very
>> useful to me.
>>
>> Cheers,
>> Katherine
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> --
> =====================================================
> KODAMA Chihiro, Ph.D.
> Japan Agency for Marine-Earth Science and Technology
> Research Institute for Global Change
> PHONE: +81-45-778-5623 (direct)
> E-Mail: kodamac at jamstec.go.jp
> URL: http://kodama.fubuki.info/
> =====================================================
>
>
>
> ------------------------------
>
> Message: 6
> Date: Mon, 21 May 2012 12:52:18 -0700
> From: rjmera <rjmera at coas.oregonstate.edu>
> Subject: [gradsusr] Area average only for California
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <4FBA9CF2.1050009 at coas.oregonstate.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Hi, all. I've looked at some of the previous posts on this subject but I
> wasn't sure if it covered what I wanted to do. What I'm trying to do is
> to get area average for an irregular polygon. To be specific, the state
> of California, for example. This would exclude other states in the
> vicinity and the ocean. How should I proceed?
>
> Robert
>
> I'm using this version of grads:
>
> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1
> Copyright (c) 1988-2010 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> GrADS comes with ABSOLUTELY NO WARRANTY
> See file COPYRIGHT for more information
>
> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf
> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile
> Issue 'q config' command for more detailed configuration information
> Loading User Defined Extensions table
> </home/server/pi/homes/rjmera/grads/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/gex/udxt>
> ... ok.
>
>
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 27, Issue 37
> ****************************************
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wrf_plotex.gs
Type: application/octet-stream
Size: 1239 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120522/ed6fcadf/attachment-0003.obj
More information about the gradsusr
mailing list