grads常用命令

第一个命令：

这样将得到一个动画显示的序列，现在输入：

claer

set lon -90

set lat -90 90

set lev 1000 100

set t 1

d t

d u

该情况下我们设置y（纬度）和z（高度）变化，所以我们得到一个竖直切面图。我们显示了两个变量，他们重叠在一起。你可以让任意多的变量重叠在一起，只要你不使用clear命令就行了。

另一个例子，这种情况X和T是变化的（哈莫图）：

c

set lon -180 0

set lat 40

set lev 500

set t 1 5

d z

现在已经知道怎样使用选定的部分数据绘图了，下面我们学习怎样操作数据。设置维数环境为Z,Y变化：

clear

set lon -180 0

set lat 0 90

set lev 500

set 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方向作平均并求距平：

clear

d z-ave(z,x=1,x=72)

也可以做时间差分：

clear

d z(t=2)-z(t=1)

完整规范的变量名是：name.file(dim |-|=va lue,…) 如果我们打开了两个文件，也许一个是模式输出，另一个是分析，我们应该区分用如下方法二者：

display z.2-z.1

另一个内置的函数通过有线差分计算水平涡度相关

clear

d hcurl(u,v)

还有另外一个计算数值方向的质量积分：

clear

d vint(ps,q,275)

这儿我们计算了可降水量（单位mm）现在来讨论控制图形输出的话题。目前为止，我们一直让GrADS 选择默认的画线间隔。我们可以这样自己设定：

set cint 30

d z

也可以控制线的颜色：clear

set ccolor 3

d z

可以使用交替的方式显示数据：

clear

set gxout shaded

d hcurl(u,v)

这看起来不很平滑，可以使用三次平滑：

clear

set csmooth on

d hcurl(u,v)

我们也可以使用不同的图形样式

set gxout contour

set ccolor 0

set cint 30

d z

使用注解：draw title 500mb Heights and Vorticity

风矢量：

clear

set gxout vector

d u,v

这里我们显示两个表达式，第一个是U分量矢，第二个是V分量矢。我们也可以使用第三个字段使图像彩色化

d u;v;q

或者：

d u;v;hcurl(u,v)

你也可以显示任何区域以显示伪矢量：

clear

d mag(u,v);q*10000

此处U分量是风速和V分量是湿度。我们也可以查看流线（并且着色）

clear

set gxout stream

d u;v;hcurl(u,v)

或者我们可以查看实际的格点值

clear

set gxout grid

d u

我们可能希望改变地图背景：

clear

set lon -110 -70

set lat 30 45

set mpdset nam

set digsize 0.2

set dignum 2

d u

改变投影：

set lon -140 -40

set lat 15 80

set mpvals -120 -75 25 65

set mpproj nps

set gxout contour

set cint 30

d z

我们让grads可以操作纬度从140W到40W之间，经度从15N到80N之间的数据。但是我们告诉它使用从120W到75W和25N到65N极坐标投影绘图。多余的部分被投影程序修剪掉了。

该简单教程到此结束。此时，你可能希望深入查看数据集，或者你可能希望浏览GrADS文档并且尝试此处描述的其它选项。英文原文并没有完全按原文翻译。

reinit 初始化、关闭文件

q/query file 查询文件,xyz格点数，时间，变量，（变量层<小于z>，）

q ctlinfo 路径名字缺色值xyz格点数，初始值，steps

open 打开文件

d/display (temp) 绘图输出

c/clear 清除

draw title (name)给图添加名字,如果加上反斜杠,名字分两行

保存文件

enable print >

print

disable print

格式转换,gv32,打开文件.gmf

quit 退出系统

维数环境设置

set lat纬度（负南纬）/lon经度/lev（海平面之上百帕以下位势米）/time val1

set time 12zsep

格点坐标

set x/y/z/t val1 如，set x 10 30 set t 1 40

改变缺省图形处处类型

set gxout graphics_type

set z 40（不用此命令，默认输出第一层）

set gxout shaded（二维填色图）（用于等温图）

添加温度与颜色的对照图

cbarn/cbar

d temp