GrADS命令大全

合集下载

【grads】ncdump的使用总结

【grads】ncdump的使⽤总结1、参数ncdumpncdump [-c] [-h] [-v var1,...] [-b lang] [-f lang] [-l len] [-n name] [-d f_digits[,d_digits]] fileWhere:-cShow the values of coordinate variables (variables that are also dimensions) as well as the declarations of all dimensions, variables, and attribute values. Data values of non-coordinate variables are not included in the output. This is the most suitable option to use for a brief look at the structure and contents of a netCDF file.-hShow only the header information in the output, that is the declarations of dimensions, variables, and attributes but no data values for any variables. The output is identical to using the -c option except that the values of coordinate variables are not included. (At most one of -c or -h options may be present.)-v var1,...,varnThe output will include data values for the specified variables, in addition to the declarations of all dimensions, variables, and attributes. One or more variables must be specified by name in the comma-delimited list following this option. The list must be a single argument to the command, hence cannot contain blanks or other white space characters. The named variables must be valid netCDF variables in the input-file. The default, without this option and in the absence of the -c or -h options, is to include data values for all variables in the output.-b langA brief annotation in the form of a CDL comment (text beginning with the characters ``//'') will be included in the data section of the output for each `row' of data, to help identify data values for multidimensional variables. If lang begins with C or c, then C language conventions will be used (zero-based indices, last dimension varying fastest). If lang begins with F or f, then Fortran language conventions will be used (one-based indices, first dimension varying fastest). In either case, the data will be presented in the same order; only the annotations will differ. This option is useful for browsing through large volumes of multidimensional data.-f langFull annotations in the form of trailing CDL comments (text beginning with the characters ``//'') for every data value (except individual characters in character arrays) will be included in the data section. If lang begins with C or c, then C language conventions will be used (zero-based indices, last dimension varying fastest). If lang begins with F or f, then Fortran language conventions will be used (one-based indices, first dimension varying fastest). In either case, the data will be presented in the same order; only the annotations will differ. This option may be useful for piping data into other filters, since each data value appears on a separate line, fully identified.-l lenChanges the default maximum line length (80) used in formatting lists of non-character data values.-n nameCDL requires a name for a netCDF data set, for use by ncgen -b in generating a default netCDF file name. By default, ncdump constructs this name from the last component of the pathname of the input netCDF file by stripping off any extension it has. Use the -n option to specify a different name. Although the output file name used by ncgen -b can be specified, it may be wise to have ncdump change the default name to avoid inadvertantly overwriting a valuable netCDF file when using ncdump, editing the resulting CDL file, and using ncgen -b to generate a new netCDF file from the edited CDL file.-d float_digits[,double_digits]Specifies default number of significant digits to use in displaying floating-point or double precision data values for variables that don't have a `C_format' attribute. Floating-point data will be displayed with float_digits significant digits. If double_digits is also specified, double-precision values will be displayed with that many significant digits. If a variable has a `C_format' attribute, that overrides any specified floating-point default. In the absence of any -d specifications, floating-point and double-precision data are displayed with 7 and 15 significant digits respectively. CDL files can be made smaller if less precision is required. If both floating-point and double-presision precisions are specified, the two values must appear separated by a comma (no blanks) as a single argument to the command. If you really want every last bit of precision from the netCDF file represented in the CDL file for all possible floating- point values, you will have to specify this with -d 9,17.Usage Notesncdump generates an ASCII representation of a specified netCDF file on standard output. The ASCII representation is in a form called CDL (``network Common Data form Language'') that can be viewed, edited, or serve as input to ncgen. ncgen is a companion program that can generate a binary netCDF file from a CDL file. Hence ncgen and ncdump can be used as inverses to transform the data representation between binary and ASCII representations. See ncgen for a description of CDL and netCDF representations.ncdump defines a default format used for each type of netCDF data, but this can be changed if a `C_format' attribute is defined for a netCDF variable. In this case, ncdump will use the `C_format' attribute to format each value. For example, if floating-point data for the netCDF variable Z is known to be accurate to only three significant digits, it would be appropriate touse the variable attributeZ:C_format = "%.3g"ncdump may also be used as a simple browser for netCDF data files, to display the dimension names and sizes; variable names, types, and shapes; attribute names and values; and optionally, the values of data for all variables or selected variables in a netCDF file.ExamplesLook at the structure of the data in the netCDF file foo.nc:ncdump -c foo.ncProduce an annotated CDL version of the structure and data in the netCDF file foo.nc, using C-style indexing for the annotations:ncdump -b c foo.nc > foo.cdlOutput data for only the variables uwind and vwind from the netCDF file foo.nc, and show the floating-point data with only three significant digits of precision:ncdump -v uwind,vwind -d 3 foo.ncProduce a fully-annotated (one data value per line) listing of the data for the variable omega, using Fortran conventions for indices, and changing the netCDF dataset name in the resulting CDL file to omega:ncdump -v omega -f fortran -n omega foo.nc > Z.cdl2、注意①它是在dos底下⽤的，如果在gradsnc中⽤，请在命令前加感叹号。

GrADS命令大全

set rbcols c1 c2<c3…cn>
设置彩虹色颜色序列值（缺省：9，14，4，11，5，13，3，10，7，12，8，2，6）
set rbcols <auto>
打开系统缺省设置的彩虹色设置
set rbrange low high
彩虹色颜色对应等值线开始和结结束的范围
set grid on|off|<style>
清干净图形显示区并重新初始化，但1）不关闭任何文件，2）不取消变量定义，3）不改变set display设置
option＝events; graphics; hbuff; norset
reinit
与reset相似，但关闭所有文件并取消所有定义
quit
退出
!shell-command
在GrADS命令交互模式下运行外部系统命令
用三次样条或线性插值的方式计算细网格的值
set cint value
设置固定等直线间隔值
set cmax value
设置最大等直线的值，超过的不再画出
set cmin value
设置最小等直线的值，小于的不再画出
set clevs lev1 lev2…
按lev1 lev2…指定值绘等直线
set ccols col1 col2…
按顺序执行fname文件中的命令组
clear <option>
c <option>
清干净图形显示区，但不关闭当前数据文件可带参数：命令可简写为“c”并可带以下参数：
events
清空事件缓冲区
graphics
清空图形缓冲区，但不清空组件缓冲区
hbuff
在双缓冲区模式下，清空显示缓冲区

grads中文说明手册

维数环境是 GrADS 的一个重要概念，GrADS 视每一个物理变量场为一个四维数据集（4D data set）。包括空间三维和时间一维，也可固定其中的一维或几维以获得实际的低于四维的数据子集。GrADS 中设置维数用以说明或指定随后的分析或图形操作时参加操作的原始数据集的维数范围，即通过设定工作数据的起止点数、取点频度（间隔）来设定最后工作数据场的数组成份。该工作数据集可以是整个原始数据场，也可以是原始数据场的一部分，详见§2，§6。所谓维数环境是对格点数据而言的。维数环境的定义可在两种空间坐标上进行。一种是地球坐标（world coordinate），以经纬度为度量单位；另一种是格点坐标（grid coordinate），以网格点数为度量单位。二者一般对应于同一个绘图坐标，都对应到网格点上。在维数环境表达式中 x，y，z，t 与 lon，lat，lev，time 是分别对应于两套坐标的专用维数变量，含义固定，如 x 与 lon 都指西到东指向的（缺省方向）水平坐标，y 与 lat 都指南到北指向的（缺省方向）水平坐标，z 与 lev 都指从地面到高空的（缺省方向）垂直坐标，t 与 time 都是时序坐标，不过 t 用的是格点时次序号，而 time 用的是格林威治标准时的 GrADS 绝对表达格式。详见§2，§4。
GrADS 当前版本为 1.8s11，可用于如下平台：DEC、Intel/LINUX、SUN、Macintosh OSX、 SGI、IBM/AIX、MS Windows。
GrADS 软件包所包含的文件：
gradsc
GrADS Classic - readline, printim only
gradsdods GrADS with readline, printim, LATS, read/write NetCDF, DODS

IDL、GRAD、NCL绘图命令对应关系一览表

IDL、NCL、GRADS、MATLAB绘图命令对应关系一览表1.grads的数据文件与ncl的什么文件对应？grads只支持按照一定顺序存储的二进制数据文件，后缀名以.grd或者.dat或者.grb结束。

这种数据必须以时间为最外层，然后每个变量按照向量形式存储，每个变量由外向内的存储顺序是高度（或等压面）==》纬度==》经度。

这必须注意，否则画图容易出现一堆一堆的乱线条，这就说明你的数据没有按照grads的要求存储。

再看看ncl，可以说ncl支持绝大多数各种数据的读写，包括netcdf，hdf,以及二进制数据甚至ASCII码(如果说是十进制数据或许你会更熟悉)，前两种数据一般都有头文件，不能用C 语言或者Fortran读取，都需要插件才可以读取，ncl可以直接读取，matlab中也可以直接读取NETCDF格式(.nc)的数据。

grads中可以读取按照说明存取的NC数据，这种数据必须又正确的时间说明，也就是说时间必须是真实的，有些模式模拟出来的数据grads的sdfopen 命令是打不开的，因为一般模式都是nonleap run，都是平年，没有设定闰年，造成了时间说明不真实，grads就会报错。

那么grads'如何使用NC数据呢？所以建议使用ncl转换数据，将NC数据，hdf数据或者十进制数据转换成grd数据，供grads使用。

这样说明是在是太空洞了，那么下面我举个例子吧。

eg1) 使用grads将netcdf数据转换成grd数据'reinit'var.1='air';var.2='hgt';var.3='uwnd'var.4='vwnd';var.5='omega';var.6='shum'j=6while (j<=6)'set fwrite/disk3/users/Rao_Jian/ERA-Interim-daily/entropy/'var.j'.daily.1979-2010.grd' 'set gxout fwrite'i=1979while (i<=2010)'sdfopen/disk3/users/lbq/ERA-Interim-daily/pressure/'var.j'.interim.'i'.nc'tt=1if(i=1980|i=1984|i=1988|i=1992|i=1996|i=2000|i=2004|i=2008)while (tt<=366)'set t 'ttzz=1while (zz<=37)'set x 1 240''set y 1 121''set z 'zz'd 'var.j''zz=zz+1endwhilett=tt+1endwhileelsewhile (tt<=365)'set t 'ttzz=1zz=1while (zz<=37)'set x 1 240''set y 1 121''set z 'zz'd 'var.j''zz=zz+1endwhilett=tt+1endwhileendifi=i+1name='/disk3/users/lbq/ERA-Interim-daily/pressure/'var.j'.interim.'i'.nc''close 1'endwhile'disable fwrite'j=j+1endwhileeg2)使用ncl将数据输出成二进制数据，grads可以使用(截取部分)patho = "/disk3/users/Rao_Jian/Hadley/"system("rm "+patho+"ev_ts.grd")system("rm "+patho+"ev_ts2.grd")system("rm "+patho+"ev_ts.txt")system("rm"+patho+"ev_ts2.txt")do nt=0,719fbindirwrite(patho+"ev_ts.grd",ev_ts(:,nt));写成二进制数据，供grads使用end dofbindirwrite(patho+"ev_ts2.grd",ev_ts(time|:,evn|:))asciiwrite(patho+"ev_ts.txt",ev_ts(time|:,evn|:))；写成十进制数据，可以贴到EXCEL中使用asciiwrite(patho+"ev_ts2.txt",ev_ts);print(ev_ts(0,:))此外ncl中还有其他的读写函数，如fbinrecread,fbinrecwrite,fbinwrite,fbinreadncl读取netcdf3/4、hdf4、grib-1/2也是通过文件操作的eg3)ncl支持直接读取格式文件的读法fi = addfile(pathi+"HadISST_sst.nc","r")sst0 =fi->sst(960:1679,:,:) ; load 50 year data duringfrom 1950 to 2009注意：这类文件的后缀名一般为.nc /.hdf/ .grb/.hdfeos/.ccm2.grads中的描述文件与ncl中的什么对应描述文件（.ctl)是一个纯文本文件，我们有数据grads还是不能出图，需要一个描述文件来指定他的存储数据个数，维度（时间长度、层数和经纬度信息）。

