气象统计实习报告.docVIP

PAGE PAGE 1 气象统计实习报告 专业：大气科学 班级：xxxx级x班 学号：2012130xxxx 姓名： *实习一 求500hPa高度场气候场、距平场和均方差场 一、实习结果 1981.1距平场 1981.1 500hpa高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。 1981.1气候场 1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。 1981.1均方差场 1981.1的均方差场在欧亚大陆的南侧有最小值，在10N~20N的南侧岁纬度而降低，在其北侧随纬度升高。 二、相关的fortran程序，gs文件 （1）Fortran程序: program main parameter(nx=37,ny=17,mo=12,yr=4) real var(nx,ny,mo,4)!数据 real vars(nx,ny,mo)!4年气候态 real jp(nx,ny, mo,4)!距平 real fc(nx,ny,mo)!方差 real jfc(nx,ny,mo)!均方差 integer i,j,m,y,irec real::summ=0.0 !*********************求均方差********************* do m=1,12 do i=1,37 do j=1,17 do y=1,4 summ=summ+jp(i,j,m,y)**2 enddo fc(i,j,m)=summ/4!方差 jfc(i,j,m)=sqrt(fc(i,j,m))!均方差 summ =0.0 enddo enddo enddo !****************************************** 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 ) do m=1,12 write(11)((vars(i,j,m),i=1,nx),j=1,ny) enddo do m=1,12 write(13)((jfc(i,j,m),i=1,nx),j=1,ny) enddo do y=1,4 do m=1,12 write(12)((jp(i,j,m,y),i=1,nx),j=1,ny) enddo enddo close(11) close(12) close(13) End （2）Gs程序 距平场 open e:\anomaly.ctl set gxout shaded d xiaobao set gxout contour d xiaobao draw title 1982年1月 enable print e:\anomaly.gmf print disable print 2.气候场 open e:\climate.ctl set gxout shaded d xiaobao set gxout contour d xiaobao draw title 1982.1 enable print e:\climate.gmf print disable print 均方差场 open e:\deviation.ctl set gxout shaded d xiaolu set gxout contour d xiaolu draw title 1982.1 deviation enable print e:\deviation.gmf print disable print *实习二 计算给定数据资料的简单相关系数和自相关系数 根据下表中年平均气温和冬季平均气温的等级数据进行下列计算： 1）计算两个气温之间的简单相关系数。 2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10） 一、实习结果 计算简单相关系数 计算出相关系数为r=0.4685170 （2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10） 可以知道，年平均气温在滞后长度j=7，冬季j=4最大 二、相关的fortran程序（部分） implicit non