气象统计实习报告材料(2)

气象统计实习报告

专业：大气科学

班级：xxxx级x班

学号：2012130xxxx

姓名：

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

一、实习结果

1981.1距平场

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

1981.1气候场

1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。

1981.1均方差场

1981.1的均方差场在欧亚大陆的南侧有最小值，在10N~20N的南侧岁纬度而降低，在其北侧随纬度升高。

二、相关的fortran程序，gs文件

（1）Fortran程序:

program main

parameter(nx=37,ny=17,mo=12,yr=4)

real var(nx,ny,mo,4)!数据

real vars(nx,ny,mo)!4年气候态

real jp(nx,ny, mo,4)!距平

real fc(nx,ny,mo)!方差

real jfc(nx,ny,mo)!均方差

integer i,j,m,y,irec

real::summ=0.0

!*********************求均方差*********************

do m=1,12

do i=1,37

do j=1,17

do y=1,4

summ=summ+jp(i,j,m,y)**2

enddo

fc(i,j,m)=summ/4!方差

jfc(i,j,m)=sqrt(fc(i,j,m))!均方差

summ =0.0

enddo

enddo

enddo

!******************************************

open(11,file='d:\shixi\vars.grd',form='binary' )

open(12,file='d:\shixi\jp.grd',form='binary' )

open(13,file='d:\shixi\jfc.grd',form='binary' )

do m=1,12

write(11)((vars(i,j,m),i=1,nx),j=1,ny)

enddo

do m=1,12

write(13)((jfc(i,j,m),i=1,nx),j=1,ny)

enddo

do y=1,4

do m=1,12

write(12)((jp(i,j,m,y),i=1,nx),j=1,ny)

enddo

enddo

close(11)

close(12)

close(13)

End

（2）Gs程序

1.距平场

'open e:\anomaly.ctl'

'set gxout shaded'

'd xiaobao'

'set gxout contour'

'd xiaobao'

'draw title 1982年1月'

'enable print e:\anomaly.gmf'

'print'

'disable print'

2.气候场

'open e:\climate.ctl'

'set gxout shaded'

'd xiaobao'

'set gxout contour'

'd xiaobao'

'draw title 1982.1'

'enable print e:\climate.gmf'

'print'

'disable print'

3.均方差场

'open e:\deviation.ctl'

'set gxout shaded'

'd xiaolu'

'set gxout contour'

'd xiaolu'

'draw title 1982.1 deviation'

'enable print e:\deviation.gmf'

'print'

'disable print'

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

根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。

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

一、实习结果

（1）计算简单相关系数

计算出相关系数为r=0.4685170

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

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

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

implicit none

real x(20)!年平均气温

real y(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

integer i,t,j

real ::m=0.0

data x/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/

data y/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/

do i=1,20

jx=jx+x(i)

jy=jy+y(i)

end do

do i=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

end do

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