grads处理NCEP再分析数据

shatong

说说这几天用grads处理1°×1°的NCEP再分析数据（grib2格式）的心得和常出现的一些问题，这里和大家分享一下，希望会有帮助。
可以看看兰溪的PDF，里面有详细的步骤，但素作为初学者的我，还是研究了好久呢
1)先用g2ctl，生成fnl数据初始时刻的那个数据文件相对应的ctl文件;2)然后将ctl里面的一些东西改一下，使之成为可以批量处理fnl数据文件的ctl:3)用gribmap,生成可以批量处理多个fnl数据文件的idx文件；4）编写gs：从多个fnl数据文件(grib2格式)中读取想要的变量（vars）和所研究的区域(set lat,set lon)，生成一个grads可以读取的数据格式（比如：dat,grb格式）
中间遇到的问题：
1.ctl文件生成是0kb
答：g2ctl.exe版本有些老，兰溪的pdf教程上就阔以下载新哒
2.idx索引文件木有生成
答：gribmap.exe有点老，建议用2.1版本的grads，这个气象家园上有，我就素下的气象家园整合版的grads
3.编写gs文件时，从ctl文件中读取不出来一些变量，但有的却可以。比如我想读取vvelprs（垂直速度），这个物理量就木有，并且显示的是
display error:invalid expression    expression = vvelprs
答：是因为g2ctl生成的ctl文件中对于该量的描述太长了，你阔以把**后的内容删了，就阔以读出该变量了。
这几天处理fnl资料遇到的问题，虽然很幼稚，但希望对大家有一点点帮助吧。
上传了自己的gs文件（读取自己想要的变量，和研究的区域），供大家参考哈

river

升级一下你的g2ctl.exe就可以解决变量读不出来的问题了

visus

{:5_213:}{:5_213:}{:5_213:}{:5_213:}{:5_213:}

旺旺仙贝

{:eb328:}{:eb328:}{:eb328:}{:eb328:}{:eb328:}

不许捏我脸

非常有用，感谢~~

xiarixuepian

非常有用啊，感谢分享