请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6566|回复: 7

Grads如何画季节平均风场

[复制链接]

新浪微博达人勋

发表于 2018-7-20 10:38:43 | 显示全部楼层 |阅读模式
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)'



这是数据属性

这是数据属性

这是出图结果

这是出图结果
asianwind.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-20 10:54:18 | 显示全部楼层

回帖奖励 +1 金钱

确实没用,报错应该告诉你没找到vwnd变量了吧?
display的时候加了.2那前面ave也要记得加上啊。另外define的变量最好换个别的变量名。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-20 11:07:12 | 显示全部楼层
楼上是正解。define的变量换个变量名应该就没有问题了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-20 12:03:10 | 显示全部楼层
lqouc 发表于 2018-7-20 10:54
确实没用,报错应该告诉你没找到vwnd变量了吧?
display的时候加了.2那前面ave也要记得加上啊。另外define ...

非常感谢,我改成了define a和b;ave里改成uwnd.1和vwnd.2
这样就好了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-20 12:03:16 | 显示全部楼层
lqouc 发表于 2018-7-20 10:54
确实没用,报错应该告诉你没找到vwnd变量了吧?
display的时候加了.2那前面ave也要记得加上啊。另外define ...

非常感谢,我改成了define a和b;ave里改成uwnd.1和vwnd.2
这样就好了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-20 12:03:19 | 显示全部楼层
lqouc 发表于 2018-7-20 10:54
确实没用,报错应该告诉你没找到vwnd变量了吧?
display的时候加了.2那前面ave也要记得加上啊。另外define ...

非常感谢,我改成了define a和b;ave里改成uwnd.1和vwnd.2
这样就好了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-20 12:03:23 | 显示全部楼层
lqouc 发表于 2018-7-20 10:54
确实没用,报错应该告诉你没找到vwnd变量了吧?
display的时候加了.2那前面ave也要记得加上啊。另外define ...

非常感谢,我改成了define a和b;ave里改成uwnd.1和vwnd.2
这样就好了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-8-13 16:17:22 | 显示全部楼层
有用,感谢楼主分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表