资源描述:
《短期气候预测基础实习三》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实习三:大气遥相关一、实习目的及要求掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。要求运用资料,计算北半球1月遥相关指数,并分析它与环流和我国气候变化的关系;用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布,正确分析结果数据,完成实习报告。二、实习内容1、计算EU(欧亚)遥相关指数,输出1月份该指数年际变化的时间序列;2、计算EU遥相关指数与同期环流场(500hPa高度场或海平面气压场)的相关系数;3、计算EU遥相关指数与同期我国气温的相关系数。三、实习资料大气环流场资料:1、NCEP/NCAR1
2、948-2008年(61年)的500百帕月平均高度场资料2、资料范围为(900S-900N,00-3600E)3、网格距为2.50×2.50,纬向格点数为144,经向格点数为734、资料为GRD格式,资料从南到北、自西向东排列,每月为一个记录,按年逐月排放。站点资料:1、我国气候中心整编的160站月平均气温资料;2、全国160个台站;3、所给的资料是1月份的;4、资料为txt格式,参加资料说明。四、实习步骤1、编程计算1月份EU遥相关指数(标准化),画出指数图编写fortran程序:programmainparameter(nx=144,ny
3、=73,n=nx*ny,nt=61)realind(nt),z(n,nt),eu(nt),hh(144,73,12,61),avey(n),sy(n),r1(n),r(n)realx,avex,sxopen(10,file='d:yclimatesh3datahgt500.grd',form='binary')open(20,file='d:yclimatesh3ind.grd',form='binary')open(30,file='d:yclimatesh3eu.grd',form='binary')open(4
4、0,file='d:yclimatesh3rxy1.grd',form='binary')doit=1,61dok=1,12doj=1,73doi=1,144read(10)hh(i,j,k,it)enddo;enddo;enddo;enddodoit=1,61doj=1,73doi=1,144z((j-1)*144+i,it)=hh(i,j,1,it)enddo;enddo;enddodoit=1,61eu(it)=-z(9+58*144,it)/4+z(31+58*144,it)/2-z(59+52*144,it)/4enddow
5、rite(20)(eu(it),it=1,61)doit=1,61x=x+eu(it)enddoavex=x/real(nt)doit=1,61sx=sx+((eu(it)-avex)*(eu(it)-avex))/real(nt)enddodoit=1,61ind(it)=(eu(it)-avex)/(sqrt(sx))enddowrite(30)(ind(it),it=1,nt)第二步中的fortran程序:doi=1,navey(i)=0.0doit=1,61avey(i)=avey(i)+z(i,it)enddoavey(i)=ave
6、y(i)/real(nt)enddodoi=1,nr1(i)=0.0doit=1,61r1(i)=r1(i)+(eu(it)-avex)*(z(i,it)-avey(i))enddo;enddodoi=1,nsy(i)=0.0doit=1,61sy(i)=sy(i)+(z(i,it)-avey(i))*(z(i,it)-avey(i))enddo;enddodoi=1,nr(i)=r1(i)/sqrt(sx*real(nt)*sy(i))write(40)r(i)enddoend注:此处给出的是第一步和第二步中的fortran程序,因此,步骤
7、二中此程序不重复给出!数据描述文件:dsetd:yclimatesh3eu.grdundef-9.99E+33titleyearlymeanhgtfromtheNCEPReanalysisxdef1linear11ydef1linear11zdef1levels500tdef61linearjan19481yrvars1eu199yeareuendvars;数据执行(gs)文件:'reinit''opend:yclimatesh3eu.ctl''enableprintd:yclimatesh3eu1.gmf''sett
8、161''deu''print''disableprint';给出1月份EU遥相关指数图:1月份EU遥相关指数图给出1月份标准化的EU遥相关指数1月份标准化的EU遥