气象图形图象软件

二 图形要素设置
1. 对图形类型为等值线（gxout=contour）起作用的设置： set ccolor color [设置等值线颜色]，颜色号color为0黑、1白、2红、3绿、4兰、5青
6紫、7黄、8桔黄、9灰。clear或display即重新设定颜色 set ccolor rainbow [设置等值线用七色序列表示] set ccolor revrain [设置等值线用反序的七色序列表示] set cstyle style [设置等值线线形， style：1实线、2长虚线、3短虚线、4长短虚线、5点线， clear或display即重新设定] set cthick thckns [设置等值线线宽， 取值1-10之间的整数，缺省为4] set clab on/off/forced/string/auto [设置等值线标记方式， 一直持续到重新设置]
*.gs文件简单范例： 'reinit' 'open h:\2\model.ctl' 'set lat 20 40' 'set lon 90 120' 'set z 5' 'set t 1' 'd t'
运行*.gs文件：直接在ga-> 输入run ＊.gs文件(含路径)，run 可 省略
'reinit‘
‘set t 1‘
*[设置时间]
‘d z’
*[显示图形命令：d（或display）]
‘c’
*[清屏（或clear）]
......
‘d u;v‘
*[显示图形命令：d（或display）]
‘print’
*[将当前显示的图形保存到输出文件中）]
‘disable print’
*[关闭该输出文件]
根据需要，GrADS可画以下不同类型的图
DSET ^model.dat OPTIONS little_endian UNDEF -2.56E33 TITLE 5 Days of Sample Model Output XDEF 72 LINEAR 0.0 5.0 YDEF 46 LINEAR -90.0 4.0 ZDEF 7 LEVELS 1000 850 700 500 300 200 100 TDEF 5 LINEAR 02JAN1987 1DY VARS 8 PS 0 99 Surface Pressure U 7 99 U Winds V 7 99 V Winds Z 7 99 Geopotential Heights T 7 99 Temperature Q 5 99 Specific Humidity TS 0 99 Surface Temperature P 0 99 Precipitation ENDVARS
§1．运行环境和基本操作
DSET data-set-name（给定二进制原始数据文件的文件名（包含路径））
DTYPE （数据类型）(station / grib) （可省略）
OPTIONS <keywords> （可省略） <keywords>选项： <yrev><zrev><sequential><byteswapped><big_endian><little_endian><template> yrev: y维数方向反向 zrev: z维数方向反向 template：用于多个时间序列原始文件仅用一个描述 ......
PS 0 99 Surface Pressure
Model.ctl ENDVARS:表示数据描述文件结束
§1．运行环境和基本操作
(3).gs GrADS控制文件，用GrADS命令 run执行。这也是一个纯文本文件， 内含用描述语言 script language写成的批处理 GrADS系统设置和命令， 可集成处理 GrADS命令。
ftp://iges.org/grads/
安装
§1．运行环境和基本操作
操作界面
提示：选择风景画（横放，11*8.5）还是肖 像画（竖放， 8.5*11 ）(硬拷贝输出) 直接回车－横放 no －竖放
§1．运行环境和基本操作
以versions 1.9为例，各程序功能：
gradsc Reads GRIB, gridded binary, BUFR, GrADS station data Writes GRIB
§1．运行环境和基本操作
VARS number: 表示变量描述开始，并给出变量个数 每个变量描述记录格式如下：
abrev:由１－８个字符组成的变量名，要求以字母开头 levs：表示该变量在本数据集中含有的垂直层次数，整型数，若 该量为地表变量，仅一层，则设为０ units：预留量，暂设为常数99 description: 为一段说明该变量的字符串，最多40个字符
UNDEF value（定义缺测值） TITLE string（用字符串简略描述数据文件的内容） （可省略）
Model.ctl
§1．运行环境和基本操作
时空维数环境设置
XDEF number LINEAR start increment XDEF number LEVELS value-list
number：x方向格点数，用整数； LINEAR：格距均匀；start:起 始经度，负数表示西经； increment：格距； LEVELS：格距不均 匀，用枚举法列出各网格点对应的经度值； value-list：顺序列出 各格点的经度值，可在下一行续行。
例如： ZDEF 7 LEVELS 1000 850 700 500 300 200 100
TDEF number LINEAR start-time increment
number:时次数 start-time格式：hh:mmZddmmmyyyy
hh:两位数的小时 mm：两位数的分钟 dd：为一或二位数的日期 mmm：三个字符的月份缩写 yyyy：为两位或四位数年份（两位代表年份：1950-2049）
3.
4. 5.
三 坐标要素控制
§3．变量和函数的定义、功能
§4．台站资料的使用
台站型离散资料画图的一般原则
•1、把台站资料r.dat写成二进制文件r.grd •2、编写r.grd数据的数据说明文件r.ctl •3、用r.ctl生成r.map文件 •4、生成插值用的背景格点文件 •5、编写grid.grd的.ctl文件：grid.ctl •6、把r.grd插值到grid.grd的网格点上，并画图
注：不给出时，hh缺省为00时，mm缺省为00分，dd缺省为1号 ，月年值不能缺省。整个时间串不能有空格。
例： 14:20Z22JAN1987 12Z1JAN1990 JUN1960
increment为时间增量，格式为vvkk,其中vv为增量值，用１或２位整型数表示， kk为增量类型，有如下几种：
HDF 数据格式
Binary数据格式
站点数据格式
GRIB数据格式
数据存放：
§1．运行环境和基本操作
(2).ctl 原始数据描述文件。是一个文本文件，用于描述原始数据文件的基本 信息，包括数据集文件名、数据类型、数据结构、变量描述等。在GrADS中 至少得首先打开一个数据描述文件。
格点数据的ctl文件
GRADS是当今气象界广泛使用的一种数据处理和显示软件 系统。该软件系统通过其集成环境，可以对气象数据进行 读取、加工、图形显示和打印输出。数据可以是格点资料， 也可以是站点资料；数据格式可以是二进制，也可以是 GRIB码，从而具有操作简单、功能强大、显示快速、出图 类型多样化、图形美观等特点。正因为如此，GRADS已迅 速成为国内外气象界通用的标准图形环境之一。
无格式数据。
整个数据集是一个大的五维数据场，包括三维物理空间、一
维物理变量、一维时间变量，存放时以二维数组片的形式按照水
平、垂直、物理变量、时间序列的顺序排放。
第一维：从西到东
第二维；从南到北
GrADS数据数据类型与格式
第三维：从低层到高层
Fra Baidu bibliotek
第四维： 变量
支持的数据格式：
第五维： 时间
netCDF数据格式
气象图形图像软件
课程安排
GraDS绘图软件 Origin绘图软件 Compaq Array Viewer介绍
GrADS
Grid Analysis and Display System
GrADS主页 iges.org/grads/ LASG中文讨论区: bbs.lasg.ac/?action-forumdisplayfid-3
gradsdods Reads GRIB, gridded binary, BUFR, GrADS station data, NetCDF, OPeNDAP (fka DODS) Writes GRIB, NetCDF
stnmap Maps station data
§1．运行环境和基本操作
基本概念
三种文件： .dat: 原始数据文件 .ctl ：原始数据描述文件 .gs：GrADS控制文件
§1．运行环境和基本操作
(1).dat 二进制无格式记录的原始数据文件，其数据既可以是格
点数据，也可以是站点数据。它们是从其它气象数据(如站点气象
报、格点气象报、模式格点输出结果)转换生成的，格式是二进制
*[同时关闭所有打开的文件，并释放所有定义的对象]
‘enable print h:\2\1.gmf ’ *[与print合用，输出图形]
'open h:\2\model.ctl'
‘set lat 20 40‘
*[设置纬度]
‘set lon 90 120‘
*[设置经度]
‘set z 5‘
*[设置高度]
2. 对图形类型为等值线（gxout=contour或shaded）起作用的设置：
set cint value [设置等值线间隔。 clear或display即重新设定或回到缺省设置] set clevs lev1 lev2…… [设置特定等值线值，只画lev1 lev2……值所在的等值线， clear或d即重新设定] set cmin value [不画低于此value的等值线， c或d即重新设定] set cmax value [不画高于此value的等值线， c或d即重新设定] set black val1 val2 [不画值介于val1和val2之间的等值线， c或d即重新设定]
例：XDEF 72 LINEAR 0.0 5.0
§1．运行环境和基本操作
YDEF number mapping start <increment> YDEF number mapping <LEVELS value-list> Mapping: 映射方式：
§1．运行环境和基本操作
ZDEF number mapping < start increment> ZDEF number mapping < value-list>
§1．运行环境和基本操作
操作系统： Linux Macintosh OSX Sun OS MS Windows SGI / IRIX DEC IBM / AIX 版本： GrADS 2.0（a0、a1、a2、a3）、GrADS 1.8s11、 GrADS 1.9b4、 GrADS 1.7...
(1)二维等值线绘图 (2)二维填色等值线绘图 (3)标示网格点数值 (4)矢量箭头形式绘二维风场 (5)天气填图 (６)直方图 (７)单线图 (８)用指定颜色填充特定点 (９)风向杆形式绘二维风场
§2．基本绘图指令
一 图形类型设置
命令： set gxout graphics-type
graphics-type包括以下图形类型：
gradsnc Reads GRIB, gridded binary, BUFR, GrADS station data, NetCDF Writes GRIB, NetCDF
gradshdf Reads GRIB, gridded binary, BUFR, GrADS station data, NetCDF, HDF-SDS Writes GRIB, HDF-SDS
coutour 二维等值线绘图 shaded 二维填色等值线绘图 grid 二维场不画图，在网格中央标示网格点数值 vector 矢量箭头形式绘二维风场 stream 流线形式绘二维风场 barb 风向杆形式绘二维风场 bar 直方图 line 单线图 fgrid 二维场不画等值线，仅用指定颜色填充特定点 model 天气填图