4.GrADS软件介绍2

记录2
记录3
记录4
GrADS基本命令和基本文件：
1. 命令窗与图形窗；Landscape 与 portrait（Y或N）选项的区别
2. 基本命令和基本文件 open， display（d），clear（c） dat文件，ctl文件，gs文件
本节内容：
1. ctl文件
2. 维数环境设置 3.图形文件的保存
2 维数环境设置
• 维数说明分为两类：
• set lon|lat|lev|time var1 • set x|y|z|t var1 <var2> <var2>
• 注意：
• 两者对应于同一组数据，只是前者为地球坐标，后 者为网格坐标。地球坐标的单位分别为：水平空间 单位用“度”，经度方向缺省为由西向东，东经为 正，西经为负或用大于180度表示；纬度方向缺省为 由南向北，南纬负，北纬为正；垂直方向由下向上， 单位为“百帕”；时间用绝对时间格式，格点坐标 用网格点数直接表示。 • var1表示起始坐标，var2表示终止坐标，不出现 var2时表示该维数方向是固定维数，规定var1< var2 • 两种坐标可以混用，其内部对应于同一数组维数环 境。
GrADS软件介绍
Grid Analysis and Display System （第二节）
2) test2.dat文件：共有4个记录，每个记录长度 为3。要求：将该文件内容读出，并写入到一个新的 有格式文件中，每个记录对应一行，每行3个数字。
记录1 real,dimension(4,3) :: a
options
20130101.dat 20130102.dat 20130103.dat …. 20130131.dat 20130201.dat …. 20130228.dat 20130301.dat …. 20131231.dat
template
dset %y4%m2%d2.dat options template … … tdef …. 365 linear %id2%im2%iy4 1dy
变量 说明
变量名 层数
默认变量循环顺序：
X(lon) Y(lat) Z(lev)VarT(time)
可选关键字：dtype
Binary Grib Station Netcdf 普通二进制数据（默认） 站点数据 NC数据
可选关键字：options
yrev Zrev Sequential 365_day_calendar Template byteswapped big_endian little_endian 自北向南的数据 自高层到低层的数据 顺序存储数据 忽略闰年（2月29号） 多个文件 文件存储的字节次序与本机相反 在SUN，SGI等机器上生成的数据 在PC等机器上生成的数据
3. 图形文件的保存
Grads2.1 以上 gxprint **.png (jpg,pdf,eps…) white
Grads2.0以下 • gmf文件: enable print pirnt disable print
*.gmf
4. 图形类型设置
display： 显示变量（默认为单线图（一维情况）或者等值线图 （二维情况））
0
起始 坐标
2.5
坐标 间隔
方向： 自西向东、自南向北（不可改变） 单位： 度（经度、纬度） 正负值： 东（西）半球经度为正（负） 北（南）半球纬度为正（负）
Z方向维数说明
zdef
Z方向 维数 定义
5
共有 5层
levels
列举法 给出 坐标值
1000 850 700 500 200
方向：自下而上 单位：hPa （也可用其他单位）
用一个.ctl文件描述多个数据文件
注意事项
• 数据描述文件为文本格式文件 • 每行记录的各项以空格分开 • 注释行在第一列打“*”，注释行不能出现在 变量列表中 • 每行记录不超过80个字符 • 对某一层某一变量在某一时刻，X，Y数据点 构成了一个水平网格，该网格严格对应于 FORTRAN中的数组存放顺序，第一维总是从 西变到东，第二维从南变到北，并按X（经 度）、Y（纬度）、Z（高度层数），然后是 各种物理变量，最后是T（时次）这样的顺序 存放的
时间维说明：
tdef
时间 维数 定义
24
共有 24个 时刻
linear
线性 间隔
00:00z01jan1979
起始 时间
1mo
时间 间隔
必须为线性间隔:linear 时间必须用规定格式表达
要素变量说明:
zdef 5 levels 1000 850 700 500 200
air
5
99
排放 顺序
monthly air temperature(K)
set gxout : 设置图形类型
set gxout 可选项
• contour：二维等值线绘图（缺省） • shaded：二维填色等值线绘图 • grid：二维场不绘图，以网格形式在各网格点中央标出该 场点数值 • vector：矢量箭头形式绘二维风场（缺省） • stream：流线形式绘二维风场 • barb：风向杆形式绘二维风场 • bar：对一维场不绘单线图，而绘直方图 • line：对一维场绘直线图（缺省） • fgrid：对二维场不绘等值线图，只将特定值的格点用指 定颜色填充该格，与命令set fgval value color value color……一起使用 • grfill：二维填色，与shaded的区别在于shaded有对网格 的平滑，而grfill是按网格填色 • 如果图中要做矢量、流线或风向杆绘图时，要求显示命令 后给出用分号“；”分隔的两个分量场，前者理解为U分 量，后者为V分量。
1.ctl文件的内容：
数据文件名 数据标题 特殊格式说明 缺测值 X方向维数 Y方向维数 Z方向维数 时间维数 变量总数 大气温度 位势高度 U风分量 V风分量 变量说明结束
X, Y方向维数说明：
xdef
X方向 维数 定义
1来自百度文库4
共有 144个 格点
linear
坐标值 线性增 加