grads基本命令

set strsiz hsiz <vsiz> 设置字符串中字符的大小， 设置字符串中字符的大小，hsiz是 是 字符的宽度， 是字符的高度， 字符的宽度，vsiz是字符的高度， 是字符的高度 单位是虚页英寸。 单位是虚页英寸。
eg.
‘set string 7 c 6 0’ ‘set strsiz 0.4 0.4’ ‘draw string 6.06 5.25 trough’
第5章 grads基础绘图指令 drawstring eg.draw string 5.5 6.0 anticyclone setstring color >> string color justification thickness rotation setstrsiz hsiz hsiz vsiz eg. 'set string 'setstrsiz 0.4 0.4' 'draw string 6.06 5.25 trough' drawline x1 y1 x2 y2 x1 y1 x2 y2 set line" setline color
设置命令如下： 设置命令如下： set clip xlo xhi ylo yhi 其中， 其中，xlo，xhi，ylo，yhi是剪裁区在 ， ， ， 是剪裁区在 实页中的英寸坐标。 实页中的英寸坐标。 用户可以设置一块剪裁区（ 用户可以设置一块剪裁区（clipping area ）用于绘制一些基础图形，设置剪裁区 用于绘制一些基础图形， 这些基础画图功能只能在剪裁区内进行。 后，这些基础画图功能只能在剪裁区内进行。
2. draw line x1 y1 x2 y2
从坐标点x1， 到 ， 从坐标点 ，y1到x2，y2 画一 条直线，可以使用“ 条直线，可以使用“set line”命令来 命令来 设置线条的性质（单位：英寸）。 设置线条的性质（单位：英寸）。

set rbcols c1 c2<c3…cn>
设置彩虹色颜色序列值（缺省：9，14，4，11，5，13，3，10，7，12，8，2，6）
set rbcols <auto>
打开系统缺省设置的彩虹色设置
set rbrange low high
彩虹色颜色对应等值线开始和结结束的范围
set grid on|off|<style>
清干净图形显示区并重新初始化，但1）不关闭任何文件，2）不取消变量定义，3）不改变set display设置
option＝events; graphics; hbuff; norset
reinit
与reset相似，但关闭所有文件并取消所有定义
quit
退出
!shell-command
在GrADS命令交互模式下运行外部系统命令
用三次样条或线性插值的方式计算细网格的值
set cint value
设置固定等直线间隔值
set cmax value
设置最大等直线的值，超过的不再画出
set cmin value
设置最小等直线的值，小于的不再画出
set clevs lev1 lev2…
按lev1 lev2…指定值绘等直线
set ccols col1 col2…
按顺序执行fname文件中的命令组
clear <option>
c <option>
清干净图形显示区，但不关闭当前数据文件可带参数：命令可简写为“c”并可带以下参数：
events
清空事件缓冲区
graphics
清空图形缓冲区，但不清空组件缓冲区
hbuff
在双缓冲区模式下，清空显示缓冲区

第一个命令：这样将得到一个动画显示的序列，现在输入：claerset lon -90set lat -90 90set lev 1000 100set t 1d td u该情况下我们设置y（纬度）和z（高度）变化，所以我们得到一个竖直切面图。

我们显示了两个变量，他们重叠在一起。

你可以让任意多的变量重叠在一起，只要你不使用clear命令就行了。

另一个例子，这种情况X和T是变化的（哈莫图）：cset lon -180 0set lat 40set lev 500set t 1 5d z现在已经知道怎样使用选定的部分数据绘图了，下面我们学习怎样操作数据。

设置维数环境为Z,Y变化：clearset lon -180 0set lat 0 90set lev 500set t 1现在假设我们想显示华氏温度代替开尔文温度。

可以这样转换：display (t-273.16)*9/5 32任何由标准的运算符,-,* 和/ 组成的表达式都是允许的，并且操作数可以包含常量，变量，或者函数。

一个包含函数的例子：d sqrt(u*u v*v)有一个函数用来计算风的级数。

