[已解决]Grads中不同时次的barb图可以叠加吗

lingxiyizhi

如题，利用生成的 fnl_20130101_00_00.ctl 和Grads画0点的风场图，gs命令如下：

'reinit'
'open D:/test/201301/fnl_20130101_00_00.ctl'
'set grid off'
'set grads off'
'set gxout contour'
*设置地图
'set mproj scaled'
'set mpdset cnworld'
*设置绘图区域
'set vpage 0 11 0 8.5'
'set parea 1 10.5 1 7'
'set lat 15 53'
'set lon 73 136'
'set lev 500'
*设置等值线
'set cterp on'
'set csmooth on'
'set clevs 520 524 528 532 536 540 544 548 552 556 560 564 568 572 576 580 584 588 592'
'set cint 6'
'set cthick 6'
'set ccolor rainbow'
'd HGTprs/10'

'set gxout barb'
'set cthick 1'
'set digsize 0.05'
'set rbcols'
'd skip(UGRDprs*2.5,3);skip(VGRDprs*2.5,3);mag(UGRDprs,VGRDprs)'

'draw title 2013/01/01 00:00(BJT) 500hPa HGT'
'printim D:\test\201301\fnl_20130101_00_00BJT_hgt_wind.png white x1000 y800'
;
效果图：


本人Grads菜鸟，有几个疑问，想请教各位大虾：
1、敢问各位大虾上面的gs还缺少什么必要的命令吗？“d HGTprs/10”这个命令表示什么意思啊？
2、参考兰溪兄的资料，本人已生成2013年1月份124个时次（每天4个时次）的批处理ctl文件和idx文件，想知道怎么使用生成的这个批处理ctl文件画其他时次的图啊？？因为当用它替换上面gs中的ctl文件时，生成的还是1月1号0点的图。。。
3、不同时次的barb图可以叠加吗？如果不能叠加的话，假如想知道1月份的主导风向，这个问题可以用Grads来分析吗？

==》》》》》》根据各位大虾以及river大神的指导，画出了1月份124个时次的平均风场，gs命令（位势高度已去掉），如下：

'reinit'
'open D:/test/201301/fnl_201301.ctl'
'set grid off'
'set grads off'
'set gxout contour'
*设置地图
'set mproj scaled'
'set mpdset cnworld'
*设置绘图区域
'set vpage 0 11 0 8.5'
'set parea 1 10.5 1 7'
'set lat 15 53'
'set lon 73 136'
'set lev 500'

'set gxout barb'
'set cthick 1'
'set digsize 0.05'
'set rbcols'

'define a1=ave(UGRDprs,t=1,t=124)'
'define a2=ave(VGRDprs,t=1,t=124)'
'd skip(a1*2.5,3);skip(a2*2.5,3);mag(a1,a2)'
'define a3=ave(mag(a1,a2),t=1,t=124)'
'd a3'

'draw title 2013/01 500hPa WIND'
'printim D:\test\201301\fnl_201301wind.png white x1000 y800'
;

利用Grads2.1，出图如下：

river

建议去看grads实用手册，多看下基础

我是彤-Tina

d HGTprs/10  就是画高度场呀，/10的原因可能是原来的资料数值是4位，现在你除了10，等值线上的数值不就变成3位数了吗。

hm_others

生成idx时，ctl里对应的什么时间，idx就是什么时间
建议批量描述后生成idx或一个一个生成

lingxiyizhi

hm_others 发表于 2015-11-30 10:40
生成idx时，ctl里对应的什么时间，idx就是什么时间
建议批量描述后生成idx或一个一个生成

我是按照你说的，先生成的批量ctl文件，然后对批量ctl文件生成的对应的idx文件

lingxiyizhi

我是彤-Tina 发表于 2015-11-30 10:09
d HGTprs/10  就是画高度场呀，/10的原因可能是原来的资料数值是4位，现在你除了10，等值线上的数值不就变 ...

受教了，多谢多谢！

lingxiyizhi

路过的大虾，请留步...

Jying

批处理之后，算平均场

lingxiyizhi

Jying 发表于 2015-11-30 14:42
批处理之后，算平均场

多谢提醒！非常感谢！根据提示，我现在正在研究平均风场的画法...您手头有类似平均场的gs吗？

lingxiyizhi

river 发表于 2015-11-30 00:01
建议去看grads实用手册，多看下基础

我今天研究了一下手册，能不能麻烦你帮我看一下我画得1月份的平均风场的图和gs命令对不对啊多谢了