2020气象统计方法实习报告.docx

目录 目录 TOC \o "1-5" \h \z 实习一求500hPa高度场气候场、距平场和均方差场 3 1、 资料介绍 - 3 2．要求 - 4 3、实习结果 - 4 1）、FORTRAW程序 4 （2） 、 grads 文件- 9 （3） 、实习结果 10 实习二 计算给定数据资料的简单相关系数和自相关系数 13 1 、资料介绍 - 13 2、 要求 - 13 3、 实习结果 - 14 （1） 、 Fortran 源程序 - 14 （2） 、程序运行结果： 19 实习三 分析中国夏季降水线性趋势的分布特征 19 19资料介绍及要求： 19 TOC \o "1-5" \h \z 实习结果 - 19 (1).matlab 程序 - 19 ( 2) . 程序运行结果 21 实习四 求给定数据的一元线性回归方程 21 1、资料介绍及要求 - 21 2、实习结果 - 22 (1)、MATLAB?序 22 ( 2)、程序运行结果 24 ( 3)、结果分析 26 实习五 对给定的海温数据进行 EOF分析 27 1 、资料介绍 - 27 2、要求 - 28 3、实习结果： - 28 (1)、FORTRA?程序 28 34( 2)空间场和时间序列的 ctl 文件 32 34 3)运行结果 TOC \o "1-5" \h \z （4）分析 35 实习三（附加） 计算给定数据的 11 年滑动平均和累积距平 36 1、资料介绍 - 36 2、要求 - 37 3、实习结果 - 38 实习四（附加） 求给定数据的多元线性回归方程 40 1、说明 - 40 2、要求 - 40 3、实习结果： - 41 （1）Matlab 源程序 - 41 （2）运行结果 50 实习一 求 500hPa 高度场气候场、距平场和均方差场 1、资料介绍 有一 500hPa高度场资料，文件名 h500.dat，范围：60~150E 0~40N. 时段： 1982.1~1985.12 共48个月。水平分辨率： 2.5*2.5 ，格点数： 37*17 ! ! Opening file 2．要求 编 fortran 程序，求 500hPa 高度场的 (1)气候场； (2)距平场； 均方差场。 并能用Grads做出图形，实习报告中气候场、距平场、均方差场任意给出 两张图，图注要 清楚，即要注明是哪个时间的图形，并做简单分析。 注：h500.For给出了如何用fortran 读取ASCII码资料h500.dat. 3、实习结果 1)、FORTRA?程序 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) open(10,file='g:\gradsdata\h500.dat') write(16) ((var(i,j,iz,it),i=1,nx),j=1,ny) write(16) ((var(i,j,iz,it),i=1,nx),j=1,ny) 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 ! ! Output climate-file 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 open(12,file='g:\gradsdata\climate.grd',form='binary') !Output anomaly-file !Output anomaly-file do it=1,nt write(12)