Grads在气象中主要配合Fortran使用，Fortran负责计算处理数据，Grads负责画图，但由于Fortran能处理的数据类型有限制，因此某些时候，直接使用Grads进行简单的计算反而是一个更方便的选择。那么Grads里的数据处理该从哪些方面入手呢？一起来学习吧！

1.数据描述文件

Grads并不直接读取或处理数据文件，而是通过读取数据描述文件(ctl格式)，间接使用数据文件。

每个数据描述文件一般包含以下几项：

①被描述的数据文件名(dset)

②此数据描述文件的标题(title)

③所描述数据的类型、格式和选项(dtype、format、option)

④时间、空间维数环境设置(xdef、ydef、zdef、tdef)

⑤变量定义(vars、endvars)

先通过一个气象家园里翻到的一个例子认识一下其用法：

dset G:\FY2E_TBB_IR1_OTG_2010.AWX title 2010 twpice fy2e_tbb options yrev * yrev表示y轴反向 undef 0.0 xdef 1201 linear  45 0.1 ydef 1201 linear -60 0.1 zdef 1 linear 1000 1 tdef 288 linear 00z19JUN2010 1hr vars  1 tbb   1    -1,40,1    Geopotential Heights (m) endvars

(来自家园用户nx花儿-LJQ：bbs.06climate.com/forum.php?mod=viewthread&tid=2528)

其中可选择使用的参数有(可参考官方网址cola.gmu.edu/grads/gadoc/descriptorfile.html)：

dset  path/filename：表示需要打开的数据文件的路径及文件名

title name：此数据描述文件的标题

undef value：定义数据文件中的缺测值

options  keywords：定义了与二进制文件有关的选项。

keywords可选择的有：

sequential：输出格式为顺序记录方式，缺省时为direct直接记录格式。

yrev：Y维反向。

zrev：Z维反向。

big_endian、little_endian：改变二进制位存放顺序，不常用。

byteswapped：反序位存放，不常用。

template：同时描述多个数据文件时使用。

xdef  xnum  linear  xstart increment或

xdef  xnum  levels value_list：xnumber给定X方向格点数，取linear参数时，网格点格距均匀，xstart为起点坐标(负数表示西经)，increment为网格间距；取levels参数时，网格点格距不均匀，枚举法按顺序列出所有网格点(格点值之间用空格隔开)。

ydef ynum mapping ：定义网格点与纬度的映射关系，包括：

linear：线性映射

levels：不均匀格点，枚举表示

gaust62：高斯T62网格

gausr15：高斯R15网格

gausr20：高斯R20网格

gausr30：高斯R30网格

gausr40：高斯R40网格

zdef的设置方式与xdef相似

tdef只有linear参数，没有levels参数

vars num：数据文件中变量的数量，其下一行同时给出变量的描述：

varname levs units descripition：varname变量名称，变量层数，0表示只有一层；units单位；description，描述变量的一段字符串。

BOOK

2.运算符

Grads的主要功能并不是用于计算，因此在计算方面，相对其他语言来说缺乏了很多。

Grads的运算符仅有：+加、-   减、*   乘、/ 除四种，甚至都没有其他语言里必备的幂运算。

在逻辑运算方面，仅在版本2.1.1.b0中针对GRIDDED数据的处理引入了逻辑运算符：

符号	含义
=	等于(== 含义相同)
!=	不等于
>	大于
>=	大于或等于
<	小于
<=	小于等于
&	和( && 含义相同)
|	或(|| 含义相同)

3.函数

