<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear GrADS users,</DIV>
<DIV>I have taken the many tips from previous postings and external templates to produce a script that displays a time series of rainfall anomalies taken from the monthly GPCP dataset. The script works great as long as the geographic coordinates are relayed as a point location (ex: lat=3N, lon=6N), but, I am more interested in determining rainfall anomalies for average area regions (ex: lat=6N,lat=11N,lon=5W,lon=2E). </DIV>
<DIV> </DIV>
<DIV>I have tried to incorporate the GrADS command AAVE to calculate the area average, but I am not getting the desired results. Please see the below script for details:</DIV>
<DIV> </DIV><SPAN lang=EN>
<DIV>'reinit'</DIV>
<DIV>*open gpcp monthly data file</DIV>
<DIV>'sdfopen c:\testdata\netcdf\precip.mon.mean.nc'</DIV>
<DIV>'set grid off'</DIV>
<DIV>'set grads off'</DIV>
<DIV>*create 12 x-axis labels; one for each month</DIV>
<DIV>'set t 1 12'</DIV>
<DIV> </DIV>
<DIV>*define new variable to create monthly climatology for rainfall starting from Jan 1997<BR>'define rainclim = ave(precip, t+216, t=347, 12)'<BR>'modify rainclim seasonal'<BR>'set t 1 last'</DIV>
<DIV>*define new variable for anomaly of monthly rainfall<BR>'define anom = precip - rainclim'</DIV>
<DIV>*set time period for anomaly time series</DIV>
<DIV>'set t 336 347'</DIV>
<DIV>*define geographic point location to assess rainfall anomaly for the given time period</DIV>
<DIV>'set lat 8'</DIV>
<DIV>'set lon 1'</DIV>
<DIV>'d anom'</DIV>
<DIV> </DIV>
<DIV>I would like to replace the set lat and set lon commands with:</DIV>
<DIV>'define region = aave(anom,lat=6,lat=11,lon=355,lon=2)</DIV>
<DIV>which should average rainfall anomaly for that area. I then change 'd anom' to 'd region', but nothing happens. Can anybody spot the error? Also, is there a problem when the x dimension varies from 0 to 360 and you are interested in the region from 357 to 1 (5W to 2E), obviously the order from 357 to 1 does not obey the proper value order? Has anybody else overcome issues with the AAVE command?</DIV>
<DIV> </DIV>
<DIV>Thank you for your time.</DIV>
<DIV>sincerely,</DIV>
<DIV>teddy allen</DIV>
<DIV></SPAN><BR><BR> </DIV>
<DIV>
<DIV><STRONG>web: <A href="http://personal.ecu.edu/tla0923/HOME.html" target=_blank rel=nofollow>http://personal.ecu.edu/tla0923/HOME.html</A></STRONG></DIV>
<DIV><STRONG>"To move quickly, go alone. To move far, go</STRONG> <STRONG><EM><FONT color=#ff0000>together</FONT></EM>."</STRONG></DIV>
<DIV><STRONG>African Proverb </STRONG><STRONG> *</STRONG></DIV></DIV>
<DIV><STRONG></STRONG> </DIV></td></tr></table><br>