气象常用绘图软件grads剖析共62页文档

16
二、文本格式数据
GrADS格式
现有如下资料： 名称： U850，V850，U200，V200，H500 和TSFC 范围： 60~150E，0~40N 分辨率：2.5×2.5 时间： 1982年1月~1985年12月的逐月资料。
如何将上述资料写成GrADS下的数据格式？
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月的时段
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‟)
7
简
介
• GrADS 是Grid Analysis and Display System的缩写,它由美国马里兰大学气象 系Brian E. Doty开发。利用该系统可实现 包括格点数据和站点数据的彩色气象图 形分析和显示。具有操作简便、功能丰 富、图形美观、显示快速的特点。在国 内外气象界得到广泛应用。 • 版本：V1.5, V1.6.8, V1.7.7，V1.7.9等

设置y轴标注的颜色<粗细<大小>>
x轴标注位置，offset：英寸；side=b,下|t,上
y轴标注位置，offset：英寸；side=r,右|l,左
数据中z维用对数坐标
设置坐标轴框的颜色<粗细>
y轴数据变化范围从vlo到vhi
x轴数据变化范围从vlo到vhi
把断点数据连上/不连上
在图上方居中写标题string
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息
fwrite
返回当前fwrite命领输出文件的名称
string s
给出经度－纬度坐标值 lon lat, 返回xy坐标值
pp2xy ppx ppy
给出页面坐标 ppx ppy，返回xy坐标值
set imprun script-name
在每个display命令前自动执行一遍script-name文件中的命令组
run file-name <params> file-name <params>
set zlog on|off
设置z坐标取对数坐标,或关闭此设制
set strmden value
设置流线密度[1-10],缺省为5
set stnopts <dig3> <nodig3>
以后3位数标数值
set mdlopts noblank|blank|dig3|nodig3

维数环境是 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

文章编号:1007-6190(2004)04-0035-02用G rADS制作ncep再分析资料气象要素图杨兆礼(中国气象局广州热带海洋气象研究所,广州　510080)中图分类号:P409 文献标识码:B G rads气象图形系统是目前国内外气象界通用的标准图形环境之一。

用它来画气象要素图很方便。

本文简单介绍了如何使用G rads气象图形系统来对NCEP/NC AR再分析资料进行简单的绘图处理。

1　NCEP/NC AR再分析资料简介NCEP/NC AR再分析资料是美国国家环境预报中心(NCEP)和国家大气研究中心(NC AR)采用美国全球气候资料同化分析系统进行四维同化分析后获得的全球同化资料。

其全部历史资料来源于美国国家气候资料中心(NC DC)。

使用ncdum p命令,可以看到高空NCEP/NC AR 再分析资料的属性如表1所示。

每天共有4个时次,有73个纬度和144个经度,等压面共有8～17层(hPa),分别是:1000,925,850,700,600,500,400, 300,250,200,150,100,70,50,30,20,10hPa。

表1　NCEP高空等压面再分析资料说明(以1991年为例)文件名要素名称单位层次Air.1991.nc气温K°17Hgt1991.nc位势高度m17rhum1991.nc相对湿度%8shum1991.nc比湿kg/kg8omega1991.nc垂直速度pascal/s12uwnd1991.nc风—U分量m/s17vwnd1991.nc风—V分量m/s17 地面NCEP/NC AR再分析资料共有94个纬度和192个经度。

包括了气温、云强迫净长波通量、对流降水率、晴空长波射入通量、射入长波辐射通量、射入日辐射通量、地表热通量、冰盖判示、可能蒸发率、降水率、气压、水径流、风—U分量(10m 处)、风—V分量(10m处)等共42种要素。

2　选择绘图要素NCEP/NC AR再分析资料的文件名中,第1个单词便是我们绘图时要用到的变量名。

收稿日期:2002-11-1作者简介:高文良,男,31岁,主要从事短期气候预测等研究工作。

气象绘图软件GrADS 高级使用技巧高文良1　刘晓燕2　曾小东3(11成都高原气象研究所　成都　610072;　21阿坝州金川县气象局　62410031阿坝州马尔康县气象局　624000)摘　要:本文通过介绍气象绘图软件GrADS 使用中的一些高级技巧,分析了GrADS Script 语言中的难理解之处,并介绍了在GrADS 中做合成分析和t 检验的程序,可以对学习和使用GrADS 的科研人员起参考、帮助作用,达到事半功倍的效果,促进GrADS 软件的应用更广泛、深入,提高科研人员的工作效率。

