[gradsusr] plotting time in hovmoller pictures

[Jeffrey Duda]

Liping,
You could try templating the files using CHSUB:

http://grads.iges.org/grads/gadoc/descriptorfile.html#CHSUB

Jeff Duda

On Tue, May 18, 2010 at 5:32 AM, Liping Liu <liping_liuwang at hotmail.com>wrote:

>  Thanks for the answers to my previous question. The "parea" and "vpage"
> commands work, but I have to specify the picture size every time. I want a
> more robotic (generic) approach, if possible.
> The GrADS template data descriptor is a wonderful idea! But does it only do
> aggregating or concatenating? I want to manipulate the fields before I plot
> out the hovmollers:
> First I plot the hovmoller using data from file 1 for Aug10-15:
> 'set yflip on'
> 'set lev 1000'
> 'set tlsupp year'
> 'set t 2 48'
> 'set lat 15'
> 'set gxout shaded'
> 'set grads off'
> 'd ave(olr.1,lat=5,lat=25)'
> Then I want to add the hovmoller plot using data from file 1 and file 2 for
> Aug15-20:
> 'set t 49 80'
> 'd ave((olr.1+olr.2)/2,lat=5,lat=25)'
>
> Question: Is it possible that I write the fields (olr) from different data
> files (with different ctl files) into one data file with one ctl file?
>
> Thanks!
>
> Liping
>
>
> > Date: Wed, 12 May 2010 18:10:26 -0400
> > From: Charles.Seman at noaa.gov
> > To: gradsusr at gradsusr.org
> > Subject: Re: [gradsusr] plotting time in hovmoller pictures
> >
> > Liping Liu,
> >
> > You might also be able to use GrADS templates to "concatenate" the "olr"
> > fields from file 1 and file 2 into a single time series within GrADS as
> > follows:
> >
> > 'open olr_template.ctl'
> > 'set yflip on'
> > 'set lev 1000'
> > 'set tlsupp year'
> > 'set t 2 80'
> > 'set lat 15'
> > 'set gxout shaded'
> > 'set grads off'
> > 'd ave(olr,lat=5,lat=25)'
> >
> > where "olr" is the name of the variable in the "vars" section of the
> > GrADS template data descriptor file "olr_template.ctl" (please see the
> > examples at http://grads.iges.org/grads/gadoc/templates.html which can
> > give you some ideas for what your "olr_template.ctl" might need)...
> > briefly, a template data descriptor file for time differs from a single
> > dataset data descriptor file in the "dset" record, a required "options
> > template" record, and in the "tdef" record (in which the number of time
> > levels in the template data descriptor file is the sum of the time
> > levels for each of the component datasets).
> >
> > I hope this helps,
> > Chuck
> >
> > Liping Liu wrote:
> > > Hi, (I am a new user of GrADs)
> > >
> > > I am trying to plot two hovmollers into one figure. The problem is the
> > > time.
> > > First I plot the hovmollers using data from file 1 for Aug10 to Aug15,
> > > with the following commands:
> > >
> > > 'set yflip on'
> > > 'set lev 1000'
> > > 'set tlsupp year'
> > > 'set t 2 48'
> > > 'set lat 15'
> > > 'set gxout shaded'
> > > 'set grads off'
> > > 'd ave(olr.1,lat=5,lat=25)'
> > >
> > > So it gives me a figure for Aug10--Aug15. But, when I repeat the
> > > plotting for another period of time Aug15-Aug20:
> > >
> > > 'set t 49 80'
> > > 'd ave(olr.2,lat=5,lat=25)'
> > >
> > > It will simply overlap the whole figure, the time is not sequenced.
> > >
> > > Question: what do I do to get the pictures of Aug10-15 and Aug15-20
> > > laid out vertically (not overlapping each other)?
> > >
> > > An alternative way is to write the two sets of data into one data
> > > file. But I don't know how to do that either.
> > >
> > > Thanks for your help. Suggestions to either of the approaches are
> > > greatly appreciated!
> > >
> > > Liping
> > >
> > >
> > >
> ------------------------------------------------------------------------
> > > Hotmail has tools for the New Busy. Search, chat and e-mail from your
> > > inbox. Learn more.
> > > <
> http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_1>
>
> > >
> > >
> ------------------------------------------------------------------------
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> >
> > --
> >
> > Please note that Charles.Seman at noaa.gov should be considered my NOAA
> > email address, not cjs at gfdl.noaa.gov.
> >
> > ********************************************************************
> > Charles Seman Charles.Seman at noaa.gov
> > U.S. Department of Commerce / NOAA / OAR
> > Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
> > 201 Forrestal Road fax: (609) 987-5063
> > Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/<http://www.gfdl.noaa.gov/%7Ecjs/>
> > ********************************************************************
> >
> > "The contents of this message are mine personally and do not reflect any
> > official or unofficial position of the United States Federal Government,
> > the United States Department of Commerce, or NOAA."
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
> ------------------------------
> Hotmail is redefining busy with tools for the New Busy. Get more from your
> inbox. See how.<http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_2>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100518/71221e18/attachment-0003.html