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

目录 实习一 求 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.Fo 给出了如何用 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 fil 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) (