分享一个Grads画风场图的脚本和详细说明

幽居寒舍

给大家分享一个Grads画风场图的脚本，对于初学者来说比较管用，脚本特点在于大家可以根据自己的需要更改里面的参数，这里我做脚本里的一些详细使用说明，希望对初学者有一点帮助（我也是初学者），如有不当的地方，望大家指正:
1、求季节或任一时间跨度平均风场：
'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'
这里计算了从上年10月到当年5月的8个月平均风速，资料的时间跨度为1:838，如果你只是需要某一个时次的风速，你可以删除这一代码，如改为'set t 1'。如果你不知道1对应时间是多少，你可以打开grads，输入以下代码查看：
'sdfopen e:\grads\NCEP\uwnd.mon.mean.nc'（路径自己设置，见下图）
'q ctlinfo'
2、地图设置
'set mpdset country1 yunn guiz guangx hun'
这里使用的底图是气象家园里的，不是grads自带的，这个底图处理了国界问题，这是非常重要的改进。
country1是底图，大家可以搜一搜就找到了，放到opengrads文件夹里就可以调用了，如果还需要分省的地图，可以在后面继续添加，如yunn（云南）等等。下载后放到grads的安装文件SupportData里，E:\Grads\opengrads\Contents\Resources\SupportData
3、去掉自带风矢比例尺
系统自带的这个很难看，位置也不好处理，可以先去掉。
'set arrlab off'
4、设置画图区范围
我们经常会发现出图缺失了一部分信息，所以需要调整画图区范围，大家可以参考grads使用说明或脚本卡去理解，这里就不一一说明，设置里面的参数看看会发生什么变化有助于理解。
'set parea 0 11 1.2 7.8'
5、设置矢量箭头大小
'set cthick 4'
6、用阴影标识大于多少风速
'set cmin 4'
'set gxout shaded'
7、对风速进行颜色分类（设置彩虹色）
'set clevs 1 2 3 4 5 6 7'
8、自定义风矢大小（这一步可以省略，因为我们已经用颜色对风速进行了区分）
'drawvc ua va 1 7.63 0.3 5'
9、出图
可以根据自己的需要设置输出图片格式，如jpg、png、pdf等等
'gxprint e:/grads/picture/windp10c5.png png white'

另外，别说你没有资料，NCEP资料下载网址：
https://wenku.baidu.com/view/e8622c268e9951e79b892764.html
https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html
废话不多说了，附上里面所需要的以下脚本。

biubiubiu

{:5_214:}

biubiubiu

{:5_213:}

feifei20141010

什么情况？

cz33

怎么都被屏蔽了？

幽居寒舍

为啥把我屏蔽了

jamesWang

这个帖子好，适合初学者，顶

desertsky

谢谢，正在学习Grads...

lm8005507771

学习了！感谢楼主分享！

blnuist

回帖不增加积分，建议坚持签到，多分享！