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