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

气象统计方法实习报告

目录 实习一 求500hPa高度场气候场、距平场和均方差场3 1、资料介绍3 2．要求3 3、实习结果3 1)、FORTRAN源程序3 (2)、grads文件5 (3)、实习结果6 实习二 计算给定数据资料的简单相关系数和自相关系数11 1、资料介绍11 2、要求11 3、实习结果12 (1)、Fortran源程序12 (2)、程序运行结果：14 实习三 分析中国夏季降水线性趋势的分布特征15 1.资料介绍及要求：15 2.实习结果15 (1).matlab程序15 (2).程序运行结果16 实习四 求给定数据的一元线性回归方程17 1、资料介绍及要求17 2、实习结果18 (1)、MATLAB程序18 (2)、程序运行结果18 (3)、结果分析19 实习五 对给定的海温数据进行EOF分析21 1、资料介绍21 2、要求21 3、实习结果：21 (1)、FORTRAN源程序21 (2)空间场和时间序列的ctl文件23 (3)运行结果24 (4)分析26 实习三(附加) 计算给定数据的11年滑动平均和累积距平28 1、资料介绍28 2、要求28 3、实习结果28 实习四(附加) 求给定数据的多元线性回归方程31 1、说明31 2、要求31 3、实习结果：31 (1)Matlab源程序31 (2)运行结果35 实习一 求500hPa高度场气候场、距平场和均方差场 1、资料介绍 有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N. 时段：1982.1~1985.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_grads implicit none integer,parameter::nx=37,ny=17,nz=4,nt=12 integer i,j,iz,it real var(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt) ! Opening file open(10,file= g:\gradsdata\h500.dat ) do iz=1,nz do it=1,nt read(10,1000) read(10,3000) ((var(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo 1000 at(2i7) 2000 at(37f6.2) 3000 at(37f8.1) 4000 at(37f7.2) close(10) !Output open(16,file= g:\gradsdata\h500.grd ,= binary ) do iz=1,nz do it=1,nt write(16)((var(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo !Calculating the Climatological Field do it=1,nt do i=1,nx do j=1,ny sum=0 do iz=1,nz sum=sum+var(i,j,iz,it) enddo cl(i,j,it)=sum/4 enddo enddo enddo ! Output climate-file open(12,file= g:\gradsdata\climate.grd ,= binary ) do it=1,nt write(12)((cl(i,j,it),i=1,nx),j=1,ny) enddo !Calculating the Anomaly do iz=1,nz do it=1,nt do i=1,nx do j=1,ny jp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it) enddo enddo enddo enddo open(13,file= g:\gradsdata\anomaly.grd ,= binary ) !Output anomaly-file do iz=1,nz do it=1,nt write(13)((jp(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo !Calculating the Mean-square Deviation do it=1,nt do i=1,nx do j=1,ny sum=0 do iz=1,nz sum=sum+(jp(i,j,iz,it))**2 enddo jfc(i,j,it)=sqrt(sum/4) enddo enddo enddo !Output mean-square deviation-file open(14,file= g:\gradsdata\deviation.grd ,= binary ) do it=1,nt write(14)((jfc(i,j,it),i=1,nx),j=1,ny) enddo end (2)、grads文件 open g:\gradsdata\***.ctl (***为所求场对应的ctl文件名) set lat 0 40 set lon 60 150 set lev 500 enable print g:\gradsdata\***.gmf (***为所求场名称) i=1 while(i<=48(或 12)) set t i d h print c i=i+1 endwhile disable print ; (3)、实习结果 ①、原始场 1982年1月 1982年7月 结果分析： 冬季(