《短期气候预测实习报告三.docx》由会员分享，可在线阅读，更多相关《短期气候预测实习报告三.docx（4页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、南京信息工程大学实验（实习）报告实验（实习）名称 短期气候预测实习日期11.12得分 指导教师系二专业 二年级二班次二姓名Trichtu学号实习目的：掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。实习要求：要求运用资料，计算北半球1月遥相关指数，并分析它与环流和我国气候变化的关系；用图 形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确分析结 果数据，完成实习报告。实习内容：计算EU （欧亚）遥相关指数，输出1月份该指数年际变化的时间序列；计算EU遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数；计算EU遥相关指数与同期我国气温
2、的相关系数。实习资料：（大气环流场资料）NCEP/NCAR 1948-2008年（61年）的500百帕月平均高度场资料资料范围为（900S-900N,00-3600E）网格距为2.50X2.50,纬向格点数为144,经向格点数为73资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。（站点资料）：我国气候中心整编的160站月平均气温资料；全国160个台站，所给的资料是1月份的，资料为txt格式，参加资料说明。实习方法简介：EU遥相关指数的定义:EU = - z(55 N,20 E) + - z(55 N,75 E)-z(40 N,145 E)4,相关系数的计算：（为一幻
3、（口-刃实习步骤：编程计算1月份EU遥相关指数（标准化），画出指数图；编程计算1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形;program calculationreal ind(6l),h500(l44,73,12,6l),r(l44,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 194
4、8-2008do it=l,61do k=l,12do j=1,73do i= 1,144read(2) h500(i,j,k,it)enddoenddoenddoenddowrite(*,*),reading is ok*Icalculationdo it=l,61ind(it)=-0.25*h500(9,59,1 ,it)+0.5*h500(31,59,it)0.25*h500(59,53,l ,it) end dowrite(*,*),calculation is ok! correlationsumy=.Oyy=.Odo it=l,61sumy=sumy+ind(it)yy=yy+md
5、(it)*2enddodo i=l,144do j= 1,73sumx=.Oxx=.0xy=.0do it=l,61sumx=sumx+h500(i,j, 1 ,it)xx=xx+h500(i,j,lxy=xy+h500(i,j,l ,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=l,61w
6、rite(3) ind(it)enddo!write(*,*) r(2,3)write(4) (r(ij),i=l,144)j=l,73)write(*,*)writing is okEnd编程计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。 program EPreal a(160, 60), ind(61), r(160), aa(61), bb (61), lat (160),Ion (160)character8 id(160)real avex, avey, xx, xy, yyopen(2, file=，e:copy3datatl601.txt)open(3, fi
7、le二e:copy3ind. grd，, form=binary，)open(4,file二e:copy3datalat_lon. txt)open (5,file=，e:copy3rl60. grd，, form=，binary，) ccccccccccccccc读数据(指数、经纬度、160站温度) ccc a: 160站气温(19512010年)ind：指数序列(19482008年)read(2, *) (a(i, j), i=l, 160), j=l, 60)close (2)do it=l,60read(3) ind(it)enddodo is=l, 160read (4,*) lat
8、 (is), Ion (is)enddocccccccccccccc编程求相关年份不同sumy=. 0yy三0do it= 1, 58sumy=sumy+ind (it+3)yy=yy+ind (it+3)*2enddo do is=l, 160sumx=. 0xx二.0 xy=. 0 do it=l,58 sumx=sumx+a(is, it) xx=xx+a(is, it)*2 xy=xy+a(is, it)*ind(it+3) enddor (is)二(xy*58-siHnx*sumy) / (sqrt (58*xx-sunix*2) *sqrt (58*yy-sumy*2) end d
9、owrite (*, *)，corralation is ok ccccccccccccccccccc 写站点数据do j=l, 160 id(j)=char (j) tim=0. 0 nlev=lnflag=lwrite (5) id(j), lat (j), lon(j), tim, nlev, nflag, r (j) enddo tim=0. 0 nlev=0 nflag=lwrite(5)id(j-1), lat(j-1), Ion(j-1), tim, nlev, nflag end实习结果：1月份EU遥相关指数，画出指数图；可以看出EU指数随时间的波动，在 可72、1973、19
10、84、1997、2006出现低值,1963、1976、1981、1984、1986、1995、2004 出现高值。1月份EU遥相关指数与500hPa高度场的相关系数分布图，并绘制图形;从图中可以看出EU指数和亚种中部、北美洲北部及大西洋北部有显著正相关，与欧洲、亚 洲东部、太平洋北部及北美洲南部有负相关。这些系统基本为500hPa的永久性或半永久性 环流中心。而与南半球的相关性不显著也反映了这个指数是来自于北半球的资料。该指数反 映了北半球500hPa的环流中心的强弱。EU指数越强，则位于亚洲中部的位势高度越高， 即亚洲低压系统较弱，而北太平洋上的位势越低，则反映副高的减弱。因此此图很好的反映 出了两槽一脊的强度变化，即东亚大槽和欧洲东部的槽加深时，西伯利亚附近的脊发展。计算EU遥相关指数与1月份我国气温的相关系数分布图，并绘制图形。从图中可以看出，EC指数与我国东北南部、山东半岛即东南沿海地区的气温呈显著负相关。 与2图的全球500hPa高度场的在东亚地区的负相关吻合。另外还有一些地区，如内蒙古中 部、陕西、青海、新疆及西藏的少部分地区由显著的负相关。根据热成风原理，EC正距平 的年份，我国中南沿海会有偏北的热成风。北风可能会引起冷空气南下导致该年的冬季气温 偏低。