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 内定，文件名随用户自定，
中国科学院大气物理研究所 大气科学和地球流体力学数值模拟国家重点实验室 National Key Laboratory of Numerical Modelling for Atmospheric Sciences and Geophysical Fluid Dynamics ( LASG ) Institute of Atmospheric Physics, Chinese Academy of Sciences
gradsnc
GrADS with readline, printim, LATS, read/write NetCDF
gradshdf GrADS with readline, printim, LATS, read/write HDF, read NetCDF
gribmap
"Maps" a GRIB data set to a GrADS descriptor file
§1.3 基本操作
新版本的 GrADS 交互环境中可以使用上下左右箭头键对曾经使用的命令进行调用和编 辑，但仅限于本次启动 GrADS 交互环境后所使用的命令。
GrADS画图入门
（缺省）
命令：
display express
ga-> d strp Contouring: 0 to 100 interval 10
设置：
set *** ###
维数参数设置 等值线设置
Set set Set set lon 100 cint 15 120 指定等值线间隔 lat 103 25 指定等值线颜色 ccolor
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
• 3、图形显示类型及绘图设置：
•
2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5_to_Grads,wrf_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
data open(3,file=filename.ctl,form='formatted',status='unknown') zlev/1000,925,850,700,500,400,300,250,200,150,100/ write(3,'(a)')'dset ^'//filename.dat tim=0.0 ；ilev=11；nflag=0 open(2,file=filename.dat,form='unformatted',access='direct',recl=73*4) write(3,'(a)')'dtype station' irecord=0 write(3,'(a)')'stnmap ^sounding.map' do j=1,irec write(3,'(a)')'undef -999.0' irecord=irecord+1 write(3,'(a)')'title sounding' write(2,rec=irecord)'abc'//sta(j)(1:5),lat(j),lon(j),tim,ilev,nflag, write(3,'(a)')'tdef 1 linear 02JAN1987 1DY ' + (zlev(i),z(j,i),t(j,i),td(j,i),d(j,i),f(j,i),i=1,11) write(3,'(a)')'vars 5' enddo write(3,'(a)')'zz 11 99 '//staa//' height1' ilev=0 write(3,'(a)')'tt 11 99 '//staa//' temper1' irecord=irecord+1 write(3,'(a)')'t_td 11 99 '//staa//' t_td1' write(2,rec=irecord)'abcdefgh',lat(1),lon(1),tim,ilev,nflag, write(3,'(a)')'uu 11 99 '//staa//' u wind1' + write(3,'(a)')'vv (zlev(i),i=1,11),(z(1,i),i=1,11),(t(1,i),i=1,11), 11 99 '//staa//' v wind1' + write(3,'(a)')'endvars‘ (td(1,i),i=1,11),(d(1,i),i=1,11),(f(1,i),i=1,11) close(2) close(3)
grads-第八章
秦育婧
南京信息工程大学 大气科学学院
第八章 函数
学习目标 学会使用常见函数
第八章 函数
学习要求 掌握ave函数
函数的调用
1）通过函数名直接引用， 2）参数放在括号中用逗号分开， 3）可以嵌套调用， 4）有些函数在运算时会改变维数
环境。
常用函数
（1）ave函数 格式： ave (expr,dexpr1,dexpr2<,tincr<,flags>>) 功能：通用的求平均函数。 说明：expr是由dexpr1和dexpr2定义的维数 范围内t（expr,constant<,flag>） 功能：设置部分网格点的值取为常数
constant。 说明：所有非缺测格点处的expr值取为常
数，flag为选项，如果加上选项-a，则所有 网格点值均设定为指定的常数，如果加上选 项-u，则只把缺测格点处的expr值设定为常 数。该函数对格点和台站资料均适用。
示V风速分量，风速单位用m/s。 边界上的涡度值设定为缺测。 例如：d hcurl(u,v)
常用函数
（5）hdivg函数
格式：hdivg（uexpr,vexpr） 功能：计算水平散度。 说明：uexpr表示U风速分量，vexpr表
示V风速分量，风速单位用m/s。
常用函数
（6）skip函数
格式：skip（expr,skipx,skipy） 功能：设定样本的取样密度。 说明：skipx,skipy 数值决定X和Y方向的 取样密度（取值1可以省略不给） *该函数主要用于对矢量场的稀疏化显示。
常用函数
例： ga->open model.ctl
set lev 500 d ave(z,t=1,t=5) （显示500hPa
气象绘图软件GrADS高级使用技巧
收稿日期:2002-11-1作者简介:高文良,男,31岁,主要从事短期气候预测等研究工作。
气象绘图软件GrADS 高级使用技巧高文良1 刘晓燕2 曾小东3(11成都高原气象研究所 成都 610072; 21阿坝州金川县气象局 62410031阿坝州马尔康县气象局 624000)摘 要:本文通过介绍气象绘图软件GrADS 使用中的一些高级技巧,分析了GrADS Script 语言中的难理解之处,并介绍了在GrADS 中做合成分析和t 检验的程序,可以对学习和使用GrADS 的科研人员起参考、帮助作用,达到事半功倍的效果,促进GrADS 软件的应用更广泛、深入,提高科研人员的工作效率。
关键词:GrADS ;技巧;程序中图分类号:TP391文献标识码:C文章编号:1003-7187(2002)04-0057-031 引言气象绘图软件GrADS (Grid Analysis and Dis 2play System )是免费共享软件,可随时从互联网上下载(http :///grads/),后续版本正陆续推出[1]。
GrADS 有丰富的内部函数,可以对数据进行计算和分析处理。
它支持处理格点资料和站点资料,并且支持对GriB 码文件、特殊格式文件(如一字节整型、二字节整型、大中型机器二进制数据等)的直接读取,气象科研领域应用非常广泛[2]。
在其最新1.8SL 9版本中,GrADS 又将应用领域推进到了海洋学科,功能也得到了进一步地增强和扩展。
但在使用GrADS 的过程中,特别是使用GrADS Script 语言编程当中,有一些问题需要特别注意。
因为GrADS Script 语言是一种类似于VB Script (或MA TL AB Script )的高级语言,稍不注意就容易出现错误,且查错十分困难,这点与其他语言有较大的差别。
2 站点数据处理GrADS 中站点数据处理基于格点数据的基础之上,需先将站点数据通过Cressman 客观分析方法内插至格点上,然后再依照格点资料的处理方式对站点资料数据进行分析和处理。
Grads 使用笔记
****************************************************************
!!!!重要命令
'set tlsupp year|month'
控制时间坐标轴的标记
例如'set tlsupp year'就会去掉年的标记只留下月份。
t5=t5+1
endwhile
'reinit'
这样可以提取你想要的年数据，然后你大可运用fortran对数据进行随心所欲的处理。
***************************************************************
'set timelab on|off'
*****************************************************************
而不能仅仅'd u;v'
之————控制矢量密度与颜色
'set ccolor 1'*控制颜色
'd skip(u,2);v'*控制密度
*******************************************************************************************
'set lat 15 55'
'd hgt'
t1=t1+1
endwhile
to=to+1
endwhile
*这里必须先观点上述运行的文件，grads最多同时可以打开20个文件左右。
第十章GradsGrADS 描述语言及编程
GRADS (Grid Analysis and Display System)讲义（10）GrADS 描述语言及编程一、描述语言概述描述语言是GrADS软件包中自带的高级语言，可以用于GrADS的高级操作，更方便地使用GrADS的功能。
通过描述语言可以编写具有复杂功能的命令集（即gs文件）。
二、gs文件中的两类型语句1）描述语句：不用引号，需要经过编译的。
2）指令性语句：GrADS命令，需用引号括起来。
通过描述语句来控制指令描述语言用到指令性语句中要用引号，即，i=1‘set t ‘i’’(或者‘set t ’ i )三、描述语言程序特点1、程序通过GrADS的run命令来编译执行；run file-name options2、可以有变量、流程控制也可以执行GrADS 的命令。
3、描述语言可以1）编写函数子程序，2）制作屏幕菜单，3）可以显示动画图形，4）可以使用从GrADS的query命令中得到的信息作注释等。
3、描述语言的所有变量都是字符串型的，描述变量可以进行计算。
流程控制是通过if/else/endif和while/endwhile语句块实现。
循环过程可以用continue或break命令改变。
变量中包含的或者由表达式结果产生的字符串可以作为GrADS命令项使用。
而由GrADS命令执行后产生的文字结果可以读入变量并在描述语言程序中使用。
描述语言也支持函数子程序功能。
四、描述语言的组成1、一个描述语言程序是由许多语句组成,每个语句之间由换行号或分号分隔。
2、下面是GrADS描述语言中语句的类型：赋值语句If / Else / Endifwhile / endwhilebreak / continuefunction header / returnsay / pull3、描述表达式描述表达式是由运算对象和运算符组成。
运算对象：描述变量、子程序函数或常数；运算符：数学、逻辑或连结运算符。
下面分别介绍变量、运算符、表达式4、描述语言变量1）构成变量名是由1至8个字符组成，以字母开头，仅可以包含字母和数字，变量名区分大小写。
GrADS命令大全
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
第6章 GrADS变量和函数课件
• 当定义的变量具有某些变化的维数环境，然后再 将这些维数固定起来，这时变量将显示固定在该 维数上的值。
• 用户也可以使用局地维数环境。
第6章 GrADS变量和函数
例如：
d temp(lat=50) 以上显示命令中，变量temp在定义时Y维
数是变化的，因此lat=50的局地维数环境的设 置有效，图形将显示变量temp在500N上的资 料。又如： d temp(t=4)
第6章 GrADS变量和函数
• 例如： …… set lat 40 set lon -120 set lev 500 set t 1 5 d tloop(aave(t,lon=-140,lon=-120,lat=20,lat=40))
注意：要显示图形时，必须使得X，Y，Z三维固定，只有时间t变化。
第6章 GrADS变量和函数
2. 常用函数
（1）ave函数 • 格式： ave（expr,dexpr1,dexpr2<用的求平均函数。 • 说明：expr是由dexpr1和dexpr2定义的维数
范围内求平均的量。
第6章 GrADS变量和函数
例如： ga->open model.ctl
括号内的dimexpr：是对当前维数 环境进行的局域维数设置表达式，该 表达式仅对该变量的维数环境进行局 域修正，因而仅对固定的维有效。
第6章 GrADS变量和函数
设置形式可以采用以下两种方式：
绝对维数表达式为： X|Y|Z|T|LON|LAT|LEV|TIME = value （value为绝 对维数值）；
小时时刻的1号文件中的变量tv）
rh
（缺省的当前文件中的变量rh）
q.2（t-1,lev=850） （2号文件中相对于当前
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
GrADS的函数分两类,
一类是对格点/站点数据执行运算的,这一类我们姑且称之为分析函数;
另一类是脚本编程(gs)时使用的,这后一类我们姑且称之为脚本函数.
第一类分析函数又分为格点分析和站点分析两个子类:
1)格点分析
绝对值abs(expr)
余弦/反余弦cos(expr)acos(expr)
正弦/反正弦sin(expr)asin(expr)
正切/反正切tan(expr)atan2(expr1,expr2)[expr1/expr2]
指数exp(expr)
对数log(expr)log10(expr)
幂函数pow(expr1,expr2)expr1expr2
直角三角形的弦长/全风速mag(u_expr,v_expr)sqrt(u*u+v*v)
平方根sqrt(expr)
加权/不加权平均ave(expr,dim_expr1,dim_expr2,tincr,{|-b})mean()
加权/不加权区域平均aave(expr,{xdim1,xdim2,ydim1,ydim2|global|g}) amean()
剔除值后时间求平均tmave(undefine_or_weight_expr,expr,texpr1,texpr2)
加权/不加权区域和asum(expr,{xdim1,xdim2,ydim1,ydim2|global|g}) asumg()
质量加权垂直积分vint(ps_expr,expr,top)
普通积分gint(expr)
最大/小值Max(expr,dim1,dim2,tincr)min()
最大/小值所在格点Maxloc(expr,dim1,dim2,tincr)minloc()
空间相关/回归scorr(expr1,expr2,{xdim1,xdim2,ydim1,ydim2|global| g})sregr()
时间相关/回归tcorr(expr1,expr2,tdim1,tdim2)tregr()
常数const(expr,const,{|-u|-a})
剔除maskout(expr,mask_expr)
跳点skip(expr,skip_x,skip_y)
加权9点平滑smth9(expr)1-D3点平滑
中央差分cdiff(expr,dim)
时间循环tloop(expr)
极值位置maxloc(expr,dim1,dim2,tincr)minloc()
特征值位置findlvl(expr,expr_to_find,lev1,lev2)
气象上物理量计算的函数:
虚温,相对湿度求比湿tvrh2q(tv_expr,rh_expr)g/g
虚温,相对湿度求温度tvrh2t(tv_expr,rh_expr)degK
垂直涡度hcurl(u_expr,v_expr)
中央差求水平散度hdivg(expr1,expr2)
全风速mag(u_expr,v_expr)
2)站点分析
从格点插到站点gr2stn(grid_expr,{stn_expr|lon,lat},-a)
收集生成格点场coll2gr(cnum,{num|-u})
Cressman插值oacres(grid_expr,stn_expr,radii,1st_guess)
网格平均oabin(grid_expr,stn_expr,{-f|-c})
时间平均stnave(expr,dexpr1,dexpr2,-m count)
时间序列最小值stnmin(expr,dexpr1,dexpr2,-m count)
时间序列最大值stnmax(expr,dexpr1,dexpr2,-m count)
转成1-D格点序列s2g1d(expr)
所有数学函数abs,cos,acos,sin,asin,tan,atan2,exp,log,log10,pow,mag,sqrt
第二类,脚本函数包括有:
abs,log,log10:math_abs,math_log,math_log10
sin,cos,tan:math_sin,math_cos,math_tan
asin,acos,atan2:math_asin,math_acos,math_atan
math_sinh,math_cosh,math_tanh
math_asinh,math_acosh,math_atanh
pow,sqrt:math_pow,math_sqrt,math_exp
取整:math_nint,math_int
求余:math_fmod,math_mod
字符串处理函数:math_strlen,valnum,wrdpos,strlen
sublin,subwrd,substr
数值格式转换函数:math_format
文件操作函数:read,write,close。