Dear Grads users,<br>I am a very beginner in GrADS. <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>--------------------------------------------------------------<br>'sdfopen <a href="http://sst_mon_81_07.nc">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>
<br><br><br>