GrADS实用绘图(基3)

这是GrADS所用的 数据排列方式
♥♥♥
数据
.ctl文件
5
LASG动力论坛版
五. 数据准备(Gr)
数据文件名 数据标题 特殊格式说明 缺测值 X方向维数 Y方向维数 Z方向维数 时间维数 变量总数 地面变量SLP 地面变量2米温度 等压面各层温度 等压面各层U风分量 等压面各层V风分量 变量说明结束
16
LASG动力论坛版
五. 数据准备(Gr)
— .ctl文件语法9-3
dset %ch 和 chsub t1 t2 string 若有如下NetCDF 数据文件, pr.1851-1900.nc包含了50 年 的逐月资料(即第1-600月), pr.1901-2000.nc包含100 年(第 601-1800月), 则:
– 方向: 自西向东, 自南向北 (永远不变) – 单位: 度 纬度/经度 – 正负值: 以东为正/ 以北为正
7
LASG动力论坛版
五. 数据准备(Gr)
• Z 方向的维数说明
— .ctl文件语法3
zdef 7 levels 1000 850 700 500 300 200 100 列举法 给出 坐标值
HDF4(SDS)
20
LASG动力论坛版
五. 数据准备(Gr)
Binary格式数据
— Binary格式1 ♥♥♥
– GrADS最基本, 我们最容易准备 – 常用后缀 .dat .grd .bin ,…… – 需要 .ctl 文件一起使用 – 是无格式文件 (form=‘unformatted’) – 若用Fortran 读写
18
LASG动力论坛版
五. 数据准备(Gr)
• 语法小结
– – – –
— .ctl文件语法10
.ctl文件都包含哪些方面的内容? X, Y, Z, T 4个维数的说明有什么特点? 数据说明的关键之处在哪里? 如何说明一系列的数据文件?
• 注意事项
♥♥♥
– 每行顶格写; – 虽然大小写没关系, 我们推荐在一个关键字中都用 大写或都用小写; 如 写成 dset, 而不写成 DsEt. – 注释行请用 开始, 并顶格写; – 在 vars ~ endvars 之间不要写注释行; – 不要在GrADS中使用制表(tab)键!
19
LASG动力论坛版
五. 数据准备(Gr)
支持的数据格式
– – – – – – netCDF (自描述) HDF-SDS (卫星) GRIB Binary BUFR
— 支持的数据格式
GrADS Binary
ASCII
netCDF
GRIB
(自制)
– 站点
ASCII (台风路径) BUFR ( v1.9版本 )
— .ctl文件语法7-1
pdef关键字
－lcc | lccr －nps | sps －BILIN | file －eta.u | pse －ops
－特殊投影坐标系下的数据 Lambert投影 极射赤面投影 其它直线网格投影 US NMC eta模式 | SSM/I极射赤面投影 CSU RAMS斜极射赤面投影
zdef 7 levels 1000 850 700 500 300 200 100
• 文件路径推荐用UNIX惯例写
– d:/mydata/test.grd
10
LASG动力论坛版
五. 数据准备(Gr)
– – – – – – – – – yrev zrev sequential 365_day_calendar template byteswapped big_endian little_endian cray_32bit_ieee
3
LASG动力论坛版
五. 数据准备(Gr)
本节内容提要
– .ctl文件语法
– – – – – –
GrADS支持的数据格式 Binary数据格式 netCDF数据格式 GRIB数据格式 站点数据格式 进阶向导
4
LASG动力论坛版
五. 数据准备(Gr)
• 3-D空间数据排列方式
– 自西向东 – 从南到北 – 从低层到高层
12
LASG动力论坛版
五. 数据准备(Gr)
pdef关键字
－使用注意
— .ctl文件语法7-2
－u,v分量需与GRIB数据格式中风场变量33和34的说明格式 相同, 以便执行矢量的旋转运算
u 18 V 18 u10 0 v10 0 33,100 34,100 33,105 34,105 U-Wind Components on Pressure Levels V-Wind Components on Pressure Levels 10 Meter U Wind 10 Meter V Wind
— .ctl文件语法6
• options 是可选的关键字, 其选项有
自北向南的数据 自高层到低层的数据 顺序存贮的数据 任何年份均是365天 多个文件 文件存储的字节次序与本机相反 在SUN, SGI等机器上生成的数据 在PC, DEC等机器上生成的数据 在Cray机上生成的数据
11
LASG 数据准备(Gr) 五. 动力论坛版
pdef isize jsize { lcc | lccr } latref lonref iref jref Slat Nlat std_lon dx dy lcc 表示数据中风分量已经旋转成经纬度坐标, lccr 表示风分量仍然是lambert 正型投影,需GrADS 对其进行旋转. isize, jsize X, Y 方向的总点数; latref, lonref 参考点的经, 纬度; iref, jref 参考点的(i,j)位置; Slat, Nlat Lambert 投影m=1 时南, 北两个纬度; std_lon 标准经度; dx, dy 以M 为单位的X 和Y 方向的网格距.
－或者用vectorpairs来说明(v1.9)
vectorpairs u,v u10,v10 uflux,vflux …
－BILIN可用于netCDF格式数据
13
LASG动力论坛版
五. 数据准备(Gr)
— .ctl文件语法8
• dtype关键字 — 用于说明数据格式的类 型, 当数据为格点binary格式时省略. • dtype的可选参数
dset ^pr.%ch.nc chsub 1 600 1851-1900 chsub 601 1800 1901-2000 options template dtype netcdf tdef 1800 linear jan1851 1mo
若两个文件分别位于不同的磁盘上:
dset %ch chsub 1 600 /disk1/pr.1851-1900.nc chsub 601 1800 /disk2/pr.1901-2000.nc
– 方向: 自下而上 – 单位: hPa / 其他
可用linear 关键字替换
8
LASG动力论坛版
五. 数据准备(Gr)
• T 方向的维数说明
— .ctl文件语法4
1mo 时间 步长 3mn 6hr 1dy 1mo 1yr 分钟 小时 天 月 年
– tdef 12 linear 00:00z21jan2000
1
LASG动力论坛版
GrADS实用绘图
第三讲 数据准备
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
2
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1 — 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2 — 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
21日00:00时 2000年1月 (按情况可省略) (不可省略)
– 只有linear关键字 – 起始时间必须按照上述格式写 – 年份用4位数写
9
LASG动力论坛版
五. 数据准备(Gr)
• 要素说明
地面变量 2米温度 等压面各层 温度
— .ctl文件语法5
t2m 0 99 2m temperature( degC) temp 7 99 upper level temperature (degC)
— .ctl文件语法1
• 格点数据描述文件 .ctl 的写法
dset d:/mydata/test.grd title This is a GrADS test data options sequential yrev undef -9999.0 xdef 144 linear 0 2.5 ydef 73 linear -90 2.5 zdef 7 levels 1000 850 700 500 300 200 100 tdef 12 linear 00:00z01jan2000 1mo vars 5 slp 0 99 Sea Surface Temperature (hPa) t2m 0 99 2m temperature( degC) temp 7 99 upper level temperature (degC) uwnd 7 99 upper level U component of wind (m/s) vwnd 7 99 upper level V component of wind (m/s) endvars
1位数字表示的年代 3位数字表示的年代 2位数字表示的年 4位数字表示的年 1位或2位数字表示的月 2位数字表示的月 3位字符表示的月 1位或2位数字表示的日 2位数字表示的日 1位或2位数字表示的小时 2位数字表示的小时 3位数字表示的小时 2位或3位数字表wk.baidu.com的预报时效 3位数字表示的预报时效 2位数字表示的分钟
6
LASG动力论坛版
五. 数据准备(Gr)
• X, Y 方向的维数说明
– xdef
144 linear 坐标值 线性 增加
— .ctl文件语法2
0 起始 坐标 值 2.5 坐标 间隔 可用levels 关键字替换
Y: gausT62 gausR15 gausR20 gausR40
X(lon) 共有 方向 144个 维数 格点 定义
17
LASG动力论坛版
五. 数据准备(Gr)
%ix1 %ix3 %iy2 %iy4 %im1 %im2 %imc %id1 %d1 %id2 %ih1 %ih2 %ih3 %in2
— .ctl文件语法9-4
tdef 中起始时间的替代方式
起始的1 位数年代(v1.9) 起始的3 维数年代 (v1.9) 起始的2 位数年 起始的4 位数年 起始的1/2 位数月 起始的2 位数月(必要时前置0) 起始的3 字符月 起始的1/2 位数天(必要时前置0) 1 位数字天 起始的2 位数天 起始的1/2 位数时 起始的2 位数时 起始的3 位数时 起始的2 位数分(必要时前置0)
– – – – – grib station bufr hdfsds netcdf 站点数据 v1.9版本新增
14
LASG动力论坛版
五. 数据准备(Gr)
20070101.grd 20070102.grd … 20070131.grd 20070201.grd … 20070228.grd 20070301.grd … 20071231.grd
15
LASG动力论坛版
五. 数据准备(Gr)
文 件 名 中 时 间 格 式 说 明 符
%x1 %x3 %y2 %y4 %m1 %m2 %mc %d1 %d2 %h1 %h2 %h3 %f2 %f3 %n2
— .ctl文件语法9-2
(如: 1910, 1920, 1930) (如: 1980, 1990, 2000) (如: 1, 2, …, 11,12) (如: 01, 03, …, 12) (如: Jan, Feb, Mar, …) (如: 1, 2, …, 31) (如: 01, 02, 10, 31) (如: 1z, 2z, …, 12z) (如: 01, 02, 10, 31) (如: 120, 012) (如: 12, 24, 48, 72, 120) (如: 012, 024, 072, 120) (如: 01, 35, 60)
— .ctl文件语法9-1
• 用一个 .ctl 文件描述多个数据文件
dset %y4%m2%d2.grd options template … … tdef 365 linear %id2%im2%iy4 1dy …
.ctl
• %y4%m2%d2 是时间的替代格式符
%y4 | %iy4 4位数的年份 %m2 | %im2 2位数的月份 %d2 | %id2 2位数的日期