请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 53262|回复: 89

[分享资料] 多层次站点资料的Grads绘图

  [复制链接]

新浪微博达人勋

发表于 2013-1-19 20:09:26 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 whmooncool 于 2013-1-21 10:09 编辑

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。在使用是发现如果不进行赋值,后续工作也无法进行,会出现错误。

Grads站点资料的多层次绘图.doc

43.5 KB, 下载次数: 368, 下载积分: 金钱 -5

评分

参与人数 5金钱 +49 贡献 +14 收起 理由
meehooqq + 2 很给力!
Aires + 14 + 6 赞一个!
qxtlyf + 10 + 2
river + 8 + 1 很给力!
mofangbao + 15 + 5

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-19 22:27:25 | 显示全部楼层
学习了。期待好帖子!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-20 13:24:30 | 显示全部楼层
期待好久了都没有人整理出来,今天终于有了,非常感谢。要是再加上一个多时次的就更完美了,支持!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-1-20 15:41:03 | 显示全部楼层
这个确实还没有人完整的整理过,我的那个帖子里面也只是顺带提了一下,可是还是有很多人不会,非常感谢楼主的分享,最近放假好多人都回家了,没及时加分啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-21 09:16:53 | 显示全部楼层

恩恩,当时懒了就没有整理了,等放假回来,我整理整理一并发上来~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-21 09:23:35 | 显示全部楼层
mofangbao 发表于 2013-1-20 15:41
这个确实还没有人完整的整理过,我的那个帖子里面也只是顺带提了一下,可是还是有很多人不会,非常感谢楼主 ...

从家园下载了好多资料,也得贡献点~~不过,等待清风的肯定好久了,呵呵~~家园加油!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-23 20:27:19 | 显示全部楼层
清风 这个很好 客户四斑竹 第一 里面的trim函数可以批量读取文件的话 那就不需要用cmd那个命令来批量读取了吗? 第二 既然高度场是 经度 纬度 层次 站点 为什么不定义为4维的呢 是不是站点确定了纬度经度就都确定了呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-23 20:42:11 | 显示全部楼层
很好的资料,风场如何转化?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-3-24 15:39:43 | 显示全部楼层
wolfwang 发表于 2013-3-23 20:42
很好的资料,风场如何转化?

有说明的,根据当时的风向和风速,分别计算U和V
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-24 16:49:00 | 显示全部楼层
这个帖子很实用,很值得学习!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表