第2章GrADS数据格式

tdef 4 linear jan1994 1mo
vars……
2.3 其他格式数据资料的使用
1. NetCDF数据格式的资料（*.nc） GrADS可以直接处理这种格式的数据，
不需要另外编写数据描述文件。
使用命令： ga_> sdfopen <路径>文件名.nc
2. GRIB码数据格式的资料
GrADS可以直接读取该格式的数据，但是 须使用“grib2ctl”和“gribmap”命令先生 成该文件的数据描述文件（.ctl）和指针 文件（.idx）。
2. 数据描述文件的形式 以上述格点数据描述文件为例：
DSET D:\model.le.dat 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
ts 0 99 Surface Temperature
p
0 99 Precipitation
ENDVARS
数据描述文件中各记录的含义：
• DSET data-set-name 给定二进制原始数据文件的
文件名（包含路径）。 例如：dset D:\model.le.dat
• UNDEF value 定义缺测值。
• TITLE string
用字符串string简略描述数 据文件的内容。
• XDEF number LINEAR start increment 或XDEF number LEVELS value-list
设置X方向网格点值与经度的对应关系。
• YDEF number LINEAR start increment 或YDEF number LEVELS value-list
第2章 数据文件
• 常用数据格式的文件特点 • Binary数据格式文件的使用 • NetCDF数据格式文件的使用 • GRIB码数据格式文件的使用
在使用GrADS绘图前，必须确保 所用数据文件满足GrADS的数据格式 要求。
2.1 常用数据格式的文件特点
• Binary（二进制）数据格式：可以是 格点数据或站点数据，可以用Fortran 或者C语言读写。
本章小结
本 章 主 要 介 绍 了 GrADS 绘 图 软 件 中 常 用 数据文件的使用。要求掌握数据描述文件的建 立；Binary文件和NetCDF文件的使用。
实习一 撰写ctl文件
1、数据说明：现有一个二进制格点数据文件 d:/ERA-interim-trend.dat 数据资料空间范 围：0-359经度，90S-90N；物理量:Tr；层 次：23层， 1000 850 800 750 700 650 600 550 500 450 400 350 300 275 250 225 200 175 150 125 100 70 50 hPa；时间：2011年4 月2日00时；分辨率：1.0*1.0。要求编写 数据描述文件trend.ctl。
2：现有一个二进制格点数据文件 d:/mhy.grd 数据资料空间范围： 60-150°E，0-40°N；物理量：u、 v、sst；层次：850、 200hPa；时 段：1982.1-1985.12的逐月资料； 分辨率：2.5*2.5。要求编写数据描 述文件mhy.ctl。
*关于template使用
多个连续时间序列原始数据文件想用一个数据描 述文件统一地描述这些原始数据时采用该选项。
例如：
h199401.grd h199402.grd h199403.grd h199404.grd
. . .
dset h%y4%m2.grd title sample dtype grid undef 9999 options template xdef …… ydef……. zdef……
（xdef,ydef,zdef,tdef） • 变量定义（vars,endvars）
• DTYPE <option> option：grid 或 station
说明数据文件的数据类型。 默认设置: dtype grid
• OPTIONS <keywords>
定义数据存放方式的选项，keywords 有：
<yrev>
<zrev>
<sequential> <byteswapped>
详见P22
<big-endian>
<little-endian>
<template>
3.数据描述文件的建立
数据描述文件为纯文本格式，可 用一般的文本编辑器产生（如记事本， 写字板等）。描述文件编辑完成后， 保存时后缀名为.ctl。
数据描述文件（ctl文件）是原始 数据文件的描述文件，用以描述原始 数据集的基本信息。
数据说明
model.le.dat是二进制格点数据文件，为 一组模式输出的全球5天的数值预报结果。 有7个层次：1000、850、700、500、300、 200、100hPa，多个要素：ps、u、v、z、 t…，时间从1987年1月2日开始，间隔1天， 已知X方向间隔为5个单位，Y方向间隔为4 个单位。
关于数据文件model.le.dat的具 体说明，可以打开数据描述文件 model.le.ctl查看。
model.le.ctl内容：dset d:\model.le.dat 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
时间表示法 ：hh:mmZddmmmyyyy 时间增量表示法：vvkk
Fra Baidu bibliotek
• VARS number
表示变量描述开始，同时给出变 量个数number 。
用ENDVARS表示数据描述文件结束。
3.数据描述文件的构成
每个数据描述文件一般包含以下几项： • 被描述的数据文件名（dset） • 该数据说明文件的标题（title） • 数据类型、和选项（dtype,options） • 缺测值说明（undef） • 时间、空间维数环境设置
100
TDEF 4 LINEAR 0Z10apr1991 12hr
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
定义Y方向网格点值与纬度的映射关系。
• ZDEF number LINEAR start increment 或ZDEF number LEVELS value-list
设置垂直网格点与气压面的映射关系。
• TDEF number LINEAR start-time increment
设置网格值与时间的映射关系 。
• NetCDF网络通用数据格式：自定义 数据格式，精确性好，便于传输。
• GRIB码数据格式：是自定义的数据 格式，压缩率高，占用空间小。
2.2 Binary数据格式文件的使用
GrADS绘图不能直接使用“Binary 数据文件”，而是通过“数据描述文 件”间接使用“数据文件”。
1. 数据描述文件的概念