如何用GRADS画二维风场矢量图不同年份的合成图

lwy326

谢谢，学习了

skyinglll

楼主好像这样做不对，我今天才发现，求平均求错了，哎，比如求Ave（uwnd,t=1,t=5)，求的不是第一和第五个时刻的平均，而是第一时刻到第五时刻的平均！！！！！！！！！！心好痛！！！！！！

skyinglll

感觉心痛到无法呼吸，我以为已经做完了，结果都错了。。。。。。

940918yjy

楼主做出来了吗

风儿飘飘

时间上没有规律的话，只能逐个相加，再求平均

幽居寒舍

这是我的，求上年10月到当年5月的平均，不知道对不对，仅供参考
'reinit'
'sdfopen e:\grads\NCEP\uwnd.mon.mean.nc'
'sdfopen e:\grads\NCEP\vwnd.mon.mean.nc'
'set lon 55 150'
'set lat 5 60'
'set grads off'
'set grid off'
'set gxout vector'
'set mpdset country1 yunn guiz guangx hun'
'set xlpos -4'
'set ylpos -4'
'set parea 0 11 1.2 7.8'
'define uwnd=(ave(uwnd,t=10,t=826,12)+ave(uwnd,t=11,t=827,12))+ave

(uwnd,t=12,t=828,12)+ave(uwnd,t=13,t=829,12)+ave(uwnd,t=14,t=830,12)+ave

(uwnd,t=15,t=831,12)+ave(uwnd,t=16,t=832,12)+ave(uwnd,t=17,t=833,12))/8'
'define vwnd=(ave(vwnd,t=10,t=826,12)+ave(vwnd,t=11,t=827,12))+ave

(vwnd,t=12,t=828,12)+ave(vwnd,t=13,t=829,12)+ave(vwnd,t=14,t=830,12)+ave

(vwnd,t=15,t=831,12)+ave(vwnd,t=16,t=832,12)+ave(vwnd,t=17,t=833,12))/8'
'd uwnd.1;vwnd.2'
'run axis.gs -type L -lsize 0.15 -lfont 5 -position o -interval 10 -start 10 -end 50'
'run axis.gs -type T -lsize 0.15 -lfont 5 -position o -interval 20 -start 60 -end 140'
'gxprint e:/grads/picture/asianwind.png png white'
;