Grads求标准差，结果看不懂

余梦于心

'reinit'
'open E:\gcm\bnu\historicalts.ctl'
'set fwrite E:\gcm\bnu\03.grd'
'set gxout fwrite'
'set lat 0 40'
'set lon 70 120'
'set lev 1000'
'set gxout shaded'
'set t 1563'
i=1563
while(i<1864)
'define a=ave(ts,t='i',t='i+2')'
i=i+12
endwhile
'set lat 5 30'
'set lon 89 110'
'define r=re(a,50,linear,70,1,40,linear,0,1,bs)'
'define anom=r-ave(r,t=1,t=26)'
'define s=sqrt(ave(anom*anom,t=1,t=26))'
'd s'
'cbarn 0.7 0'
'disable fwrite'
;

oucdaqi

至于问题三，我想是因为'set gxout shaded'覆盖了'set gxout fwrite'的命令以至于出不来新的grd

余梦于心

oucdaqi 发表于 2017-4-3 09:06
至于问题三，我想是因为'set gxout shaded'覆盖了'set gxout fwrite'的命令以至于出不来新的grd

确实是这样，O(∩_∩)O谢谢

余梦于心

第一问和第二问有没有大神帮看看啊