关键词:GrADS ;技巧;程序中图分类号:TP391文献标识码:C文章编号:1003-7187(2002)04-0057-031　引言气象绘图软件GrADS (Grid Analysis and Dis 2play System )是免费共享软件,可随时从互联网上下载(http :///grads/),后续版本正陆续推出[1]。

GrADS 有丰富的内部函数,可以对数据进行计算和分析处理。

它支持处理格点资料和站点资料,并且支持对GriB 码文件、特殊格式文件(如一字节整型、二字节整型、大中型机器二进制数据等)的直接读取,气象科研领域应用非常广泛[2]。

在其最新1.8SL 9版本中,GrADS 又将应用领域推进到了海洋学科,功能也得到了进一步地增强和扩展。

但在使用GrADS 的过程中,特别是使用GrADS Script 语言编程当中,有一些问题需要特别注意。

因为GrADS Script 语言是一种类似于VB Script (或MA TL AB Script )的高级语言,稍不注意就容易出现错误,且查错十分困难,这点与其他语言有较大的差别。

2　站点数据处理GrADS 中站点数据处理基于格点数据的基础之上,需先将站点数据通过Cressman 客观分析方法内插至格点上,然后再依照格点资料的处理方式对站点资料数据进行分析和处理。

Descriptor: strike024_041714.ctl Binary: strike024_041714.dat File 2 : TC0801-StrkPro20 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat
命令：query file n
当指定颜色少于等值线条数时，最后一种颜色适用于超出的所有等值线。
等值线标注设置
set clab %.1f 等值线标记方式
set clab on（等值线标数值）/off（不标记）/forded（强制标记）/%gK（在 标记后加上字符K）/%g%%（加上字符%）/%.2f(保留2位小数)/ %03.1f（整数 部分保留3为，若不足3位，前部用0补齐）
GrADS 画图入门
Grid Analysis and Display System
国家气象中心 2008年12月7日
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
strp 0 0 strike probability strn 0 0 strike number
查询命令
命令：query ctlinfo
ga-> q ctlinfo dset strike024_041714.dat title TC0801-StrkPro14 undef 9999 xdef 101 linear 100 0.5 ydef 81 linear 0 0.5 zdef 2 levels 1000 500 tdef 1 linear 00Z19JUL2005 60mn vars 2 strp 0 0 strike probability strn 0 0 strike number

GrADS下任意方向垂直剖面的实现缪明余清平（解放军理工大学气象学院）廖德敏（93966部队气象室）引言GrADS（Grid Analysis and Display System）是当前气象学中应用最为广泛的图形图像处理软件之一，它功能强大，使用方便而深受用户的欢迎。

利用它可直接实现图像的动画显示，在三维数值模拟研究中，在固定模式高度的情况下，可较为方便地显示某一模式预报量在该层次上的水平分布；在固定经（纬）度的情况下，还可方便地显示某一模式预报量沿该经（纬）度的经（纬）向垂直剖面分布。

但在实际问题的三维数值模拟研究中，用户关心的不仅仅是模式预报量沿某一经纬向的垂直剖面分布，更多的时候用户关心的是沿某一方向的垂直剖面分布，但在GrADS下如何实现这一功能，这可能是困扰一般的GrADS用户的问题。

其实，GrADS除了强大的图形图象的分析处理能力外，还有较强的数据处理能力，我们可利用GrADS的交互性和数据处理的特点，利用GrADS脚本语言的编程技术，按GrADS 的数据格式，通过提取所要分析的模式预报量沿某一方向上的值（ascII码形式），然后再将提取ascII码值转换成GrADS支持的二进制格式，并为转化后的数据编写相应的数据描述文件（ctl文件），在GrADS下打开它就可实现某一模式预报量沿该方向的垂直剖面分布的分析。

一、实现的思想和步骤在介绍实现的思想和gs编程（GrADS Script programing）之前先说明GrADS下某些基本的常识和某些基本命令的含义和功能。

1．GrADS基本常识和基本命令的简介（1）进入GrADS后，在“ga ->”后每执行一个GrADS命令，在文本窗口都会返回1行或2行的的文本代码，返回的代码即为GrADS 的默认变量“result”值，每执行一次GrADS命令，原来的“result”值将被刷新成刚执行的GrADS命令返回码。

在gs编程中，可通过代码“say result ”看到默认变量“result”对应的返回文本代码。

#$%& !’ ($& !
#$%&%#’(&!)* %% + ,-) .
/ 资料文件的数量 0 #$%&%#’(&!!* 12%3( + 4 . 1%’%12%3( 5 6 17 8 9:;(2%3( 8 <;=> 1%’ 6 ? 6 17 8 9:;(2%3( 8 <<’> 1%’ 6 ? @ 6 17 8 9:;(2%3( 8 A4-*> 1%’ 6 ? 6 17 8 9:;(2%3( 8 A-**> 1%’ 6 ? @ 6 17 8 9:;(2%3( 8 A!**> 1%’ 6 ? 6 17 8 9:;(2%3( 8 B4-*> 1%’ 6 ? @ 6 17 8 9:;(2%3( 8 B-**> 1%’ 6 ? 6 17 8 9:;(2%3( 8 B!**> 1%’ 6 5 / 须描述资料文件的列举 0 1C ,* 2 D ,? 4 C=(2 + ,* E 2? 9:;( D 12%3( + 2 . . ,* #C2’:2A( / 依次打开各个资料文件 0 C=(2 + !*? 9:;( D 6 17 8 9:;(2%3( 8 %;;1%’> F&1 6 ? 9C&3 D 6 G:2%&H 6 ? @%##(<< D 6 1:&(#’ 6 ? &(#; D 2I! 2H!J . / 把所有须描述的资料文件数据同时放入 , 个文件 %;;1%’ 中， 以便只须 , 次转换后便可以随意调取 0 :2A3G(& D * 1C ,**K D ,? ,-) 1C !* 2 D ,? 4 &(%1 + ,* E 2， -* . %% + K . &(%1 + ,* E 2? J* . + + F&:1 + :? L . ? : D ,? 2I . ? L D ,? 2H . :2A3G(& D :2A3G(& E , M&:’( + !*? &(# D :2A3G(& . + + F&:1 + :? L . ? : D ,? 2I . ? L D ,? 2H . !* #C2’:2A( ,** #C2’:2A( / 本程序按经纬度顺序先读写第一时刻的 <;=、 再读写第 <<’、 A4-*、 A-**、 A!**、 B4-*、 B-**、 B!**， 二时刻的资料， 依次类推 > > > 0 J* 9C&3%’ + ,-9J> * . -* 9C&3%’ + %)* . (21 !> ! 数据描述文件的编写 数据描述文件 1%’> #’; 的编写格式为： 1<(’17 8 9:;(2%3( 8 %;;1%’> F&1 / 1<(’ 是给出所描述文件的文件名 0 ’:’;(M(%’$(&1%’% A21(9 N "> ""(OO I1(9J- ;:2(%&J* !> H1(9,- ;:2(%& N ,*> * !> Q1(9J ;(B(;<,*** 4-* -** !** ’1(9 ,-) ;:2(%& L%2,"4! ,3C B%&<J / 共 <;=、 ’、 A、 B J 个变量 0 <;= * * / <(%;(B( ;=&(<<A&( 0 ’ * * / ’(3=< 0 ・!"・ / 该数据描述文件的标题 0 / 缺省记录的标记 0 / I 方向共 J- 个格点， 起始纬度为 J*P， 步长为 !> -P 0 / H 方向共 ,- 个格点， 起始经度为 N ,*P， 步长为 !> -P 0 / Q 方向分为 J 层， 即 , ***、 4-*、 -**、 !** R=S 0 / 时间共 ,-) 个月， 起始时间为 ,"4! 年 , 月， 步长为 , 个月 0

GrADS下任意方向垂直剖面的实现缪明余清平（解放军理工大学气象学院）廖德敏（93966部队气象室）引言GrADS（Grid Analysis and Display System）是当前气象学中应用最为广泛的图形图像处理软件之一，它功能强大，使用方便而深受用户的欢迎。

利用它可直接实现图像的动画显示，在三维数值模拟研究中，在固定模式高度的情况下，可较为方便地显示某一模式预报量在该层次上的水平分布；在固定经（纬）度的情况下，还可方便地显示某一模式预报量沿该经（纬）度的经（纬）向垂直剖面分布。

但在实际问题的三维数值模拟研究中，用户关心的不仅仅是模式预报量沿某一经纬向的垂直剖面分布，更多的时候用户关心的是沿某一方向的垂直剖面分布，但在GrADS下如何实现这一功能，这可能是困扰一般的GrADS用户的问题。

其实，GrADS除了强大的图形图象的分析处理能力外，还有较强的数据处理能力，我们可利用GrADS的交互性和数据处理的特点，利用GrADS脚本语言的编程技术，按GrADS 的数据格式，通过提取所要分析的模式预报量沿某一方向上的值（ascII码形式），然后再将提取ascII码值转换成GrADS支持的二进制格式，并为转化后的数据编写相应的数据描述文件（ctl文件），在GrADS下打开它就可实现某一模式预报量沿该方向的垂直剖面分布的分析。

一、实现的思想和步骤在介绍实现的思想和gs编程（GrADS Script programing）之前先说明GrADS下某些基本的常识和某些基本命令的含义和功能。

1．GrADS基本常识和基本命令的简介（1）进入GrADS后，在“ga ->”后每执行一个GrADS命令，在文本窗口都会返回1行或2行的的文本代码，返回的代码即为GrADS 的默认变量“result”值，每执行一次GrADS命令，原来的“result”值将被刷新成刚执行的GrADS命令返回码。

在gs编程中，可通过代码“say result ”看到默认变量“result”对应的返回文本代码。

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 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：用户可在图中点击给定图例的放置位置。

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采用命令行输入的方式交互式地显示图形。

gxtran
Displays metafiles
stnmap
Maps station data
wgrib
See /wgrib.html
GrADS 英文主页：/grads/
§1.2 基本概念
按磁盘文件记录格式，GrADS 软件包的文件分为如下几类：
gxps，gxpsc，gxpscw
将 GrADS 内 print 命令产生的图元文件分别转换为单色（灰度）、彩色以及反白（即白底）
彩色的 postscript 格式文件，用于硬拷贝输出。语法为：
Technical Report
GrADS究所
大气科学和地球流体力学数值模拟国家重点实验室
（ LASG ）

关于本手册
本手册是在 1995 年张洪编写的《GrADS 气象图形系统使用方法》及其它相关文档的基础上 汇编而成。主要介绍了 GrADS 的基本用法、使用技巧、及与 GrADS 相关的绘图技术。本 手册将根据软件发展不定期进行修订。
1

令。详见§6。 (4) .exe 系统命令文件，是 GrADS 系统在 DOS 环境下的各执行文件。如 grads.exe 为 GrADS
图形分析和显示命令；dos4gw.exe 为 DOS 的扩展环境；gxps.exe、gxpsc.exe、gxpscw.exe 都是图元文件转换为 postscript 文件的执行文件；gxtran.exe 是图元文件转换到显示器上 显示的执行文件；gx.exe 是将图元文件转换为各种不带 ps 解释器的打印机输出的执行文 件。详见§1.3。 (5) .gmf（.met） GrADS 系统图元输出文件，格式由 GrADS 内定，文件名随用户自定，
(2) .ctl 原始数据描述文件。该文件为纯 ASCII 文件，用以描述原始数据集的基本信息，包 括数据集文件名、数据类型、数据结构、变量描述等等，详见§4。在 GrADS 环境中至 少得首先打开（open）一个数据描述文件，以便后续的操作有数据对象。

GrADS软件的优点
1）软件的专业性强，它是专门为从事气象研究 的科研和业务人员开发的，因此特别适用于气象 数据资料的分析和显示；
2）软件适用于多种操作系统，如Windows7， Linux和工作站下的Unix等； 3）软件采用命令输入和图形显示交互式的操作 方式
GrADS操作界面由命令窗口和图形显示窗口组成， 在命令窗输入命令，则在图形窗出现相应图形，这种 操作方式简单直观，容易理解和使用。
GrADS主要功能
1）可以处理四维数据，包含经度、纬度、层次和 时间； 2）可以处理多种数据格式，包括二进制数据格式 ，GRIB格式，NetCDF格式等通用数据格式； 3）可以对数据进行再加工，如fwrite语句； 4）可以显示多种类型图形； 5）可以精细设置各种绘图要素； 6）可以按多种方式存储图形文件； 7）提供多种后期处理的软件工具。
授课方式
本课程一半学时为老师讲授，另一半学 时为学生上机实习，老师在一旁指导。
每周上课边讲授边上机实习，有利于学 生及时消化掌握老师讲授的知识，使教 学效果达到最好。
总评成绩
平时成绩：40%
包括考勤、课堂作业、课后作业等
期末成绩：60%
机考
课程内容
第一章 气象绘图软件GrADS基础知识 第二章 GrADS绘图的各种设置 第三章 变量和函数的定义及功能 第四章 站点数据格式及绘制站点图 第五章 气象数据的处理 第六章 其他气象常用绘图软件简介（选讲）
多时刻多个层次多个变量的数据存放次序
数据文件的读取和写入
Readuwnd.f90
数据描述文件（即.ctl文件）（例）：
DSET ua.dat
TITLE Upper Air Data
DTYPE grid

气象图形软件应用——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环境
由GrADS.exe文件调用 产生格点资料映射文件(*.idx)
看Grib码资料，转为TXT格式文件
显示*.gmf格式的图形文件 产生台站资料的映射文件(*.map)
Grib码资料的解读程序
将*.gmf格式的图形文件转换为*.wmf 格式
10
启动和退出GrADS
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
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

●计算机与气象仪器　●收稿日期:2002-02-20气象绘图软件GrADS 使用要点高文良　陈忠明　闵文彬(四川省气象科学研究所　成都　610071)摘要:通过介绍GrADS 绘图软件,分析了其优点及难点,指出了GrADS 的使用要点。

关键词:GrADS ;绘图软件;要点中图分类号:TP30214文献标识码:B文章编号:1003-7187(2002)02-0049-021　引言在众多的绘图软件中,针对气象部门的绘图软件GrADS (Grid Analysis and Display System )是一款非常不错的免费共享软件。

它整个系统由C 语言编成,运行速度快,图形显示快速、美观大方,支持多种图形类型的输出。

GrADS 具有丰富的内部函数,可以对资料进行初步的计算和分析处理,而且操作十分简便,易于理解掌握。

它既可以处理格点资料,又可以处理站点资料。

它支持对GriB 码文件、NetCDF 文件的直接读取,并且对气象界最流行的各种数据格式都能正确分析和处理。

用户还可以利用GrADS 的描述语言(GrADS Script )很方便地进行编程、开发。

GrADS 推出后在气象界得到了广泛的使用。

2　基本概念在GrADS 中,每个原始资料数据集都需要一个数据描述文件文件与它相对应。

在数据描述文件中,包括了对原始资料数据文件的名字、数据格式、纬向格点数及间隔、经向格点数及间隔、空中层次数及间隔、时次及间隔、变量数等的说明。

数据描述文件从时间和空间上对数据集的维数作出了明确的规定。

这样,在GrADS 启动后将数据描述文件打开,就可以对整个资料数据集进行处理和分析。

所以,数据描述文件一定要准确,避免任何差错。

数据文件中原始数据的排放次序也有一定的规则。

其排放次序为首先是经度、纬度,然后是高度和物理变量,最后是时次变化。

[1]。

如果改变了数据的排放规则,必须在数据描述文件中进行特别说明。

因此可以将原始数据序列看成是一个五维(或少于五维)的数据场,将每个经向和纬向的数据面看成是一个记录,是整个数据集的一个二维切片。

参考资料



1、GrADS Tutorial 2、GrADS Reference Card 3、LASG: GrADS Mannual 4、LASG: GrADS 实用绘图讲义，funny@动 力论坛 5、动力论坛GrADS中文讨论区
常见软件


运行GrADS




打开数据： ga->open *.ctl ga->sdfopen *.nc ga->xdfopen test.ddf 查询信息：query ga->query ctlinfo ga->q dims ga->q file 关闭数据文件 ga-> close 3 ga-> close 2 Ga-> close 1 设置维数环境 ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵ 26 LASG动力论坛版 三. 运行GrADS绘图软件包 • 设置维数环境 ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵

set annot #color(1) thick(6) set frame on | off | circle set xyrev on set zlog { on | off | swap | undefine } set coslat {on | off } set tlsupp { year | month } set grid {on| off| horizontal| vertical } linestyle #color set gridln {auto | off | #color } set xlint int set ylint set xaxis starrt end int set yaxis set vrange vlo vhi set vrange2 set xlpos offset(0) { b | t } set ylpos…{l|r} set xlopts #color(1) thick(4) size(0.12) set clopts #col(-1) thick(-1) size(0.09) set ylopts set ylevs set xlevs x1 x2 x3… set ylabs set xlabs b1|b2|b3|… set ylab set xflip { on | off } set yflip set xlab { on | off | aut draw ylab draw title title_string