d mag(u,v)另一个内建函数计算平均值：clear d ave(a,t=1,t=5)这种情况我们可以计算5天的平均。

我们也可以从数据中移除平均值（距平值）：d z-ave(z,t=1,t=5)也可以在x方向作平均并求距平：cleard z-ave(z,x=1,x=72)也可以做时间差分：cleard z(t=2)-z(t=1)完整规范的变量名是：name.file(dim |-|=va lue,…) 如果我们打开了两个文件，也许一个是模式输出，另一个是分析，我们应该区分用如下方法二者：display z.2-z.1另一个内置的函数通过有线差分计算水平涡度相关cleard hcurl(u,v)还有另外一个计算数值方向的质量积分：cleard vint(ps,q,275)这儿我们计算了可降水量（单位mm）现在来讨论控制图形输出的话题。

每间隔val条等直线作标注
set cthick thickness
设置等直线粗细,thickness=1~10
set csmooth on|off|linear
用三次样条或线性插值的方式计算细网格的值
set cint value
设置固定等直线间隔值
set cmax value
设置最大等直线的值，超过的不再画出
GrADS快速参考手册
(GrADS Version 1.7 beta 7 ) 编：Meier-fleischer@dkrz.de 译： 朱 禾
命令行参数grads -options
-c‘命令组’
启动GrADS后执行该命令组
-b
以批处理方式运行
-l
以landscape 模式运行
-p
以portrait 模式运行
是否标注台站号
set wxcols c1 c2 c3 c4 c5 c6
设置气象符号底颜色
坐标轴及边框设定参数
设置x轴标注的范围和<增量>
设置y轴标注的范围和<增量>
按lab1 lab2…标注x轴
按lab1 lab2…标注y轴
设置x轴标注的间隔
设置y轴标注的间隔
x轴y轴颠倒（缺省：x轴水平；y轴垂直）
modify varname<time type>
定义变量为气候序列平均，时间类型如下：
seasonal
月平均
diurnal
一日内多时刻平均
query <option>
q <option>
查询命令及其参数：
files
列出所有已打开的数据描述文件
file n

GrADS绘图软件实用手册2002年1月目录第一章 GrADS绘图软件概述1.GrADS绘图软件简介2.Internet上的GrADS资源2.1GrADS在Internet上的主页2.2 windows环境下GrADS资源3.GrADS绘图软件的安装（windows环境）3.1在windows环境下安装GrADS软件包3.2X server 的安装第二章 GrADS绘图模板1.GrADS示例演示1.1 启动GrADS1.2 退出GrADS1.3 示例演示GrADS命令的使用2.GrADS绘图模板3.GrADS模板的高级应用3.1GrADS描述语言3.2GrADS高级模板的应用第三章 GrADS数据格式1.格点数据描述文件1.1 数据描述文件各项解释1.2 生成model.le.dat和model.le.ctl文件的程序代码片段2.站点数据的格式附录1.如何精确控制图形输出的尺寸—Landscape纸型2.台站资料的显示3.Linux环境下的安装第一章 GrADS绘图软件概述1GrADS绘图软件简介The Grid Analysis and Display System(GrADS) 是一套应用广泛、使用方便的科学数据绘图软件包。

其主要特点：z GrADS属于自由软件，可以从Internet上免费获得。

z可运行于各种Windows 和Unix工作平台。

z GrADS可用于4D数据的分析。

既经度、纬度、层（气压层、高度层等）和时间/xyzt 4维。

数据可以是格点化的数据或离散点数据。

GrADS特别适用于气象类数据的分析。

但也完全可以用于更广泛类型的数据分析。

z GrADS有多种显示方式：等值线、流线、矢量图、风矢量图、站点填图、折线图、直方图等多种两维图形。

z可处理多种数据格式的数据。

GRIB、NetCDF、HDF-SDS等通用数据格式和系统自定义的一种二进制数据格式。

z采用命令行输入的方式交互式地显示图形。

