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命令大全
设置彩虹色颜色序列值（缺省：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命令大全
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
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绘图的各种设置
draw xlab string draw ylab string 分别在水平坐标轴或垂直坐标轴左侧写字符串
string，作为x轴或y轴的说明。
set annot color <thickness> 设置上述图注所用的颜色和线宽，缺省为白色，
线宽为6，该命令同时设置了坐标轴线以及刻 度的颜色和线宽，坐标刻度和标尺的线宽为图 注设置的线宽再减1。
例：draw title Geopotential Hight (October 2010)
set xyrev on 交换水平和垂oj 设置当前地图投影方式。
常用proj取值包括： latlon：缺省设置，用固定的投影角进行Lat/lon投影
； scaled：用不固定的投影角进行latlon投影，地图比例
失效； nps | sps：北半球|南半球极地投影； robinson：robinson投影，x：-180 180；y：-90 90 off：不画出地图，坐标轴也不代表lat/lon。
（缺省时即按数据描述文件的维数进行绘图。）
特别地， 当数据文件为全球多层次多时次资料时，例如全球
144*73*12的多时次资料，缺省时维数为x、y方向为 （0-360E，90S-90N），z方向和t为第一个点。
GrADS最擅长画2-D图，即全球数据的缺省维数。
维数说明分为两类：
set lon | lat | lev | time val1 <val2> set x | y | z | t val1 <val2> 两者对应于同一组数据，前者是地球坐标
set grid on|off|value|horizontal|vertical <style> <color>
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中文说明
目录前言 (5)第一章启动GrADS 及打开文件 (5)一.启动GrADS (5)二.打开文件 (6)第二章维数环境 (6)set lat|lon|lev|time val1<val2> (6)set x|y|z|t vall <val2 > (6)第三章显示命令及动画 (7)一、显示命令 (7)二、动画 (7)set loopdim x︱y|z|t (7)set looping on|off (7)第四章图形输出类型 (8)对格点数据： (8)bar: 直方图 (8)barb:风向杆形式绘二维风场 (8)contour:二维等值现图 (8)errbar:单线图及误差分布 (8)frgid:用指定颜色填充二维格点场 (8)fwrite:图形不在屏幕上显示，而是将输出结果存入一个由“set fwrite 文件名”所指定的文件中。
(8)gid:以网格形式在各网点中央标出该点数值 (8)line:单线图 (8)linefill:两单曲线之间填色 (8)scatter:散点图 (8)shaded:二维填色图 (8)stat:输出有关的统计量 (8)stream:流线形式绘二维风场 (8)vector:向量间箭头形式绘二维风场 (8)对站点数据： (8)barb:在各站点绘风向标 (8)findstn:搜索最近的站点（详见描述语言部分） (8)model:以天气图形式将天气观测各分量填放在站点四周 (8)value:在各站点标值 (8)wxsym:绘wx天气符号 (8)Draw wxsym symbol x y size<color<thicks>> (9)对于model,命令为set gxout model ,然后再添图，形式为displayu;v;t;d;slp;delta;cld;wx;vis (9)第五章 GrADS资料格式 (9)格点资料资料描述文件 (9)以下是一个资料描述文件的例子： (9)１、DSET data-set-name (11)２、TITLE string (11)３、UNDEF value (11)4、OPTIONS <keywords> (12)5、XDEF number<LINEAR start increment>或XDEF number<LEVELS value-list> (12)6、YDEF number mapping start<increment>或YDEF number<LEVELS value-list> (13)７、ZDEF number mapping <start increment>或ZDEF number<value-list> (13)８、TDEF number LINEAR start-time increment (14)９、VARS number (14)二. 台站资料资料描述文件 (15)第六章变量名和表达式 (15)一、变量名和表达式 (15)二、表达式 (16)第七章函数 (16)一. aave函数 (16)二．abs函数 (16)三．acos函数 (17)四、asin函数 (17)五、atan2函数 (17)六、ave函数 (17)七、cdiff函数 (17)八、const函数 (18)九、cos函数 (18)十、exp函数 (18)十一、gr2stn函数 (19)十二、hcurl函数 (19)十三、hdivg函数 (19)十四、log函数 (19)十五、log10函数 (19)十六、mag函数 (19)十七、maskout函数 (20)十八、oacres函数 (20)十九、pow函数 (20)二十、sin函数 (20)二十一、skip函数 (20)二十二、smth9函数 (21)二十三、sqrt函数 (21)二十四、stnave函数 (21)二十五、stnmin函数 (21)二十六、stnmax函数 (21)二十七、tan函数 (22)二十八、tloop函数 (22)二十九、tvrh2q函数 (22)三十、tvrh2t函数 (22)三十一、vint函数 (23)第八章图形要素设置 (23)对于图形类型为contour起作用的设置 (23)set ccolor color——设置等值线颜色 (23)set ccolor rainbow——设定等值线颜色用彩虹表示 (23)set cstyle style——设定等值线线型 (23)Set cthick thckns——设定等值线线宽hckns(取值范围1～20) (23)Set cterp on｜off——设置样条差值光滑开关 (23)Set clab on｜off｜forced｜string｜auto——控制等值线的标记方式。
GrADS第4章 各类参数设置
'set gxout line' 'set cstyle 1' 'set cmark 0' 'set ccolor 2' set lfcols col1 col2 'd v.4' 'set gxout line' d line1;line2 当line1<line2时，用颜色1（col1） 'set cstyle 1' 'set cmark 0' 当line1>line2时，用颜色2（col2） 'set ccolor 1' 'd v.5 ' 'set gxout linefill' 'set lfcols 4 10' 'd v.4; v.5'
不固定 不固定
不固定
一维曲线图
4.2 图形类型的设置
当维数环境确定后，缺省情况下， 一维变量输出的图形为单线图，二 维变量为等值线图，改变缺省图形 输出类型的命令为： set gxout graphics_type
1. 格点数据
contour: 二维等值线图（缺省设置） set gxout contour shaded: 二维填色图(通常与contour并用) 先画填色图，再画等值线图；否则等值线会 被覆盖。 set gxout shaded d hgt set gxout contour d hgt
3. 维数环境的定义
（1）一种是地球坐标（world coordinate）
以经纬度为度量单位。 形式如下： set lon|lat|lev|time val1 <val2> （2）一种是格点坐标（gridcoordinate） 以网格点数为度量单位。 形式如下： set x | y | z | t val1 <val2>
GRADS软件的常用GS文件
while(i<=593)
aa=read('e:eddy-statisticdens-westpacific-1.txt')
aa1=sublin(aa,2)
lat1=subwrd(aa1,3)
lon1=subwrd(aa1,2)
data1=subwrd(aa1,8)
之字符：'draw string 9.09 3.03 `52900141' *then you can see bold roman numbers 2900141
之符号：'draw string xxx xxx `3#' then you can see a arrow,remember this the symbol must come together with the number "3"
aa1=sublin(aa,2)
lat=subwrd(aa1,5)
lon=subwrd(aa1,4)
data=subwrd(aa1,6)
'q w2xy 'lon' 'lat
x = subwrd(result,3)
y = subwrd(result,6)
y = subwrd(result,6)
经纬度对应图上位置——之2
lon='xxx.x xxx.x xxx.x'
lat='xx.x xx.x xx.x'
j=1
cn=999
while(j<cn)
'q w2xy '%subwrd(lon,i)%' '%subwrd(lat,i)
第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第4章各类参数设置.
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout grid du
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout fgrid set fgvals -1 4 Set fgvals -3 5 du
3. 维数环境的定义
（1）一种是地球坐标（world coordinate） 以经纬度为度量单位。
形式如下： set lon|lat|lev|time val1 <val2>
（2）一种是格点坐标（gridcoordinate） 以网格点数为度量单位。 形式如下：
set x | y | z | t val1 <val2>
1. 维数环境的概念
GrADS系统视每一个物理变量（VAR） 场为一个四维数据集（4D data set）， 即包括三维空间（x，y，z）和一维时间 （t）。也可固定其中的一维或几维以获 得实际的低于四维的数据子集。
2.设置的作用
GrADS 中设置维数用以说明或指定随 后的分析或图形操作时参加操作的原 始数据集的维数范围。该工作数据集 可以是整个原始数据场，也可以是原 始数据场的一部分。
或 set lev 1000 850
set t 1 3
或 set time 01jan1985 03jan1985
注意: grads只能显示两维或一维不固定的图！
当所有维数都固定时，得到的是一个单 值数据点
维数环境的设置一直保持到再次设定时 都有效。要查看当前维数情况，可以在 运行环境中输入命令： ga-> q dims
第4章 各类参数设置
第十章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常用命令
第一个命令：这样将得到一个动画显示的序列，现在输入：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常用指令
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.绘制路径图。
GrADS命令大全
set imprun script-name
在每个display命令前自动执行一遍script-name文件中的命令组
run file-name <params> file-name <params>
按顺序执行file-name文件中的命令组
exec fname <arg0,…,arg9>
GrADS快速参考手册
(GrADS Version 1.7 beta7 )编：Meier-fleischer@dkrz.de译：朱禾
命令行参数grads -options
-c‘命令组’
启动GrADS后执行该命令组
-b
以批处理方式运行
-l
以landscape模式运行
-p
以portrait模式运行
-C
以两位数记年时，自动设定世纪值
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>
11
椭圆0
绘图参数设置
set clip xlo xhi ylo yhi
设置剪裁区x范围y范围（英寸）
设置等值线颜色，index=0,1,2…
set cstyle style
设置线型,style=0,1,2…
set cmark mark_index
GrADS使用简介
GrADS使用简介谢逸乔云亭1 前言GrADS(Grid Analysis and Display System)是美国马里兰大学气象系Brian E.Doty 开发的气象数据分析与显示软件。
GrADS提供一个全32位交互操作的气象格点数据与站点数据的分析与显示环境,具有强大的气象数据分析能力、丰富的地图投影坐标,以及彩色动画等功能,因此迅速成为国内外气象数据显示的标准平台之一。
近几年,随着PC机的普及、性能的不断提高,GrADS推出了基于PC机的DOS版本, 为需要用到GrADS,却不具备UNIX图形工作站的广大气象工作者带来了方便。
2 软件安装与硬件配置微机版本的GrADS,运行于386或以上的PC机上,也可以在WIN95或WIN NT的DOS环境下运行。
安装GrADS时,将整个目录拷贝到硬盘的某个目录下(例如:C:\GrADS),并根据显示器的分辨率类型,用Set命令在DOS下设置GAVIDEO参数。
有以下几种:Set gavideo=vga VGA,16 color,640*480Set gavideo=ega EGA,16 color,640*350Set gavideo=vga256 VGA,16 color,640*480Set gavideo=svga SVGA,256 color,800*600Set gavideo=xvga XVGA,256 color,1024*768如果显示器和图形卡(显示卡)不支持16色以上的VGA模式,不可强行将参数设成VGA,16色或以上,否则会损坏显示器。
完成以后,在系统的Config.sys文件中加入:shell=c:\/p/e:20483 GrADS文件组成GrADS包括以下几种类型的文件:(1)二进制无格式的原始气象数据文件(*.dat)其数据既可以是格点数据,也可以是站点数据。
它们是从其它气象数据(如站点气象报、格点气象报、模式格点输出结果)转换生成的,格式是二进制无格式数据。
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 在这个例子中，赋值语句相当于：
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一个命令：这样将得到一个动画显示的序列，现在输入：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 选择默认的画线间隔。
我们可以这样自己设定：set cint 30d z也可以控制线的颜色：clearset ccolor 3d z可以使用交替的方式显示数据：clearset gxout shadedd hcurl(u,v)这看起来不很平滑，可以使用三次平滑：clearset csmooth ond hcurl(u,v)我们也可以使用不同的图形样式set gxout contourset ccolor 0set cint 30d z使用注解：draw title 500mb Heights and Vorticity风矢量：clearset gxout vectord u,v这里我们显示两个表达式，第一个是U分量矢，第二个是V分量矢。
我们也可以使用第三个字段使图像彩色化d u;v;q或者：d u;v;hcurl(u,v)你也可以显示任何区域以显示伪矢量：cleard mag(u,v);q*10000此处U分量是风速和V分量是湿度。
我们也可以查看流线（并且着色）clearset gxout streamd u;v;hcurl(u,v)或者我们可以查看实际的格点值clearset gxout gridd u我们可能希望改变地图背景：clearset lon -110 -70set lat 30 45set mpdset namset digsize 0.2set dignum 2d u改变投影：set lon -140 -40set lat 15 80set mpvals -120 -75 25 65set mpproj npsset gxout contourset cint 30d z我们让grads可以操作纬度从140W到40W之间，经度从15N到80N之间的数据。
但是我们告诉它使用从120W到75W和25N到65N极坐标投影绘图。
多余的部分被投影程序修剪掉了。
该简单教程到此结束。
此时，你可能希望深入查看数据集，或者你可能希望浏览GrADS文档并且尝试此处描述的其它选项。
英文原文并没有完全按原文翻译。
reinit 初始化、关闭文件q/query file 查询文件,xyz格点数，时间，变量，（变量层<小于z>，）q ctlinfo 路径名字缺色值xyz格点数，初始值，stepsopen 打开文件d/display (temp) 绘图输出c/clear 清除draw title (name)给图添加名字,如果加上反斜杠,名字分两行保存文件enable print <routine<file.gmf>>printdisable print格式转换,gv32,打开文件.gmfquit 退出系统维数环境设置set lat纬度（负南纬）/lon经度/lev（海平面之上百帕以下位势米）/time val1<val2>set time 12zsep格点坐标set x/y/z/t val1<val2> 如，set x 10 30 set t 1 40改变缺省图形处处类型set gxout graphics_typeset z 40（不用此命令，默认输出第一层）set gxout shaded（二维填色图）（用于等温图）添加温度与颜色的对照图cbarn/cbard tempvector矢量形式绘二维风场set gxout vectord u;vd skip(u,4,4);skip(v,4,4)每四个格点取一个计算风速大小d sqrt(u*u*,v*v*)d mag(u,v)chapter 2 图形要素设置set ccolor <color number>设置图形颜色set cstyle <style number>设置等值线线型set cthick <thickness number>设置等值线线宽set cterp on/off 设置样条差值光滑开关off可使添色图与等值线的边缘重合set clab on/off/forced/dtring/autoon快速标记，缺省值水平位置offforce强迫标记所有等值线clab %g<需要的字符如%gK,%g%%>set clopts color<thickness><size> 默认值为-1,-1,0.09set cint value 设置等值线间隔,clear之后重设set clevs lev1 lev2设置特定的等值线,即lev1,2 的等值线显示set ccols col1 col2设置等值线的颜色,各数值对应于set clevs命令对应的个等值线set cmin/cmax value 画出不低于/高于value的等值线set black val1 val2 画出介于v1 v2之间的等值线clear之后均重置set stream value 设置流线密度,value值在1到10之间,5为缺省set rbcols clo1 clo2 <clo3> 设置新的彩虹序列,set color 用set rgb命令定义set rbcols auto启用内定的彩虹颜色set rbrange low high设置彩虹颜色序列对应的等值线范围,介于最高最低值之间的等值线按彩虹排列,之外的同色(两种)对于line类型的设置set cmark cmarker设置线条上的标记符号0无标记set wrang vlo vhi 设置坐标轴通常指y轴标尺的取值范围set missconn on/offvectorset arrscl size 设置矢量箭头长度0.5-1.0set arrowhead sizeset zlog on/off z方向取对数尺度的开关,重设后才改变set xaxis/yaxis start end <incr>设置坐标轴的起始值和结束值,以incr作为增量set grid on/off 设置是否显示网格线set annot color<thickness> 设置坐标轴的线宽和颜色set xyrev on 交换水平竖直坐标轴set xflip/yflip on 水平/竖直坐标轴维数转向地图投影设置set mproj proj设置地图投影方式nps/sps北/南半球基地投影set map color style thicknessset mpdraw on/off off不绘地图背景,但是地图标尺起作用chapter 5页面控制set vpage xmin xmax ymin ymax 设置一个给定的虚页set vpage offset parea xmin xmas ymin ymax 在给定位置绘图set parea offchapter 6 基础绘图指令draw string x y string在xy位置上写一串字符string性质由set string 和set strsiz命令指令，如”drawstring 2 7 sunjia”set strsiz hsiz<vsiz>设置字符大小,hsiz高度vsiz宽度set string colordraw line x1 y1 x2 y2set line color<style <thickness>设置线条draw rec xlo ylo xhi yhi 画一个矩形draw recf xlo ylo xhi yhi 在矩形中添色draw mark marktype x y size 在xy位置处画一个标记为marker的标记chapter 7 函数ave()求平均值aave(expr,x=x1,x=x2,y=y1,y=y2)求面积平均abs(expr)绝对值cdiff(expr,dim)中央差分函数hcurl(u,v)涡度mag(u,v) 求风速chapter 9 定义变量define varname=vpr。