关于多时刻，多变量的站点资料处理(肯定有用)

司空摘月

关于多时刻，多变量的站点资料处理

        本人也是在看了气象家园上面的资料学习后，总结出来的。有获得就得有回报啊，O(∩_∩)O~家园上关于多时刻多变量的介绍不是很多，而我又是要处理这样的数据，所以写出来供大家参考，有问题可以提出来一起交流啊。废话不多说，开始啦：

一：首先跟大家说下基本的流程：

1.      产生一个含有所有数据名称的文档（我的是obs.txt）

2.      将观测资料编译成二进制文件（cwd.f90→obs.grd）

3.      编写读取obs.grd的ctl文件（obs.ctl）

4.      生成格点的背景场，为插值做准备（grid.f90→grid.grd）

5.      编写读取grid.grd的ctl文件（grid.ctl）

6.      编写gs文件画图

二：具体操作：

1.我的数据格式见图1（我处理的是北京的数据，总共是2天的数据，每一天数据的每个变量有217个数，前面2行是一个说明我们在读取时直接跳过）

图1

1.      产生obs.txt

图1是一个操作说明，图2是我的操作。

图2

图3

2        将观测资料编译成二进制文件（cwd.f90→obs.grd）

图4是我的程序，有人是用动态数组处理的，我会在后面给个网址大家看看。（关于程序里面具体的一些说明，见附录的文件，这些文件讲的很详细，需要大家最好看看，受益匪浅）

图4

3        编写读取obs.grd的ctl文件（obs.ctl）

图5是对应的ctl文件，仔细看的同学会发现有个

stnmap F:\2011072627\obs.map

这是为二进制数据生成站点映射文件，这个必须有。所以会有图6的操作，打开grads在里面输入图6的操作，注意路径中必须是用 /  。

都操作完之后，就可以话站点数据图了（自己可以试试），但是没有插值到格点上，所以接下来是匹配相应的格点。

图5

图6

4        生成格点的背景场，为插值做准备（grid.f90→grid.grd）

图7是生成格点背景场的程序

图7

5        编写读取grid.grd的ctl文件（grid.ctl）

图8为对应的文件

图8

6        编写gs文件画图

图9是写的gs文件，注意需要打开2个ctl

'open F:\2011072627\grid.ctl'

'open F:\2011072627\obs.ctl'

图10为所画的一个图

      

图9

        

图10

好啦，到这里已经基本上是把站点资料图画出来啦，但是每个fortra和ctl、gs的意义我没有具体说明，我在后面附上一些前辈们写的文档可以供大家参考，一定要看啊。

三 思考：

• 关于最后一个gs文件里面有个oacres，这是一个插值函数，至于如何很好运用这个函数我还不是很了解，高手指点
• 虽然插值是一个不错的方法，但是误差也存在不少，所以我觉得可以不插值。把站点的数据大小直接显示出来，如图11所示，可以通过不同颜色来表示其大小或者圆圈的大小表示值的大小。（可惜我现在不会这个，所以还需要会的同学给予指导啊）
  

图11

动态数组：http://bbs.06climate.com/forum.php?mod=viewthread&tid=13422

                                  很有用的网址：http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903&extra=page%3D2

http://bbs.06climate.com/forum.php?mod=viewthread&tid=3079

图片我这边复制不上去，所以我把word文档放在附件了

司空摘月

图片没有复制的上去，所以大家要看的话下载一下那个word 文档吧

lqouc

司空摘月 发表于 2013-5-15 11:15

                               
登录/注册后可看大图

图片没有复制的上去，所以大家要看的话下载一下那个word 文档吧

图片要先上传再贴

lqouc

司空摘月 发表于 2013-5-15 11:15

                               
登录/注册后可看大图

图片没有复制的上去，所以大家要看的话下载一下那个word 文档吧

你说多变量多时刻的介绍不多，就你贴出来的清风的教程上面就有吧。清风的教程上把所有的程序也都给了，你卖的贡献不低，好歹说说你的内容做了什么改进，或者有什么补充。

river

这还挺贵的呀，那就回个帖算了。

mofangbao

先谢谢楼主的辛苦劳动啦！
不过，我已经去除我那两个帖子的附件啦，你可以直接把帖子地址贴上来的~

喵~

谢楼主分享。正需要这些方法。被插值折磨的死去活来的了。

司空摘月

lqouc 发表于 2013-5-15 11:24

                               
登录/注册后可看大图

你说多变量多时刻的介绍不多，就你贴出来的清风的教程上面就有吧。清风的教程上把所有的程序也都给了，你 ...

他那个里面讲的比较系统，我写的只是针对多时刻多变量的这一种，需要系统学习的话还是参照清风的比较好啊。我也是第一次发帖，贡献值就随便写了一个，如果贵了还请见谅，需要的话 我可以直接发给你。

司空摘月

mofangbao 发表于 2013-5-15 13:24

                               
登录/注册后可看大图

先谢谢楼主的辛苦劳动啦！
不过，我已经去除我那两个帖子的附件啦，你可以直接把帖子地址贴上来的~

好的，主要还是学习的你的材料，第一次发帖，有些细节没注意，下次会注意的O(∩_∩)O~

mofangbao

建议楼主调低贡献值啊，要不然可能不太好卖哦~薄利多销才是王道啊