气象统计实习报告.docVIP

PAGE

PAGE4

气象统计实习报告

专业：大气科学

班级：xxxx级x班

学号：2012130xxxx

姓名：

*实习一求500hPa高度场气候场、距平场和均方差场

一、实习结果

1981.1距平场

1981.1500hpa高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。

setgxoutcontour

dxiaobao

drawtitle1982.1

enableprinte:\climate.gmf

print

disableprint

均方差场

opene:\deviation.ctl

setgxoutshaded

dxiaolu

setgxoutcontour

dxiaolu

drawtitle1982.1deviation

enableprinte:\deviation.gmf

print

disableprint

*实习二计算给定数据资料的简单相关系数和自相关系数

根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：

1）计算两个气温之间的简单相关系数。

2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）

一、实习结果

计算简单相关系数

计算出相关系数为r=0.4685170

（2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）

可以知道，年平均气温在滞后长度j=7，冬季j=4最大

二、相关的fortran程序（部分）

implicitnone

realx(20)!年平均气温

realy(20)!冬季平均气温

real::zx(10)=(/0,0,0,0,0,0,0,0,0,0/)

real::rzx(10)=(/0,0,0,0,0,0,0,0,0,0/)

real::s=0.0!协方差

real::jx=0.0

real::jy=0.0

real::jfx=0.0

real::jfy=0.0,r

integeri,t,j

real::m=0.0

datax/3.40,3.30,3.20,2.90,3.40,2.80,3.60,3.00,2.80,3.00,

3.10,3.00,2.90,2.70,3.50,3.20,3.10,2.80,2.90,2.90/

datay/3.24,3.14,3.26,2.38,3.32,2.71,2.84,3.94,2.75,1.83,

2.80,2.81,2.63,3.20,3.60,3.40,3.07,1.87,2.63,2.47/

doi=1,20

jx=jx+x(i)

jy=jy+y(i)

enddo

doi=1,20

s=s+(x(i)-jx/20)*(y(i)-jy/20)

jfx=jfx+(x(i)-jx/20)**2

jfy=jfy+(y(i)-jy/20)**2

enddo

r=s/20/sqrt(jfx/20*jfy/20)

print*,r=,r

dot=1,10

doi=1,20-t

zx(t)=zx(t)+(x(i)-jx/20)*(x(i+t)-jx/20)

enddo

rzx(t)=zx(t)/(20-t)/(jfx/20)

print*,t,rzx(t)

if(abs(rzx(t))m)then

m=abs(rzx(t))

j=t

endif

enddo

print*,年平均温度的自相关系数绝对值最大的滞后时间长度

print*,j

实习三计算给定数据的落后交叉相关系数和偏相关系数

根据下表北京冬季（12月~2月）气温资料计算：12月气温与1月和2月气温的落后交叉相关系数（滞后长度τ最大取10）和偏相关系数。在实习报告中给出程序。

年份

12月

1月

2月

1951

1.0

-2.7

-4.3

1952

-5.3