短期气候预测实习报告一

实验（实习）名称 短期气候预测实习 日期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月平均高度场，绘平均图；

1mH(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月的高度场纬偏值，绘纬偏图

1mHLAk(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纬偏场

北大西洋和欧洲大部分地区由强的正距平中心，东亚由一个较小的正距平中心。北大西洋和北极圈地区及青藏高原上有强的副局评中心。而南半球表现为低纬地区的四个正距平中心和高纬地区一个位于海上的负距平中心。