S7 Graph 的学习心得
关于Graph 中常用的命令
（1）命令S；当步为活动步时，使输出置位为 1 状态并保持。

（2）命令R：当步为活动步时，使输出复位为0 状态并保持。

（3）命令N：当步为活动步时，输出被置为1；该步变为不活动
步时，输出被复位为0。

（4）命令L：用来产生宽度受限的脉冲，相当于脉冲定时器。

（5）命令CALL：用来调用块，当该步为活动步时，调用命令中
指定的块。

（6）命令D：使某一动作的执行延时，延时时间在该命令右下方
的方框中设置。

在"直接"模式用鼠标右键点击动作框，在弹出的菜单中选择插
入动作行。

这种编程的重点就是要能够认清楚监控和互锁条件的区别与应用，如果能够准确的认识和掌握，其实这样的顺序编程也就会了一大半了，没有什么困难的。

顺序控制器中的条件
1．转换条件
2．互锁条件：如果互锁条件的逻辑满足，受互锁控制的动作
被执行。

3．监控条件：如果监控条件的逻辑运算满足，表示有干扰事件
V1 发生。

顺序控制器不会转换到下一步，保持当前步为活动步。

如果监控条件的逻辑运算不满足，表示没有干扰，如果转换条件满足，转换到下一步。

只有活动步被监控。

动作的分类：
∙标准动作
o与interlock 有关或者是无关
∙基于事件的动作
o o 与interlock 有关或者是无关与步的激活与非激活有关
∙计数器，定时器以及算术指令的动作。

grads基本命令d ave(olr,t+0,t+359,12) ave表示求平均，olr 为变量， t+0表示起始点为当前时次，即00z01jul1974；t+359表示终止时次为当前时次后推359个时次，即00z01jul2003；12表示每12个时次取一次数据；打开nc文件：ga->sdfopen查询信息：ga->q filega->q dimga->q ctlinfo写标题：ga-> draw title清除图面：ga-> c同时打开两个文件： ga->open model.ctlga->open model.le.ctl删除所有设置重新回到刚进入GrADS状态：ga->reinit删除open命令后的所有设置：ga->reset由风场导出涡度场：ga->d hcurl(u,v)以分色图形方式输出：ga->set gxout shaded以等值线方式输出（缺省方式）：ga->set gxout contour画图例：ga->cbarn 1 0 （1:相对长短，>1放大；0:水平；1:垂直。

）以箭头方式表示矢量场: ga->gxout vector风矢量场：ga->d u; v; q （显示矢量时，d x分量；y分量<；标量>。

“<>”内的部分只起标颜色的作用。

）以流线方式表示矢量场：ga->set gxout stream风流线场：ga->d u; v; q以WMO风标方式表示矢量场： ga->set gxout barb风标：ga->d u; v; q直接输出网格点数值：ga->set gxout grid设置保留小数位数：ga->set dignum 0设置数字大小：ga->set digsize 0.1ga->set mpdraw on 如为off，不画地图背景（非经纬度数据需此项）ga->set poli on 如为off不画国界省界等。

目录前言 (5)第一章启动GrADS 及打开文件 (5)一.启动GrADS (5)二.打开文件 (6)第二章维数环境 (6)set lat|lon|lev|time val1<val2> (6)set x|y|z|t vall <val2 > (6)第三章显示命令及动画 (7)一、显示命令 (7)二、动画 (7)set loopdim x︱y|z|t (7)set looping on|off (7)第四章图形输出类型 (8)对格点数据： (8)bar: 直方图 (8)barb:风向杆形式绘二维风场 (8)contour:二维等值现图 (8)errbar:单线图及误差分布 (8)frgid:用指定颜色填充二维格点场 (8)fwrite:图形不在屏幕上显示，而是将输出结果存入一个由“set fwrite 文件名”所指定的文件中。