(完整版)GrADS绘图软件使用手册

GrADS绘图软件实用手册2002年1月目录第一章GrADS绘图软件概述1.GrADS绘图软件简介2.Internet上的GrADS资源2.1GrADS在Internet上的主页2.2windows环境下GrADS资源3.GrADS绘图软件的安装（windows环境）3.1在windows环境下安装GrADS软件包X server 的安装第二章GrADS绘图模板1.GrADS示例演示1.1 启动GrADS1.2 退出GrADS1.3 示例演示GrADS命令的使用2.GrADS绘图模板3.GrADS模板的高级应用GrADS描述语言GrADS高级模板的应用第三章GrADS数据格式1.格点数据描述文件1.1 数据描述文件各项解释1.2 生成model.le.dat和model.le.ctl文件的程序代码片段2.站点数据的格式附录1.如何精确控制图形输出的尺寸—Landscape纸型2.台站资料的显示3.Linux环境下的安装第二章GrADS绘图软件概述1GrADS绘图软件简介The Grid Analysis and Display System(GrADS) 是一套应用广泛、使用方便的科学数据绘图软件包。

其主要特点：●GrADS属于自由软件，可以从Internet上免费获得。

●可运行于各种Windows 和Unix工作平台。

●GrADS可用于4D数据的分析。

