matlab 累积距平,气象统计方法实习报告

《气象统计方法实习报告》由会员分享，可在线阅读，更多相关《气象统计方法实习报告(36页珍藏版)》请在人人文库网上搜索。

1、目录实习一 求500hPa高度场气候场、距平场和均方差场31、资料介绍32要求33、实习结果31)、FORTRAN源程序3(2)、grads文件5(3)、实习结果6实习二 计算给定数据资料的简单相关系数和自相关系数111、资料介绍112、要求113、实习结果12(1)、Fortran源程序12(2)、程序运行结果：14实习三 分析中国夏季降水线性趋势的分布特征151.资料介绍及要求：152.实习结果15(1).matlab程序15(2).程序运行结果16实习四 求给定数据的一元线性回归方程171、资料介绍及要求172、实习结果18(1)、MATLAB程序18(2)、程序运行结果18(3)、结果。

2、分析19实习五 对给定的海温数据进行EOF分析211、资料介绍212、要求213、实习结果：21(1)、FORTRAN源程序21(2)空间场和时间序列的ctl文件23(3)运行结果24(4)分析26实习三(附加) 计算给定数据的11年滑动平均和累积距平281、资料介绍282、要求283、实习结果28实习四(附加) 求给定数据的多元线性回归方程311、说明312、要求313、实习结果：31(1)Matlab源程序31(2)运行结果35实习一 求500hPa高度场气候场、距平场和均方差场 1、资料介绍有一500hPa高度场资料，文件名h500.dat，范围：60150E，040N.时段：1982.。

3、11985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17。2要求编fortran程序，求500hPa高度场的(1)气候场；(2)距平场；(3)均方差场。并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.3、实习结果 1)、FORTRAN源程序program ex_gradsimplicit noneinteger,parameter:nx=37,ny=17,nz=4,nt=12integer i,j,iz,itrea。

4、l var(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt)! Opening fileopen(10,file=g:gradsdatah500.dat)do iz=1,nzdo it=1,ntread(10,1000)read(10,3000) (var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo1000 format(2i7)2000 format(37f6.2)3000 format(37f8.1)4000 format(37f7.2)close(10)!Outputopen(16,fil。

5、e=g:gradsdatah500.grd,form=binary)do iz=1,nzdo it=1,ntwrite(16)(var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!Calculating the Climatological Fielddo it=1,ntdo i=1,nxdo j=1,nysum=0do iz=1,nzsum=sum+var(i,j,iz,it)enddocl(i,j,it)=sum/4enddoenddoenddo! Output climate-fileopen(12,file=g:gradsdataclimate.grd,f。

6、orm=binary)do it=1,ntwrite(12)(cl(i,j,it),i=1,nx),j=1,ny)enddo!Calculating the Anomalydo iz=1,nzdo it=1,ntdo i=1,nxdo j=1,nyjp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it)enddoenddoenddoenddoopen(13,file=g:gradsdataanomaly.grd,form=binary)!Output anomaly-filedo iz=1,nzdo it=1,ntwrite(13)(jp(i,j,iz,it),i=1,n。

7、x),j=1,ny)enddoenddo!Calculating the Mean-square Deviationdo it=1,ntdo i=1,nxdo j=1,nysum=0do iz=1,nzsum=sum+(jp(i,j,iz,it)*2enddojfc(i,j,it)=sqrt(sum/4)enddoenddoenddo!Output mean-square deviation-fileopen(14,file=g:gradsdatadeviation.grd,form=binary)do it=1,ntwrite(14)(jfc(i,j,it),i=1,nx),j=1,ny)e。

8、nddo end(2)、grads文件open g:gradsdata*.ctl (*为所求场对应的ctl文件名)set lat 0 40set lon 60 150set lev 500enable print g:gradsdata*.gmf (*为所求场名称)i=1while(imax_y) THENmax_y=rxy_ty(ty)k=tyEND IFEND DOPRINT (全年平均气温绝对值最大自相关系数rx