用GrADS制作ncep再分析资料气象要素图

3. 3
选择时次 NCEP/ NCAR 再分析资料每天有 4 个时次 , 而
时次选择的用法是用 set t 时次 , 其中, 1 月 1 日 00 时为 1, 06 时为 2, 12 时为 3, 18 时为 4, 1 月 2 日
( 上接 42 页) 2. 9 日照的预审 参照各定时云量、 能见度、 地面温度 , 检查日照 是否有漏计。一般云量少 , 能见度在 3. 0 km 以上 , 应有日照。 2. 10 封面、 封底的预审 封面封底有些数据每月都相同 , 这就造成部分 预审员的粗心大意 , 结果往往在最 简单的地方出 错。 ( 1) 概况、 备注栏内的字母与数字都应在全角 状态下输入( 特别注意标点符号的输入 ) 。 ( 2) 备注栏的日期应在半角状态下输两位数 , 不足十位补 0; 备注内容不能有非正常换行, 输入 的内容一定要连续; 若备注太长, 可以分段写 , 但每 段均要对应有日期。 ( 3) 如有搬迁值班室、 整改观测场等, 日期最好 能备注整改结束日期 , 而且对观测场高度、 气压表 高度等有关数据注意修改正确。
图 3 1991 年 1 月 1 日 00 时 250 hPa 气温剖面图
( 20 ~ 35 N 、 120 ~ 135 E)
draw t itle air NCEP DATA at 1000 hPa in 19919- 6 printim p. gif white 其中, 第 1 行 sdfopen air. 91. nc 用于打开所用 的 ncep 要素文件 , 每 2~ 7 行是写上图的标题 air NCEP DATA at 1000 hPa in 1991- 9- 6 , 第 8 行是 生成名叫 p. gif 的图形文件。
报表的预审 , 除了掌握一定的经验 , 还要在每 一次预审时都认真、 细心、 耐心 , 并多设想产生错误 记录的可能性 , 只要擅于借鉴以往的错误 , 避免错 误的再生, 日常值班保证记录的质量, 预审质量就 能保持高水平。
2004. 4 文章编号 : 1007- 6190( 2004) 04- 0035- 02
广东气象
35
用 GrADS 制作 ncep 再分析资料气象要素图
杨兆礼
( 中国气象局广州热带海洋气象研究所 , 广州 510080)
中图分类号 : P409
文献标识码 : B
Grads 气象图形系统是目前国内外气象界通用 的标准图形环境之一。用它来画气象要素图很方 便。本文简单 介绍了如何使用 Grads 气象图形系 统来对 NCEP/ NCAR 再分析资料进 行简单的绘图 处理。 1 NCEP/ NCAR 再分析资料简介 NCEP/ NCAR 再分析资料是美国国家环境预报 中心( NCEP) 和国家大气研究 中心 ( NCAR) 采用美 国全球气候资料同化分析系统进行四维同化分析 后获得的全球同化资料。其全部历史资料来源于 美国国家气候资料中心( NCDC) 。 使用 ncdump 命令 , 可以看到高空 NCEP/ NCAR 再分析资料的属性如表 1 所示。每天共有 4 个时 次, 有 73 个纬度和 144 个经度, 等压面共有 8~ 17 层( hPa) , 分别是: 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10hPa 。
NCEP/ NCAR 再分析资料的文件名中 , 第 1 个 单词便是我们绘图时要用到的变量名。例如其中 的 air. 91. nc, 91 代表了年份, 而 air 便是变量名 , 我 们可用命令 d air 来绘该变量的要素图。
图 2 1991 年 1 月 1 日 00 时的高空气温剖面图
( 10 N 、 60 W~ 60 E, 其它默认)
此外, 还可以固定纬度不变 , 只选择经度范围, 此时生成的 是沿经度变化 的剖面图, 如图 2 所示 ( set lon 10) 。
地面 NCEP/ NCAR 再分析资料共有 94 个纬度 和 192 个经度。包括了气温、 云强迫净长波通量、 对流降水率、 晴空长波射入通量、 射入长波辐射通 量、 射入日辐射通量、 地表热通量、 冰盖判示、 可能 蒸发率、 降 水率、 气 压、 水径 流、 风 U 分量 ( 10m 处) 、 风 V 分量 ( 10m 处) 等共 42 种要素。 2 选择绘图要素
( 4) 对缺测等异常记录的备注, 要检查是否与 报表里异常记录的时间、 处理结果等一致。 ( 5) 有关受灾情况应记录在纪要栏的重要天气 现象栏。 3 报表的修改 报表预审完后, 还应对 D 文件重新查对一遍, 确认所有疑误记录已纠正再上报, 有时由于急于求 成或其它原因 , 查出了疑误却没修改就上传 , 造成 错情的情况也屡见不鲜。 4 结 语
文件名 Air. 1991. nc Hgt1991. nc rhum1991. nc shum1991. nc omega1991. nc uwnd1991. nc vwnd1991. nc
图 1 1991 年 1 月 1 日 00 时 850hPa 气温图
( 20~ 35 N、 120~ 135 E)
表1 NCEP 高空等压面再分析资料说明 ( 以 1991 年为例 ) 要素名称 气温 位势高度 相对湿度 比湿 垂直速度 风 U 分量 风 V 分量 单位 K m % kg/ kg pascal/ s m/ s m/ s 层次 17 17 8 8 12 17 17
3 3. 1
选择绘图范围
经纬度选择 我们可以使用 set lat 20 35 的格式来选择 20 ~ 35 N 的范围 ( 默认情况为 90 S~ 90 N) , 相应地, 经度范围为 set lon 120 135 。此时用 air. 91. nc 文 件为例生成的 850 hPa 和要素图如图 1 所示。
36
杨兆礼 : 用 GrADS 制作 ncep 再分析资 料气象要素图
2004. 4
3. 2
等压面选择
00 时为 5, Байду номын сангаас次类推。例如 , 9 月 6 日 18 时为 1000,
类似地, 我们可以使用 set lev 250 来选择 250 则时次选择为 set t 1000 。 hPa 面上的要素图 , 如图 3 所示。 4 gs 文件的用法 gs 文件是 grads 的控制文件 , 用 grads 命令 run 执 行。 下面 举 个 简 单 例子 说 明 gs 文 件的 用 法 ( ncep. gs) : sdfopen air. 91. nc set lev 925 set t 1000 set lat - 60 80 set lon 50 320 d air