GrADS绘图与编程(值得收藏)
GrADS基础绘图
自由图形、符号
设置线属性-> 画线-> 画符号-> 画空心|实心方框-> 画多边形-> 画多边形-> 画天气符号-> 设置字体-> 写标题-> 设置字符串属性-> 设置字符串大小-> 写字符串->
‘set line #color #style #thick’ ‘draw line x1 y1 x2 y2’ ‘draw mark marktype x y size’ ‘draw { rec | recf } xlo ylo xhi yhi’ ‘draw ployf x1 y1 x2 y2 ……xn yn’ ‘draw mappoly lon1 lat1 lon2 lat2 … lon1 lat1’ ‘draw wxsym #symbol x y size #color #thick’ ‘set font [0~5]0’ ‘draw title title_string’ ‘set string #color { tl | tc| tr| l | c| r| bl | bc| br} ‘set strsiz width height’ ‘draw string x y string’
常用绘图软件，及GrADS优缺点
matlab画图和grads画图的比较 总体上来说matlab画图比grads好 matlab的优点： （1），matlab是用对象属性来控制图形中的每一个部分，因此在对图片进行操作是相当容 易。 （2），它画二维曲线图时，画的图非常光滑，视觉效果很好。 （3），matlab画二维等值线图时，投影方式相当丰富，基本覆盖各行业所需。它的底图 （地图）比较漂亮，而且还可以用海拔高度函数是不同海拔用不同的颜色表示出来。 （4），它可以画三维图，图形非常漂亮。grads在这点目前是办不到的. Matlab的缺点：
《GrADS基础绘》课件
GrADS基本语法及操作
GrADS基础语法介 绍
掌握GrADS的基本语法规则 和函数，包括变量定义、数 据读取和运算。
GrADS文件格式及 读取
学习不同类型的文件格式 （如NetCDF）在GrADS中的 读取和处理方法。
GrADS数据处理技 巧
探索利用GrADS进行数据处 理和分析的高级技巧，包括 绘图、统计和过滤。
总结
GrADS的优缺点评价
综合评价GrADS的优点和不足，帮助您了解 该工具是否适合自己的数据可视化需求。
GrADS使用心得交流
分享使用GrADS的一些心得体会和实用技巧， 促进经验交流和学习。
GrADS与其他软件集成
GrADS和Python的 数据交互
探索如何结合GrADS和 Python的强大功能，进行数 据交互和自动化处理。
GrADS与Matlab的 数据联动
了解如何在GrADS和Matlab 之间无缝传递数据，实现更 多分析和可视化功能。
GrADS在GIS领域的 应用
发现如何将GrADS与地理信 息系统（GIS）工具结合，实 现空间数据的分析和可视化。
GrADS案例Leabharlann Baidu析
1
GrADS数据处理应用案例
2
分享一些使用GrADS进行数据处理和
分析的实际案例，帮助您更好地应用
GrADS。
grads中文说明手册
(3) .gs GrADS 控制文件，用 GrADS 命令 run 执行之。这也是一个纯文本文件，内含用描 述语言（script language）写成的批处理 GrADS 系统设置和命令，可集成处理 GrADS 命
版权
本手册版权属中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验 室（LASG）。本手册可免费下载及使用。
本手册的维护与更新
本手册由王军负责维护与更新，如果有什么意见和建议，请与他联系，或在论坛上讨论。 Email : wjun@lasg.iap.ac.cn
有用的链接
GrADS Homepage ： GrADS Script Library：
Technical Report
GrADS使用手册
LASG 编
中国科学院大气物理研究所
大气科学和地球流体力学数值模拟国家重点实验室
（ LASG ）
http://www.lasg.ac.cn
关于本手册
本手册是在 1995 年张洪编写的《GrADS 气象图形系统使用方法》及其它相关文档的基础上 汇编而成。主要介绍了 GrADS 的基本用法、使用技巧、及与 GrADS 相关的绘图技术。本 手册将根据软件发展不定期进行修订。
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 ↵
– 运行GrADS安装程序,建议安装到盘区根目录下
• 如 C:\PcGrADS 或 D:\PcGrADS • 安装后, PcGrADS目录下的子目录主要有
– dat/ – lib/ – win32/ win32e/ 字体, 底图, 要素单位udunits.dat 批处理脚本 执行文件(X-Win32版本) 执行文件(win32e版本)
19
LASG动力论坛版
二. GrADS绘图软件包的安装
• 备份
– dat\ 目录的内容 – lib\ 目录的内容 – udft\ 目录的内容
最具扩展性的内容都保存在这几个目录下, 这也是分目 录存放不同文件的好处.
• 搬家
– 将整个PcGrADS\目录完整拷贝 – 按照前述的环境变量设置方法设置好环境变量
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 sst
open(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件
GrADS绘图基础
Grads网络资源 Grads网络资源
GrADS主页和中文讨论区 GrADS主页和中文讨论区 – http://grads.iges.org/grads/ 或ftp:// grads.iges.org/grads (软件下载/ 联机文档/ email讨论区) 软件下载/ 联机文档/ email讨论区 讨论区) – http://bbs.lasg.ac.cn/cgi-bin/forum/forums.cgi?forum=3 http://bbs.lasg.ac.cn/cgi(中文手册下载/ 中文讨论区/ 脚本和其它资源) 中文手册下载/ 中文讨论区/ 脚本和其它资源) GrADS实用文档 GrADS实用文档 – LASG编《GrADS实用手册》 LASG编 GrADS实用手册 实用手册》
显示图形 ga-> display u gaga-> d u; v.2 gaga-> d u; v.2; mag(u,v.2) ga 保存图形 ga-> enable print gaa.gmf ga-> print gaga-> disable print ga 清屏 ga-> clear gaga-> c ga重置缺省环境, 重置缺省环境, 不关闭 数据文件 ga-> reset ga重新初始化, 重新初始化, 关闭 reinit
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 MATLABIDV (Integrated Data Viewer) Origin surfer
grads-第九章
内部函数
4、read (name) 、
返还码是： 0 - 正常 1 - 打开文件错误 2 - 文件结束 8 - 打开的文件用于写入 9 - 输入、输出错误
内部函数
5、write (name, record <, append>) 、 写函数，当第一次调用’write’函数时，会生 成一个文件名为’name’的文件，并处于打开 的写入状态。如果这个文件早已存在，这个 文件会损坏。如果使用’append’选项，写入 的内容就会接到原来文件的后面。 返还码是： 0 - 正常 1 - 打开文件错误 8 - 打开的文件用于读取。
GrADS绘图与编程 绘图与编程
秦育婧
南京信息工程大学 大气科学学院
第九章 脚本语言及编程
学习目标 学会使用脚本语言编程
第九章 脚本语言及编程
学习要求 掌握流程控制、内部函数等
脚本语言概述
脚本语言（ script language ）是GrADS 软件包中自带的高级语言，可以用于 GrADS的高级操作，更 方便地使用 GrADS的功能。 通过脚本语言可以编写具有复杂功能的 命令集（即gs文件）。
gs文件中的两类型语句 文件中的两类型语句
1）脚本语句：不用引号，需要经过编译的。 2）指令性语句：GrADS命令 命令，需用单引 命令 单引 号括起来。 脚本语言用到指令性语句中要用引号 i=1 ‘set t ‘i’’ (或者‘set t ’i )
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 sst
open(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件
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 启动GrADS
1.2 退出GrADS
1.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
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 sst
open(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件
GrADS精致绘图说讲解
GrADS 精致绘图说
【转】GrADS 精致绘图说(一)
分类:技术文摘
2007.2.13 21:37 作者：步行者| 评论：0 | 阅读：2872
[程序见GrADS Script Library]
1. basemap.gs:
basemap L | O | U
在低分辨率海岸廓线范围内用颜色覆盖陆地／海洋。适用于各种投影方式，需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 ymid
cbarc 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 -p
cbar_line -x X -y Y -c color -m mark -l linestyle -t text -p
GrADS实用绘图--脚本编程
{exec| run} batch.gs run script.gs
使用脚本编程规则
八. 批处理和脚本语言(S) — 变量/数组/运算符
GrADS 脚 本 程 序
LASG动力论坛版
5
变量
定义变量 全局变量 abbrev.#file(dimexpr,dimexpr,...) 逻辑变量 abbrev.#file(stid=‘...’) 定义变量 ‘define va = vwnd 2 ’ 预定义变量 变量类型 预定义变量 lat, lon, lev, time (向量/矩阵变量) 变量引用方式 t=1 _global = 1 TRUE = 1; FALSE = 0 result, rc “随意/不随意” ♥♥♥
绝对值, 对数 三角函数 反三角函数 双曲函数 反双曲函数 幂函数 平方根, e的幂 最接近整数, 截去小数 余数
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(expr) math_nint(expr), math_int(expr) 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)
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 sst
open(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件
GrADS画图入门
命令：open filename2.ctl
ga-> open strike024_041720 Scanning description file: strike024_041720 Data file strike024_041720.dat is open as file 2
查询命令
命令：query files
• 3、图形显示类型及绘图设置：
（1）一维：直方、折现、散点、误差图....（excel 2003） （2）二维：标量场——等值线、阴影图 矢量场——风矢量、流线图 站点—— （3） draw line/string/title/mark……
标量场分析
1. 等值线
2. 阴影图 1. 等值线
设置 set gxout contour
设置 set gxout shaded
（缺省）
命令：
display express
ga-> d strp Contouring: 0 to 100 interval 10
设置：
set *** ###
维数参数设置 等值线设置
Set set Set set lon 100 cint 15 120 指定等值线间隔 lat 103 25 指定等值线颜色 ccolor
命令：query file n
ga-> q file 2 File 2 : TC0801-StrkPro024 2008041720 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat Type = Gridded Xsize = 101 Ysize = 81 Zsize = 2 Tsize = 1 Number of Variables = 2 strp 0 0 strike probability strn 0 0 strike number
GrADS实用绘图--站点图 函数使用
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
9
LASG动力论坛版
六. 图形显示(D)2
(x,y)的最近站点 set gxout findstn d stnexpr; x; y 注意:
1. 维数只能X,Y变化 2. x, y 为虚页坐标
— 站点绘图3
风羽序列图 set gxout tserbarb d uexpr; vexpr 注意:
1. 只能时间维变化 2. 变量注明站号stid
20
LASG动力论坛版
七. 数据分析(A)
本节内容提要
• 数学函数 • 统计函数 • 网格处理函数 • 滤波平滑函数 • 差分函数 • 物理量计算函数 • 站点数据专用函数
— 函数的使用
21
LASG动力论坛版
七. 数据分析(A)
功 能
abs(expr) cos(expr) sin(expr) tan(expr) exp(expr) log(expr) pow(expr1, expr2)
GrADS学习资料：第7章 GrADS 描述语言及编程(1)
描述语言函数总是得到一个单一的字符串 结果，但可以有一个或多个字符串自变量。
它的调用格式为： name (arg, arg, arg, ……)
如果函数没有参数，仍必须写出一对空括 号，例如：name（）。
通过函数子程序的定义语句在你的描述语 言程序中构造自己的函数： function name(variable, variable, ...)
如果字符串‘string’太短，结果是一个空字符串。 ‘word’必须是整数。
例如： x=subwrd(how are you,3) 结果是将字符串“how are you”中的第3个单词 “you”作为新的字符串赋值给变量x。
sublin (string, line) 功能：这个函数可以从字符串’string’中截取第
变量“rc”被赋予 GrADS 命令执行的返回码（总是一个整型数值）
'open e:\grads\1.ctl' 'set lat 0 90' 'set lon -180 0' 'set lev 500' 'enable print e:\htry.gmf' i=1 while(i<=5) 'set t 'i 'd z' 'print' 'c' i=i+1 endwhile 'disable print'
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
文本格式数据源
Open(9,file=„data.grd‟,form=„unformatted,access=„direct‟,recl=ii*jj*4)
Irec=1 Do 200 iy=1,4 Do 100 m=1,12 Read(1,1000) Read(1,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j),I=1,ii),j=1,jj) Irec=irec+1
以下详细说明数据描述文件中各记录的含义： 1、DSET data-set-name 给定二进制原始数据文件的文件名（可包含路径），若数据 文件与描述文件在同一路径下，可用省缺路径符号“^”代表。
23
2、title string 用字符串string简略描述数据文件的内容。 3、undef value 定义缺测值，GrADS在运算操作时和图形操作时将忽略这些格点。 4、options <keywords> <keywords>可以是 <yrev><zrev><sequential><byteswapped><big_endian><little_endian ><template> 其中：sequential表示数据是以顺序无格式形式存放，每个记录为一 个(x,y)场;yrev表示 y维数方向与 ydef中说明的方向相反，即为由 北到南； zrev表示 z维数方向与 zdef中说明的方向相反，即为由 上到下；byteswapped表示二进制数据的位存放顺序取反序； big_endian和 little_endian不同机器之间二进制位存放顺序的自动 改变； template表示 多个时间序列原始数据文件共用一个数据描 述文件。例如：
7
简
介
• GrADS 是Grid Analysis and Display System的缩写,它由美国马里兰大学气象 系Brian E. Doty开发。利用该系统可实现 包括格点数据和站点数据的彩色气象图 形分析和显示。具有操作简便、功能丰 富、图形美观、显示快速的特点。在国 内外气象界得到广泛应用。 • 版本：V1.5, V1.6.8, V1.7.7，V1.7.9等
16
二、文本格式数据
GrADS格式
现有如下资料： 名称： U850，V850，U200，V200，H500 和TSFC 范围： 60~150E，0~40N 分辨率：2.5×2.5 时间： 1982年1月~1985年12月的逐月资料。
如何将上述资料写成GrADS下的数据格式？
17
顺序：(x,y)
U850
无格式，直接记录
Read(3,1000) Read(3,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(2,1000) Read(2,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(4,1000) Read(4,2000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(5,1000) Read(5,3000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1 Read(6,1000) Read(6,4000) ((var(I,j),I=1,ii),j=1,jj) Write(9,rec=irec) ((var(I,j)I=1,ii),j=1,jj) Irec=irec+1
24
有1may92.dat,2may92.dat,…..等数据文件，需用一个共同的描述文件， 这时的CTL文件可写为： Dset %d1%mc%y2.dat ….. Options template …… Tdef 72 linear 0z1may92 1 dy 其中：%y2:2位数年；%y4:4位数年；%m1:1或2位数月；%m2:2位数 月；%mc:3个字符的月份英文缩写;%d1：1或2位数天；%d2：2位 数天；%h1：1或2位数小时；%h2：2位数小时 5、xdef number linear start inc 或 xdef number levels <value-list> 设置网格点与经度的对应关系。其中number是x方向网格点数， linear或 levels表明网格点映射类型， start起始经度， inc格距大小， value-list表示X方向各格点的列表。
所以，固定的维数环境和变化的维数环境相结合，就构成 了当前的维数环境。
14
(1)x,y,z,t均固定时，得到一个单值数据点 (2)x,y,z,t中只有一维变化时，得到一维曲线 (3)二维发生变化时，对应于二维剖面图， 如X-Y平面图，X-Z，Y-Z，X-T，Y-T， Z-T剖面图 (4)当三维发生变化时，GrADS以动画方式 显示二维切片。
8
安
Fra Baidu bibliotek
装
1、将系统软件拷入硬盘某一子目录 下，如：c:\grads
2、修改autoexec.bat文件，加上 path c:\grads;%path% set gaddir=c:\grads set gascrp=c:\grads
9
GrADS 软件包的组成
• • • • • • • • Grads.exe Dos4gw.exe Gribmap.exe Gribscan.exe Gxtran.exe Stnmap.exe Wgrib.exe Gv.exe
x (或lon) y (或lat) z (或lev) t (或time)
从西向东的水平坐标 从南到北的水平坐标 从地面到高空的垂直坐标 时序坐标。
如 ： set lon 40 160 ；set lat 0 60 定义了水平变化范围 set lat 30；set lon 0 180 定义了沿30N的纬向变化范围 set time jan81 dec94 定义了从81年1月到94年12月的时段
关键词：如dset, undef, title,xdef ，ydef 等
被描述的数据文件名
缺测标记 标题
起始值
确定维数环境 间 隔 指定变量名 层次数
22
数据描述文件一般包含以下几项： （1）被描述的数据文件名（DSET） （2）该数据说明文件的标题（TITLE） （3）数据类型、格式和选项（DTYPE ，FORMAT，OPTION） （4）时间、空间维数环境设置（XDEF，YDEF，ZDEF，TDEF） （5）变量定义（VARS，ENDVARS）
11 8.5
11
所需预备知识
DOS和Windows基础 FORTRAN编程 WORD基础
学习要求
强调上机编程实践
12
第二章 GrADS 数据格式
• 目标 (1)能将文本格式数据转为GrADS格式
(2)能写出数据描述文件(*.ctl)
重要概念
一、维数环境的概念
1、含义：GrADS视每一个变量(VAR)场为一个四维数据集， 即包括三维空间(x,y,z)和一维时间(t)。 2、作用：说明和指定随后的分析或图形操作时参加的原 始数据集的维数范围。 3、定义方法： Set lat|lon|lev|time val1 <val2> 地球坐标 Set x|y|z|t val1 <val2> 格点坐标 13
系统的核心文件,由此进入GrADS环境
由GrADS.exe文件调用 产生格点资料映射文件(*.idx)
看Grib码资料，转为TXT格式文件
显示*.gmf格式的图形文件 产生台站资料的映射文件(*.map)
Grib码资料的解读程序
将*.gmf格式的图形文件转换为*.wmf 格式
10
启动和退出GrADS
GrADS绘图与编程
图例1 1月份500hPa位势高度场
2
图例2 500hPa高度场（1、4、7、10月）
3
图例3 亚洲季风区7月850hPa流场图
4
图例4 SST和Sea Level分布图
5
图例5（130E，25N）处的感热和潜热的时间演变
6
第一章
★ ★ ★ ★ ★ ★
概述
简介 GrADS的安装 GrADS软件的组成 GrADS的启动和退出 所需预备知识 学习本课程的要求
z
(X,Y) (X,Y)
v
t
• T=1 • T=2 • T=3
U V H T
U200
分(z)
时(v)
5月20日(t)
秒(x,y)
18
FORTRAN源程序
Parameter(ii=37,jj=17) Real var(ii,jj) Open(1,file=„u850.dat‟) Open(2,file=„v850.dat‟) Open(3,file=„u200.dat‟) Open(4,file=„v850.dat‟) Open(5,file=„h500.dat‟) Open(6,file=„tsfc.dat‟)
25
6、ydef number mapping start inc 或 ydef number levels values-list 设置Y方向格点与纬度的映射关系。其中number为Y方向的格点数， mapping表示映射方式，有： Linear:线性映射 Gausr15：高斯R15纬度，Gausr20,gausr30,gausr40等 对于线性映射linear，start 为起始纬度，inc为Y方向格距。 对于高斯映射，start为第一高斯网格数。 对于levels映射， value-list为Y方向取值表 7、zdef number linear start inc 或 zdef number levels value-list 设置气压面与垂直网格点的映射关系。如： Zdef 10 linear 1000 –100 Zdef 10 levels 1000,925,850,700,600,500,400,300,200,100
15
说出下列维数环境的含义
Set lon 60 150 Set lat 0 30 Set t 7 Set lev 200
Set lat 30 Set lon 60 130 Set t 1 12 Set z 1
Set lon 120 Set lat 30 Set z 1 12 Set t 1
Set lon 60 150 Set lat 0 30 Set t 1 12 Set lev 200
100 continue 200 continue 1000 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) end
三、数据描述文件
Data.ctl 文件内容：
Dset data.grd Undef –9.99e+33 Title Ncep/Ncar reanalysis project Xdef 37 linear 60 2.5 Ydef 17 linear 0 2.5 Zdef 2 levels 850 200 Tdef 48 linear jan1982 1mo Vars 4 U 2 99 u wind (m/s) V 2 99 v wind(m/s) H 1 99 H500 T 1 99 Tsfc data endvars
在DOS下键入： 既可进入；键入: 即可退出grads,返回DOS。 grads quit
GrADS在初始化绘图环境之前，将提示用户选择风景画(Landscape)或肖像画 (Portrait)形式。风景画形式的大小为11×8.5英寸，肖像画为 8.5 ×11英寸。 启动GrADS系统时，有以下选择： 8.5 -b 以批处理形式运行GrADS。 -l 以风景画形式运行GrADS。 11 -p 以肖像画形式运行GrADS。 -c 进入GrADS后首先执行随后跟在RUN命令后的描述文件，如： grads –c “run profile.gs” 上述选项可结合起来使用，如： grads –blc “run batch.gs”