短期气候预测实习报告一
南京信息工程大学实验（实习）报告
实验（实习）名称 短期气候预测实习 日期10.29得分 指导教师 系 ~ 专业 ~ 年级 ~ 班次~ 姓名 Trichtu 学号 ~ 实习目的：
1掌握大气环流基本状况的表征－－平均、距平、纬偏图
2大气环流基本状况分析
3熟悉了解：资料格式FORTRAN语言
4 GRADS绘图（数据描述文件及执行程序
实习内容：
(1)计算1948-2007年（60年）1月的平均高度场，绘制环流平均图；
(2)计算2008年1月的高度距平，绘制高度距平场图（相对于1948-2007年共60年的平均）；
(3)计算2008年1月的高度场纬偏值，绘制环流纬偏图。
实习资料：
NCEP/NCAR 1948－2008年1～12月的500hPa月平均高度场资料 范围（90°S-90°N,0-360°E)
网格距2.5°×2.5°,纬向格点数144，经向格点数73
GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 实习方法：
①计算（1948-2007）1月平均高度场，绘平均图；
1m(i,j)=∑Hk(i,j)mk=1
(i,j)为经、纬格点，m为样本长度，k为年份
②计算2008年1月的高度距平，绘距平图；
m1 HAHk(i,j)∑k'(i,j)=Hk'(i,j)-mk=1
(i,j)为经、纬格点，m为纬圈格点数，k为年份
③计算2008年1月的高度场纬偏值，绘纬偏图
1m
HLAk(i,j)=Hk(i,j)-∑Hk(i,j) mi=1
(i,j)为经、纬格点，m为样本长度，k为年份
实习步骤
①编写计算1948-2007年（共60年）1月500hPa高度平均程序；
②编写计算2008年1月500hPa高度场距平程序（相对于1948-2007年共60年的平均）程序；
③编写2008年1月500hPa高度场纬向平均和偏差程序；
④根据①②③结果数据编写GRADS数据描述文件及执行文件，并绘图。 program EX1
real a(144,73,12,61),ave1(144,73),ave7(144,73),a2008(144,73)
real aalt(73),p2008(144,73)
integer i,j,imo,it
open(2,file='hgt500.grd',form='binary') !补充正确路径
open(4,file='ave60.grd',form='binary') !补充数据输出路径
do it=1,61
do imo=1,12
do j=1,73
do i=1,144
read(2)a(i,j,imo,it)
enddo;enddo;enddo;enddo
cccccccccccccc 求平均
do i=1,144
do j=1,73
do it=1,60
ave1(i,j)=ave1(i,j)+a(i,j,1,it)
ave7(i,j)=ave7(i,j)+a(i,j,7,it)
enddo
ave1(i,j)=ave1(i,j)/60.0
ave7(i,j)=ave7(i,j)/60.0`
enddo;enddo
do j=1,73
do i=1,144
write(4)ave1(i,j)
enddo;enddo
do j=1,73
do i=1,144
write(4)ave7(i,j)
enddo;enddo
CCCCCC2008 距平
do j=1,73
do i=1,144
a2008(i,j)=a(i,j,1,61)-ave1(i,j)
end do
end do
open(5,file='a2008.grd',form='binary')
do j=1,73
do i=1,144
write(5) a2008(i,j)
enddo
enddo
CCCCCCC偏距平
do j=1,73
do i=1,144
aalt(j)=aalt(j)+a2008(i,j)/144.0
end do
end do
do j=1,73
do i=1,144
partial(i,j)=a2008(i,j)-aalt(j)
end do
end do
open(6,file='p2008.grd',form='binary')
do j=1,73
do i=1,144
write(6) p2008(i,j)
enddo
enddo
end
Ctl文件：（其他的类似，不做重复）
dset e:\copy\a2008.grd
undef -9.99E+33
title monthly mean hgt sta value
xdef 144 linear 0.000 2.500
ydef 73 linear -90.000 2.500
zdef 1 levels 500
tdef 1 linear jan1948 1mo
vars 1
ab 1 -999 hgt abnormal 2008
endvars
;
画图gs文件：（其他类似，不做重复）
'reinit'
'open e:\copy\ave60.ctl'
'open e:\copy\a2008.ctl'
'open e:\copy\partial.ctl'
'enable print e:\copy\ave60.gmf'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.13'
'set gxout shaded'
'd ave1.1'
'set gxout contour'
'd ave1.1'
'draw title Jan average 500hPa height'
'print'
'c'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.13'
'set gxout shaded'
'd ave7.1'
'set gxout contour'
'd ave7.1'
'draw title Jul average 500hPa height'
'print'
'c'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.13'
'set gxout shaded'
'd ab.2'
'set gxout contour'
'd ab.2'
'draw title 2008 Jan 500hPa anormal height'
'print'
'c'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.13'
'set gxout shaded'
'd pa.3'
'set gxout contour'
'd pa.3'
'draw title 2008 Jan 500hPa partial height'
'print'
'c'
'disable print'
'reinit'
;
实习结果分析
(1)计算1948-2007年（60年）1月的平均高度场，绘制环流平均图
60年冬季全球500hPa环流平均场
从全球500hPa的高度场中看出，位势梯度在低纬地区较小，在中高纬地区较大。由于北半球的海陆分布特点，北半球的位势呈现冬季三槽三脊的特点可以明显的表现出来。而南半球受陆地较小，海洋热惯性较大的特点，分布较为平直。
60年夏季全球500hPa环流平均场
从全球位势梯度分布来看，南大北少，低纬小，中高纬大。最高的位势位于加勒比海地区。由于太阳直射北半球，北半球的受热增多，气块膨胀，平均位势高度均大于南半球。槽脊分布来看，北半球由弱的四槽四脊。
(2)计算2008年1月的高度距平，绘制高度距平场图（相对于1948-2007年共60年的平均）
2008年1月500hPa距平场
500hPa，北大西洋和欧洲大部分地区由强的正距平中心，东亚由一个较小的正距平中心。北大西洋和北极圈地区有强的副局评中心，青藏高原有弱的距平中心。而南半球表现为低纬地区的四个正距平中心和高纬地区一个位于海上的负距平中心。
(3)计算2008年1月的高度场纬偏值，绘制环流纬偏图
2008年1月500hPa纬偏场
北大西洋和欧洲大部分地区由强的正距平中心，东亚由一个较小的正距平中心。北大西洋和北极圈地区及青藏高原上有强的副局评中心。而南半球表现为低纬地区的四个正距平中心和高纬地区一个位于海上的负距平中心。