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中文说明手册
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分量<；标量>。
“<>”内的部分只起标颜色的作用。
）以流线方式表示矢量场：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不画国界省界等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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
列出第n个数据描述文件的信息，n=1,2,3……
set zlog on|off
设置z坐标取对数坐标,或关闭此设制
set strmden value
设置流线密度[1-10],缺省为5
set stnopts <dig3> <nodig3>
以后3位数标数值
set mdlopts noblank|blank|dig3|nodig3
以3位数标模式值
set stid on|off
set gxout graphic-type
graphic-type可取以下参数：
bar
柱状图形模式
model
不规则分布测站模式
rb
WMO风标表示的风矢量图
scatter
不规则分布数据模式
contour
等值线图模式
shaded
分色图
errbar
误差图模式
stat
按WMO规范显示测站观测数据
fgrid
测试s代表的字符串的长度
defval vl i j
为第（i,j）网格点重新赋值, i,j可以取小数
udft
显示用户自定义函数列表
lats
显示GrADS－LATS界面信息
xy2w v1 v2
给出xy坐标值(v1,v2)，返回经度－纬度坐标值
xy2gr v1 v2
给出xy坐标值(v1,v2)，返回网格坐标值
地图投影方式
设置地图投影方式，关键字proj取：
等比矩形经纬度投影(缺省)
不等比矩形经纬度投影
北半球极射投影
南半球极射投影
robinson投影
orthogr投影
mollweide投影
兰勃托投影
关闭投影
设置投影关闭|或颜色（或－1，由系统决定）、线型、粗细, type=0~255或*-任意类型投影
为极射投影设置标准经度 lnmn，lnmx和标准纬度ltmn，ltmx
短虚线
4
一长一短虚线
5
点线
6
点划线
7
点点划线
GrADS缺省定义的12种标注和它实际代表的意义
0
不作标记
1
十字
2
圆圈○
3
圆圈●
4
方块□
5
方块■
6
X
7
菱形◇
8
三角△
9
不作标记
10
11
椭圆0
绘图参数设置
set clip xlo xhi ylo yhi
设置剪裁区x范围y范围（英寸）
设置等值线颜色，index =0,1,2…
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
set lfcols col1 col2
为gxout linefill模式设置颜色
set black off|val1 val2
设置不画在val1 val2之间的等直线/关闭
set rbcols c1 c2<c3…cn>
设置彩虹色颜色序列值（缺省：9，14，4，11，5，13，3，10，7，12，8，2，6）
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息
fwrite
返回当前fwrite命领输出文件的名称
string s
z 方向网格坐标值为vl1或从vl1到 vl2范围
页面控制
set vpage off
关闭虚页设置，此时物理页等于虚页
set vpage xmin xmax ymin ymax
设置虚页范围，单位： 英寸
set parea xmin xmax ymin ymax
在虚页内，设置绘图页面大小（英寸）
图形显示模式控制
在（x，y）画大小为size的标记，marktype＝1～11
画多边形
在（x，y）画大小为size，颜色为color，粗细为thickness的气象标记，symbol＝1～41
设置矩形剪裁区
基本字符串绘图命令
设置字符串绘图的颜色<对齐><粗细><旋转>，对齐＝l、c、r、tl、tc、tr、bl、bc、br，旋转＝＋度/－度，顺时针/逆时针
hbuff
在双缓冲区模式下，清空显示缓冲区
reset <option>
清干净图形显示区并重新初始化，但1）不关闭任何文件，2）不取消变量定义，3）不改变set display 设置
option＝events; graphics; hbuff; norset
reinit
与reset相似，但关闭所有文件并取消所有定义
设置层取val1或从val1到 val2范围
set t val1 <val2>
设置时间取val1或从val1到 val2范围
set x vl1 <vl2>
x 方向网格坐标值为vl1或从vl1到 vl2范围
set y vl1 <vl2>
y 方向网格坐标值为vl1或从vl1到 vl2范围
set z vl1 <vl2>
-C
以两位数记年时，自动设定世纪值
一般性设置
列出帮助选项
set grads on|off
打开/关闭 显示“GrADS” 标志
设置显示模式<和色彩方式>为：
grey|greyscale
用灰度方式
color <black|white>
以彩虹模式绘图
set frame <option>
设置图形边框绘图方式为：
quit
退出
!shell-command
在GrADS命令交互模式下运行外部系统命令
维数参数设置
set lon val1 <val2>
设置经度取val1或从val1到 val2范围
set lat val1 <val2>
设置纬度取val1或从val1到 val2范围
set lev val1 <val2>
标数字时，指定数字大小（英寸）
set arrlab on|off
设置是否标注矢量图题注（箭头）
set arrsclsize <magnitude>
矢量图题注箭头放大缩小(<1)比例
set arrowhead size
矢量图题注箭头尺寸（英寸）
set fgvals v1 c1 <v2 c2>…
为gxout fgrid模式设置数值、颜色对应列表
x轴数值由大到小（缺省：由小到大）
y轴数值由大到小（缺省：由小到大）
x轴标注|不标注|自动|指定字符串
y轴标注|不标注|自动|指定字符串
指定x轴按lab1 |lab2 |…标注
指定y轴按lab1 |lab2 |…标注
x轴中心下方写字串string
y轴中心左侧写字串string
设置x轴标注的颜色<粗细<大小>>
set cmin value
设置最小等直线的值，小于的不再画出
set clevs lev1 lev2…
按lev1 lev2…指定值绘等直线
set ccols col1 col2…
按指定值绘等直线时，为每一条等直线指定颜色