请教一下大家grib资料怎么画时间序列图

水台

问题如上，就是横坐标是时间，纵坐标是位势高度，该怎样处理呢？

edwardli

不管是用NetCDF还是grib(1 or 2)还是其他的什么，画时间序列图最关键的是合并文件，写描述文件，生成控制文件。
我用GFS还有T639还有NC都画了，方法如下：
第一步：合并文件，cmd下  copy /b *.grib1 XXXX.grib1
第二步：描述文件，cmd下 grib1或者nc资料  grib2ctl XXXX.grib1 > YYYY.ctl
                                           grib2资料（GFS）g2ctl XXXX.grib2 > YYYY.ctl
              **描述文件中比较重要的是：时间纬度与间隔，比如T639的240时效的资料，60小时以内的是3hr一个，60~120的是6hr，根据你要整出来的时间轴大小、用途确定、更改**
第三步：控制文件：cmd下 gribmap -v -i YYYY.ctl >ZZZZ.txt    ——这个.txt有两个作用，保存信息+确保不刷屏
第四步：gs文件： 前面的就不说了。。。'open YYYY.ctl'.....
                                                                'set lat 40'
                                                                'set lon 120'
                                                                'set z 1 9'  ——你要画的层次  比如这个是T639的是1000hPa~100hPa
                                                                 ‘set t 1 21’ ——T639的72小时，6hr一个
                                                                不是正式格式  比如要 湿度 垂直速度 风。。。包括图层属性：颜色、色阶、色标等等具体的就不写了 可以参考GrADs的东东。保存成Slice.gs
第五步：画图：cmd下 grads -cl Slice.gs    或者grads中打开

我的万用gs文件 Slice.gs (964 Bytes, 下载次数: 309)

2012-3-15 20:49 上传

点击文件名下载附件
以下是我画的。。。

NCEP—Re1的2012年3月1日至5日分析

GFS的2012031412起报的时序图

T639的

maplesnow82

edwardli 发表于 2012-3-15 20:50

                               
登录/注册后可看大图

不管是用NetCDF还是grib(1 or 2)还是其他的什么，画时间序列图最关键的是合并文件，写描述文件，生成控制文 ...

edwardi，答复真的是很详细啊，佩服！

edwardli

呵呵  大前天晚上一宿在研究这个，昨天一宿研究c#调用。头发都掉光了~~~

edwardli

@兰溪之水  客气了 还是要谢谢大神你的，看了你的FNL处理 我才鼓捣通这个的~~~

mofangbao

edwardli 发表于 2012-3-15 20:50

                               
登录/注册后可看大图

不管是用NetCDF还是grib(1 or 2)还是其他的什么，画时间序列图最关键的是合并文件，写描述文件，生成控制文 ...

真是很不错，一定要给分

edwardli

mofangbao 发表于 2012-3-16 08:56

                               
登录/注册后可看大图

真是很不错，一定要给分

谢谢了  最近在搞CMAcast，数据一坨，多数数据都有。回头整理好了发到论坛里大家多指教。

水台

感谢edwardli的指教，很详细，回头学习学习，有不懂的还要再请教你

浅水

太强大了 我最近研究这个 都没找到突破口啊 感激啊

浅水

合并了t639得grib文件 生成了ctl文件 但是文件里的时间描述怎么改啊？请教不改就画不出时间演变图