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

气象统计实习报告 专业：大气科学 班级： xxxx 级 x 班 学号： 2012130xxxx 实习一 求500hPa高度场气候场、距平场和均方差场 、实习结果 1981.1距平场 1981.1 500hpa高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。 1981.1气候场 1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。 I982J deviation 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,m o,4)! 数据 real vars( nx, ny,mo)!4 年气候态 real jp( nx,ny, mo,4)! 距平 real fc(nx,n y,mo)! 方差 real jfc(n x, ny,mo)! 均方差 in teger i,j,m,y,irec real::summ=0.0 I *********************求均方差********************* I ********************* 求均方差 ********************* do m=1,12 do i=1,37 do j=1,17 do y=1,4summ=summ+jp(i,j,m,y)**2 r r 、 一 、/?- r r 、 一 、/?- 均方差 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 气候场 open e:\climate.ctl set gxout shaded d xiaobao set gxout con tour d xiaobao draw title 1982.1 en able print e:\climate.gmf print disable print 均方差场 ope n e:\deviati on. ctl set gxout shaded d xiaolu set gxout con tour d xiaolu draw title 1982.1 deviatio n en able print e:\deviati on. gmf print disable print *实习二计算给定数据资料的简单相关系数和自相关系数 根据下表中年平均气温和冬季平均气温的等级数据进行下列计算： 1） 计算两个气温之间的简单相关系数。 2） 分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。 （滞后长度T最大取10） 中国1970-1989年年平卿冬季网關 197M979 140 no 120 290 140 280 360 3.00 180 3.00 1980-1989 3.W 100 2.90 2.70 3,50 3.20 □ 2.80 2.90 2.90 縛平帖温 197M979 3.24 3.14 3.26 2.38 332 2.71 2.84 3.94 2.75 1.832加 2.81 2.63 320 3.60 3.07 1.87 2.63 2.47 、实习结果 (1)计算简单相关系数 ■ *C：\05er$\AdministratQr,\De!Ektop\z