短期实习报告三.docVIP

实习三：大气遥相关 实验目的 掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。 实习要求 要求运用资料，计算北半球1月遥相关指数，并分析它与环流和我国气候变化的关系；用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确分析结果数据，完成实习报告。 实习内容： （1）计算EU（欧亚）遥相关指数，输出1月份该指数年际变化的时间序列； （2）计算EU遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数； （3）计算EU遥相关指数与同期我国气温的相关系数。 四、实验步骤 1、熟悉资料 （1）大气环流场资料： NCEP/NCAR 1948－2008年（61年）的500百帕月平均高度场资料 资料范围为（900S-900N,00-3600E) 网格距为2.50×2.50,纬向格点数为144，经向格点数为73 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 （2）站点资料： 我国气候中心整编的160站月平均气温资料； 全国160个台站； 所给的资料是1月份的； 资料为txt格式，参加资料说明。 2、方法简介： EU遥相关指数的定义: 相关系数的计算： 编写程序 （1）编程计算1月份EU遥相关指数（标准化），画出指数图； （2）编程计算1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形； （3）编程计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。 五、实验结果与分析 ①程序(实验一和实验二的程序内容，实验三与之类似，就不全部将程序全部展开赘述了） program ex3 parameter nt=61 real t(160,nt),eu(61),h(144,73,61),hsum(144,73),have(144,73),euave,rh(144,73),r(144,73),r2(160) real a(144,73,12,61),ave1(144,73),rup(144,73),rh_2(144,73),reu_2,eu2,reu,tave(160),tsum(160) real rh2(160),rh2_2(160),rup2(160) real lat(160),lon(160),tim character*8 id(160) integer nlev,nflag !读数据 open(3,file=e:\duanqi\shixi3\data\t1601.txt) open(2,file=e:\duanqi\shixi3\data\hgt500.grd,form=binary) open(4,file=e:\duanqi\shixi3\data\rheu.grd,form=binary) !h和eu的相关 open(5,file=e:\duanqi\shixi3\data\eu.grd,form=binary) !eu指数 open(6,file=e:\duanqi\shixi3\data\rteu.grd,status=replace,form=binary) !t和eu的相关 open(7,file=e:\duanqi\shixi3\data\lat_lon.txt) !读160站温度 read(3,*)((t(i,j),i=1,160),j=1,nt) do i=1,160 do j=1,nt t(i,j)=t(i,j)/10.0 enddo enddo !读经纬度 do i=1,160 read(7,*)lat(i),lon(i) enddo !读高度场 do it=1,61 do imo=1,12 do j=1,73 do i=1,144 read(2)a(i,j,imo,it) enddo;enddo;enddo;enddo !计算EU指数 do it=1,61 eu(it)=-0.25*a(9,59,1,it)+0.5*a(31,59,1,it)-0.25*a(59,53,1,it) enddo !!!!计算EU指数与高度场的相关系数 !h--高度场 hsum--高度场和 have--高度场平均值 eu--EU指数,eusum,euave类似 !1,提取1月份高度场 do it=1,61 do j=1,73 do i=1,144 h(i,j,it)=a(i,j,1,it) enddo enddo enddo !2,计算高度场和EU指数的平均值 do j=1,73 do i=1,144 do it=1,61 hsum(i,j)=hs