[gradsusr] problem with an extract script

Ismaila Diallo ismailadiallo64 at yahoo.fr
Tue Jan 4 15:48:29 EST 2011


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


      



More information about the gradsusr mailing list