既经度、纬度、层（气压层、高度层等）和时间/xyzt 4维。

数据可以是格点化的数据或离散点数据。

GrADS特别适用于气象类数据的分析。

但也完全可以用于更广泛类型的数据分析。

●GrADS有多种显示方式：等值线、流线、矢量图、风矢量图、站点填图、折线图、直方图等多种两维图形。

●可处理多种数据格式的数据。

GRIB、NetCDF、HDF-SDS等通用数据格式和系统自定义的一种二进制数据格式。

●采用命令行输入的方式交互式地显示图形。

并有多种命令对数据进行再加工。

grads常用命令

第一个命令：这样将得到一个动画显示的序列，现在输入：claerset lon -90set lat -90 90set lev 1000 100set t 1d td u该情况下我们设置y（纬度）和z（高度）变化，所以我们得到一个竖直切面图。

我们显示了两个变量，他们重叠在一起。

你可以让任意多的变量重叠在一起，只要你不使用clear命令就行了。

另一个例子，这种情况X和T是变化的（哈莫图）：cset lon -180 0set lat 40set lev 500set t 1 5d z现在已经知道怎样使用选定的部分数据绘图了，下面我们学习怎样操作数据。

设置维数环境为Z,Y变化：clearset lon -180 0set lat 0 90set lev 500set 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方向作平均并求距平：cleard z-ave(z,x=1,x=72)也可以做时间差分：cleard z(t=2)-z(t=1)完整规范的变量名是：name.file(dim |-|=va lue,…) 如果我们打开了两个文件，也许一个是模式输出，另一个是分析，我们应该区分用如下方法二者：display z.2-z.1另一个内置的函数通过有线差分计算水平涡度相关cleard hcurl(u,v)还有另外一个计算数值方向的质量积分：cleard vint(ps,q,275)这儿我们计算了可降水量（单位mm）现在来讨论控制图形输出的话题。