(8)gid:以网格形式在各网点中央标出该点数值 (8)line:单线图 (8)linefill:两单曲线之间填色 (8)scatter:散点图 (8)shaded:二维填色图 (8)stat:输出有关的统计量 (8)stream:流线形式绘二维风场 (8)vector:向量间箭头形式绘二维风场 (8)对站点数据： (8)barb:在各站点绘风向标 (8)findstn:搜索最近的站点（详见描述语言部分） (8)model:以天气图形式将天气观测各分量填放在站点四周 (8)value:在各站点标值 (8)wxsym:绘wx天气符号 (8)Draw wxsym symbol x y size<color<thicks>> (9)对于model,命令为set gxout model ,然后再添图，形式为displayu;v;t;d;slp;delta;cld;wx;vis (9)第五章 GrADS资料格式 (9)格点资料资料描述文件 (9)以下是一个资料描述文件的例子： (9)１、DSET data-set-name (11)２、TITLE string (11)３、UNDEF value (11)4、OPTIONS <keywords> (12)5、XDEF number<LINEAR start increment>或XDEF number<LEVELS value-list> (12)6、YDEF number mapping start<increment>或YDEF number<LEVELS value-list> (13)７、ZDEF number mapping <start increment>或ZDEF number<value-list> (13)８、TDEF number LINEAR start-time increment (14)９、VARS number (14)二. 台站资料资料描述文件 (15)第六章变量名和表达式 (15)一、变量名和表达式 (15)二、表达式 (16)第七章函数 (16)一. aave函数 (16)二．abs函数 (16)三．acos函数 (17)四、asin函数 (17)五、atan2函数 (17)六、ave函数 (17)七、cdiff函数 (17)八、const函数 (18)九、cos函数 (18)十、exp函数 (18)十一、gr2stn函数 (19)十二、hcurl函数 (19)十三、hdivg函数 (19)十四、log函数 (19)十五、log10函数 (19)十六、mag函数 (19)十七、maskout函数 (20)十八、oacres函数 (20)十九、pow函数 (20)二十、sin函数 (20)二十一、skip函数 (20)二十二、smth9函数 (21)二十三、sqrt函数 (21)二十四、stnave函数 (21)二十五、stnmin函数 (21)二十六、stnmax函数 (21)二十七、tan函数 (22)二十八、tloop函数 (22)二十九、tvrh2q函数 (22)三十、tvrh2t函数 (22)三十一、vint函数 (23)第八章图形要素设置 (23)对于图形类型为contour起作用的设置 (23)set ccolor color——设置等值线颜色 (23)set ccolor rainbow——设定等值线颜色用彩虹表示 (23)set cstyle style——设定等值线线型 (23)Set cthick thckns——设定等值线线宽hckns(取值范围1～20) (23)Set cterp on｜off——设置样条差值光滑开关 (23)Set clab on｜off｜forced｜string｜auto——控制等值线的标记方式。

第一讲 GrADS简介一、G r A D S的应用领域及其功能*GrADS的全称：“The Grid Analysis and Display System”*应用领域：可在UNIX工作站以及个人微机上进行地球科学领域的数据资料分析和绘图*功能:对数据进行访问、分析和绘图1.可以根据需要绘制单线图、直方图、等值线图、填色等值线图、流线图、矢量图、站点模型图等各类图形.2.用描述语言编程,达到理想的绘图效果3.可以把在GrADS中绘制的图形以文件的形式保存起来,以备对其进行编辑4.调用GrADS的内部函数，能够对数据进行某些特定的计算，然后输出计算结果二、GrADS中常用的基本概念*数据格式：GrADS能够识别的数据为二进制无格式直接或顺序记录格式，该种格式数据的生成可以通过Fortran语言编程来实现。

*数据类型：格点数据（NCEP/NCAR的再分析资料）、站点数据（站点实测资料）、Grib数据（NMC产品）。

*维数环境：GrADS的操作对象为4维的数据集（4D data set）,包括空间三维（纬度、经度、高度）和时间一维。

