[gradsusr] Help invoking user function
Davide Sacchetti
davide.sacchetti at arpal.gov.it
Tue May 22 11:42:10 EDT 2012
I understand
now your script works
Pay attention on variables naming: 'result' is a reserved name: you
should not use it as a function input argument name
I made some change to the script just for testing on my system (minimal
changes) and I changed "time format" to an easier one
It should help ...
bye bye
Davide
On Tue, 2012-05-22 at 10:31 -0400, Zilore Mumba wrote:
> 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
> > ****************************************
>
> _______________________________________________
> 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
-------------- next part --------------
*
#'c'
#'reinit'
#'open wrfout_d01_2012-05-19_00.ctl'
'set display color white'
'c'
'set grads off'
'set mpdset hires'
'set mpt 0 4 1 6'
'set mpt 1 1 1 6'
'set frame on'
'set gxout shaded'
'set grid on 3 1'
************************
'set t 1'
BaseDate=ValidD()
**************************************
'set t 6'
ValidDat=ValidD()
*'run grads_date.gs'
* say ValidDat
#'d rainc'
'd mslp-995'
'set clopts 0 1 0.09'
*'set clskip 3'
'set clevs 0 2 4 8 10 12 14 16 18 20 22 24 26 28 30'
'set ccolor 0'
'set cmin 1'
'set gxout contour'
#'d rainc'
'd mslp-995'
'Draw title WRF Rainfall T+06H\Valid to 06Z on 'ValidDat' From 00Z on 'BaseDate''
'run cbarc.gs'
#'enable print raint+06'
#'print'
#'disable print'
#'!gxps -ic raint+06 -o raint+06.ps'
#'!ps2pdf13 raint+06.ps raint+06.pdf'
#'!rm raint+06'
#'!rm raint+06.ps'
*
function ValidD()
'q time'
word=subwrd(result,3)
hour=substr(word,1,2)
day=substr(word,4,2)
month=substr(word,6,3)
year=substr(word,9,4)
date=day'-'month'-'year
return date
More information about the gradsusr
mailing list