[gradsusr] plotting time in hovmoller pictures

[Liping Liu]

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/
> ********************************************************************
> 
> "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.
http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_2
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100518/791850dd/attachment-0003.html