中国科学院大气物理研究所
大气科学和地球流体力学数值模拟国家重点实验室
National Key Laboratory of Numerical Modelling for Atmospheric
Sciences and Geophysical Fluid Dynamics ( LASG )
Institute of Atmospheric Physics, Chinese Academy of Sciences
Technical Report
GrADS使用手册
LASG 编
中国科学院大气物理研究所
大气科学和地球流体力学数值模拟国家重点实验室
( LASG )
http://www.lasg.ac.cn
关于本手册
本手册是在 1995 年张洪编写的《GrADS 气象图形系统使用方法》及其它相关文档的基础上
汇编而成。主要介绍了 GrADS 的基本用法、使用技巧、及与 GrADS 相关的绘图技术。本
手册将根据软件发展不定期进行修订。
版权
本手册版权属中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验
室(LASG)。本手册可免费下载及使用。
本手册的维护与更新
本手册由王军负责维护与更新,如果有什么意见和建议,请与他联系,或在论坛上讨论。
Email : wjun@lasg.iap.ac.cn
有用的链接
GrADS Homepage :
GrADS Script Library:
http://grads.iges.org/grads/
http://grads.iges.org/grads/gadoc/library.html
http://www.lasg.ac.cn/atmos/gadoc/library.html (本地)
http://web.lasg.ac.cn/grads/index.html
http://www.lasg.ac.cn/cgi-bin/forum/forums.cgi?forum=3
本手册下载:
GrADS 中文讨论区:
致谢
感谢张洪编译了《GrADS 气象图形系统使用方法》。
感谢 funny 对本手册耐心细致的修订工作,以及在论坛上给予大家的热情帮助。
GrADS 使用手册
§1.1 GrADS 简介
第一章 基本概念与基本操作
GrADS(Grid Analysis and Display System)提供一个全 32 位的交互操作的气象格点数
据和站点数据的分析和显示环境。该系统具有气象数据分析功能强、地图投影坐标丰富、高
级编程语言使用容易、图形显示快速,并具有彩色动画功能等特点,因此迅速成为国外气象
数据显示的标准平台之一。
GrADS 当前版本为 1.8s11,可用于如下平台:DEC、Intel/LINUX、SUN、Macintosh OSX、
SGI、IBM/AIX、MS Windows。
GrADS 软件包所包含的文件:
gradsc
gradsdods
gradsnc
gradshdf
gribmap
gribscan
gxps
gxeps
gxtran
stnmap
wgrib
GrADS Classic - readline, printim only
GrADS with readline, printim, LATS, read/write NetCDF, DODS
GrADS with readline, printim, LATS, read/write NetCDF
GrADS with readline, printim, LATS, read/write HDF, read NetCDF
"Maps" a GRIB data set to a GrADS descriptor file
Extracts grid and variable info from a GRIB data set
Converts GrADS metafiles to Postscript (older version)
Converts GrADS metafiles to Postscript (Level 1 or Level 2)
Displays metafiles
Maps station data
See http://wesley.wwb.noaa.gov/wgrib.html
GrADS 英文主页:http://grads.iges.org/grads/
§1.2 基本概念
按磁盘文件记录格式,GrADS 软件包的文件分为如下几类:
(1) .dat 二进制无格式记录的原始气象数据文件,其数据既可以是格点数据,也可以是站
点数据。它们是从其它气象数据(如站点气象报、格点气象报、模式格点输出结果)转
换生成的。对格点数据而言,其格式为二进制无格式直接或顺序记录格式。总之对格点
数据,整个数据集是一个大的五维数据场,包括三维物理空间、一维物理变量、一维时
间变量,存放时以二维数组片的形式按水平、垂直、物理变量、时间序列的顺序排放,
详情参见§4。
(2) .ctl 原始数据描述文件。该文件为纯 ASCII 文件,用以描述原始数据集的基本信息,包
括数据集文件名、数据类型、数据结构、变量描述等等,详见§4。在 GrADS 环境中至
少得首先打开(open)一个数据描述文件,以便后续的操作有数据对象。
(3) .gs GrADS 控制文件,用 GrADS 命令 run 执行之。这也是一个纯文本文件,内含用描
述语言(script language)写成的批处理 GrADS 系统设置和命令,可集成处理 GrADS 命
1
http://www.lasg.ac.cn
令。详见§6。
(4) .exe 系统命令文件,是 GrADS 系统在 DOS 环境下的各执行文件。如 grads.exe 为 GrADS
图形分析和显示命令;dos4gw.exe 为 DOS 的扩展环境;gxps.exe、gxpsc.exe、gxpscw.exe
都是图元文件转换为 postscript 文件的执行文件;gxtran.exe 是图元文件转换到显示器上
显示的执行文件;gx.exe 是将图元文件转换为各种不带 ps 解释器的打印机输出的执行文
件。详见§1.3。
(5) .gmf(.met) GrADS 系统图元输出文件,格式由 GrADS 内定,文件名随用户自定,
其内容为屏幕显示图形的二进制图元数据,用于产生图形的硬拷贝输出。在 windows
平台,用 gv.exe 或 gv32.exe 可以查看此图元文件,并可将其另存为.wmf 格式
的图形文件。详见§1.3,§5。
(6) .ps Postscript 格式文件,其内容为 ASCII 码形式的 Postscript 语言格式的图形数据,
它是图元文件*.gmf(*.met)经 gx.exe、gxpsc、gxpscw 转换生成的,可用于 ps 打印机
的直接硬拷贝输出,也可被其它应用软件调用,只要该软件识别 ps 格式数据。详见附
录二“PS 文件详解”。
(7) .exc 直接执行批处理文件,其内容为 GrADS 交互环境下所打命令的直接集成,按记录
存放在一个 ASCII 码文件中,在 GrADS 环境下用 exec 命令执行之。
启动 GrADS 后首先需打开至少一个数据描述文件,按打开文件的次序系统自动给打开
文件编号,第一个打开的文件为 1 号文件,以后顺排。以后的维数环境设置和图形操作都是
针对当时的缺省文件进行的,除非变量名下标标出文件编号,自动的缺省文件为 1 号文件,
改变当前缺省文件序号用 set dfile # 命令来设定第#号文件为当前缺省文件。
维数环境是 GrADS 的一个重要概念,GrADS 视每一个物理变量场为一个四维数据集
(4D data set)。包括空间三维和时间一维,也可固定其中的一维或几维以获得实际的低于四
维的数据子集。GrADS 中设置维数用以说明或指定随后的分析或图形操作时参加操作的原
始数据集的维数范围,即通过设定工作数据的起止点数、取点频度(间隔)来设定最后工作
数据场的数组成份。该工作数据集可以是整个原始数据场,也可以是原始数据场的一部分,
详见§2,§6。所谓维数环境是对格点数据而言的。维数环境的定义可在两种空间坐标上进
行。一种是地球坐标(world coordinate),以经纬度为度量单位;另一种是格点坐标(grid
coordinate),以网格点数为度量单位。二者一般对应于同一个绘图坐标,都对应到网格点上。
在维数环境表达式中 x,y,z,t 与 lon,lat,lev,time 是分别对应于两套坐标的专用维数变
量,含义固定,如 x 与 lon 都指西到东指向的(缺省方向)水平坐标,y 与 lat 都指南到北指
向的(缺省方向)水平坐标,z 与 lev 都指从地面到高空的(缺省方向)垂直坐标,t 与 time
都是时序坐标,不过 t 用的是格点时次序号,而 time 用的是格林威治标准时的 GrADS 绝对
表达格式。详见§2,§4。
当所有维数都固定时,得到的是一个单值数据点;如果只有一维变化,得到的时一维数
据线,屏幕显示时为一条曲线;二维发生变化时对应于二维切片(slice),屏幕显示时缺省
表达为二维平面图,也可显示为一维曲线的动画序列;三维发生变化时 GrADS 解释为一个
二维切片的序列,屏幕显示时须设定一维作为动画维,以动画方式显示;四维变化就须指定
两维为固定或用动画方式才能显示。总之图形输出只能以二维或一维方式表达多维变量。
GrADS 定义的绘图区域为横放或竖放两种矩形区域(缺省为横放),它们缺省的硬拷贝
输出区域大小分别为 11x8.5 和 8.5x11 英寸。所谓横放即所谓风景画形式(Landscape),简
2
GrADS 使用手册
称 L,竖放即所谓肖像画形式(Portrait),简称 P,两者都是通常的 A4 纸大小,所以我们在
计算机屏幕上使用的工作窗口最好也按上述长宽比例设置。GrADS 启动时首先提示用横放
还是竖放形式,缺省为横放,如想竖放则键入 no 后回车,随后窗口工作区就从原先的字符
窗口状态进入图形窗口状态。需要说明的是现在键盘输入显示和执行后的回应信息显示在字
符窗口,而图形显示在图形窗口,两者未分离,而是重叠显示。不过通过 print 命令产生的
图形硬拷贝文件中不包含所有字符窗口的信息。GrADS 的命令提示符为“ga>”。可以利用
frame.gs 或 stack.gs 来半自动分离字符和图形窗口。
GrADS 的绘图工作区分三个层次,一层是实际页(real page),即硬拷贝的 A4 纸大小,
单位为英寸(注意横放或竖放);一层是虚拟页(virtual page),单位也是英寸,缺省时虚页
等同于实页;第三层是在虚页中指定绘图区域,其单位用的虚页中的虚英寸,即缺省时等同
于实际英寸,当设置虚页后按比例度量。注意第三层所定义的区域只包含图形,不包括坐标
轴、标题等附属信息的位置,即要预留出附属信息的区域。
§1.3 基本操作
新版本的 GrADS 交互环境中可以使用上下左右箭头键对曾经使用的命令进行调用和编
辑,但仅限于本次启动 GrADS 交互环境后所使用的命令。
在 DOS 或 Unix shell 环境下的命令:
grads
启动 GrADS 系统,在交互式方式或批处理方式下分析或显示气象数据,并可生成硬拷贝图
元文件。其语法为:
grads [ -lpbc "run 描述文件名" ]
其中 l 表示横放,p 表示竖放,输入该选项后启动 GrADS 时系统将不再提示选择方向,b 表
示批处理,屏幕不再显示图形结果,直接完成全部操作,c 表示进入 GrADS 环境后首先执
行随后跟在 run 命令后的描述文件,该文件由用户设定,类型为*.gs。几个选项可组合使用,
也可键入 grads 直接交互操作。
gxps,gxpsc,gxpscw
将 GrADS 内 print 命令产生的图元文件分别转换为单色(灰度)、彩色以及反白(即白底)
彩色的 postscript 格式文件,用于硬拷贝输出。语法为:
gxps -i 图元文件名 -o ps 文件名
gxpsc -i 图元文件名 -o ps 文件名
gxpscw -i 图元文件名 -o ps 文件名
也可交互使用。
stnmap
用于将观测站点数据与图形显示时的站点图建立链接以加速显示。语法为:
stnmap -i 站点数据描述文件
也可直接键入 stnmap 交互式操作。
gs -sDevice=vga ps 文件名
3
http://www.lasg.ac.cn
用于显式或打印 ps 文件。
gribmap 形成 GRIB 数据文件的 index 文件(*.idx),用法:
~ [-i fname.ctl] [-v] [-t0] [-0] [-fhr] [-sxxx] [-hxxx]
-v:详细输出信息,
-0:跳过预报记录,仅取诊断记录,
-t0:仅取基准时间与*.ctl 文件中初始时刻一致的记录(如 0,12,24,…,72hr 的预报结果)
-fhr:仅取预报时效为 hr 的记录,
-sxxx:搜索'GRIB'字符串的有效字节长度为 xxx,
-hxxx:跳过 xxx 字节开始搜索,对 NMC 格式的 GRIB 数据,用-hnmc。
gribscan 抽取 GRIB 数据文件的格点信息,用法:
~ [-i ifname] [-o ofname] [-file options] [-processing options] [-display options]
ifname:输入的 GRIB 文件, ofname:无后缀的输出数据文件名,
file options:输出文件格式,
-og:GRIB,
-oa:ASCII(%8g),
-of:32-bit(64-bit CRAY)流式实数格式,
processing options:
-snnn:文件 GRIB 信息的最大字节数,缺省为 500,
-spnnn:选择第 nnn 号要素,
-slnnn:选择第 nnn(hPa)层的要素场,
-stnnn:选择τ=nnn 的要素场,
-hnnn:固定文件头长度 nnn,NMC 为-h81。
display options:
-q:快速输出 gribmap 所需信息,
-q1:1 行的快速输出,
-d:逗号分隔模式,
-v:诊断用详细输出,
-bd:Binary 数据块信息,
-gv:使用 NMC GRIB 要素表输出要素代号,标题和单位,
-gd:由 grid defn sec 输出信息,
-S:无标准输出。
gxtran 显示 GRADS 输出图形文件(单祯或动画),用法:
~ [ -airg ]
-i fname:GRADS 图形文件名,
-a:无需用户击回车键,自动放映文件中的各祯图,
-r:反转背景色,
-g geom:设置 X window 的尺寸 geom(WWW×HHH+X+Y)。
ncdump 在屏幕上输出指定 netCDF 文件的描述,用法:
~ [-c] [-h] [-v var1,...] [-b lang] [-f lang] [-l len] [-n name] [-d f_digits[,d_digits]] file
ncgen 生成 netCDF 文件,用法:
4
GrADS 使用手册
~ [-b] [-c] [-f] [-n] [-o output_file] input_file
-b:生成(binary) netCDF 文件,文件名后缀为.nc,
-c:在屏幕上输出生成 netCDF 文件的 C 语言代码,
-f:在屏幕上输出生成 netCDF 文件的 Fortran 语言代码,
-o outputfile:指定输出的 netCDF 文件名(暗示使用了-b 项),
-n:同-b,输出文件使用废弃的.cdf 后缀,本选项支持向上的兼容性。
gxgif [-i *.gmf] [-o *.gif],将 GrADS 输出的图形文件转换成 gif 图形格式。
§1.4 在 GrADS 环境内的基本命令(在 ga>提示符后输入):
open 数据描述文件名
打开数据文件。
set 各类选项
设置各种环境参数。
display(或 d) 表达式
对表达式处理后进行屏幕图形显示。
clear(或 c)
清屏命令,清除字符窗口和图形窗口的内容。
quit
退出 GrADS 环境,回到 DOS shell 或 Unix shell。
query(或 q)
系统环境设置的查询命令。
define 临时变量名=表达式
在交互方式下定义临时变量场,用于随后的操作。
modify 临时变量名 时间序列类型
对自定义的临时变量声明为气候型的,用于时次代换。
draw 选项
低级绘图指令,直接进行所指定图形元素操作,如绘一些字符串、线条、各种标记等。
enable print 磁盘图元文件名
打开磁盘文件(没有时创建,已存在时刷新),用于存放随后 print 命令转换生成的当前屏幕
图形的图元数据。
print
5
http://www.lasg.ac.cn
将当前图形窗口重的图形转化为图元数据,存放在先前 enable 命令指定的磁盘文件中,每
执行一次 print 即向该文件中附加一幅窗口图形。
disable print
关闭图形硬拷贝输出转换。
exec 批处理命令集文件名<变量 0,变量 1,…>
在 GrADS 环境中集中执行文件名(形式为*.exe)中定义的操作。变量 0……9 用于参量替
换。
run 命令描述文件名
执行文件名(形式为*.gs)中定义的操作。
swap on|off
打开双缓冲区,在动画显示时用于文件交换。
collect cnum expr
采集站点或时间序列数据作为资料阵列用于绘剖面图。与 coll2gr 函数联合使用。其中:
cnum:采集资料序号(0~31), expr:任意表达式,设为'free'表示空阵列并释放内存。
printim filename options
在 GrADS v1.8 以上版本有效,可在批处理文件中使用。其中:
filename:输出的目标文件名,若已经存在,则将覆盖;
options 有多个选项时可以任意次序排列。选项如下:
gif:输出 GIF 格式文件(缺省为 PNG 格式),
black:采用黑色背景(缺省为当前的 display 设置),
white:采用白色背景(缺省为当前的 display 设置),
xNNN:水平尺寸为 NNN 个象素,
例: 输出 1000x800 的 PNG 图象:
yNNN:垂直尺寸为 NNN 个象素。
printim out.png x1000 y800
输出 800x600 白色背景的 GIF 图象:
printim gifimage.out gif x800 y600 white
wi filename.fmt
将 GrADS 绘图窗口中的图形直接存成文件,批处理模式下不可用。注意使用本命令时
图形窗口上不应有其它窗口迭放。fmt 代表图形文件的格式,可为多种:
AVS,BIE,BMP,BMP24,CGM,CMYK,DCX,DIB,EPS,EPS2,EPSF,EPSI,
FAX,FIG,FITS,FPX,GIF,GIF87,GRAY,GRADATION,GRANITE,HDF,HISTOGRAM,
HTML,JBIG,JPEG,ICO,LABEL,MAP,MIFF,MNG,MONO,MPEG,MTV,
NETSCAPE,NULL,PBM,PCD,PCL,PCX,PDF,PGM,PICT,PLASMA,PNG,
PNM,PPM,PS,PS2,RAD,RGB,RGBA,RLA,RLE,SGI, SUN,SHTML, TEXT,
TGA,TIFF,TIFF24,TILE,UIL,VICAR,VID,VIFF,X,XC,XBM,XPM,XWD,……
outxwd filename
将 GrADS 绘图窗口中的图形直接存成 xwd (X window dump)文件,批处理模式下不可
用。xwd 格式与 GIF 格式更为匹配,用 ImageMagick 转换更快且效果更佳。注意使用本命
6