- 积分
- 244
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-6-29
- 最后登录
- 1970-1-1
|
GrADS
系统平台: |
|
问题截图: |
- |
问题概况: |
如何使用Grads,利用NCEP资料U、V分量画季节平均风场 |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
3 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 幽居寒舍 于 2018-7-20 12:08 编辑
求大神帮忙看看我的代码哪儿出了问题。我想画上年10月-当年5月平均风场,并标识>5m/s的区域。
资料时间跨度为1-838(1948:1:1:0-2017:10:1:0)。
近地面资料。
图是能出来,但是明显感觉我求平均的公式没有起作用,这个求平均的公式对吗?
下面是我的代码:
'reinit'
'sdfopen e:\grads\NCEP\uwnd.mon.mean.nc'
'sdfopen e:\grads\NCEP\vwnd.mon.mean.nc'
'set grads off'
'set grid off'
'set lon 55 150'
'set lat 5 60'
'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'
'set mpdset country1 yunn guiz guangx hun'
'set xlpos -4'
'set ylpos -4'
'set parea 0 11 1.2 7.8'
'set cthick 5'
'set cmin 5'
'set gxout shaded'
'd mag(uwnd.1,vwnd.2)'
'set gxout vector'
'd uwnd.1;vwnd.2;mag(uwnd.1,vwnd.2)'
'cbarn'
'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'
;
遵照大神的指导,以下是代码里被修改的部分,如下:
'define a=(ave(uwnd.1,t=10,t=826,12)+ave(uwnd.1,t=11,t=827,12)+ave(uwnd.1,t=12,t=828,12)+ave(uwnd.1,t=13,t=829,12)
+ave(uwnd.1,t=14,t=830,12)+ave(uwnd.1,t=15,t=831,12)+ave(uwnd.1,t=16,t=832,12)+ave(uwnd.1,t=17,t=833,12))/8'
'define b=(ave(vwnd.2,t=10,t=826,12)+ave(vwnd.2,t=11,t=827,12)+ave(vwnd.2,t=12,t=828,12)+ave(vwnd.2,t=13,t=829,12)
+ave(vwnd.2,t=14,t=830,12)+ave(vwnd.2,t=15,t=831,12)+ave(vwnd.2,t=16,t=832,12)+ave(vwnd.2,t=17,t=833,12))/8'
'd mag(a,b)'
'd a;b;mag(a,b)'
|
|