GrADS第5章 基础绘图指令
GrADS基础绘图
绘图速度较慢，当数据文件较大时，读入很容易内存溢出 远程处理比较麻烦。
grads的优点; （1），最多能读入20个数据文件，数据读取迅速，出图速度快，但是要写描述文件和gs文 件比较麻烦。 （2），在画二维等值线图时，图形非常平滑，这点是我到现在还没大规模用matlab的主要 原因。 （3），它包含了很多气象上所需要的函数，在gs文件里可以直接调用，因此在用它时也省 了不少事。
Grads在win/Linux下的下载、安装、运行
程序：
Windows：grads-1.8sl11-win32e.exe Linux：grads-bin-1.9b4-linuxRH9.tar.gz grads-1.8sl11-linux.tar.gz data.tar.Z
环境变量：
Windows：
Panoply OpenDX (formerly IBM Data Explorer) NCL (NCAR Command Language) ncBrowse (netCDF File Browser) HDF (Hierarchical Data Format) interface GrADS (Grid Analysis and Display System) FERRET NetCDF Toolbox for MATLAB-5 IDV (Integrated Data Viewer) Origin surfer
命令顺序 grads –b|l|p|c open/run set display draw print quit/reinit
显示类型
GrADS实用绘图--基本操作
• 关闭数据文件
ga-> close 3 ↵ ga-> close 2 ↵ ga-> close 1 ↵
文件打开从1 开始依次编号, 文件关闭时先 关闭序号最大 的文件
26
LASG动力论坛版
三. 运行GrADS绘图软件包
• 设置维数环境
ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵
♥♥♥
18
LASG动力论坛版
二. GrADS绘图软件包的安装
• Windows下的安装(2)
– 设置环境变量 ♥♥♥ “我的电脑”右键菜单 属性 高级 环境变量
• 新建变量 GASDIR c:\PcGrADS\dat GASCRP c:\PcGrADS\lib GAUDFT c:\PcGrADS\udft\udft.list • 修改PATH变量, 增加路径 c:\PcGrADS\win32;c:\PcGrADS\lib;c:\PcGrADS\dat;
8
LASG动力论坛版
一. 引
• GrADS实用文档
言
– LASG编 《GrADS实用手册》 ♥♥♥ 中文, 较全面, 可用来自学 – GrADS reference card 英文, 绘图命令和函数速查 – GrADS scripting language reference card 英文, 脚本语言语法速查 – GrADS Documentation Index 联机英文索引, 最全面
12
LASG动力论坛版
二. GrADS绘图软件包的安装
• 操作系统
GrADS绘图学习技巧与实例
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
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图形输出类型GRADS (Grid Analysis and Display System)讲义（5）图形输出类型内容提要⼀、图形类型总介⼆、各类图形详细介绍⼀、图形类型总介set gxout graphics_type图形类型⼆、各类图形详细介绍主要介绍格点数据的数据类型，结合具体实例1、line: 单线图（缺省）2、bar: 直⽅图3、linefill: 两单曲线之间填⾊⽤法：set gxout linefilld a;b4、barb: 风向杆形式绘⼆维风场5、contour: ⼆维等值线图6、shaded: ⼆维填⾊图7、errbar: 单线图及误差分布8、fgrid: ⽤指定颜⾊填充⼆维格点场⽤法： set gxout fgridset fgvals Val1 col1 val2 col29、grid: 以⽹格形式在各⽹格点中央标出该点数值10、scatter: 散点图 d u;v11、stat: 输出有关的统计量12、print: 输出ASCII码值12、stream: 流线形式绘⼆维风场13、vector: ⽮量箭头形式绘⼆维风场14、fwrite: 图形不在屏幕上显⽰,⽽是将输出结果存⼊⼀个由“set fwrite ⽂件名”所指定的⽂件中⽤法：’open data.ctl’* 要输出的数据资料⽂件名a.grd‘set fwrite d:\li\a.grd’‘set gxout fwrite’* 设定维数环境‘set lev 500’;‘set t 1 30’‘set lon 120 270’‘set lat –10 10’‘d u’*关闭⽂件a.grd‘disable fwrite’如果资料为NC格式资料，写法如下 ‘reinit’‘sdfopen d:\llp\u.nc’其余同上想要了解nc资料的信息：‘sdfopen d:\llp\u.nc’‘q file’‘set x 1’‘set x 37’注意：（1）当图形输出类型为vector、stream和barb 时, 在显⽰命令后需给出两种数据，两者间⽤分号“;”隔开，前者理解为U分量，后者理解为V分量；（2）当图形输出类型为errbar、scatter、linefill时，也须⽤分号将两种数据分开，如：display u ; vdisplay ave(u,t=1,t=10) ; ave(v,t=1,t=10)（3）对于vector 和 stream类型的图形, 可以给出第三种数据，如：display u;v;hcurl(u,v)display u;v;mag(u,v)⾼级篇输出ASCII格式数据： File=‘uwind.txt’‘set gxout print’‘set prnopts %g 8 1’‘d u’Rc=write(file,result) Rc=close(file)站点资料以后讨论对站点数据：barb: 在各站点绘风向标findstn: 搜索最进的站点(详见描述语⾔部分) model: 以天⽓图形式将天⽓观测各分量填放在站点四周value: 在各站点标值（缺省）wxsym: 绘wx天⽓符号（3）对于wxsym类型的图形, 站点坐标的每⼀个数值都假定为⼀种wx天⽓符号的代码，台风和热带风暴的符号包含在wx天⽓符号集中。
GrADS绘图学习技巧与实例
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
GrADS绘图基础
命令顺序 grads –b|l|p|c open/run set display draw print quit/reinit
显示类型
bar 直方图 line 折线图 errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill 按网格填色(马赛克图) 按网格填色(马赛克图) fgrid 按set fgvals指定值填网格颜色 fgvals指定值填网格颜色 grid 网格填值 findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值 contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 stat 输出统计数据 print 十进制屏幕输出(set prnopts) 十进制屏幕输出(set fwrite 输出二进制文件(set fwrite) 输出二进制文件(set tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt) 在站点上标记天气符号(set
Grads在win/Linux下的下载、安装、 Grads在wi程序：
Windows：grads-1.8sl11Windows：grads-1.8sl11-win32e.exe Linux：grads-bin-1.9b4Linux：grads-bin-1.9b4-linuxRH9.tar.gz grads-1.8sl11grads-1.8sl11-linux.tar.gz data.tar.Z
Linux： Linux： 在.cshrc 文件中添加: 文件中添加: # 执行文件路径：setenv PATH /usr/local/grads/bin:${PATH} 执行文件路径： # 字体和底图文件路径：setenv GADDIR /usr/local/grads/dat 字体和底图文件路径： # 脚本库的路径：setenv GASCRP /usr/local/grads/lib 脚本库的路径：
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绘图软件使用手册
GrADS绘图软件实用手册2002年1月目录第一章 GrADS绘图软件概述1.GrADS绘图软件简介2.Internet上的GrADS资源2.1GrADS在Internet上的主页2.2 windows环境下GrADS资源3.GrADS绘图软件的安装（windows环境）3.1在windows环境下安装GrADS软件包3.2X server 的安装第二章 GrADS绘图模板1.GrADS示例演示1.1 启动GrADS1.2 退出GrADS1.3 示例演示GrADS命令的使用2.GrADS绘图模板3.GrADS模板的高级应用3.1GrADS描述语言3.2GrADS高级模板的应用第三章 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) 是一套应用广泛、使用方便的科学数据绘图软件包。
其主要特点：z GrADS属于自由软件，可以从Internet上免费获得。
z可运行于各种Windows 和Unix工作平台。
z GrADS可用于4D数据的分析。
既经度、纬度、层（气压层、高度层等）和时间/xyzt 4维。
数据可以是格点化的数据或离散点数据。
GrADS特别适用于气象类数据的分析。
但也完全可以用于更广泛类型的数据分析。
z GrADS有多种显示方式：等值线、流线、矢量图、风矢量图、站点填图、折线图、直方图等多种两维图形。
z可处理多种数据格式的数据。
GRIB、NetCDF、HDF-SDS等通用数据格式和系统自定义的一种二进制数据格式。
z采用命令行输入的方式交互式地显示图形。
GrADS实用绘图(基5)
LASG动力论坛版
8
res = functionname(args) ;
* res为函数值
math_abs(expr), math_log(expr), math_log10(expr) math_sin(radian), math_cos(radian), math_tan(radian) math_asin(expr), math_acos(expr), math_atan(expr) math_sinh(radian), math_cosh(radian), math_tanh(radian) math_asinh(expr), math_acosh(expr), math_atanh(expr) math_pow(expr, exponent) math_sqrt(expr), math_exp(exexpr) math_fmod(expr1,expr2), math_mod(expr) math_strlen(string) strlen(string) valnum(string) wrdpos(string,int) sublin(string,n) Subwrd(string,n) substr(string,start,length)
res = write(filename, record, { | append } ) res: 状态值 =0 OK =1 打开错误 =8 文件为read打开 res = close(filename) res: 状态值 =0 OK =1 文件未打开
LASG动力论坛版
11
八. 批处理和脚本语言(S) —控制语句
GrADS绘图学习技巧与实例(阿木)
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！（笔者：阿木）1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
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精细绘图
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绘图软件的使用
气象图形软件应用——G r A D S绘图软件的使用成都信息工程学院大气科学系巩远发主要内容一、GrADS的一般介绍二、GrADS的数据结构详解三、GrADS数据数据处理和绘图详解四、GrADS绘图编程（脚本）语言第一讲 GrADS的一般介绍什么是GrADS相关的关于GrADS的网络资源启动和退出GrADS最基本的GrADS操作（命令）GrADS使用文件的分类GrADS的主要功能介绍1. 什么是GrADS？G r A D S——Grid data Analysis and Display System 格点数据分析和显示系统GrADS（Grid data Analysis and Display System）是一个全32 位的交互操作的气象格点数据和站点数据的分析和显示环境。
系统具有气象数据分析功能强、地图投影坐标丰富、高级编程语言使用容易、图形显示快速，并具有彩色动画功能等特点。
目前是国内外气象数据显示的标准平台之一。
2. 相关的关于GrADS的网络资源WWW Sites/grads/head.html/home.htmlGrADS论坛动力论坛专业绘图软件格点资料下载http://www.ecmwf.int//cdc/reanalysis/3. 启动和退出GrADSgrads 启动GrADSga-> GrADS的系统提示符help给出GrADS的主要命令quit 退出GrADS4. 最基本的GrADS操作Open :打开一个网格点或站点数据文件供GrADS使用(sdfopen: 打开netCDF数据)Query (或仅用q) :查看当前GrADS环境的各种参数Display (或仅用d) :显示(绘制)一幅GrADS “表达式”的图形Set : 设置绘图的操作方式(包括绘什么图，在那里绘图，怎样绘图……)Clear (或仅用c) : 清除当前显示的图形5. GrADS使用文件的分类数据文件: 二进制无格式记录的原始数据，可以是格点的，也可以是站点的。
第二讲GrADS绘图的各种设置
对图形类型为等值线、填色等 值线、矢量图、流线图起作 用的设置：
set strmden value 设置流线密度，值value为1到10，
缺省为5。 例: set strmden 7
set rgb cnum red green blue 设置自定义颜色号 cnum取值范围为16-99. red、green和blue分别表示该颜色号
符串中反斜杠表示起新行。
draw xlab string draw ylab string 分别在水平坐标轴或垂直坐标轴左侧写字符串
string，作为x轴或y轴的说明。
set annot color <thickness> 设置上述图注所用的颜色和线宽，缺省为白色，
线宽为6，该命令同时设置了坐标轴线以及刻 度的颜色和线宽，坐标刻度和标尺的线宽为图 注设置的线宽再减1。
数值 vector: 矢量箭头形式绘二维风场 stream: 流线形式绘二维风场 barb: 风向杆形式绘二维风场 bar: 对一维场不绘单线图，而绘直方图 line: 对一维场绘单线图（缺省） 例：set gxout shaded
另外三种图形类型设置：
set gxout fgrid 用指定颜色填充二维格点场 与命令set fgvals val col合用 例：set gxout fgrid
定等值线的颜色。
set cmin value 不画低于此value的值。 set cmax value 不画高于此value的值。 set black val1 val2 不画值介于val1和val2之间的等值线
。 例：set cint 20 或 set clevs -20 -10 -5 0 5 10 20
set grid on|off|value|horizontal|vertical <style> <color>
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
6. draw title sting
在图形顶部写字符串为标题 。
7. draw wxsym symbol x y size <color<thickness>>
在指定的位置（x、y）画出指定的天 气符号symbol，颜色为color，大小为 size，线宽为thickness。
5.2 设置剪裁区
旋转是以justification定义的点为中心，反时针旋转
justification设置字符串的位置， 描述了与“draw string” 命令中x、y点位置上的相互关系
'set string 2 bl 8' 'set strsiz 0.3 0.3' 'draw string 6 4.58 G '
设置命令如下： set clip xlo xhi ylo yhi 其中，xlo，xhi，ylo，yhi是剪裁区在实页中 的英寸坐标。 用户可以设置一块剪裁区（clipping area ） 用于绘制一些基础图形，设置剪裁区后，这 些基础画图功能只能在剪裁区内进行。
在（x、y）点画一个类型为marktype大小为 size的符号，marktype取值为： 1 - 十字线 2 - 空心圆圈 3 - 实心色圆圈 4 - 空心长方形 5 - 实心长方形
'set gxout stream' 'd u;v' 'draw mark 3 4 4 0.5'
第5章
基础绘图指令
GrADS系统提供了一系列基础绘图 指令来控制和实现一些基础的绘图功能， 比如绘制线条、符号、字符串等。
y string 在x、y点位置上写字符串string， x和y是 虚页上的英寸坐标。
注：x、y点为字符串左下角的位置！ 可以使用“set string” 和 “set strsiz”命令 来设置字符串的性质。
set
strsiz hsiz <vsiz>
设置字符串中字符的大小，hsiz是字符 的宽度，vsiz是字符的高度，单位是虚
页英寸。
set
string color <justification <thickness <rotation>>> 设置字符串string 的性质。 color 设置字符串颜色； justification设置字符串的位置； thickness设置字符串的粗细； rotation设置字符串旋转的角度。
4. draw recf xlo ylo xhi yhi
以（xlo, ylo）、（xhi, ylo）、（xhi, yhi） 和（xlo, yhi）四点为顶点，画一个填色 的长方形。
'set line 2 6 8' 'draw recf 3 3 6 6'
5. draw mark marktype x y size
3. draw rec xlo ylo xhi yhi
以（xlo, ylo）、（xhi, ylo）、（xhi, yhi） 和（xlo, yhi）四点为顶点，画一个不填 色的长方形。
set
line color <style> <thickness> 用于设置线条性质。
'set line 2 6 8' 'draw rec 3 3 6 6'
2. draw line x1 y1 x2 y2
从坐标点x1，y1到x2，y2 画一条直 线，可以使用“set line”命令来设置 线条的性质（单位：英寸）。
set
line color <style> <thickness> 用于设置线条性质。
'set line 2 4 6' 'draw line 3 3 6 6'