短期气候预测基础实习一

实习一：大气环流状况的表征

一、实习目的及要求

1. 掌握大气环流基本状况的表征－－平均、距平、纬偏图

2. 大气环流基本状况分析

3. 熟悉了解：资料格式

FORTRAN 语言

GRADS 绘图（数据描述文件及执行程序）

二、实习内容

1. 计算1948-2007年（60年）1月的平均高度场，绘制环流平均图；

2. 计算2008年1月的高度距平，绘制高度距平场图（相对于1948-2007年共60年的平均）；

3. 计算2008年1月的高度场纬偏值，绘制环流纬偏图。

三、实习步骤

1、实习资料

NCEP/NCAR 1948－2008年1～12月的500百帕月平均高度场资料

范围：（90°S-90°N,0-360°E)

网格距：2.5°×2.5°,纬向格点数144，经向格点数73

资料为GRD 格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 2、实习方法

① 计算（1948-2007）1月平均高度场，绘平均图；

(i,j)为经、纬格点 m 为样本长度，k 为年份 ② 计算2008年1月的高度距平，绘距平图；

(i,j)为经、纬格点 m 为样本长度，k 为年份 ③ 计算2008年1月的高度场纬偏值，绘纬偏图

(i,j)为经、纬格点 m 为纬圈格点数，k 为年份

3、实习步骤

（1）编写计算1948-2007年（共60年）1月500hPa 高度平均程序；

gs 文件：

'reinit'

'open d:\nyclimate\sh1\hgt500.ctl'

'enable print d:\nyclimate\sh1\have.gmf' 'set lev 500'

),(1),(1

j i H m j i H m

k k ∑==),(1),(),(1j i H m j i H j i HA m

k k k k ∑=-=),(1),(),(1

j i H m j i H j i HLA m

i k k k ∑=-=

'set lon 0 360'

'set t 1'

'define have=ave(hgt,t+0,t=720,12)'

'd have'

'print'

'disable print'

;

1月份500hPa高度场平均北半球投影的gs文件：'reinit'

'open d:\nyclimate\sh1\hgt500.ctl'

'enable print d:\nyclimate\sh1\have1.gmf' 'set frame circle'

'set mproj nps'

'set mpdest mres'

'set poli on'

'set map 1 1 0.6'

'set lev 500'

'set lat 0 90'

'set lon 0 360'

'set t 1'

'define have=ave(hgt,t+0,t=720,12)'

'set gxout shaded'

'set cmin 5880'

'd have'

'set ccolor rainbow'

'set gxout contour'

'd have'

'draw title Jan average 500hPa height'

'print'

'disable print'

;

2008年1月份500hPa高度场北半球投影的gs文件：'reinit'

'open d:\nyclimate\sh1\hgt500.ctl'

'enable print d:\nyclimate\sh1\sh2.gmf'

'set frame circle'

'set mproj nps'

'set mpdest mres'

'set poli on'

'set map 1 1 0.6'

'set lev 500'

'set lon 0 360'

'set t 721'

'set time jan2008'

'set gxout shaded'

'set cmin 5880'

'd hgt'

'set ccolor rainbow'

'set gxout contour'

'd hgt'

'draw title Jan 2008 500hPa height'

'print'

'disable print'

;

（2）编写计算2008年1月500hPa高度场距平程序（相对于1948-2007年共60年的平均）程序:

gs文件：

'reinit'

'open d:\nyclimate\sh1\hgt500.ctl'

'enable print d:\nyclimate\sh1\uave2.gmf'

'set lev 500'

'set lat -90 90'

'set lon 0 360'

'set t 1'

'define have=ave(hgt,t+0,t=720,12)'

'modify have seasonal'

'set t 721'

'hano=hgt-have'

'set time jan2008'

'set gxout shaded'

'set cmax 0'

'd hano'

'set gxout contour'

'd hano'

'print'

'disable print'

;

（3）编写2008年1月500hPa高度场纬向平均和偏差程序；

gs文件:

'reinit'

'open d:\nyclimate\sh1\hgt500.ctl'

'enable print d:\nyclimate\sh1\uave.gmf'

'set lev 500'