气象统计实习报告.docx

1、气象统计实习报告实习一 求 500hPa高度场气候场、距平场和均方差场一、实习内容有一 500hPa高度场资料， 文件名 h500.dat， 范围： 60150E， 040N.时段： 1982.11985.12共 48 个月。水平分辨率： 2.5*2.5，格点数：37*17。编 fortran 程序，求 500hPa高度场的( 1)气候场；( 2)距平场；( 3)均方差场。并能用 Grads做出图形 ,实习报告中气候场、 距平场、 均方差场任意给出 两张图 ，图注要清楚，即要注明是哪个时间的图形，并做简单分析。二、 fortran 程序program mainimplicit noneinte

2、ger ii,jj,mo,yrparameter(ii=37,jj=17,mo=12,yr=4)real var(ii,jj,mo,yr)real climate(ii,jj,mo),juping(ii,jj,mo,yr),junfang(ii,jj,mo)integer i,j,m,iyopen(5,file=d:1h500.dat)open(1,file=d:1h500aver.grd,form=binary)open(2,file=d:1h500juping.grd,form=binary)open(3,file=d:1h500junfang.grd,form=binary)do iy=

3、1,4 !读数据资料 h500do m=1,12read(5,1000)read(5,3000) (var(i,j,m,iy),i=1,ii),j=1,jj)End doEnd doformat(2i7)FORMAT(37F6.2)format(37f8.1)format(37f7.2)do m=1,mo !气候场do j=1,jjdo i=1,iiclimate(i,j,m)=0do iy=1,yrclimate(i,j,m)=climate(i,j,m)+var(i,j,m,iy)end doclimate(i,j,m)=climate(i,j,m)/ 4end doend doend d

4、owrite(1)(climate(i,j,m),i=1,ii),j=1,jj),m=1,mo)do i=1,ii ! 距平场do j=1,jjdo m=1,modo iy=1,yrjuping(i,j,m,iy)=var(i,j,m,iy)-climate(i,j,m)end doend doend doend dowrite(2)(juping(i,j,m,iy),i=1,ii),j=1,jj),m=1,mo),iy=1,yr)do m=1,mo ! 均方差场do i=1,iido j=1,jjjunfang(i,j,m)=0do iy=1,yrjunfang(i,j,m)=junfang

5、(i,j,m)+juping(i,j,m,iy)*2end dojunfang(i,j,m)=sqrt(junfang(i,j,m)/yr)end doend doend dowrite(3)(junfang(i,j,m),i=1,ii),j=1,jj),m=1,mo)End1000200030004000三、图形显示1 、气候场：2、距平场：T=1982.7分析：对于 7 月份而言， 1982 年的高度场偏低，大部分地区为负距平。实习二 计算给定数据资料的简单相关系数和自相关系数一、实习内容根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。2）分别

6、找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度 最大取 10）二、 fortran 程序program mainimplicit noneinteger i,j,maxireal:k=0,sumx=0,sumy=0,xa,ya,sx,sy,r1,y(20),x(20),rx(10),ry(10)real:data1(20)=(/ 3.4,3.3,3.2,2.9,3.4,2.8,3.6,3.0,2.8,3.0,3.1,3.0,2.9,2.7,3.5,3.2,3.1,2.8,2.9,2.9/) real:data2(20)=(/ 3.24,3.14,3.26,2.38,3.32,2

7、.71,2.84,3.94,2.75,1.83,2.8,2.81,2.63,3.2,3.6,3.4,3.07, 1.87,2.63,2.47/)do i=1,20sumx=sumx+data1(i)sumy=sumy+data2(i)end doxa=sumx/20;ya=sumy/20sumx=0;sumy=0do i=1,20sumx=sumx+(data1(i)-xa)*2 ； sumy=sumy+(data2(i)-ya)*2 end dosx=sqrt(sumx/20);sy=sqrt(sumy/20)sumx=0;sumy=0do i=1,20x(i)=(data1(i)-xa)/

