气象统计实习报告 解读.docx

气象统计实习报告

专业：大气科学

班级：xxxx级x班学号：2012130xxxx姓名：

1

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

一、实习结果

1981.1距平场

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

2

1981.1气候场

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

1981.1均方差场

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

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

Fortran程序:

programmain

parameter(nx=37,ny=17,mo=12,yr=4)realvar(nx,ny,mo,4)!数据

realvars(nx,ny,mo)!4年气候态realjp(nx,ny,mo,4)!距平realfc(nx,ny,mo)!方差

realjfc(nx,ny,mo)!均方差integeri,j,m,y,irecreal::summ=0.0

!*********************求均方差*********************dom=1,12

doi=1,37

doj=1,17

doy=1,4

3

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

fc(i,j,m)=summ/4!方差jfc(i,j,m)=sqrt(fc(i,j,m))!均方差summ=0.0

enddoenddo

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')dom=1,12

enddo

dom=1,12

enddo

doy=1,4

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

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

dom=1,12

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

enddoenddoclose(11)close(12)close(13)End

Gs程序

1.距平场

'opene:\anomaly.ctl''setgxoutshaded'

'dxiaobao'

'setgxoutcontour''dxiaobao'

'drawtitle1982年1月'

'enableprinte:\anomaly.gmf''print'

'disableprint'2.气候场

'opene:\climate.ctl'

'setgxoutshaded''dxiaobao'

4

'setgxoutcontour''dxiaobao'

'drawtitle1982.1'

'enableprinte:\climate.gmf''print'

'disableprint'

3.均方差场

'opene:\deviation.ctl''setgxoutshaded'

'dxiaolu'

'setgxoutcontour''dxiaolu'

'drawtitle1982.1deviation''enableprinte:\deviation.gmf''print'

'disableprint'

*实习二计算给定数据资料的简单相关系数和自相关系数根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。

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

（滞后长度τ最大取10）

5

一、实习结果

计算简单相关系数

计算出相关系数为r=0.4685170

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

τ最大取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!协