短期气候预测实习报告三

南京信息工程大学实验（实习）报告实验（实习）名称短期气候预测实习日期11.12得分指导教师系 ~专业~年级~班次~姓名Trichtu学号~实习目的：掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。实习要求：要求运用资料，计算北半球1月遥相关指数，并分析它与环流和我国气候变化的关系；用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确分析结果数据，完成实习报告。实习内容：计算EU（欧亚）遥相关指数，输出1月份该指数年际变化的时间序列；计算EU遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数；计算EU遥相关指数与同期我国气温的相关系数。实习资料：（大气环流场资料）NCEP/NCAR1948－2008年（61年）的500百帕月平均高度场资料资料范围为（900S-900N,00-3600E)网格距为2.50×2.50,纬向格点数为144，经向格点数为73资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。(站点资料)：我国气候中心整编的160站月平均气温资料；全国160个台站，所给的资料是1月份的，资料为txt格式，参加资料说明。实习方法简介：EU遥相关指数的定义:相关系数的计算：实习步骤：编程计算1月份EU遥相关指数（标准化），画出指数图；编程计算1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形；programcalculationrealind(61),h500(144,73,12,61),r(144,73)realsumx,sumy,xx,xy,yyopen(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')!readdata1948-2008doit=1,61dok=1,12doj=1,73doi=1,144read(2)h500(i,j,k,it)enddoenddoenddoenddowrite(*,*)'readingisok'!calculationdoit=1,61ind(it)=-0.25*h500(9,59,1,it)+0.5*h500(31,59,1,it)-0.25*h500(59,53,1,it)enddowrite(*,*)'calculationisok'!correlationsumy=.0yy=.0doit=1,61sumy=sumy+ind(it)yy=yy+ind(it)**2enddodoi=1,144doj=1,73sumx=.0xx=.0xy=.0doit=1,61sumx=sumx+h500(i,j,1,it)xx=xx+h500(i,j,1,it)**2xy=xy+h500(i,j,1,it)*ind(it)enddo!write(*,*)61*yy-sumy**2,61*xx-sumx**2r(i,j)=(xy*61-sumx*sumy)/(sqrt(61*xx-sumx**2)*sqrt(61*yy-sumy**2))enddoenddowrite(*,*)'corralationisok'!writingdoit=1,61write(3)ind(it)enddo!write(*,*)r(2,3)write(4)((r(i,j),i=1,144),j=1,73)write(*,*)'writingisok'End编程计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。programEP reala(160,60),ind(61),r(160),aa(61),bb(61),lat(160),lon(160) character*8id(160)realavex,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.grd',form='binary')ccccccccccccccc读数据（指数、经纬度、160站温度）ccca:160站气温(1951~2010年)ind：指数序列(1948~2008年)read(2,*)((a(i,j),i=1,160),j=1,60) close(2) doit=1,60 read(3)ind(it) enddodois=1,160 read(4,*)lat(is),lon(is) enddocccccccccccccc编程求相关年份不同sumy=.0yy=.0doit=1,58sumy=sumy+ind(it+3)yy=yy+ind(it+3)**2enddodois=1,160sumx=.0xx=.0xy=.0doit=1,58sumx=sumx+a(is,it)xx=xx+a(is,it)**2xy=xy+a(is,it)*ind(it+3)enddor(is)=(xy*58-sumx*sumy)/(sqrt(58*xx-sumx**2)*sqrt(58*yy-sumy**2))enddowrite(*,*)'corralationisok'ccccccccccccccccccc写站点数据 doj=1,160 id(j)=char(j) tim=0.0 nlev=1 nflag=1 write(5)id(j),lat(j),lon(j),tim,nlev,nflag,r(j) enddo tim=0.0 nlev=0 nflag=1 write(5)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflag end实习结果:1月份EU遥相关指数，画出指数图；可以看出EU指数随时间的波动，，在1972、1973、1984、1997、2006出现低值，1963、1976、1981、1984、1986、1995、2004出现高值。1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形；从图中可以看出EU指数和亚种中部、北美洲北部及大西洋北部有显著正相关，与欧洲、亚洲东部、太平洋北部及北美洲南部有负相关。这些系统基本为500hPa的永久性或半永久性环流中心。而与南半球的相关性不显著也反映了这个指数是来自于北半球的资料。该指数反映了北半球500hPa的环流中心的强弱。EU指数越强，则位于亚洲中部的位势高度越高，即亚洲低压系统较弱，而北太平洋上的位势越低，则反映副高的减弱。因此此图很好的反映出了两槽一脊的强度变化，即东亚大槽和欧洲东部的槽加深时，西伯利亚附近的脊发展。计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。