Grads的程序脚本是由.gs文件表示的，可以通过编写脚本文件来对数据文件进行计算、绘图等操作，具体可参考官方文档(cola.gmu.edu/grads/gadoc/script.html#intro)。

在计算数据的过程中，你可能用到的函数有：

函数	功能
abs(expr)	绝对值
cos(expr)	余弦
acos(expr)	反余弦
sin(expr)	正弦
asin(expr)	反正弦
tan(expr)	正切
atan2(expr1, expr2)	反正切
exp(expr)	指数
log(expr)	对数
pow(expr1, expr2)	幂函数
sqrt(expr)	平方根
ave (expr, dim_expr1,dim_expr2, tincr, { |-b } )	加权平均
max (expr, dim1, dim2, tincr)	最大值
maxloc(expr, dim1,dim2, tincr)	最大值所在格点
scorr(expr1,expr2, {xdim1, xdim2, ydim1, ydim2| global | g})	空间相关
smth9 ( expr )	加权9点平滑
cdiff ( expr, dim )	中央差分

简单的介绍就到这里，看一百本书不如自己亲手实践一下，文末为大家准备了新手学习练手的材料，欢迎各位自取练习~

GRADS初级入门的文件，希望能提供帮助：by：做个霸气的木头

bbs.06climate.com/forum.php?mod=viewthread&tid=7483

Grads关于变量和表达式的讲义：by：lqy671bbs.06climate.com/forum.php?mod=viewthread&tid=17186

Grads函数：by：wxxxren

bbs.06climate.com/forum.php?mod=viewthread&tid=7329

drawvc 函数自定义参考矢量箭头画矢量图：by：兰溪之水

bbs.06climate.com/forum.php?mod=viewthread&tid=15303

drawfigstr 函数给图标号：by：兰溪之水

bbs.06climate.com/forum.php?mod=viewthread&tid=15203

一元线性回归：by：1092995269

bbs.06climate.com/forum.php?mod=viewthread&tid=19633

ltrend命令：求线性趋势：by：刷牙

bbs.06climate.com/forum.php?mod=viewthread&tid=31103

用ltrend 命令，快速求出数据的趋势线！或者给数据去掉倾向！：by：zl.87625

bbs.06climate.com/forum.php?mod=viewthread&tid=15576

tcorr用法：by：叩容

bbs.06climate.com/forum.php?mod=viewthread&tid=13326

标准差场，相关系数及相关系数场的程序示例：by：Aires

bbs.06climate.com/forum.php?mod=viewthread&tid=13110

超前滞后相关程序：by：Aires

bbs.06climate.com/forum.php?mod=viewthread&tid=14527

Grads打点信度检验的绘制方法：by：游子

bbs.06climate.com/forum.php?mod=viewthread&tid=18475

[求助]还是置信度打点的问题：

bbs.06climate.com/forum.php?mod=viewthread&tid=29909

[求助]Grads中可以定义数组吗？：

bbs.06climate.com/forum.php?mod=viewthread&tid=21262

总温度平流：by：马小邦

bbs.06climate.com/forum.php?mod=viewthread&tid=56117

Grads计算温度平流：by：shayingying

bbs.06climate.com/forum.php?mod=viewthread&tid=35764

直接用Grads计算温度平流和涡度平流的gs：by：Michael_Qichao

bbs.06climate.com/forum.php?mod=viewthread&tid=21680

流线叠加温度平流的计算gs：by：qiubabbs.06climate.com/forum.php?mod=viewthread&tid=94843

温度平流GS：by：巫毒公主

bbs.06climate.com/forum.php?mod=viewthread&tid=2463

[求助]温度平流随时间的剖面图：

bbs.06climate.com/forum.php?mod=viewthread&tid=12727

[求助]关于绘制某测站温度平流随时间变化的垂直剖面图问题

bbs.06climate.com/forum.php?mod=viewthread&tid=94866

Grads中计算水平螺旋度和垂直螺旋度的gs程序：

bbs.06climate.com/forum.php?mod=viewthread&tid=2396

Grads计算水汽通量和水汽通量散度的gs：

bbs.06climate.com/forum.php?mod=viewthread&tid=22268

整层水汽通量散度计算及其单位：by：捷然不同

bbs.06climate.com/forum.php?mod=viewthread&tid=4920

整层水汽通量的计算(根据自己的1961-2014夏季平均)necp资料：by：帅帅f

bbs.06climate.com/forum.php?mod=viewthread&tid=92300

Grads单层水汽通量和水汽通量散度的gs：by：姿姿

bbs.06climate.com/forum.php?mod=viewthread&tid=68183

利用Grads2.1制作散度、垂直速度、水汽通量散度剖面图的gs：by：qiuba

bbs.06climate.com/forum.php?mod=viewthread&tid=61654

FNL风场、温度、相对湿度资料做水汽通量和散度gs：by：思浓思雨

bbs.06climate.com/forum.php?mod=viewthread&tid=53686

计算水汽通量的终极心得及问题：by：泰达巴神

bbs.06climate.com/forum.php?mod=viewthread&tid=39141

解决水汽通量和水汽通量散度量级的问题：by：278803532

bbs.06climate.com/forum.php?mod=viewthread&tid=24090

分享水汽通量合成分析脚本：by：陌小默的

bbs.06climate.com/forum.php?mod=viewthread&tid=7290

水汽通量及水汽通量散度程序：by：wuxiao

bbs.06climate.com/forum.php?mod=viewthread&tid=36576

Grads计算水汽通量和散度：by：帅帅f

bbs.06climate.com/forum.php?mod=viewthread&tid=92277

[求助]水汽通量散度的量级问题：

bbs.06climate.com/forum.php?mod=viewthread&tid=1697

Grads计算湿位涡全套程序：by：mofangbaobbs.06climate.com/forum.php?mod=viewthread&tid=9884

湿位涡剖面图gs示例：by：Mouse

bbs.06climate.com/forum.php?mod=viewthread&tid=1864

湿位涡的GS程序：by：蓝月亮bbs.06climate.com/forum.php?mod=viewthread&tid=2395

广义位温和广义湿位涡：by：平流层的萝卜

bbs.06climate.com/forum.php?mod=viewthread&tid=17828

位温位涡计算及垂直剖面图：by：wu460748616

bbs.06climate.com/forum.php?mod=viewthread&tid=96016

分享用fnl资料计算等熵位涡的gs：by：风荷1

bbs.06climate.com/forum.php?mod=viewthread&tid=13259

假相当位温gs！个人总结修改增值版：by：ljh110011

bbs.06climate.com/forum.php?mod=viewthread&tid=45075

计算假相当位温的gs文件：by：liantai4567

bbs.06climate.com/forum.php?mod=viewthread&tid=2927

假相当位温Θse计算—GS：by：zsmnantong

bbs.06climate.com/forum.php?mod=viewthread&tid=1913

发一个假相当位温的剖面图gs文件，供大家交流：by：carlye2002

bbs.06climate.com/forum.php?mod=viewthread&tid=2449

Grads绘降水等值线图，高空形势图，湿位涡和假相当位温图：by：20111367047

bbs.06climate.com/forum.php?mod=viewthread&tid=35318

全总指数TT.gs：by：movlan

bbs.06climate.com/forum.php?mod=viewthread&tid=98316

锋生函数.gs：by：movlan

bbs.06climate.com/forum.php?mod=viewthread&tid=97740

锋生函数小程序：by：★翼★

bbs.06climate.com/forum.php?mod=viewthread&tid=36986

分享10个时间-高度剖面图及gs(温度平流、涡度平流、水汽通量散度、垂直速度......)：by：qiuba

bbs.06climate.com/forum.php?mod=viewthread&tid=94917

当初因为毕业论文而写的一些Grads脚本文件~  ：by：占戈示申

bbs.06climate.com/forum.php?mod=viewthread&tid=7986

做毕业论文的脚本，LP指数里面含有垂直风切变、水平风切变和温度平流的计算：by：扫眉才子

bbs.06climate.com/forum.php?mod=viewthread&tid=24184

分享一个理查森数计算的gs：by：平流层的萝卜

bbs.06climate.com/forum.php?mod=viewthread&tid=13173

Grads中斜温图的绘制(Skew-T)：by：传说中的谁

bbs.06climate.com/forum.php?mod=viewthread&tid=1104

关于露点温度计算的探讨(已知干球温度和相对湿度)：by：river

bbs.06climate.com/forum.php?mod=viewthread&tid=51560

分享自己写的P坐标涡度方程Grads脚本：by：桔梗花开

bbs.06climate.com/forum.php?mod=viewthread&tid=57579

GRADS 计算大气可降水量：by：无敌泡

bbs.06climate.com/forum.php?mod=viewthread&tid=29347

K指数脚本gs：by：gzdsg

bbs.06climate.com/forum.php?mod=viewthread&tid=90916

K指数计算方法：by：传说中的谁

bbs.06climate.com/forum.php?mod=viewthread&tid=1854

[求助]怎么把垂直速度单位由pa/s转为m/s

bbs.06climate.com/forum.php?mod=viewthread&tid=14293

[求助]关于大气热源Q1的计算：by：哈库拉玛塔塔000

bbs.06climate.com/forum.php?mod=viewthread&tid=43376

[求助]垂直风切变是如何计算的：

bbs.06climate.com/forum.php?mod=viewthread&tid=5608

本期内容到此结束！意犹未尽？

欢迎到气象家园 Grads板块，一起讨论吧！

http://bbs.06climate.com/

编辑：星星伴月