[gradsusr] Fw: problem with an extract script

[Ismaila Diallo]

Hi,
muhamed,you are right for your suggestion (how must i do the loop in this 
script)???????So the data time is more than 10 years and i need only a period of 
10 years it's why i set the tmax at 3920. 1) in each year i need only the june 
july august september data (JJAS).2) the data covert a big domain and i need 
data for a region (in the big domain).3) with this average i must get one value 
for each year (this value must be the 


average of JJAS for the year).

Until now i try to change this script but something is wrong again. 


the script is:

'reinit'
'sdfopen /home/lpaosf/RT.1996-2007_DAY.AVG.NC'
t=518
tmax=3926
'set fwrite test.dat'
'set gxout fwrite'
while (t<=3926)
   'set t 't
   'define jun=ave(rt,t=t+0,t='t+30')'
   'set t 't+30
   'define jul=ave(rt,t=t+0,t='t+31')'
   'set t 't+61
   'define aug=ave(rt,t=t+0,t='t+31')'
   'set t 't+92
   'define sep=ave(rt,t=t+0,t='t=t+30')'
   'define jjas=(jun+jul+aug+sep)/4'
   'define aa=ave(jjas,lon='-8.5',lon='8.5')'
   'define bb=ave(aa,lat='10',lat='16')'
   'd bb'
  t=t+365
endwhile
'disable fwrite'

so again i think in your great help to resolve this problem.

Sincerely
ismaila