求助，我想在850hPa叠加风场和降水，在滤掉风场的时候，为什么出来的风场图是t=1的呢

askira2011

'reinit'
'sdfopen D:/postgraduate/ERAI/0501-0731uv.nc'
'open D:\postgraduate\ERAI\orog.ctl'
'sdfopen E:/cdo/daily.0520-0709.nc'

'set dfile 3'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set map 1 1 2'
'set lon 60 150'
'set lat 10 60'
'set z 1'
'set gxout shaded'
'set cmin 5'
'run 200colors.gs'
'd ave(r,t=1,t=9)'
'run cbarn.gs'

'set dfile 1'
'set dfile 2'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set map 1 1 2'
'set lon 60 150'
'set lat 10 60'
'set lev 850'
'topo=orog.2(t=1)'
'set gxout vector'
'set arrscl 0.5 10'
'define u1=maskout(u.1,1500-topo)'
'define v1=maskout(v.1,1500-topo)'
'uave=ave(u1,t=77,t=112)'
'vave=ave(v1,t=77,t=112)'

'd  uave;vave'

'q w2xy 94.33 29.67'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.08'

'draw title 05.20-05.28 850hPa wind rainfall'

'printim D:/postgraduate/ERAI/2u850.png white'

'reinit'
;

askira2011

求哪位大神指教一下，谢谢啦！

风之牧语

同求答案啊！

alameda

楼主有答案了吗

askira2011

alameda 发表于 2014-10-16 09:07
楼主有答案了吗

把set dfile 2删掉就可以了，原来这么简单。