grads基本命令

grads基本命令d ave(olr,t+0,t+359,12) ave表示求平均，olr 为变量， t+0表示起始点为当前时次，即00z01jul1974；t+359表示终止时次为当前时次后推359个时次，即00z01jul2003；12表示每12个时次取一次数据；打开nc文件：ga->sdfopen查询信息：ga->q filega->q dimga->q ctlinfo写标题：ga-> draw title清除图面：ga-> c同时打开两个文件： ga->open model.ctlga->open model.le.ctl删除所有设置重新回到刚进入GrADS状态：ga->reinit删除open命令后的所有设置：ga->reset由风场导出涡度场：ga->d hcurl(u,v)以分色图形方式输出：ga->set gxout shaded以等值线方式输出（缺省方式）：ga->set gxout contour画图例：ga->cbarn 1 0 （1:相对长短，>1放大；0:水平；1:垂直。

）以箭头方式表示矢量场: ga->gxout vector风矢量场：ga->d u; v; q （显示矢量时，d x分量；y分量<；标量>。

“<>”内的部分只起标颜色的作用。

）以流线方式表示矢量场：ga->set gxout stream风流线场：ga->d u; v; q以WMO风标方式表示矢量场： ga->set gxout barb风标：ga->d u; v; q直接输出网格点数值：ga->set gxout grid设置保留小数位数：ga->set dignum 0设置数字大小：ga->set digsize 0.1ga->set mpdraw on 如为off，不画地图背景（非经纬度数据需此项）ga->set poli on 如为off不画国界省界等。

grads基本命令

d ave(olr,t+0,t+359,12) ave表示求平均，olr 为变量， t+0表示起始点为当前时次，即00z01jul1974；t+359表示终止时次为当前时次后推359个时次，即00z01jul2003；12表示每12个时次取一次数据；打开nc文件：ga->sdfopen查询信息：ga->q filega->q dimga->q ctlinfo写标题：ga-> draw title清除图面：ga-> c同时打开两个文件： ga->open model.ctlga->open model.le.ctl删除所有设置重新回到刚进入GrADS状态：ga->reinit删除open命令后的所有设置：ga->reset由风场导出涡度场：ga->d hcurl(u,v)以分色图形方式输出：ga->set gxout shaded以等值线方式输出（缺省方式）：ga->set gxout contour画图例：ga->cbarn 1 0 （1:相对长短，>1放大；0:水平；1:垂直。

）以箭头方式表示矢量场: ga->gxout vector风矢量场：ga->d u; v; q （显示矢量时，d x分量；y分量<；标量>。

“<>”内的部分只起标颜色的作用。

第3章 grads基本操作命令

关于数据文件model.le.dat的具的具关于数据文件体说明，可以打开model.le.ctl文件体说明，可以打开文件查看。查看。
2. 示例演示
启动GrADS后，选择图形输出窗口后启动的方式（或），然后在命令窗口ga-> ），然后在命令窗口的方式（L或P），然后在命令窗口提示符后输入命令，提示符后输入命令，通过图形输出窗口可以观察输出结果。可以观察输出结果。
open <路径数据描述文件名路径>数据描述文件名路径
open命令用于打开命令用于打开GrADS的数命令用于打开的数据文件，启动GrADS后首先需要打开据文件，启动后首先需要打开至少一个数据描述文件命令如下：数据描述文件，至少一个数据描述文件，命令如下： open filename
define 临时变量名表达式临时变量名=表达式
define命令用于定义新的变命令用于定义新的变量，所定义的新变量可以用于随后的表达式中。后的表达式中。
modify 临时变量名时间序列的类型<seasonal diurnal> 类型 modify命令可将自定义的变命令可将自定义的变量声明为气候值，量声明为气候值，用于后面的时次代换。次代换。
（以下set命令是设置维数环境，关于维数环境的介绍详见第4章）
（显示该变量）（清屏）
set lon -180 0 dz c set lat 0 90 dz quit
（重新设置维数环境）
（退出GrADS系统）
实际操作后，实际操作后，在图形输出窗口观察输出结果。口观察输出结果。比较不同的维数环境设置下，数环境设置下，输出结果有何不关于set命令的应用将在第命令的应用将在第4 同？关于命令的应用将在第章具体介绍。章具体介绍。

grads常用命令

我们显示了两个变量，他们重叠在一起。

你可以让任意多的变量重叠在一起，只要你不使用clear命令就行了。

设置维数环境为Z,Y变化：clearset lon -180 0set lat 0 90set lev 500set 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天的平均。

grads常用GS程序打包