8、sx； y(i)=(data2(i)-ya)/syend dodo i=1,20sumx=sumx+x(i)*y(i)end dor1=sumx/20;sumx=0print *, 简单相关系数为 ,r1do j=1,10do i=1,20-jsumx=sumx+x(i)*x(i+j) ； sumy=sumy+y(i)*y(i+j)end dorx(j)=abs(sumx/(20-j) ； ry(j)=abs(sumy/(20-j)sumx=0;sumy=0end domaxi=1do i=1,10if(rx(i)=rx(maxi) thenmaxi=iend ifend doprint *,

9、 年平均气温对应的 j 为 ,maximaxi=1do i=1,10if(ry(i)=ry(maxi) thenmaxi=iend ifend doprint *, 冬季平均气温对应的 j为 ,maxiend三、运行结果显示实习三 计算给定数据的落后交叉相关系数和偏相关系数一、实习内容根据下表 北京冬季（ 12 月 2 月）气温资料 计算： 12 月气温与 1 月和 2 月气温的落后交叉相关系数 （滞后长度 最大取 10） 和偏相关系数。 在实习报告中给出程序。年份12月1月2月19511.0-2.7-4.31952-5.3-5.9-3.51953-2.0-3.4-0.81954-5.7-4.

10、7-1.11955-0.9-3.8-3.11956-5.7-5.3-5.91957-2.1-5.0-1.619580.6-4.30.21959-1.7-5.72.01960-3.6-3.61.31961-3.0-3.1-0.819620.1-3.9-1.11963-2.6-3.0-5.21964-1.4-4.9-1.71965-3.9-5.7-2.51966-4.7-4.8-3.31967-6.0-5.6-4.91968-1.7-6.4-5.11969-3.4-5.6-2.01970-3.1-4.2-2.91971-3.8-4.9-3.91972-2.0-4.1-2.41973-1.7-4.2

11、-2.01974-3.6-3.3-2.01975-2.7-3.70.11976-2.4-7.6-2.21977-0.9-3.5-2.31978-2.7-4.2-0.51979-1.6-4.5-2.91980-3.9-4.8-1.4二、 fortran 程序program mainimplicit noneinteger iy,t,j,ireal:var(30,3),s12(10),s13(10),s1=0,s2=0,s3=0,a1,a2,a3,sum1=0,sum2=0,sum3=0,r12(10),r13(10),p12,p13,t12,t13,t23open(18,file=d:grads

12、shiyan3.txt) ! 读入数据read(18,*) (var(iy,t),t=1,3),iy=1,30)do iy=1,30 !求均值和标准差s1=s1+var(iy,1)s2=s2+var(iy,1)s3=s3+var(iy,1)end doa1=s1/30； a2=s2/30； a3=s3/30do iy=1,30sum1=sum1+(var(iy,1)-a1)*2sum2=sum2+(var(iy,2)-a2)*2sum3=sum3+(var(iy,3)-a3)*2end dos1=sqrt(sum1/30);s2=sqrt(sum2/30);s3=sqrt(sum3/30);s

13、um1=0;sum2=0;sum3=0do j=1,10 !求 12月和 1 月落后交叉相关系数 r12 以及 12 月和 2 月系数 r13do i=1,30-jsum1=sum1+(var(i,1)-a1)*(var(i+j,2)-a2)sum2=sum2+(var(i,1)-a1)*(var(i+j,3)-a3)end dor12(j)=sum1/(30-j)/s1 /s2； r13(j)=sum2/(30-j)/s1 /s3sum1=0;sum2=0end dosum1=0;sum2=0;sum3=0print *,12 月和 1 月落后交叉相关系数 r12(j) !输出 12 月和 1 月落后交叉相关系数 r12(j)以及 12 月和 2 月系数 r13(j)do j=1,10print *,j=,j, 时 ,r12=,r12(j)end doprint *print *,12 月和 2 月落后交叉相关系数 r13(j)do j=1,10print *,j=,j, 时 ,r13=,r13(j)end doprint *do i=1,30 !求 12 月和 1 月简单相关系数 t12,求 12月和 2 月 t13,以及 1 月和 2 月 t23sum1=sum1+(var(i,1)