最新可编辑word文档最新可编辑word文档最新可编辑word文档气象统计实习报告专业：大气科学班级：xxxx级x班学号：2012130xxxx姓名：*实习一求500hPa高度场气候场、距平场和均方差场一、实习结果1981.1距平场1981.1500hpa高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。1981.1气候场1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。1981.1均方差场1981.1的均方差场在欧亚大陆的南侧有最小值，在10N~20N的南侧岁纬度而降低，在其北侧随纬度升高。二、相关的fortran程序，gs文件（1）Fortran程序:programmainparameter(nx=37,ny=17,mo=12,yr=4)realvar(nx,ny,mo,4)!数据realvars(nx,ny,mo)!4年气候态realjp(nx,ny,mo,4)!距平realfc(nx,ny,mo)!方差realjfc(nx,ny,mo)!均方差integeri,j,m,y,irecreal::summ=0.0!*********************求均方差*********************dom=1,12doi=1,37doj=1,17doy=1,4summ=summ+jp(i,j,m,y)**2enddofc(i,j,m)=summ/4!方差jfc(i,j,m)=sqrt(fc(i,j,m))!均方差summ=0.0enddoenddoenddo!******************************************open(11,file='d:\shixi\vars.grd',form='binary')open(12,file='d:\shixi\jp.grd',form='binary')open(13,file='d:\shixi\jfc.grd',form='binary')dom=1,12write(11)((vars(i,j,m),i=1,nx),j=1,ny)enddodom=1,12write(13)((jfc(i,j,m),i=1,nx),j=1,ny)enddodoy=1,4dom=1,12write(12)((jp(i,j,m,y),i=1,nx),j=1,ny)enddoenddoclose(11)close(12)close(13)End（2）Gs程序距平场'opene:\anomaly.ctl''setgxoutshaded''dxiaobao''setgxoutcontour''dxiaobao''drawtitle1982年1月''enableprinte:\anomaly.gmf''print''disableprint'2.气候场'opene:\climate.ctl''setgxoutshaded''dxiaobao''setgxoutcontour''dxiaobao''drawtitle1982.1''enableprinte:\climate.gmf''print''disableprint'均方差场'opene:\deviation.ctl''setgxoutshaded''dxiaolu''setgxoutcontour''dxiaolu''drawtitle1982.1deviation''enableprinte:\deviation.gmf''print''disableprint'*实习二计算给定数据资料的简单相关系数和自相关系数根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）一、实习结果计算简单相关系数计算出相关系数为r=0.4685170（2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）可以知道，年平均气温在滞后长度j=7，冬季j=4最大二、相关的fortran程序（部分）implicitnonerealx(20)!年平均气温realy(20)!冬季平均气温real::zx(10)=(/0,0,0,0,0,0,0,0,0,0/)real::rzx(10)=(/0,0,0,0,0,0,0,0,0,0/)real::s=0.0!协方差real::jx=0.0real::jy=0.0real::jfx=0.0rea