南京信息工程大学GrADS-复习资料.doc

目前有很多流行的绘图软件，如MICAPS、Grapher、Excel、GrADS、Surfer、NCL、MATLAB等，不同的软件有其不同的特点和优势。GrADS不仅能够绘制一维和二维图形、进行简单统计计算、物理量的计算，还可进行数据格式的转换。当前气象学研究中许多格点资料（如二进制，GRIB码，Netcdf，HDF-SDS资料）均需用GrADS软件进行处理，该软件因此也成为气象、海洋专业科研工作者和高校学生做研究论文使用最广泛的的绘图软件之一，并成为本科生和函授研究生必须选修的一门课程。 1.1 GrADS作图的一般步骤 GrADS作图一般包括四个步骤，即打开数据描述文件、输入GrADS内部命令、保存图形、退出GrADS。 1.2 几个重要的GrADS命令 打开文件的命令（open/sdfopen/xdfopen） 主要用于打开数据描述文件或者数据文件，命令形式为 ga-open data.ctl ga-sdfopen u.nc ga-xdfopen wind.ddf 2）查询命令（query或者q） 主要用于查询数据资料的时空结构、变量等信息。 ga-q file 2 ga-q ctlinfo 1 3）关闭数据文件的命令（close） ga-close 2 ga-close 1 注意要先关闭后打开的文件。 4）设置维数环境（set） ga-set lon 60 160 ga-set lat 20 50 ga-set time Jul1983 ga-set lev 200 上面的lon、lat、time、lev为地图坐标，也可以换为格点坐标（x\y\z\t）表示。 5）显示和清屏，初始化命令 clear(或者c)为清屏命令，display(或者d)为显示命令。 初始化环境（reinit/reset） 注意Reinit 和reset的区别，前者是关闭所有文件，重新初始化，后者则是重置初始环境，但不关闭数据文件，reinit常作为Grads命令批处理文件中的第一个命令。 6）图形类型的设置（set gxout） 等值线：contour；阴影：shaded；矢量图：vector，d u;v 流线图：stream d u；v；单线图：line；直方图：bar …. (ccolor，cstyle，cthick，clevs，ccols，cmin，cmax，cmark) 7）定义新变量的命令（define/modify） ga-define uave=aave(u,t=1,t=50) ga-modify uave seasonal(uave被定义为气候变量) 8）GrADS内定的变量 GrADS内定的变量有lon,lat,time,lev,result,rc。巧用这些内定变量可以让绘图变得简单。比如绘制赤道的命令如下， ga〉set clevs 0 gad lat Result是一个常用的重要的内定变量，GrADS命令的返回值、query（q）命令的查询结果均放在该变量中，比如 gad orog gaQ w2xy 60 3000 gaX = 2.35 y= 4.21 “X = 2.35 y = 4.21” 这一结果就放在result变量中，2.35、4.21分别为result变量中的第3个、第6个词。可内部函数subwrd从result变量中提取出这两个值，如 X1=subwrd(result,3) X2=sbuwrd(reselt,4) 然后可以进一步在该点确定的坐标处做标记。 9）基础绘图指令 draw命令均用在d之后，用于对图形进行标注或者做标记、连线等。 （draw string/line/rec/recf/mark/title/wxsym） 10）图形存储格式 gmf；png；gif； 11) 函数，变量名 ave/aave；mag（u，v）；skip； z.3 (lev=500)；tv(t-1,lev=300) 12）描述语言及编程 If/else/endif；if/endif；while/endwhile 2 GrADS绘图需要的三类文件 要使用GrADS进行绘图，需要准备三类型文件，即数据文件、数据描述文件(.ctl)、控制文件(.gs)文件。下面对三类型文件进行说明。 2.1数据文件 GrADS能够处理的数据文件常见格式有netCDF(自带描述)、HDF-SDS(卫星资料)、GRIB、Binary(用户自己生成)、站点资料等。NCEP/NCAR和ECWMF资料大多为netCDF格式，这种格式的数据资料自带描述信息，在GrADS中用sdfopen 命令打开。Binary是GrADS最基本的数据格式，常用文件后缀为.dat、.grd、.bin，需要和数据描述(.ctl)