grads处理多个ctl文件和nc文件
2011-10-10 21:03:59| 分类：grads学习 | 标签： |举报 |字号大中小 订阅
下载LOFTER我的照片书 |
用grads处理多个相同格式的数据时若单个单个处理非常麻烦，当文件非常多的时候是单个处理是不实际的。下面介绍一种方法;
第一步，在这种情况下可以重新写一个ctl描述文件，其文件变量都和已知的ctl相同，若原来的n文件只是时间不同，那么新描述文件的时间维数是所有原文件的时间的和。同样，若其他维数不同时也用同样的方法处理。
第二步，在第一行之后添加一行： options template 表示多个时间序列原始数据文件想用一个描述文件统一地描述。这些原数据的原文件名由dset定义的形势命名文件名。
第三步，修改dset 的文件名。原路径不变，把文件名用％表示。其中：
％y2 代表两位数年
％y4 代表四位数年
％m1 代表一位或者两位数的月
％m2 代表两位数月（用0补齐1位数）
％mc 3个字符月份的缩写
％d1 1或2位天
％d2 两位天
％h1 1或者2位时
％h2 2位时
例如：
原文件其中之一的文件名为gdas2006050812f00，且所有文件只有天和时的变化
那么新描述文件的文件名为：gdas200605%d2%h2f00
另外如果源文件里有index项的话，需要修改其idx的文件名，假设改成fnl.idx。并用在dos下用gribmap函数生成一个新的idx文件。gribmap -e -i fnl.ctl（加绝对路径）
open fnl.ctl就可以打开所有文件。
********************************************************************************************************************************** 若想要提取从1951-2006年56年nc文件中的某些数据，一个一个处理非常麻烦，这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。
收稿日期:2002-11-1作者简介:高文良,男,31岁,主要从事短期气候预测等研究工作。气象绘图软件GrADS高级使用技巧
高文良1 刘晓燕2 曾小东3
(11成都高原气象研究所 成都 610072; 21阿坝州金川县气象局 624100
31阿坝州马尔康县气象局 624000)
摘 要:本文通过介绍气象绘图软件GrADS使用中
的一些高级技巧,分析了GrADSScript语言中的难
理解之处,并介绍了在GrADS中做合成分析和t检
验的程序,可以对学习和使用GrADS的科研人员起
参考、帮助作用,达到事半功倍的效果,促进GrADS
软件的应用更广泛、深入,提高科研人员的工作效率。
关键词:GrADS;技巧;程序
中图分类号:TP391文献标识码:C
文章编号:1003-7187(2002)04-0057-03
1 引言
气象绘图软件GrADS(GridAnalysisandDis2
playSystem)是免费共享软件,可随时从互联网上下
载(/grads/),后续版本正陆续
推出[1]。GrADS有丰富的内部函数,可以对数据进
行计算和分析处理。它支持处理格点资料和站点资
料,并且支持对GriB码文件、特殊格式文件(如一字
节整型、二字节整型、大中型机器二进制数据等)的直
接读取,气象科研领域应用非常广泛[2]。在其最新
1.8SL9版本中,GrADS又将应用领域推进到了海洋学科,功能也得到了进一步地增强和扩展。
但在使用GrADS的过程中,特别是使用GrADS
Script语言编程当中,有一些问题需要特别注意。因
为GrADSScript语言是一种类似于VBScript(或
MATLABScript)的高级语言,稍不注意就容易出现
错误,且查错十分困难,这点与其他语言有较大的差
别。
2 站点数据处理
GrADS中站点数据处理基于格点数据的基础之
上,需先将站点数据通过Cressman客观分析方法内
插至格点上,然后再依照格点资料的处理方式对站点
1 GrADS 快速参考手册
(GrADS Version 1.7 beta 7 ) 编：***********************译： 朱 禾
命令行参数grads -options
-c‘命令组’ 启动GrADS后执行该命令组
-b 以批处理方式运行
-l 以landscape 模式运行
-p 以portrait 模式运行
-C 以两位数记年时，自动设定世纪值
一般性设置
help 列出帮助选项
set grads on|off 打开/关闭 显示“GrADS” 标志
set display > 设置显示模式为：
grey|greyscale 用灰度方式
color 以彩虹模式绘图
set frame 设置图形边框绘图方式为：
on 以一个矩形包围所有图形
off 不画边框
circle 对于其它投影方式画矩形边框，而极射投影画圆形边框
set background index 设置绘图区背景颜色
display expression
d expression 显示expression代表数据的图形
open control_filename 打开一个数据描述文件，可指明路径
close file# 关闭数据描述文件,file#=1,2,3…
set dfile number 改变当前ctl文件为number指定描述文件,number=1,2,…
define var=expression
var=expression 定义名为var的新变量,
undefined var 取消var定义的变量
modify varname 定义变量为气候序列平均，时间类型如下：
seasonal 月平均
diurnal 一日内多时刻平均
query
q 查询命令及其参数：
files 列出所有已打开的数据描述文件
file n 列出第n个数据描述文件的信息，n=1,2,3……
1 GrADS下任意方向垂直剖面的实现
缪明 余清平
（解放军理工大学气象学院）
廖德敏
（93966部队气象室）
引言
GrADS（Grid Analysis and Display System）是当前气象学中应用最为广泛的图形图像处理软件之一，它功能强大，使用方便而深受用户的欢迎。利用它可直接实现图像的动画显示，在三维数值模拟研究中，在固定模式高度的情况下，可较为方便地显示某一模式预报量在该层次上的水平分布；在固定经（纬）度的情况下，还可方便地显示某一模式预报量沿该经（纬）度的经（纬）向垂直剖面分布。但在实际问题的三维数值模拟研究中，用户关心的不仅仅是模式预报量沿某一经纬向的垂直剖面分布，更多的时候用户关心的是沿某一方向的垂直剖面分布，但在GrADS下如何实现这一功能，这可能是困扰一般的GrADS用户的问题。其实，GrADS除了强大的图形图象的分析处理能力外，还有较强的数据处理能力，我们可利用GrADS的交互性和数据处理的特点，利用GrADS脚本语言的编程技术，按GrADS的数据格式，通过提取所要分析的模式预报量沿某一方向上的值（ascII码形式），然后再将提取ascII码值转换成GrADS支持的二进 2 制格式，并为转化后的数据编写相应的数据描述文件（ctl文件），在GrADS下打开它就可实现某一模式预报量沿该方向的垂直剖面分布的分析。
一、 实现的思想和步骤
在介绍实现的思想和gs编程（GrADS Script programing）之前先说明GrADS下某些基本的常识和某些基本命令的含义和功能。
1． GrADS基本常识和基本命令的简介
（1）进入GrADS后，在“ga ->”后每执行一个GrADS命令，在文本窗口都会返回1行或2行的的文本代码，返回的代码即为GrADS的默认变量“result”值，每执行一次GrADS命令，原来的“result”值将被刷新成刚执行的GrADS命令返回码。在gs编程中，可通过代码“say result ”看到默认变量“result”对应的返回文本代码。
