grads 使用fnl grib1 数据 t做时间循环 利用while 绘制33个时次的图形

天堂ミ亚澜

二爷名声在外 发表于 2014-7-1 20:00
再发这个就算你灌水……

这。。厉害~~

天堂ミ亚澜

ziki 发表于 2014-7-1 21:35
依照fnl的数据命名规则，看你打开的ctl这个fnl是2014年6月16日00时的一个时次的数据？
这样可以循环t么？ ...

可以循环 我打开的是第一个文件，用索引文件加ctl.

天堂ミ亚澜

1092995269 发表于 2014-7-2 09:15
楼上正解，，文件只包含一个时次的各种物理量资料

fnl 资料 dos 处理过了。

ziki

天堂ミ亚澜 发表于 2014-7-4 20:03
可以循环 我打开的是第一个文件，用索引文件加ctl.

嗯…那很有可能就是有人说到的那个set t那儿的问题…
其实可以把ctl和报错的截图贴上来的，这样大家可能更明白一点，因为可能都没遇到过，光看gs里的语句可能也不知道哪儿有问题

ziki

首先，我是女的……我的头像难道给人感觉我是男的么……

Q1：“要画出33个时次的水汽通量叠加风场 500pa环流形式。 我想用while 循环一次性画出来并且自动更改图片名称及图像的题标。但是总是提示while 循环有问题。”

这个问题我用你发给我的文件测试过了，（话说下次麻烦把gs一起打包给我好嘛……还得到这个贴子里来找你之前的gs语句……不过你要是本来打算让我按你的要求给你重写一个当我没说……）也是报错endwhile无效，改了@sun_shine_Xia 之前说的set t那儿空一格的问题也不行，感觉应该是enable print这个不能在循环里面？具体为什么我也不清楚，反正是把这个挪出来就可以正常用了。其实没必要开那么多个gmf啊，在一个gmf里面看起来也不麻烦吧，反正标题都有的。改过的gs如下，你把路径和循环次数改了就可以了。嗯另外你说是要画水汽通量和500hPa环流形势，不麻烦的，自己研究下吧，不懂的多使用搜索功能~
'reinit'
'open d:\tddownload\fnl.ctl'
'set mpdset cnhimap'
'enable print d:\tddownload\fnl_UV.gmf'

t=1
while (t<=3)
'set grads off'
'set gxout barb'
'set lat 34 50'
'set lon 72 95'
'set lev 500'
'set t 't''
'd ugrdprs;vgrdprs'
'draw title UV_'t''
'print'
'c'
t=t+1
endwhile

'disable print'
'reinit'

Q2：“想画风场、涡度、垂直速度的剖面图，由于不了解fnl变量的具体意思，总是画出来不对。”

关于变量的意思，请参考兰溪的这个帖子，讲的很清楚~
http://bbs.06climate.com/forum.php?mod=viewthread&tid=6226

天堂ミ亚澜

ziki 发表于 2014-7-8 11:26
首先，我是女的……我的头像难道给人感觉我是男的么……

Q1：“要画出33个时次的水汽通量叠加 ...

这，有点小尴尬。感谢学姐之指点。因为我之前不做预报，最近调到气象台实习才开始做这方面的东西。最近收获很大，昨晚上已经解决了剖面图的问题。希望能跟学姐有更多交流。。

ziki

天堂ミ亚澜 发表于 2014-7-8 18:57
这，有点小尴尬。感谢学姐之指点。因为我之前不做预报，最近调到气象台实习才开始做这方面的东西。最近收 ...

在气象台的话，很多业务方面需要分析的东西micaps可以很好的实现，虽然数据网格没那么细~

天堂ミ亚澜

ziki 发表于 2014-7-8 20:59
在气象台的话，很多业务方面需要分析的东西micaps可以很好的实现，虽然数据网格没那么细~

但是精度不高，论文要求达不到。得出文章呀~~