[gradsusr] Help invoking user function
Zilore Mumba
zmumba at gmail.com
Tue May 22 10:31:21 EDT 2012
Davide,
Sorry to be a bit of a bother, my problem is that my function does not
seem to take the result from 'set t 6' at all'. I appreciate the
solution of removing :'s, this is what I was going to figure out next.
I would appreciate if you could help me with a code snippet for the
loop, am a bit lost.
Zilore Mumba
On Tue, May 22, 2012 at 6:15 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: Help invoking user function (Davide Sacchetti)
> 2. stream function in pressure level (resmi ea)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 22 May 2012 11:41:21 +0200
> From: Davide Sacchetti <davide.sacchetti at arpal.gov.it>
> Subject: Re: [gradsusr] Help invoking user function
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <1337679681.1012.116.camel at dvd>
> Content-Type: text/plain; charset="UTF-8"
>
> you cannot rely on the ':' position of 'set t 6' result ...
>
> I mean: the answer could be:
>
> Time values set: 2012:5:1:9 2012:5:1:9
>
> or
>
> Time values set: 2012:11:10:19 2012:11:10:19 ...
>
> A solution is to transform ':' in blanks, then using subwrd:
>
> loop over result chars (subwrd(result,i,1))
> if subwrd(result,i,1) = ':' set blank, else let it unchanged
> after the loop year = subwrd 1, month ...
>
> hoping it helps a bit ...
> Davide
>
> On Tue, 2012-05-22 at 01:06 -0400, Zilore Mumba wrote:
>> 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
>> > ****************************************
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> --> Attenzione cambio indirizzo: davide.sacchetti at arpal.gov.it
> --
> Davide Sacchetti
> Centro Funzionale Meteo Idrologico di Protezione Civile della Regione Liguria
> ARPAL Unit? Tecnica Complessa di livello Regionale
> V.le Brigare Partigiane 2 16121 Genova (I)
> tel: +39 010 6437535 fax: +39 010 6437520
> mail: davide.sacchetti at arpal.gov.it web: www.meteoliguria.it
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Tue, 22 May 2012 15:45:43 +0530
> From: resmi ea <reshmiea at gmail.com>
> Subject: [gradsusr] stream function in pressure level
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CAHda7jbWyDg1Kzzb0DUPCc6TomAuT9yTK7JGO0N7v8beHCax-A at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear users
> i have calculated stream function using open grads from 1000-10 hPa level
> using NCEP/NCAR reanalysis data ex.1dec2009
> The grads script give the grd file containing three variable like psi,upsi
> and vpsi.
> Here i have attached the figure. i find mistake in the figure...the
> programme doesnot undergone integration part, at each level it calculate
> three variable and finally dispalys.
> i want to know that whether it is correct or not. how can i modify the
> programme
> 'reinit'
> 'sdfopen G:\windows\ncep\uwnd\uwnd.2009.nc'
> 'sdfopen G:\windows\ncep\uwnd\vwnd.2009.nc'
> 'set lon 0 357.5'
> 'set lat -90 90'
> prompt 'key-in file name for psi:'
> pull fname
> 'set fwrite 'fname
> tt=1
> while (tt<=1)
> 'set t 'tt
> zz=1
> while (zz<=17)
> 'set z 'zz
> 'define psi= fish_psi(uwnd,vwnd.2)'
> 'define one=1+0*lat'
> 'define upsi=mvadv(one,psi)'
> 'define vpsi=-muadv(one,psi)'
> 'set gxout fwrite'
> 'd psi/1e7'
> 'd upsi'
> 'd vpsi'
> zz=zz+1
> endwhile
> tt=tt+1
> endwhile
> 'disable fwrite'
>
>
>
> --
> Resmi E A
> Research Scholar
> Cochin university of science and technology
> Cochin-16
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120522/83e6b86f/attachment.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: 1dec2008psi.jpeg
> Type: image/jpeg
> Size: 15795 bytes
> Desc: not available
> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120522/83e6b86f/attachment.jpeg
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 27, Issue 39
> ****************************************
More information about the gradsusr
mailing list