资源描述：
《气象统计方法实习报告.》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、36气象统计方法课程实践内容2013目录实习一求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第36页共36页36气
2、象统计方法课程实践内容2013（1）、MATLAB程序18（2）、程序运行结果18（3）、结果分析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第36页共36页36气象统计方法课程实践内容2013实习一求50
3、0hPa高度场气候场、距平场和均方差场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.第36页共36页36气象统计方法课程实践内容201
4、33、实习结果1）、FORTRAN源程序programex_gradsimplicitnoneinteger,parameter::nx=37,ny=17,nz=4,nt=12integeri,j,iz,itrealvar(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt)!Openingfileopen(10,file='g:gradsdatah500.dat')doiz=1,nzdoit=1,ntread(10,1000)read(10,3000)((var(i,j,iz,it),i=1,nx),j=1,n
5、y)enddoenddo1000format(2i7)2000format(37f6.2)3000format(37f8.1)4000format(37f7.2)close(10)!Outputopen(16,file='g:gradsdatah500.grd',form='binary')doiz=1,nzdoit=1,ntwrite(16)((var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!CalculatingtheClimatologicalFielddoit=1,ntdoi=1,nxdoj=1,nysum=0doiz=1,nzsum=
6、sum+var(i,j,iz,it)enddocl(i,j,it)=sum/4enddoenddoenddo第36页共36页36气象统计方法课程实践内容2013!Outputclimate-fileopen(12,file='g:gradsdataclimate.grd',form='binary')doit=1,ntwrite(12)((cl(i,j,it),i=1,nx),j=1,ny)enddo!CalculatingtheAnomalydoiz=1,nzdoit=1,ntdoi=1,nxdoj=1,nyjp(i,j,iz,it)=var(i,j,iz,it)-cl(i,
7、j,it)enddoenddoenddoenddoopen(13,file='g:gradsdataanomaly.grd',form='binary')!Outputanomaly-filedoiz=1,nzdoit=1,ntwrite(13)((jp(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!CalculatingtheMean-squareDeviationdoit=1,ntdoi=1,nxdoj=1,nysum=0doiz=
