短期气候预测实习报告一

南京信息工程大学实验（实习）报告

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

(i,j)为经、纬格点，m 为样本长度，k 为年份

②计算2008年1月的高度距平，绘距平图；

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

③计算2008年1月的高度场纬偏值，绘纬偏图

),(1),(1

j i H m j i H m k k ∑==11(,)(,)(,)m k k k k HA i j H i j H i j m ''==-∑

(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

),(1),(),(1

j i H m j i H j i HLA m

i k k k ∑=-=

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'