气象统计方法实习报告.docVIP

目录 实习一 求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 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 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) close(10) !Output open(16,file=g:\gradsdata\h500.grd,form=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,form=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