grads软件系列学习第3章基本操作命令

第3章基本操作命令

GrADS系统的运行方式，是采用命令行输入的方式交互式地显示图形。3.1基本命令

启动GrADS系统后，在ga->提示符后输入命令。在交互式环境内常用基本命令如下：

open、

set、

display、

clear……

下面逐一介绍。

open<路径>数据描述文件名

open命令用于打开GrADS的数据文件，启动GrADS后首先需要打开至少一个数据描述文件，命令如下：

openfilename主要用于打开数据描述文件或者数据文件,命令形式为ga->opendata.ctlga->sdfopenu.ncga->xdfopenwind.ddf关闭数据文件(close)ga->close2ga->close1注意要先关闭后打开的文件。set各类选项

set命令用于设置各种环境参数，包括维数环境、图形类型、图形要素、屏幕显示等。如设置维数环境(set)ga->setlon60160ga->setlat2050ga->settimeJul1983ga->setlev200lon、lat、time、lev为地图坐标,也可以换为格点坐标(x\y\z\t)表示。设定的维数环境决定了随后参与绘图或者计算的数据范围。

display(或者d)为显示命令

display命令是对表达式处理后进行屏幕图形显示。最简单的表达式是变量名的缩写。如变量名为v则显示命令:dv

clear(或c)

clear是清屏命令，清除图形窗口的内容。

quit

quit命令用于退出GrADS系统。

query(或q)选项

query是系统环境设置的查询命令。主要用于查询数据资料的时空结构、变量等信息。ga->qfile2

如querydefine：可知道定义了那些变量

dims：当前的维数环境

filen：查询第n号描述文件的内容

files：打开n个文件的次序

gxinfo：用在d之后，告诉用户图的一些信息

shades：用在d之后，告诉用户一些阴影的信息

time：时间设置信息定义新变量的命令(define/modify)

define命令用于定义新的变量，所定义的新变量可以用于随后的表达式中。ga->defineuave=aave(u,t=1,t=50)ga->modifyuaveseasonal(uave被定义为气候变量)初始化环境(Reinit/reset)注意Reinit和reset的区别,前者是关闭所有文件,重新初始化,后者则是重置初始环境,但不关闭数据文件,reinit常作为GrADS命令批处理文件中的第一个命令。基础绘图指令(draw)

draw命令均用在d之后,用于对图形进行标注字符串或者做标记、连线等run*.gs

run命令用于执行文件*.gs中定义的操作。

GrADS中几个重要gs文件(1)生成图形(.gmf)文件,并在(85°E,30°N)处标字符G′reinit′′opend:\llp\data.ctl′′enableprintd:\llp\h500.gmf′′dh′′qw2xy8530′x1=subwrd(result,3);y1=sbuwrd(result,6)′drawstring′x1′′y1′G′′print′′disableprint′draw命令都要用在display之后。第二种方法（printim命令）：形式：printim<路径>filenameoptionprintim命令在GrADSv1.8以上版本有效，可以在批处理文件中使用。例：输出1000x800的PNG图象：

printimout.pngx1000y800输出800x600白色背景的GIF图象：printimgifimage.outgifx800y600white(2)netCDF数据的读写并生成新的数据文件′reinit′′sdfopenh.nc′′opendata.ctl′′setgxoutfwrite′′setfwrited:\llp\h500.grd′*(设定数据时空范围)′settimejan1982dec1983′′setlon70100′′setlat2040′′setlev500′′dh′′disablefwrite′此例中,若资料为netCDF格式,就用sdfopen所在命令行;否则用open命令行。(3)绘制阴影图形,0线加粗′reinit′′opendata.ctl′′setgxoutshaded′′setcsmoothon′′dh′′setgxoutcontour′′dh′′setclevls0′′setcthick8′′dh′

当阴影图形和等值线图形在一起画,要先画阴影后画等直线,这一例实现了3幅图的叠加,阴影图、等值线图、单独加粗的0线图形。(4)熟悉Define命令和modify命令的使用技巧。特别注意,使用define命令前所设定的维数环境是新变量的维数环境,而不是原变量的维数环境,但是,新的GrADS版本中,define命令可以省略。′opendata.ctl′*(设定了气候变量hcli的维数环境)′sett112′′setlev500′*定义新变量hcli′′definehcli=ave(h,t+0,t=48,12)′′modifyhcliseasonal′*(设定了异常变量hano的维数环境)′sett148′′definehano=h-hcli′′settimejul1983′′dhano′(5)熟悉GrADS脚本语言中的循环语句“While.....endwhile”,下面一段程序可以连续显示多个时次500hPa高度场图形。′opendata.ctl′t=1while(t<=48)′c′′sett′t′′′dh′pulldummyt=t+1endwhile1.数据说明

nc.dat是二进制格点数据文件，为一组模式输出的全球5天的数值预报结果。有1个层次：1000，多个要素：h85,t85,u85…，时间从1948年1月开始，间隔1月，已知X方向间隔为144个单位，Y方向间隔为73个单位。

关于数据文件nc.dat的具体说明，