短期气候预测实习报告三.docVIP

南京信息工程大学实验（实习）报告 实验（实习）名称 短期气候预测实习 日期11.12得分 指导教师 系 ~专业 ~ 年级 ~ 班次~姓名Trichtu 学号~ 实习目的： 掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。 实习要求： 要求运用资料，计算北半球1月遥相关指数，并分析它与环流和我国气候变化的关系；用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确分析结果数据，完成实习报告。 实习内容： 计算EU（欧亚）遥相关指数，输出1月份该指数年际变化的时间序列； 计算EU遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数； 计算EU遥相关指数与同期我国气温的相关系数。 实习资料： （大气环流场资料） NCEP/NCAR 1948－2008年（61年）的500百帕月平均高度场资料 资料范围为（900S-900N,00-3600E) 网格距为2.50×2.50,纬向格点数为144，经向格点数为73 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 (站点资料)： 我国气候中心整编的160站月平均气温资料； 全国160个台站，所给的资料是1月份的，资料为txt格式，参加资料说明。 实习方法简介： EU遥相关指数的定义: 相关系数的计算： 实习步骤： 编程计算1月份EU遥相关指数（标准化），画出指数图； 编程计算1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形； program calculation real ind(61),h500(144,73,12,61),r(144,73) real sumx,sumy,xx,xy,yy open(2,file=e:\copy3\data\hgt500.grd,form=binary) open(3,file=e:\copy3\ind.grd,form=binary) open(4,file=e:\copy3\corr.grd,form=binary) !read data 1948-2008 do it=1,61 do k=1,12 do j=1,73 do i=1,144 read(2) h500(i,j,k,it) enddo enddo enddo enddo write(*,*)reading is ok !calculation do it=1,61 ind(it)=-0.25*h500(9,59,1,it)+0.5*h500(31,59,1,it)-0.25*h500(59,53,1,it) end do write(*,*)calculation is ok !correlation sumy=.0 yy=.0 do it=1,61 sumy=sumy+ind(it) yy=yy+ind(it)**2 enddo do i=1,144 do j=1,73 sumx=.0 xx=.0 xy=.0 do it=1,61 sumx=sumx+h500(i,j,1,it) xx=xx+h500(i,j,1,it)**2 xy=xy+h500(i,j,1,it)*ind(it) enddo !write(*,*) 61*yy-sumy**2,61*xx-sumx**2 r(i,j)=(xy*61-sumx*sumy)/(sqrt(61*xx-sumx**2)*sqrt(61*yy-sumy**2)) end do end do write(*,*)corralation is ok !writing do it=1,61 write(3) ind(it) enddo !write(*,*) r(2,3) write(4) ((r(i,j),i=1,144),j=1,73) write(*,*)writing is ok End 编程计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。 program EP real a(160,60),ind(61),r(160),aa(61),bb(61),lat(160),lon(160) character*8 id(160) real avex,avey,xx,xy,yy open(2,file=e:\copy3\data\t1601.txt) open(3,file=e:\copy3\ind.grd,form=binary) open(4,file=e:\copy3\data\lat_lon.txt) open(5,file=e:\copy3\r160.gr