多层次站点资料的Grads绘图

whmooncool

Grads站点资料的多层次绘图

看了很多关于如何使用grads和fortran转换站点资料的帖子，大多数都集中在地面或者单层的，对于高空的多个层次很少有提到，而且大多数有些模糊，不好理解，我把自己做的一些整理传上来，供大家参考参考。欢迎大家指出理解和程序的错误。

一，首先是数据的说明

选用micaps高空观测资料850hpa,700hpa,600 hpa,500hpa，200hpa，100 hpa六个层次作为例子。为了方便，时间层次就选择1个时间。

micaps高空观测资料的说明如下：

diamond 2 10年07月28日20时200百帕高空观测

10 07 28 20  200  419

因此定义mt=419,nz=6

数据每行的排列分别是：站号，经度，纬度，高度……..具体详见micaps的使用手册，这里不再赘述。根据这些定义变量。

二，fortran读取站点资料

注：

1.根据micaps的高空资料定义相应的变量，由于850-100这6个层次的站点是相同的，站号，经纬度等只需定义成一维的，如lat(mt)，而lev(nz)表示有6个层次，也为一维；高度场等又有mt个站的资料又有nz个层次，因此定义为2维的。

2.由于micaps资料分别在不同的文件内，文件名重复，我修改为100_2820,200_2820,500_2820…….然后使用trim函数顺序读取资料。

   3.为了更好的说明地面和高空资料的读取，将地形高度high(mt)设置为地面的量，nlev定义为7，即高空六层，地面一层

   4.以输出地面一个变量，高空三个变量为例。先进行站点的循环，再进行高空层次的循环。

   5.前面我们定义了lev(nz)，现在给lev(nz)赋值，使用select case.,分别赋值为850,700,600,500,200,100。在使用是发现如果不进行赋值，后续工作也无法进行，会出现错误。

kongfeng0824

学习了。期待好帖子！

river

期待好久了都没有人整理出来，今天终于有了，非常感谢。要是再加上一个多时次的就更完美了，支持！

mofangbao

这个确实还没有人完整的整理过，我的那个帖子里面也只是顺带提了一下，可是还是有很多人不会，非常感谢楼主的分享，最近放假好多人都回家了，没及时加分啊

whmooncool

river 发表于 2013-1-20 13:24

                               
登录/注册后可看大图

期待好久了都没有人整理出来，今天终于有了，非常感谢。要是再加上一个多时次的就更完美了，支持！

恩恩，当时懒了就没有整理了，等放假回来，我整理整理一并发上来~~

whmooncool

mofangbao 发表于 2013-1-20 15:41

                               
登录/注册后可看大图

这个确实还没有人完整的整理过，我的那个帖子里面也只是顺带提了一下，可是还是有很多人不会，非常感谢楼主 ...

从家园下载了好多资料，也得贡献点~~不过，等待清风的肯定好久了，呵呵~~家园加油！！

jeffggg

清风 这个很好 客户四斑竹 第一 里面的trim函数可以批量读取文件的话 那就不需要用cmd那个命令来批量读取了吗？ 第二 既然高度场是 经度 纬度 层次 站点 为什么不定义为4维的呢 是不是站点确定了纬度经度就都确定了呢？

wolfwang

很好的资料，风场如何转化？

whmooncool

wolfwang 发表于 2013-3-23 20:42

                               
登录/注册后可看大图

很好的资料，风场如何转化？

有说明的，根据当时的风向和风速，分别计算U和V

托尼托尼乔巴

这个帖子很实用，很值得学习！