Ncep资料作距平

小亭子ng

打开两个文件用Grad对ncep资料进行作差，怎么做啊？我想用2013年6-8月92天逐日资料减去35年的夏季高温平均。我定义了一个a作为平均值，想用2013年的资料来减去它，怎么把它用在下面文件打开的地方呢？新手编的比较乱，请多多指教

'reinit'
'set gxout fwrite'
'set fwrite H:\data1\NCEP-1\juping.grd'
'open H:\data1\NCEP-1\1979-20136-8.ctl'
'set x 1 192'
'set y 1 94'
tt=1
'define a=ave(tmax,t='tt+0',t='tt+34')'

'sdfopen H:\data1\NCEP-1\daily\Tmax\tmax.2m.gauss.2013.nc'
tt=152
t0=1
'set x 1 192'
'set y 1 94’
while(tt<=243)
'd tmax-a'
'set t 't0''
t0=t0+1
endwhile
'disable fwrite'
;