可以固定其中的一维或者几维以获得低于四维的数据子集。

此概念是对于格点资料而言的。

维数环境的定义可以在两种坐标上进行。

一种是地球坐标（world coordinate），以经纬度为度量单位；一种是格点坐标（grid coordinate），以网格点数为度量单位。

*几种文件类型：*.dat―数据文件*.ctl—原始数据描述文件*.gs―GrADS控制文件，用命令run执行之*.exe―GrADS在DOS环境下的各种执行文件三、启动和退出GrADS*启动GrADS的两种方式：1.在dos环境下直接输入grads命令即：切换到MS-DOS方式，进入到F:\pcgrads\msdos子目录下输入：grads回车2.从windows桌面上的“我的电脑”进入F:\pcgrads\子目录，然后双击g.exe图标*注意：在启动GrADS时，系统会问你以何种方式进入。

地图投影设置：
set mproj proj 设置当前地图投影方式。
常用proj取值包括： latlon：缺省设置，用固定的投影角进行Lat/lon投影
； scaled：用不固定的投影角进行latlon投影，地图比例
失效； nps | sps：北半球|南半球极地投影； robinson：robinson投影，x：-180 180；y：-90 90 off：不画出地图，坐标轴也不代表lat/lon。
对图形为矢量图起作用的 设置：
set arrscl size <magnitude> 设置矢量箭头的长度为size（
英寸）， 选项magnitude的值为矢量的
大小（例10m/s）。
set arrowhead size 设置箭头大小，缺省为0.05，
取0不画头，取负与矢量值 成比例。
例：set arrowhead 0.06
，后者是网格坐标。 例：set lon -180 0
set lat 0 90 set lev 500
set t 1 换成set time 00:00z01Jan2001
图形类型设置：
当维数环境确定后，缺省情况下，一维变量输出的图形为单线 图，二维变量为等值线图，若想进行特殊指定，键入命令：
set gxout graphics-type（图形类型） 对格点数据： contour: 二维数据等值线图（缺省） shaded: 二维填色等值线图 grid: 二维场不绘图，以网格形式在各网格点中央标出该场点
数值 vector: 矢量箭头形式绘二维风场 stream: 流线形式绘二维风场 barb: 风向杆形式绘二维风场 bar: 对一维场不绘单线图，而绘直方图 line: 对一维场绘单线图（缺省） 例：set gxout shaded

关于数据文件model.le.dat的具 的具 关于数据文件 体说明，可以打开model.le.ctl文件 体说明，可以打开 文件 查看。 查看。
2. 示例演示
启动GrADS后，选择图形输出窗口 后 启动 的方式（ 或 ），然后在命令窗口ga-> ），然后在命令窗口 的方式（L或P），然后在命令窗口 提示符后输入命令， 提示符后输入命令，通过图形输出窗口 可以观察输出结果。 可以观察输出结果。
open <路径 数据描述文件名 路径>数据描述文件名 路径
open命令用于打开 命令用于打开GrADS的数 命令用于打开 的数 据文件，启动GrADS后首先需要打开 据文件，启动 后首先需要打开 至少一个数据描述文件 命令如下： 数据描述文件， 至少一个数据描述文件，命令如下： open filename
define 临时变量名 表达式 临时变量名=表达式
define命令用于定义新的变 命令用于定义新的变 量，所定义的新变量可以用于随 后的表达式中。 后的表达式中。
modify 临时变量名 时间序列的 类型<seasonal diurnal> 类型 modify命令可将自定义的变 命令可将自定义的变 量声明为气候值， 量声明为气候值，用于后面的时 次代换。 次代换。
（以下set命令是设置维数环境， 关于维数环境的介绍详见第4章）
（显示该变量） （清屏）
set lon -180 0 dz c set lat 0 90 dz quit
（重新设置维数环境）
（退出GrADS系统）
实际操作后， 实际操作后，在图形输出窗 口观察输出结果。 口观察输出结果。比较不同的维 数环境设置下， 数环境设置下，输出结果有何不 关于set命令的应用将在第 命令的应用将在第4 同？关于 命令的应用将在第 章具体介绍。 章具体介绍。

