《气象统计方法实习报告》由会员分享，可在线阅读，更多相关《气象统计方法实习报告（49页珍藏版）》请在金锄头文库上搜索。
1、目录实习一 求 500hPa 高度场气候场、距平场和均方差场 31、资料介绍 32要求 43、实习结果 41）、FORTRAN 源程序 4（2）、grads 文件 9（3）、实习结果 10实习二 计算给定数据资料的简单相关系数和自相关系数 131、资料介绍 132、要求 133、实习结果 14（1）、Fortran 源程序 14（2）、程序运行结果： 19实习三 分析中国夏季降水线性趋势的分布特征 19191. 资料介绍及要求：2. 实习结果 19(1) .matlab 程序19(2) .程序运行结果20实习四 求给定数据的一元线性回归方程 201、资料介绍及要求 202、实习结果 22(1)
2、 、MATLAB 程序22(2) 、程序运行结果 23( 3 ) 、结果分析 25实习五 对给定的海温数据进行EOF分析261、资料介绍 262、要求 273、实习结果：27(1) 、FORTRAN 源程序27(2) 空间场和时间序列的ctl文件31323)运行结果（4）分析 33实习三（附加） 计算给定数据的 11 年滑动平均和累积距平 351、资料介绍 352、要求 353、实习结果 36实习四（附加） 求给定数据的多元线性回归方程 381、说明 382、要求 383、实习结果： 39（1）Matlab 源程序 39（2）运行结果 48实习一 求 500hPa 高度场气候场、距平场和均方差
3、场1、资料介绍有一 500hPa高度场资料，文件名，范围：60150E, 040N.时段：共 48 个月。水平分辨率：*，格点数：37*17。2要求编fortran程序，求500hPa高度场的（1）气候场；（2）距平场；（3）均方差场。并能用Grads做出图形，实习报告中气候场、距平场、均方差场任意给出两张图，图注要 清楚，即要注明是哪个时间的图形，并做简单分析。注：给出了如何用fortran读取ASCII码资料.3、实习结果1）、FORTRAN 源程序program ex_gradsimplicit noneinteger，parameter:nx=37，ny=17，nz=4，nt=12in
4、teger i，j，iz，itreal var（nx，ny，nz，nt），cl（nx，ny，nt），sum，jp（nx，ny，nz，nt），jfc（nx，ny，nt）open(10，file=g:gradsdata)do iz=1,nzdo it=1,ntread(10,1000)read(10,3000) (var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo1000format(2i7)2000format3000format4000formatclose(10)!Outputopen(16,file=g:gradsdata,form=binary)do iz=
5、1,nzdo it=1,ntwrite(16) (var(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!Calculating the Climatological Fielddo it=1,ntdo i=1,nxdo j=1,ny sum=0 do iz=1,nzsum=sum+var(i,j,iz,it)enddo cl(i,j,it)=sum/4 enddoenddoenddoopen(12,file=g:gradsdata,form=binary)do it=1,ntwrite(12) (cl(i,j,it),i=1,nx),j=1,ny)enddo!Cal
6、culating the Anomalydo iz=1,nzdo it=1,ntdo i=1,nxdo j=1,ny jp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it) enddoenddoenddoenddoopen(13,file=g:gradsdata,form=binary)do iz=1,nzdo it=1,ntwrite(13) (jp(i,j,iz,it),i=1,nx),j=1,ny)enddoenddo!Calculating the Mean-square Deviationdo it=1,ntdo i=1,nxdo j=1,ny sum=0 d
7、o iz=1,nz sum=sum+(jp(i,j,iz,it)*2 enddo jfc(i,j,it)=sqrt(sum/4) enddoenddoenddo!Output mean-square deviation-fileopen(14,file=g:gradsdata,form=binary)do it=1,ntwrite(14) (jfc(i,j,it),i=1,nx),j=1,ny)enddoend、grads文件open g:gradsda ta* .c tl(*为所求场对应的ctl文件名)set lat 0 40set lon 60 150set lev 500enable p
8、rint g:gradsdata*.gmf (*为所求场名称)i=1wh il e(i=48(或 12)set t id hprint， ，ci=i+1endwhiledisable print7（3）、实习结果 、原始场1982年1月1982年7月结果分析：冬季（此处以1月为代表）等高线分布整体平缓，表明高度场分布相对均匀，且北部接 近极地位势高度低，赤道地区位势高度高，这与太阳直射点在1月在南半球，使北半球整体 的辐射吸收随纬度增加而减小有关。北半球的气块受热随纬度递增而递减，因而膨胀率递 减，故南方接近赤道地区的气体膨胀大，位势高，而北方近极地地区气体膨胀小，位势低夏季（此处以7月为代表
9、），大洋上空出现副热带高压（588线位置），在东亚地区存在较为明显的位势高度槽，即东亚大槽。夏季在30 N以北的地区位于西风带中，从图中可看 出明显的西风带长波特征。除东亚大槽外，在中亚地区也存在一长波槽，这些槽线发生长 波调整时，会在部分地区发生较剧烈的天气变化。此外，从图中可以看出，1982年7月副热 带高压脊线的平均位置位于25 N,125 E附近，我国华东地区位于副高北侧西南气流控制 下，西南气流为水汽输送的主要通道，该地区发生降水较为频繁。从图中还可以看出冬夏季的位势高度分布存在明显差异，这与太阳直射点的年纪变化 密切相关。 、气候场2月8月结果分析：气候场是多年数据中同时段的平均值
10、序列，表征了区域内多年平均的位势高度变化。从气候场图形可以看出多年平均的500hPa高度场中，冬季（此处以多年平均的2月气候 场为代表）等高线较平直，大陆上等高线较稀疏，而海洋上等高线较密集，这表明大陆上 空的位势高度变化率较海洋小。此外，冬季等位势高度线分布平直，还说明冬季的天气变 化显着。夏季（此处以多年平均的8月气候场为代表）在太平洋上有副热带高压，副高是深厚的 系统，所以可以一直延伸到500hPa高度。东亚地区存在明显的西风带长波槽，即东亚大槽, 东亚大槽的移动和变化配合副热带高压和夏季风的影响，会使我国大部分地区天气造成强 烈变化，如形成大范围降水或强对流天气。同时，孟加拉湾处存在一
11、低涡。由于高原的阻 挡作用，这一系统对我国的影响并不显着。图中副高脊线8月的多年平均位置位于30 N以 北，且东部长波槽位于110 E附近，故由气候场分析，华北地区位于长波槽前，又槽线受到副高的阻挡作用，因而华北地区容易形成降水。 、距平场1982年1月1982年7月结果分析：距平场指示了位势高度的震荡趋势，因距平的平均值为0，则大于0的值表明位势高度 偏高，小于0的值表示位势高度偏低。从图中可看出冬季（此处以1982年1月距平场做代表）在大陆位势高度为正距平，而在 大洋则存在明显的负距平。则由距平场的性质得，冬季在大洋上位势高度偏高，在大陆上 位势高度偏低。其原因是，海水的比热大于陆地，则冬
12、季海洋温度比陆地高，所以海洋上 气块膨胀更多，位势高度更高。夏季（此处以1982年7月距平场做代表）相反，大洋上位势高度偏低，而陆地偏高。在 70 E90 E附近位势偏高的原因在于夏季青藏高原接受太阳辐射使之在对流层中层形成 热源，位势高度因比大陆其他位置高。由此可见，位势高度的变化不仅与海陆差异有关， 同时也与地形有关。在海洋上副热带高压所在的位置存在证据平值。 、均方差场6月12月 结果分析：均方差场反应同一时段内的位势高度变化幅度的大小。由图可以看出，整体位势高度在大陆上的变化幅度比海洋小， 且海洋上冬季的变化幅 度比夏季大，而陆地上相反，冬季的变化幅度比夏季小。因为陆地的比热小，所以陆
13、地在 夏季白天与夜间的温差大于冬季，对应的高度场震荡就比冬季剧烈。海洋上的位势高度变 化幅度的影响方面温差为次要因素，其主要受到副热带高压，西风带长波槽脊影响，位势 高度根据天气系统的移动而变化，所以震荡幅度较大。除受到天气系统影响外，海洋上的 位势高度场还受到洋流等因素的影响。实习二 计算给定数据资料的简单相关系数和自相关系数1、资料介绍根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度T最大取10）2、要求实习报告中附出简单相关系数或自相关系数程序。答案：r=年平均气温在滞后长
14、度j=3、冬季序列在j=4最大。3、实习结果(1)、Fortran 源程序PROGRAM EXAMIMPLICIT NONEINTEGER,PARAMETER:N=20INTEGER i,j,k,ty,tw,tywREAL:avr_y=0,avr_w=0,sy=0,sw=0,rxy=0,max_y=0,max_w=0,max_yw=0REAL y(N),w(N)DATAy/,)/，DATAw/,REAL syy(N),sww(N),r(N),rty(N),rtw(N),rtyw(N),rxy_ty(N),rxy_tw(N),rxy_tyw(N) !求两数组平均值DO i=1,Navr_y=avr_y+y(i)avr_w=avr_w+w(i)END DOavr_y=avr_y/Navr_w=avr_w/N!简单相关系数DO j=1,Nsyy(j)=(y(j)-avr_y)*2 sy=sy+syy(j)sww(j)=(w(j)-avr_w)*2sw=sw+sww(j)END DOsy=sqrt(sy/N)sw=sqrt(sw/N)DO j=1,Nr(j)=(y(
