气象统计实习报告 (2)解读.docxVIP

气象统计实习报告 专业：大气科学 班级：xxxx 级 x 班学号：2012130xxxx 姓名： 1 *实习一 求 500hPa 高度场气候场、距平场和均方差场 一、实习结果 1981.1 距平场 1981.1 500hpa 高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。 2 1981.1 气候场 1982.1 的气候场呈现明显的沿纬度的上升而下降的趋势。 1981.1 均方差场 1981.1 的均方差场在欧亚大陆的南侧有最小值，在10N~20N 的南侧岁纬度而降低，在其北侧随纬度升高。 二、相关的 fortran 程序，gs 文件 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 3 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 enddo do m=1,12 enddo do y=1,4 write(11)((vars(i,j,m),i=1,nx),j=1,ny) write(13)((jfc(i,j,m),i=1,nx),j=1,ny) 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 Gs 程序 1.距平场 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 4 set gxout contour d xiaobao draw title 1982.1 enable print e:\climate.gmf print disable print 3.均方差场 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） 5 一、实习结果 计算简单相关系数 计算出相关系数为 r=0.4685170 分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度 τ最大取 10） 可以知道，年平均气温在滞后长度 j=7，冬季 j=4 最大 二、相关的 fortran 程序（部分） implicit none real x(20)!年平均气温real y(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.0 real ::jy=0.0 real ::jfx=0.0 real ::jfy=0.0,r inte