'set clab forced'
'set clab %gm/s'
*等值线标记数值的最后标上m/s*
*'set cint 1'
*等值线显示间隔为2*
'set clskip 2'
*等值线标记间隔为2*
'set cterp off'
*让添色图和等值线图的边缘重合*
*'d u'
*0线加粗*
'q w2xy jan1982 0'
x1=subwrd(result,3)
y1=subwrd(result,6)
'q w2xy dec1985 0'
x2=subwrd(result,3)
y2=subwrd(result,6)
'draw line 'x1' 'y1' 'x2' 'y2''
'set ylopts 1 1 0.15'
*设置X、Y的颜色(1)、线宽(4)、字符大小(0.12)*
*'set xyrev on'
'set z 1'
'set x 1'
'set y 1'
'set t 1 48'
*设置仅取时间为变量*
'set cmark 3'
*设置线条上的标记符号（3——实心圆）*
'd a'
'print'
pull dummy

GrADS常用指令

1.ctl文件的编写dset d:\test\test1.dat 数据文件路径名和文件名title monthly precipitation data 数据标题options 365_day_calendar 特殊格式说明365_day_calendar/template(多个文件) undef -9999.0 缺测值，若有多个则用fortran改写成同一个缺测值xdef 144 linear 0 2.5 x维数格点数/排列方式/起始值/间隔*linear线性/levels列举ydef 73 linear -90 2.5 y维数zdef 5 linear 1000 850 700 500 200 z维数tdef 24 linear 00z01Jan1979 1mo 时间维数*mo月yr年vars 1 变量数precip 1 99 *precipitation data 变量名/层数/排列顺序/说明endvars*变量循环顺序：x经度y维度z高度v变量t时间*注释行第一列用‘*’，注释行不能出现在变量列表中2.维数环境设置set lon/lat/lev/time var1 <var2> *实际值set x/y/z/t var1 <var2> *格点数值*两种坐标可以混用3.图形文件的保存gxprint/printim **.png/jpg/pdf/eps…. whiteenable print **.gmfprintdisableprint4.图形类型的设置set gxout typecontour 二位等值线图*默认shaded 二维填色等值线图vector 矢量箭头二维风场图*d u;v时为默认bar 直方图line 折线图*默认…………………5.图形要素设置gxout=contour or shadedset cint value 设置等值线间隔set clevs lev1,lev2,……设置特定等值线set ccols col1,col2,……设置对应于set clevs的等值线的颜色cbarn size 0/1 x y 设置色标gxout=line or contourset ccolor color 设置线条颜色set cstyle style设置线条样式1实线2长虚线3短虚线……set cthick thickns 设置线条粗细1~10set cmark marker 设置数据点标记0无标记*line时有效gxout=barset bargap val 设置直方条间隔*0~100百分比，0为无间隔set barbase val/bottom/top 设置直方条起始值gxout=vectorset arrscl size <magnitude> 设置箭头长度set arrowhead size 设置箭头大小，缺省为0.056.gs文件的编写将命令用单引号括起来7.变量define varname=expression 定义变量‘d air(t=3)’‘d air(t=3,z=5)’‘d air(lev=1000)-air(lev=500)’设置变量维数8.平均函数ave(expr,dimexpr1,dimexpr2,<timeint>) 一维平均aave(expr,dimexpr1,dimexpr2, dimexpr3,dimexpr4,) 二维平均*不能对z、t *需将维度设置成点9.地图投影设置set mproj proj 设置地图投影方式latlon 缺省nps 北半球极地投影sps 南半球极地投影set map color style thickness 设置背景地图线条颜色、线型、宽度10.绘图区域设置set parea xmin xmax ymin ymax11.坐标要素控制set frame on/off/circle 绘制图形边框/不绘制/绘制圆形边框set grads on/off 标记grads及绘图时间draw title <string> 绘制标题draw xlab/ylab <string> 绘制x轴/y轴图注set xyrev on/off xy轴互换set zlog on/off z方向用对数坐标set xlint/ylint value 设置x/y轴坐标刻度间隔set vrange val1 val2 设置y轴范围set vrange2 val1 val2 设置x轴范围set xlopts/ylopts color thickness size x/y轴刻度值的颜色、线宽、大小*缺省1,4,0.12 set tlsupp year 去掉年份显示12.颜色的定义‘define_colors’使用内部默认颜色标号13.基本绘图命令set line color style thickness 设置颜色、线型、线宽draw line x1 y1 x2 y2 绘制直线draw rec/recf x1 y1 x2 y2 绘制不填色/填色矩形draw mark type x y size 绘制标记*1叉2空心圆3实心圆4空心方框5实心方框draw string x y string 绘制字符串set strsiz hsiz <vsiz> 设置字符大小14.文件操作(1)输出二进制变量设置维数环境‘set gxout fwrite’‘set fwrite ***.dat’‘d var’(2)读写有格式文件result=read(filename)result 第一行：返回码（0 OK，1打开错误2文件结束8文件状态为写入9I/O错误）result 第二行：读取的文件内容Text = sublin(result,2) 将result相应行的内容赋值给Textaa=subwrd(Text,#n) 将Text的第n个值赋值给aaresult=write(filename, text, <append>)15.维数查询q w2xy lon lat 将经纬度坐标转换为图形窗口坐标返回值：X=*** Y=*** q xy2w X Y 将图形窗口坐标转为经纬度坐标返回值：Lon=*** Lat=***q w2gr lon lat 将经纬度坐标转换为格点坐标返回值：Xdim=*** Ydim=***q gr2w xdim ydim 将格点坐标转换为经纬度坐标返回值：Lon=*** Lat=***q xy2gr X Y 将图形窗口坐标转换为格点坐标返回值：Xdim=*** Ydim=***q gr2xy xdim ydim 将格点坐标转换为经纬度坐标返回值：X=*** Y=***16.气候平均设置季节变量‘set t 1 12’‘define climate=ave(var,t+0,t=tmax,12)’‘modify climate seasonal’滑动平均‘set t 1 tmax’‘define running_mean=ave(var,t-n,t+n)’17.站点资料绘图(1)转换生成数据文件及格点文件数据文件的顺序为先循环站号再循环时间，每条记录包含站号-纬度-经度-时间-nlev-flag-数据，每循环完一个时间后要加一个不写数据且nlev=0的记录(2)编写ctldtype stationstnmap ***.map**变量层数为0(3)生成.map映射文件!stnmap –i */*/*.ctl(4)插值‘open ***\grid.ctl’‘open ***dat.ctl’‘define varname=oacres(grid,var.2)’(5)绘图‘cnbasemap varname’18.绘制路径图。

第4章 GrADS各类参数设置

1. 对于图形类型为contour起作用的设置
（1）set ccolor color - 设置等值线颜色。
其中，color为颜色号：0 - 黑， 1 - 白， 2 - 红，3 - 绿，4 - 兰， 5 - 青，6 - 洋红，7 - 黄，8 - 橘黄，15 - 灰。该设置在“clear”或“display”命令后即重新设定颜色。
（8）set clskip number – 设置间隔几根(number-1)等值线标示数值。
2. 对于图形类型为contour或 shaded起作用的设置
（1）set cint value - 设置等值线间隔，该设置在“clear”或“display”命令后即重新设置。（2）set clevs lev1 lev2 ...- 设置特定的等值线值，只画lev1，lev2，.......值所在的等值线，用于不等间隔绘图，该设置在“clear”或“display”命令后即重新设置。
（6）set cmax value –设置不画高于此value值的等值线。该设置在 “clear”或“display”命令后即重新设置。（7）set cblack val1 val2–设置不画介于val1和val2之间的等值线。该设置在“clear”或“display”命令后即重新设置。
（4）set csmooth on|off - 设置是否将网格值重新插值。如取“on”，则在绘等值线图前用三次插值将现网格值插到更精细网格上，重新设置才改变本次设置。
（5）set cmin value –设置不画低于此value值的等值线。该设置在 “clear”或“display”命令后即重新设置。
数据描述文件
dset d:\model.le.dat UNDEF -2.56E33 TITLE 5 Days of Sample Model Output XDEF 72 LINEAR 0.0 5.0 YDEF 46 LINEAR -90.0 4.0 ZDEF 7 LEVELS 1000 850 700 500 300 200 100 TDEF 5 LINEAR 02JAN1987 1DY vars 8 ps 0 99 Surface Pressure u 7 99 U Winds v 7 99 V Winds z 7 99 Geopotential Heights t 7 99 Temperature q 5 99 Specific Humidity ts 0 99 Surface Temperature p 0 99 Precipitation ENDVARS

S7 GRAPH指令

所有指令介绍13.1常用指令及无联锁指令在运行过程中执行动作动作注释地址范围事指令地址位置时间常数件N Q,I,M,D*M.N当前步激活时，信0.0-65535.7号置1，步失效时置0。

（无关互锁）步激活时，信号置1 S Q,I,M,D*M.N并保持置1，直到下0.0-65535.7次复位（无关互锁）R Q,I,M,D*M.N步激活时，信号值0.0-65535.70，并保持。

（无关互锁）0.0-65535.7D Q,I,M,D*M.N T#<const>当前步激活，延迟T#TIME时间后置1，当前步激活时间小于T#TIME则不置1，当前步失效后置复位。

（无关互锁）步激活时，信号值L Q,I,M,D* M.NT#<const> 当前步激活时，持续 T#TIME 时间置1，时间过后置 0，当前步失效后复位。

（无关互锁）0.0-65535.7CALLFB,FC,SFB,S FCBlock nu m ber只要当前步激活，就调用其他功能块。

（无关互锁）NC Q,I,M,D* M.N 当前步激活时，在0.0-65535.7满足互锁条件下，信号置 1，不满足互锁条件下置 0。

当前步失效后复位。

步激活时，在满足互锁条件下，信号 SC Q,I,M,D* M.N0.0-65535.7置 1 并保持置 1，直到下次复位.RC Q,I,M,D* M.N0.0-65535.70，并保持。

（无关互锁）DCQ,I,M,D*M.NT#<const> 当前步激活，延迟T#TIME 时间后置1，当前步激活时间小于 T#TIME 则不置 1，当前步失效后置复位。

（无关互锁）LCQ,I,M,D*M.NT#<const> 当前步激活时，持续 T#TIME 时间置 1，时间过后置 0，当前步失效后复位。

（无关互锁）0.0-65535.70.0-65535.7CALLCFB,FC,SFB,S FCBlock nu m ber只要当前步激活，就调用其他功能顺控器 FB 块引脚参数说明：FB 参数（上升沿有效）内部变量（静态数据区名称）块。

GrADS快速参考手册教学内容

w2xy lon lat
给出经度－纬度坐标值(lon,lat)，返回xy坐标值
w2gr lon lat
给出经度－纬度坐标值(lon,lat)，返回网格坐标值
gr2w v1 v2
给出网格坐标值(v1,v2)，返回经度－纬度坐标值
gr2xy v1 v2
给出网格坐标值(v1,v2)，返回xy坐标值
ll2xy lon lat
-C
以两位数记年时，自动设定世纪值
一般性设置
列出帮助选项
set grads on|off
打开/关闭显示“GrADS” 标志
设置显示模式<和色彩方式>为：
grey|greyscale
用灰度方式
color <black|white>
以彩虹模式绘图
set frame <option>
设置图形边框绘图方式为：
以一个矩形包围所有图形
off
不画边框
circle
对于其它投影方式画矩形边框，而极射投影画圆形边框
set background index
设置绘图区背景颜色
display expression
d expression
显示expression代表数据的图形
open control_filename
打开一个数据描述文件，可指明路径
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息

Grads精细绘图

GrADS 精致绘图说（一）1. basemap.gs:basemap L | O | U <fill_color> <out_color> <hi/lo>在低分辨率海岸廓线范围内用颜色覆盖陆地／海洋。

适用于各种投影方式，需lpoly.asc, lpoly_hires.asc, lpoly_US.asc, opoly.asc, opoly_hires.asc文件。

其中：L(l)：覆盖陆地，O(o)：覆盖海洋，U(u)：覆盖20N-50N的墨西哥和加拿大领土（低分辨率，适用美国），fill_color：填充色号，缺省为15，out_color：廓线的颜色号，缺省为15，hi/lo：高分辨率('set mpdset hires'，仅对15N-53N, 130W-60W区域)/低分辨率。

2. cbar.gs、cbarn.gs、cbarc.gs、cbar_l.gs和cbar_line.gs:cbarn sf vert xmid ymidcbarc center_x center_y back_color绘制'set gxout shaded'图形的填色标尺。

sf：色标尺寸，1为全尺寸，0.5为半尺寸；vert：0为水平，1为垂直；xmid,ymid：色标中心点的位置。

cbar_l -x X -y Y -n number -t text -pcbar_line -x X -y Y -c color -m mark -l linestyle -t text -p加'set gxout line'的图例说明。

其中：-x,-y：图中x和y的位置，-n：线条的数目（最多可为10条），-t：文字说明的内容（最多10条，需双引号括起），-c：线和标记的颜色，-m：定义标记；-l：定义线型，-p：用户可在图中点击给定图例的放置位置。

3. colors.gs:为雪盖(颜色序号40~45)、降水(颜色序号50~59)及温度(颜色序号64~85)资料的shaded图设置填充色。

GrADS实用绘图--图形显示

字符大小
set digsize #inch
特征值马赛克图 set gxout fgrid 指定特征值/颜色
set fgvals v1 #c1 v2 #c2 v3 #c3 …
格线
set gridln auto ~ ~ off ~ ~ #color
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件在实际空间上的结构
高度
经度
UNIT=20 的数据文件在实际空间上的结构
高度
纬度
纬度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
• 网格坐标系 x, y, z, t ~ 如同数组下标 • 实际坐标系 lon, lat, lev, time • 缺省维数环境
15
四. 图形显示(D)
filled bar linefill
— 1-D序列图示例
line
line
outline bar linefill
‘set t 3 9’ ‘define red=lat-lat’ ‘set t 1 12’ … ‘d red’
set xlevs x1 x2 x3… set ylevs y1 y2 y3…
(李晓峰, 2005)
20
LASG动力论坛版
四. 图形显示(D)1
网格填值图 set gxout grid 小数位数
set dignum number
— 2-D单变量图2
马赛克图 set gxout grfill 颜色

几个有用的GRADS函数

几个有用的GRADS函数monmask多年时间序列中显示每年的相同月的资料<run> monmask.gs 5 define mmask for May (5)d (expr)*mmask plot (expr) only for Maycbar_line: 显示图例exp: cbar_line -x num(x方向位置) -y num(y方向位置) -c num(线条颜色) -l num(线条类型) -n num(线条数，max=10) -m num(标记) -t "(文本)" -p (加上-p，就不要-x, -y, 鼠标点到哪里就画在那)cbar: 对shaded图形有效，做图例cbarn: 同上，可加参数，exp: cbarn 1(图例粗细) 1(图例竖放，横放=0) <5(图例位置，x方向)> <5(图例位置，y方向)>cbarc: 同上，默认图例在右上角，可加参数，exp: cbarc <5(图例位置，x方向)> <5(图例位置，y方向)>connect_the_dots: 用鼠标画线，点击图的下方退出，无参数font: 设置字体，参数0-5，默认font=0, ROMAN字体为font=1loop: 循环作图，exp: loop x(or y，z，t，lat，lon，lev) min(起始值) max(最终值) step(步长)page: 选择画图区域，usage:run page top plot on top half of screenrun page bottom plot on bottom half of screenrun page left plot on left half of screenrun page right plot on right half of screenrun page q1 plot on quadrent 1run page q2 plot on quadrent 2run page q3 plot on quadrent 3run page q4 plot on quadrent 4run page xy XY (x,y,X,Y) are single digitsxy = location, XY = sizeonly 1 blank between xy and XYquadrents:-----------| q1 | q2 |-----------| q3 | q4 |-----------Matrix notation:xy is the location of the pageXY is number of subdivisionsThe idea here is that the page is divided into a grid of X by Yrectangles. Plotting is done in the x-y rectangle. In thefollowing example, the page is divided into two columns (X=2)and 3 rows (Y=2).-----------| 13 | 23 |-----------| 12 | 22 | xy values for XY=23-----------| 11 | 21 |-----------To draw into the top left rectangle, you "run page 13 23".To get the bottom right rectangle, you "run page 21 23".The following commands are equivalennt:"run page q1" "run page 12 22""run page q2" "run page 22 22""run page top" "run page 12 12""run page bottom" "run page 11 12"Limitations: X and Y must have values from 1 to 9.关掉多图形显示：page offsave: 保存图形，exp: save <路径>文件名.gmfwrdlin: 用经纬度画直线，exp: wrdlin minlon(起始经度) minlat(起始纬度) maxlon(最后经度) maxlat(最后纬度)wrdmrk: 在经纬度上标记，exp: wrdmrk num(标记类型) lon(经度) lat(纬度) siz(大小)wrdrec: 在图上画矩形，exp: wrdrec minlon(起始经度) minlat(起始纬度)[矩形左下角坐标] maxlon(最后经度) maxlat(最后纬度)[矩形右上角坐标]wrdrecf: 掩盖矩形内的图形，exp: wrdrec minlon(起始经度) minlat(起始纬度)[矩形左下角坐标] maxlon(最后经度) maxlat(最后纬度)[矩形右上角坐标]wrdstr: 在图上写字符串，exp: wrdstr lon(经度) lat(纬度) "(字符串)"map: 设置地图投影，nps (north-pole stereographic) lat='20 90' lon='-270 90'nps2 (north-pole stereographic) lat='20 90' lon='-180 180'nps3 (north-pole stereographic) lat='20 90' lon='-60 300'sps (sorth-pole stereographic) lat='-90 -20' lon='-270 90'lola (lat-lon) (0..360)'lola2 (lat-lon) (-180..180)'usa 'latlon' lat='24 52' lon='-127 -65'usa2 'nps' lat='15 80' lon='-150 -45'n_amer 'nps' lat='5 90' lon='-270 90's_amer 'latlon' lat='-60 20' lon='-90 -30'africa 'latlon' lat='-40 50' lon='-20 60'europe 'nps' lat='5 90' lon='-180 180'euro2 'nps' lat='5 90' lon='-180 180'asia 'latlon' lat='0 80' lon='40 170'aust 'latlon' lat='-50 0' lon='100 180'robinson 'robinson' lat='-90 90' lon='-180 180'c_pac 'latlon' lat='-45 45' lon='120 290'n_pac 'nps' lat='0 90' lon='90 270'custom [lon0 dlon lat0 dlat] (lat0,lon0=left-bottom corner)。

GrADS学习资料：第7章 GrADS 描述语言及编程(1)

在i、j变量被替代前后，复合变量名的整个长度不能超过16个字符。
在GrADS描述语言中不能设置过多的变量，因此，不能利用复合变量产生一个较大的数组，例如：
i=1 while (i<=10000) var.i = i i=i+1 endwhile
显示1-3月的图和标题
t.1='Jan';t.2='feb';t.3='mar' i=1 while(i<=3) 'set lev 850' 'set t ‘i’' 'set lon 70 140' 'set lat 10 60' 'set gxout vector' 'd v;v.2' 'draw title uv for 't.i'' 'print' 'c' i=i+1 endwhile
while语句
while （表达式）描述语言语句描述语言语句 … … endwhile
当表达式结果为真，即：是一个不为0 的字符，循环就继续执行。
lon1=0 lon2=357.5 lat1=-90 lat2=90
iyear=1948 while(iyear<=2002) say ''iyear'' pull dummy m=1 while(m<=12)
复合描述变量：可以在描述语言程序中构造数组。复合变量名可以是一个由“·”分隔的多个变量名的复合体。例如： varname.i.j
如果i和j也是变量，那么i和j的内容将会被i和j的字符串值所替代。 i = 10 j=3 varname.i.j = 343 在这个例子中，赋值语句相当于：