set strmden [1~10]5
矢量长度比例
set arrscl inch value
1杠=10 knots!
display u*2.5;v*2.5
23
LASG动力论坛版
四. 图形显示(D)1
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
14
LASG动力论坛版
四. 图形显示(D)1
柱间距
set bargap percent
— 1-D单变量序列图
颜 色
set ccolor #color
基 线
set barbase top ~ ~ bottom ~ ~ value
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
画布/ 虚页/ 绘图区大小
— 格点绘图
命令顺序 set display draw print
– 图形类型
等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
set xlabs b1|b2|b3|… set ylabs b1|b2|b3|…
16
LASG动力论坛版
四. 图形显示(D)1
— 1-D双变量序列图
散点图 set gxout scatter d var1;var2 双线填色图 set gxout linefill d line1;line2 颜 色 set lfcols #c1 #c2
set ccolor rainbow ~ ~ revrain ~ ~ #color

1.ctl文件的编写dset d:\test\test1.dat 数据文件路径名和文件名title monthly precipitation data 数据标题options 365_day_calendar 特殊格式说明365_day_calendar/template(多个文件) undef -9999.0 缺测值，若有多个则用fortran改写成同一个缺测值xdef 144 linear 0 2.5 x维数格点数/排列方式/起始值/间隔*linear线性/levels列举ydef 73 linear -90 2.5 y维数zdef 5 linear 1000 850 700 500 200 z维数tdef 24 linear 00z01Jan1979 1mo 时间维数*mo月yr年vars 1 变量数precip 1 99 *precipitation data 变量名/层数/排列顺序/说明endvars*变量循环顺序：x经度y维度z高度v变量t时间*注释行第一列用‘*’，注释行不能出现在变量列表中2.维数环境设置set lon/lat/lev/time var1 <var2> *实际值set x/y/z/t var1 <var2> *格点数值*两种坐标可以混用3.图形文件的保存gxprint/printim **.png/jpg/pdf/eps…. whiteenable print **.gmfprintdisableprint4.图形类型的设置set gxout typecontour 二位等值线图*默认shaded 二维填色等值线图vector 矢量箭头二维风场图*d u;v时为默认bar 直方图line 折线图*默认…………………5.图形要素设置gxout=contour or shadedset cint value 设置等值线间隔set clevs lev1,lev2,……设置特定等值线set ccols col1,col2,……设置对应于set clevs的等值线的颜色cbarn size 0/1 x y 设置色标gxout=line or contourset ccolor color 设置线条颜色set cstyle style设置线条样式1实线2长虚线3短虚线……set cthick thickns 设置线条粗细1~10set cmark marker 设置数据点标记0无标记*line时有效gxout=barset bargap val 设置直方条间隔*0~100百分比，0为无间隔set barbase val/bottom/top 设置直方条起始值gxout=vectorset arrscl size <magnitude> 设置箭头长度set arrowhead size 设置箭头大小，缺省为0.056.gs文件的编写将命令用单引号括起来7.变量define varname=expression 定义变量‘d air(t=3)’‘d air(t=3,z=5)’‘d air(lev=1000)-air(lev=500)’设置变量维数8.平均函数ave(expr,dimexpr1,dimexpr2,<timeint>) 一维平均aave(expr,dimexpr1,dimexpr2, dimexpr3,dimexpr4,) 二维平均*不能对z、t *需将维度设置成点9.地图投影设置set mproj proj 设置地图投影方式latlon 缺省nps 北半球极地投影sps 南半球极地投影set map color style thickness 设置背景地图线条颜色、线型、宽度10.绘图区域设置set parea xmin xmax ymin ymax11.坐标要素控制set frame on/off/circle 绘制图形边框/不绘制/绘制圆形边框set grads on/off 标记grads及绘图时间draw title <string> 绘制标题draw xlab/ylab <string> 绘制x轴/y轴图注set xyrev on/off xy轴互换set zlog on/off z方向用对数坐标set xlint/ylint value 设置x/y轴坐标刻度间隔set vrange val1 val2 设置y轴范围set vrange2 val1 val2 设置x轴范围set xlopts/ylopts color thickness size x/y轴刻度值的颜色、线宽、大小*缺省1,4,0.12 set tlsupp year 去掉年份显示12.颜色的定义‘define_colors’使用内部默认颜色标号13.基本绘图命令set line color style thickness 设置颜色、线型、线宽draw line x1 y1 x2 y2 绘制直线draw rec/recf x1 y1 x2 y2 绘制不填色/填色矩形draw mark type x y size 绘制标记*1叉2空心圆3实心圆4空心方框5实心方框draw string x y string 绘制字符串set strsiz hsiz <vsiz> 设置字符大小14.文件操作(1)输出二进制变量设置维数环境‘set gxout fwrite’‘set fwrite ***.dat’‘d var’(2)读写有格式文件result=read(filename)result 第一行：返回码（0 OK，1打开错误2文件结束8文件状态为写入9I/O错误）result 第二行：读取的文件内容Text = sublin(result,2) 将result相应行的内容赋值给Textaa=subwrd(Text,#n) 将Text的第n个值赋值给aaresult=write(filename, text, <append>)15.维数查询q w2xy lon lat 将经纬度坐标转换为图形窗口坐标返回值：X=*** Y=*** q xy2w X Y 将图形窗口坐标转为经纬度坐标返回值：Lon=*** Lat=***q w2gr lon lat 将经纬度坐标转换为格点坐标返回值：Xdim=*** Ydim=***q gr2w xdim ydim 将格点坐标转换为经纬度坐标返回值：Lon=*** Lat=***q xy2gr X Y 将图形窗口坐标转换为格点坐标返回值：Xdim=*** Ydim=***q gr2xy xdim ydim 将格点坐标转换为经纬度坐标返回值：X=*** Y=***16.气候平均设置季节变量‘set t 1 12’‘define climate=ave(var,t+0,t=tmax,12)’‘modify climate seasonal’滑动平均‘set t 1 tmax’‘define running_mean=ave(var,t-n,t+n)’17.站点资料绘图(1)转换生成数据文件及格点文件数据文件的顺序为先循环站号再循环时间，每条记录包含站号-纬度-经度-时间-nlev-flag-数据，每循环完一个时间后要加一个不写数据且nlev=0的记录(2)编写ctldtype stationstnmap ***.map**变量层数为0(3)生成.map映射文件!stnmap –i */*/*.ctl(4)插值‘open ***\grid.ctl’‘open ***dat.ctl’‘define varname=oacres(grid,var.2)’(5)绘图‘cnbasemap varname’18.绘制路径图。

w2xy lon lat
给出经度－纬度坐标值(lon,lat)，返回xy坐标值
w2gr lon lat
给出经度－纬度坐标值(lon,lat)，返回网格坐标值
gr2w v1 v2
给出网格坐标值(v1,v2)，返回经度－纬度坐标值
gr2xy v1 v2
给出网格坐标值(v1,v2)，返回xy坐标值
ll2xy lon lat
-C
以两位数记年时，自动设定世纪值
一般性设置
列出帮助选项
set grads on|off
打开/关闭 显示“GrADS” 标志
设置显示模式<和色彩方式>为：
grey|greyscale
用灰度方式
color <black|white>
以彩虹模式绘图
set frame <option>
设置图形边框绘图方式为：
以一个矩形包围所有图形
off
不画边框
circle
对于其它投影方式画矩形边框，而极射投影画圆形边框
set background index
设置绘图区背景颜色
display expression
d expression
显示expression代表数据的图形
open control_filename
打开一个数据描述文件，可指明路径
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息

gxtran
Displays metafiles
stnmap
Maps station data
wgrib
See /wgrib.html
GrADS 英文主页：/grads/
§1.2 基本概念
按磁盘文件记录格式，GrADS 软件包的文件分为如下几类：
gxps，gxpsc，gxpscw
将 GrADS 内 print 命令产生的图元文件分别转换为单色（灰度）、彩色以及反白（即白底）
彩色的 postscript 格式文件，用于硬拷贝输出。语法为：
Technical Report
GrADS究所
大气科学和地球流体力学数值模拟国家重点实验室
（ LASG ）

关于本手册
本手册是在 1995 年张洪编写的《GrADS 气象图形系统使用方法》及其它相关文档的基础上 汇编而成。主要介绍了 GrADS 的基本用法、使用技巧、及与 GrADS 相关的绘图技术。本 手册将根据软件发展不定期进行修订。
1

令。详见§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 内定，文件名随用户自定，
(2) .ctl 原始数据描述文件。该文件为纯 ASCII 文件，用以描述原始数据集的基本信息，包 括数据集文件名、数据类型、数据结构、变量描述等等，详见§4。在 GrADS 环境中至 少得首先打开（open）一个数据描述文件，以便后续的操作有数据对象。

ga-> q files File 1 : TC0801-StrkPro024 2008041714 Descriptor: strike024_041714.ctl Binary: strike024_041714.dat File 2 : TC0801-StrkPro024 2008041720 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat
命令：query dims
ga-> q dims Default file number is: 1 X is varying Lon = 100 to 150 X = 1 to 101 Y is varying Lat = 0 to 40 Y = 1 to 81 Z is fixed Lev = 1000 Z = 1 T is fixed Time = 00Z19JUL2005 T = 1
命令：open filename2.ctl
ga-> open strike024_041720 Scanning description file: strike024_041720 Data file strike024_041720.dat is open as file 2
查询命令
命令：query files
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
• 3、图形显示类型及绘图设置：
set cmax 90 不画超过90（不含）以上的等值线 set clevs 495 523 534 556 560 564 572 584 588 set ccols 1 2 3 4 5 当指定颜色少于等值线条数时，最后一种颜色适用于超出的所有等值线。

气象图形软件应用——G r A D S绘图软件的使用成都信息工程学院大气科学系巩远发主要内容一、GrADS的一般介绍二、GrADS的数据结构详解三、GrADS数据数据处理和绘图详解四、GrADS绘图编程（脚本）语言第一讲 GrADS的一般介绍什么是GrADS相关的关于GrADS的网络资源启动和退出GrADS最基本的GrADS操作（命令）GrADS使用文件的分类GrADS的主要功能介绍1. 什么是GrADS？G r A D S——Grid data Analysis and Display System 格点数据分析和显示系统GrADS（Grid data Analysis and Display System）是一个全32 位的交互操作的气象格点数据和站点数据的分析和显示环境。

系统具有气象数据分析功能强、地图投影坐标丰富、高级编程语言使用容易、图形显示快速，并具有彩色动画功能等特点。

目前是国内外气象数据显示的标准平台之一。

2. 相关的关于GrADS的网络资源WWW Sites/grads/head.html/home.htmlGrADS论坛动力论坛专业绘图软件格点资料下载http://www.ecmwf.int//cdc/reanalysis/3. 启动和退出GrADSgrads 启动GrADSga-> GrADS的系统提示符help给出GrADS的主要命令quit 退出GrADS4. 最基本的GrADS操作Open :打开一个网格点或站点数据文件供GrADS使用(sdfopen: 打开netCDF数据)Query (或仅用q) :查看当前GrADS环境的各种参数Display (或仅用d) :显示(绘制)一幅GrADS “表达式”的图形Set : 设置绘图的操作方式(包括绘什么图，在那里绘图，怎样绘图……)Clear (或仅用c) : 清除当前显示的图形5. GrADS使用文件的分类数据文件: 二进制无格式记录的原始数据，可以是格点的，也可以是站点的。