内容发布更新时间 : 2024/5/17 19:07:00星期一 下面是文章的全部内容请认真阅读。

win32下安装有GrADS的所有可执行程序；dat目录下有运行GrADS所需的字库和地图数据文件；doc目录下存放帮助和说明文件；lib目录下存放了一些有趣的模板。 GrADS常用命令列表 Grads 经典GrADS，支持Grib,NetCDF,IEEE格式的数据。 GradsC 压缩模式的GrADS，占内存最小。 GradsGUI （假）图形界面的。 GradsHDF 与经典GrADS相当，多增加了HDF-SDS格式。 GradsNC 与经典GrADS相同。 gv32 Windows下看图（metafile格式）及格式转换工具。 rxvt 在Windows下打开一个命名行窗口。

第三章 GrADS绘图模板

初学者可以通过“初学指南”既参考GrADS软件的“Tutorial”联机帮助文件开始学习。

GrADS并不直接使用“数据文件”，而是通过“描述文件”间接使用“数据文件”。GrADS中“打开一个数据文件”既是指打开一个数据描述

文件。

关于GrADS数据格式和用户如何生成该格式的数据，将有下面章节讲述；关于GRIB、NETCDF等数据格式的使用也将有专门论述。

12

以下以第一组数据说明GrADS的使用。在使用数据前，先将数据，Model.ctl和Model.dat文件存于硬盘上的某一目录下，如C:\\OpenGrADS\\Contents\\Resources\\SampleDatasets。首先来了解一下数据的内容。Model.le.dat是一组模式输出的全球5天数值预报结果，包括了多个要素、多层，按经纬度网格存放的数据。以下是Model.ctl清单：

ps 0 1 Surface pressure [hPa] ua 7 33 Eastward wind [m/s] va 7 34 Northward wind [m/s] zg 7 7 Geopotential height [m] ta 7 11 Air Temperature [K]

hus 7 51 Specific humidity [kg/kg]

ts 0 11 Surface (2m) air temperature [K]

pr 0 59 Total precipitation rate [kg/(m^2*s)]

1. GrADS示例演示

目的：

? 学习GrADS的基本使用方法。 ? 介绍GrADS的常用命令。

1.1 启动GrADS

? 打开一个MS-DOS窗口。 ? 在MS-DOS窗口输C:\\OpenGrADS\\Contents\\Resources\\SampleDatasets

入命令：cd

13

? 再输入命令: grads， 显示如下提示：

首先显示GrADS的一般信息，最底行提示是用“Landscape”模式（11X8.5英寸）或portrait（8.5 X11英寸—GrADS中长度和大小单位都用英寸）模式显示图形输出窗口。输入L回车或回车是“Landscape”模式；输入P回车是用portrait模式。之后进入GrADS的命令交互模式，等待用户输入命令。此时你的桌面应如下图所示的特征。若想启动GrADS直接进入命令交互模式，输入命令：grads –l，直接进入“Landscape”模式或grads –p，直接进入portrait模式。要获得GrADS命令行帮助，请运行grads –help。

grads -b 以批处理方式运行，不显示图形输出 grads -l 以“Landscape”模式运行。 grads -p 以“portrait”模式运行。

grads -c 进入GrADS后，执行批处理命令文件，显示图形输出。 grads -cl 以上选项可联合使用。

14

左上角一个窗口是图形显示区，GrADS的所有图形输出结果在此窗口下显示；下部为原打开的MS-DOS窗口，此时正在运行GrADS，提示显示“ga->”，说明正处于GrADS命令等待状态。GrADS的所有命令都只能通过该窗口输出，所有文字信息也都由此窗口输出。GrADS命令采用以下格式：

ga->命令 <参数<参数， ，>， ，>

输入的参数不包括“< >”符号，“< >”内的部份可以省略。所为的命令交互模式 既是在GrADS命令提示符下，一步步输入各种GrADS命令产生各种图形。

1.2 退出GrADS

ga->quit

15

1.3 示例演示GrADS命令的使用。 例1 ga->open model.le.ctl 进入的第一步，打开一个描述文件（扩展名ctl可省略），GrADS并不直接使用“数据文件”，而是通过“描述文件”间接使用“数据文件”。 ga->q file 显示GrADS打开数据的内容，显示结果应与描述文件说明一致。(可以直接运行ga->q显示q 命令所代的参数及多项功能解释。) ga->d ps 显示地面气压（1987.1.2.0hr）如上图。“d”是“显示”命令，“ps”是model.le.ctl文件中定义的变量名称，代表地面气压。 例2

紧接上例输入以下命令：

ga->c 清除图面。如不清除图面，GrADS后续显示的

图形将与已存在的图形产生叠加。

ga->set lat 40 ga->set lon -90 ga->set lev 500

ga->set t 1 GrADS视所有物理量为四维数据，以

lon|lat|lev|time—称为“世界坐标“或x|y|z|t—称为“网格坐标”。两种方式表示。缺省采用第一种方式，水平（lon|lat）自动认为是经度/纬度坐标。

以x|y|z|t方式表示维数时，每一维都是一组从1开始

的序列数，对应网格的序号。本例中：x从1到73（代表从00到3600共72个格点，第73个格点即第1个格点）；y从1到46（代表从-900到900共46个格点）;z从1到7（从1000 hPa开始，共7层）;t从1到5（5天）运行“open model.le.ctl”命令后，lon|lat或x|y是可变的（lon:00~3600; lat:-900~900）。而垂直和时间维取固定值，lev=1000hPa或z=1; time=1987.1.2.0hr或t=1。

“set 维数参数 数值1 <数值2>”命令改变当前

维数设定值。当取“数值1 数值2”时，表示该维是变化的；而只取一个值时，表示该维取固定值，对于网格坐标可以是分数值，系统会自动作内插。维数参数一直保持到再次设定时都有效。要察看当前维数情况，运行：

ga->q dims

ga->d z 显示位于(400N,900W)，500hPa层，1987.1.2.0hr一点的位

势高度。

16