短期气候预测实习报告三（6页）_1

1、短期气候预测实习报告三（6页）短期气候预测实_报告三（6页） 短期气候预测实_报告三 南京信息工程大学实验（实_）报告实验（实_）名称 短期气候预测实_ 日期 11.12 得分 指导教师 系 专业 年级 班次 姓名 Trichtu 学号 实_目的：掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。实_要求：要求运用资料，计算北半球 1 月遥相关指数，并分析它与环流和我国气候变化的关系； 用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确 分析结果数据，完成实_报告。实_内容：计算EU （欧亚）遥相关指数，输出 1月份该指数年际变化的时间序列；计算E

2、U遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数； 计 算EU遥相关指数与同期我国气温的相关系数。实_资料：（大气环流场资料）NCEP/NCAR 1948 xxxx年（61 年）的 500百帕月平均高度场资料资料范围为（ 900S-900N,00-3600E）网格距为2.50 X 2.50,纬向格点数为144,经向格点数为73资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。（ 站点资料 ） ：我国气候中心整编的 160站月平均气温资料；全国 160个台站,所给的资料是 1 月份的,资料为 txt 格式,参加资料说明。实_方法简介：实_步骤：编程

3、计算1月份EU遥相关指数（标准化），画出指数图；编程计算1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形; program calculationreal ind(61),h500(144,73,12,61),r(144,73)real sumx,sumy,xx,xy,yyopen(2,file=e:copy3datahgt500.grd,form=binary)open(3,file=e:copy3ind.grd,form=binary)open(4,file=e:copy3corr.grd,form=binary)!read data 1948-xxxxdo it=1,

4、61do k=1,12do j=1,73do i=1,144read(2) h500(i,j,k,it)enddoenddoenddoenddowrite(*,*)reading is ok!calculationdo it=1,61ind(it)=-0.25*h500(9,59,1,it)+0.5*h500(31,59,1,it)-0.25*h500(59,53,1,it)end dowrite(*,*)calculation is ok!correlationsumy=.0yy=.0do it=1,61sumy=sumy+ind(it)yy=yy+ind(it)*2enddodo i=1,

5、144do j=1,73sumx=.0xx=.0xy=.0do it=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) end doend dowrite(*,*)corralation is ok!writingdo it=1,61write(3) ind(it)enddo!w

6、rite(*,*) r(2,3)write(4) (r(i,j),i=1,144),j=1,73)write(*,*)writing is okEnd编程计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形program EPreal a(160,60),ind(61),r(160),aa(61),bb(61),lat(160),lon(160)character*8 id(160)real avex,avey,xx,xy,yyopen(2,file=e:copy3datat1601.txt) open(3,file=e:copy3ind.grd,form=binary) open(

7、4,file=e:copy3datalat_lon.txt) open(5,file=e:copy3r160.grd,form=binary)ccccccccccccccc 读数据(指数、经纬度、 160 站温度)ccc a:160 站气温 (1951xxxx年) ind ：指数序列 (1948xxxx 年)read(2,*)(a(i,j),i=1,160),j=1,60)close(2)do it=1,60read(3) ind(it)enddodo is=1,160read(4,*) lat(is),lon(is)enddocccccccccccccc 编程求相关 年份不同sumy=.0y

8、y=.0 do it=1,58sumy=sumy+ind(it+3) yy=yy+ind(it+3)*2enddodo is=1,160sumx=.0xx=.0xy=.0do it=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) end do write(*,*)corralation is okccccccccccccccccccc 写站点数据do j=1,160id(j)=cha

9、r(j)tim=0.0nlev=1nflag=1write(5)id(j),lat(j),lon(j),tim,nlev,nflag,r(j)enddotim=0.0nlev=0nflag=1write(5)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagend实_结果 :1 月份 EU 遥相关指数，画出指数图；可以看出EU指数随时间的波动，在 1972、1973、1984、1997、xxxx出现低值， 1963、 1976、 1981、 1984、 1986、 1995、 201 9出现高值。1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形；从图中可以看出EU指数和亚种中部、北美洲北部及大西洋北部有显著正相关，与欧洲、亚洲东部、太平洋北部及北美洲南部有负相关。这些系统基本为500hPa的永久性或半永久性环流中心。而与南半球的相关性不显著也反映了这个指数是来自于北半球的资料。该指数反映了北半球500hPa的环流中心的强弱。EU指数越强，则位于亚洲中部的位势高 度越高，即亚洲低压系统较弱，而北太平洋上的位势越低，则反映副高的减弱。因此此图 很好的反映出了两槽一脊的强度变化，即东亚大槽和欧洲东部的槽加深时，西伯利亚附近 的脊发展。计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。从图中可以看出，E