气象统计实习报告2.docx
- 文档编号:4808958
- 上传时间:2022-12-09
- 格式:DOCX
- 页数:16
- 大小:296.68KB
气象统计实习报告2.docx
《气象统计实习报告2.docx》由会员分享,可在线阅读,更多相关《气象统计实习报告2.docx(16页珍藏版)》请在冰豆网上搜索。
气象统计实习报告2
气象统计实习报告-
(2)
气象统计实习报告
专业:
大气科学
班级:
xxxx级x班
学号:
2012130xxxx
姓名:
*实习一求500hPa高度场气候场、距平场和均方差场
一、实习结果
1981.1500hpa高度场在欧亚大陆为正距平,在印度洋和太平洋为负距平。
1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。
1981.1的均方差场在欧亚大陆的南侧有最小值,在10N~20N的南侧岁纬度而降低,在其北侧随纬度升高。
二、相关的fortran程序,gs文件
(1)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,irec
!
*********************求均方差*********************
dom=1,12
doi=1,37
doj=1,17
doy=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))!
均方差
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')
dom=1,12
write(11)((vars(i,j,m),i=1,nx),j=1,ny)
enddo
dom=1,12
write(13)((jfc(i,j,m),i=1,nx),j=1,ny)
enddo
doy=1,4
dom=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.距平场
'opene:
\anomaly.ctl'
'setgxoutshaded'
'dxiaobao'
'setgxoutcontour'
'dxiaobao'
'drawtitle1982年1月'
'enableprinte:
\anomaly.gmf'
'print'
'disableprint'
2.气候场
'opene:
\climate.ctl'
'setgxoutshaded'
'dxiaobao'
'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)
一、实习结果
(1)计算简单相关系数
(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:
:
jfy=0.0,r
integeri,t,j
datax/3.40,3.30,3.20,2.90,3.40,2.80,3.60,3.00,2.80,3.00,&
datay/3.24,3.14,3.26,2.38,3.32,2.71,2.84,3.94,2.75,1.83,&
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
1952
1953
1954
1955
1956
1957
1958
1959
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
一、实验结果
偏相关系数
计算出的12月气温与1月气温的偏相关系数为0.327,12月气温与2月气温的偏相关系数为0.290。
二、相关的fortran程序
Funtionarea(x,y)
implicitnone
realx(30)
realy(30)
real:
:
jfy=0.0,r
integeri,
doi=1,30
jx=jx+x(i)
jy=jy+y(i)
enddo
doi=1,30
s=s+(x(i)-jx/30)*(y(i)-jy/30)
jfx=jfx+(x(i)-jx/30)**2
jfy=jfy+(y(i)-jy/30)**2
enddo
r=s/30/sqrt(jfx/30*jfy/30)
print*,"r=",r
endfuntionarea
*实习四求给定数据的一元线性回归方程
利用下表数据,以环流指标为预报因子,气温为预报量,计算气温和环流指标之间的一元线性回归方程,并对回归方程进行检验。
年份
气温T
环流指标
1951
32
1952
25
1953
20
1954
26
1955
27
1956
24
1957
28
1958
0
24
1959
15
1960
16
1961
24
1962
30
1963
22
1964
30
1965
24
1966
33
1967
26
1968
20
1969
32
1970
35
一、实习结果
(1)用excel制作的气温-环流的医院线性回归方程
回归方程为:
(2)回归方程的检验
检验结果:
F=20.18>Fα,回归方程显著
二、Fortran程序(部分)
(1)回归方程的检验
!
回归方程显著性检验
!
计算两数组的距平及均方差
doi=1,m
var(i)=0
doj=1,nx
diff(j,i)=dat(j,i)-ave(i)
var(i)=var(i)+diff(j,i)**2
enddo
var(i)=sqrt(var(i)/nx)
enddo
!
计算协方差
i=1;t=0
doj=1,20
t=t+diff(j,i+1)*diff(j,i)
enddo
E=t/20
!
计算相关系数
r=E/(var
(1)*var
(2))
print*,'therelativevalueris:
',r
F=r**2*(nx-2)/(1-r**2)
实习五求给定数据的多元线性回归方程
说明:
x1-x4为四个预报因子,y为预报量;样本个数n=13
要求:
选取预报因子1、2、4,求预报量的标准化回归方程。
i
1
2
3
4
5
6
7
8
9
10
11
12
13
x1
7
1
11
11
7
11
3
1
2
21
1
11
10
x2
26
29
56
31
52
55
71
31
54
47
40
66
68
x3
6
15
8
8
6
9
17
22
18
4
23
9
8
x4
60
52
20
47
33
22
6
44
22
26
34
12
12
y
一、用excel制作的医院线性回归方程
标准化变量回归方程:
*实习七计算给定数据的11年滑动平均和累积距平
利用数据ma.dat,编写11点滑动平均的程序,ma.for给出了阅读资料的fortran程序。
数据在文件夹中单独给出。
要求:
实习报告中附出程序,并给出原数据和滑动后数据的图形(1张图)和累积距平数据图形(1张图)
一、实习结果
累积距平数据图形
滑动后数据
二、部分fortran程序:
programma
!
dimensionx(1000),x1(1000),nny1(1000)
!
write(*,10)
!
10format(5x,'n=?
ih=?
nyear=?
')
!
Read(*,*)n,ih,nyear
!
**********************************************
!
*n:
samplesizeofthetimeseries*
!
*ih:
movinglength*
!
*nyear:
firstyearoftheseries*
!
*x(n):
oroginaltimeseries*
!
*x1(n-ih+1):
movedseries*
!
**********************************************
integeri,n,ih,nyear
parameter(n=85,ih=11,nyear=1922)
30realx(n),x1(n-ih+1)
open(2,file='d:
\ma.dat')
read(2,*)(x(i),i=1,n)
close
(2)
doi=1,n-ih+1
x1(i)=sum(x(i:
i+ih-1))/ih
enddo
print*,x1
open(10,file='d:
\moveaverage.txt')
write(10,'(f10.6)')(x1(i),i=1,n-ih+1)
close(10)
end
*实习八对给定的海温数据进行EOF分析
给出海表温度距平数据资料sstpx.grd,以及相应的数据描述文件sstpx.ctl,对其进行EOF分析,资料的时空范围可以根据sstpx.ctl获知。
数据在文件夹中单独给出,距平或者标准化距平处理后再进行EOF。
Zhunsst.for给出了如何读取资料,
Ssteof.for为对距平或者标准化距平处理后的资料进行EOF分析。
要求:
实习报告中给出第一特征向量及其时间系数,并分析其时空特征。
由上图可看出,1957、1973、1985、1989年正异常较大,此时有ElNino出
现。
厄尔尼诺现象是发生在热带太平洋海温异常增暖的一种现象,大范围热带
太平洋增暖,会造成一些地区干旱,另一些地区又降雨过多的异常气候现象。
1956、1974、1988年负异常比较大,海面温度变低,有LaNina出现。
拉尼娜是指赤道太平洋东部和中部海面温度持续异常偏冷的现象(与厄尔尼
诺现象正好相反),也伴随着全球性气候混乱。
由图分析,ElNino和LaNina可能以一定的周期交替出现。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 气象 统计 实习 报告