Dear Madam,<br>Thanks for your reply.<br>In fact I sent my question yesterday to the GrADS user group, but didn't get any reply so far.<br>I still could not resolve it.<br><br>I found that there is no keyword search option in the Grads Forum archive. Could you please tell me if there is any such?<br>
<br>The data file link is given below:<br><pre><a href="http://apdrc.soest.hawaii.edu/dchart/index.html?cid=AAAAJA@@" target="_blank">http://apdrc.soest.hawaii.edu/dchart/index.html?cid=AAAAJA@@</a><br><br></pre>Please also see my yesterday's email. if you have any further suggestion, I will highly appreciate.<br>
--------------------------------------------------------------<br><br>I have a multi-year SST data. I want to get a time series monthly and annual mean plot of this data<br>in a specified region in the same figure, keeping same x-axis (time) and y-axis (SST).<br>
<br>Please see below the code I am having, I also attach the figure I am getting out of it.<br>The annual mean value is found to be constant for all the 25 years I considered.<br>Definitely some mistake I am making in the code, which I am not able to identify.<br>
<br>Can somebody please help me the correct way of doing it?<br><br>------------------------------<div id=":1ar" class="ii gt">--------------------------------<br>'sdfopen <a href="http://sst_mon_81_07.nc/" target="_blank">sst_mon_81_07.nc</a>'<br>
<br>'set x 1 360'<br>
'set y 1 180'<br>'set z 1'<br>'set t 1 311'<br><br>'set lon 30'<br>'set lat 5'<br>'set t 2 301'<br>'define sstarea=aave(sst.1,lon=30,lon=78,lat=5,lat=31)'<br>'set gxout line';'set cstyle 1';'set cmark 3';'set ccolor 1'<br>
'd tloop(sstarea)'<br><br><br>'set x 1 360'<br>'set y 1 180'<br>'set z 1'<br>'set t 1 311'<br><br>ts = 1<br>while ( ts <= 25 )<br>'set t 'ts * 12<br>'define sstave = ave(sst,t-10,t+1)'<br>
ts = ts + 1<br>endwhile<br><br>'set lon 30'<br>'set lat 5'<br>'set t 1 25'<br>'set gxout line';'set cstyle 2';'set cmark 5';'set ccolor 3'<br>'d tloop(aave(sstave,lon=30,lon=78,lat=5,lat=31))'<br>
<br><br>
</div><br><br><div class="gmail_quote">On Mon, May 24, 2010 at 6:08 PM, Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div style="">Dear Sudev, <div><br></div><div><div><div class="im"><div>On May 21, 2010, at 11:25 AM, sudev das M P wrote:</div><br><blockquote type="cite">Dear Grads users,<br> I am a beginner in GrADS. <br> I have a multi-year time series data. I want to get the annual mean in this data.<br>
<br> Please see that from the following I can get monthly time series plot, and it is working fine.<br> <br> 'set lon 30'<br> 'set lat 5'<br> <br> 'set t 2 302'<br></blockquote></div>Why is this not 'set t 1 302' ? </div>
<div><div class="im"><br><blockquote type="cite"> <br> *------------------------<br> *Monthly time series<br> *------------------------<br> 'define sstmonth=aave(sst.1,lon=30,<div>lon=78,lat=5,lat=31)'<br></div></blockquote>
<div><br></div><blockquote type="cite"><div> 'set gxout line'<br> 'set cmark 4'<br> 'set grads off'<br> 'set ccolor 9'<br> 'd sstmonth'<br> <br> In the same plot, I want to have one more curve showing the annual mean.<br>
I mean 26 annual mean values for the 26 years.<br> <br> I don't know how to do this.<br></div></blockquote></div><div>Try this after you have defined sstmon:</div><div><br></div><div>'set t 1 12' </div><div>
'sstann = ave(sstmonth,t+0,t=302,12)'</div><div>'modify sstann seasonal'</div><div><br></div>--Jennifer</div><div><br><blockquote type="cite"><div class="im"><div> <br> I will highly appreciate if someone can help me in this.<br>
<br> Thanks,<br><font color="#888888"> <br> Sudev</font></div></div><div class="im"> _______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br></div></blockquote></div><font color="#888888"><br><div> <span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="">
<span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a></div><div><br>
</div><br></span></span></div></span> </div><br></font></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br>