一、实习目的及要求
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、实习步骤
d201
（1）编写计算1948-2007年（共60年）1月500hPa高度平均程序；
gs文件：
'reinit'
'open d:\nyclimate\l'
'enable print d:\nyclimate\sh1\have.gmf'
'set lev 500'
'set lat -90 90'
'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\l'
'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\l'
'enable print d:\nyclimate\f'
'set frame circle'
'set mproj nps'
'set mpdest m西堤红山res'
'set poli on'
'set map 1 1 0.6'
'set lev 500'
'set lat 0 90'
'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\l'
'enable print d:\nyclimate\f'
'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\l'
'enable print d:\nyclimate\f'
'set lev 500'
'set lat -90 90'
'set lon 0 360'
二元二次方程的解法
'set t 721'
'define uave=ave(hgt,lon=0,lon=360,-b)'
