[gradsusr] problem with an extract script

Ismaila Diallo ismailadiallo64 at yahoo.fr
Tue Jan 4 10:02:26 EST 2011


Dear Grads User's

i've daily data for a big region. Now i try to extract since this file the data 
available for an area and this extract data must be only the data average of 4 
months (june july august and september. so that i perform this script: 


'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'

but with this script i get some errors and i don't understand what's wrong in my 
script.
so could someone help me to resolve this difficult.
Believe in your great Help.
Sincerely
Ismaila 



      



More information about the gradsusr mailing list