立即注册 登录
气象家园 返回首页

比利的个人空间 http://bbs.06climate.com/?80614 [收藏] [复制] [分享] [RSS]

日志

NC数据转grd

已有 98 次阅读2019-1-26 21:54

下面给出使用Grads读取nc格式文件的一个范例脚本文件(.gs文件):

 

   以位势高度场为例(红色字体为解释说明语句;蓝色字体为输入的命令,也就是gs文件的内容)

 

‘reinit’

‘sdfopen d:xxxxhgt.mon.mean.nc ‘  (打开需要的nc文件,注意路径)

   

‘set gxout fwrite’

‘set fwrite d:xxxxxxhgt.bin’     (指定要输出的二进制文件名和路径)

 

‘set lon 0 357.5’

‘set lat -90 90 ‘                  (设定写出区域范围,这是全球范围,第一句是必需的,如果是用默认,则会在x方向输出为145个点,要注意)

 

 

‘set lev 1000’                       (设定要输出的层次)

‘set t 1 640’                        (设定时间范围)

 

‘d hgt’                           (输出)

 

‘reinit’

 

这样就生成了需要的二进制文件,其对应的ctl文件需要另外给出。

 

以上是输出单层的二进制文件,把多层的数据放在同一个二进制文件中,就要使用循环,根据grads规定的数据存放次序(依次为经度、纬度、高度、变量、时次),如下:

 

‘reinit’

‘sdfopen d:xxxxhgt.mon.mean.nc’    (打开需要的nc文件,注意路径)

   

‘set gxout fwrite’

‘set fwrite d:xxxxxxhgt.bin’         (指定要输出的二进制文件名和路径)

 

‘set lon 0 357.5’

‘set lat -90 90’                     (设定写出区域范围,这是全球范围,第一句是必需的,如果用默认set lon 0 360,则会在x方向输出为145个点,要注意)

 

t=a                             (a为起始时次)

while(t<=b)                    (开始循环,b为结束时次)

‘set t’ t                           (设定时次)

 

‘set lev 1000’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 

‘set lev 925’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 

‘set lev 850’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 

‘set lev 700’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 

‘set lev 600’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 

‘set lev 500’                       (设定要输出的层次)

‘d hgt ‘                          (输出)

 。

。                        (可以添加任意的层次)

。                  (由于高度分布不连续而且数量有限,所以用这种方法)

 

 

t=t 1

endwhile                           (结束循环)

‘reinit’

 



全部作者的其他最新日志

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

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

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

返回顶部