预测因子地选择短期气候预测实习报告材料

1、实用标准文案短期气候预测实习报告实习 4 预测因子的选择文档大全实用标准文案南京信息工程大学短期气候预测实习报告实验名称 预测因子的选择 日期 2016 年 11 月 22 日 得分 指导教师系 大气科学 专业 大气科学 年级 班次姓名 学号实习目的： 掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流 异常在短期气候预测中物理机制的认识 。实习要求： 熟悉资料和方法程序（提供部分子程序） ；使用 GRADS 对结果进 行绘图输出和文字分析，完成实习报告。实习资料： NCEP/NCAR 再分析资料 1948 2012 年 1 12 月的 500 百帕月平均高度 场资料范围（ 90 S-

2、90 N,0-360 E）网格距 2.5 2.5 ,纬向格点数 144 ，经向格点数 73资料为 GRD 格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月 排放。国家气候中心整编的 6 、7、8 月降水量资料（时间段： 1951-2010 年，资料的格式参见 readme.txt 文件）1951-2010 年雨型分类表1 0 0 一类雨型0 1 0 二类雨型0 0 1 三类雨型文档大全实用标准文案实习方法：降水距平百分率：Ri RRPi i 100%iRRi 某年夏 季降水量R为1971 2000年夏季降水多年平均值 ?实习步骤：用 FORTRAN 编写相关程序（程序见附件 3 和附

3、件 4 ）利用绘图软件绘制 1951 2010 年夏季三类雨型年合成图；绘制前期冬季 500hPa 高度场距平合成图，指出可能存在的遥相关型环流特征。实习结果及讨论：图 1 1951 2005 年夏季一类雨型年合成图从图 1 可以看出一类雨型主要多雨带位于黄河流域及其以北地区， 江淮流域大范 围少雨，梅雨偏弱并常有较明显的伏旱，江南南部至华南一般为次要多雨区。文档大全实用标准文案图 2 1951 2005 年夏季二类雨型年合成图从图 2 可以看出二类雨型主要多雨带位于黄河至长江之间， 雨带中心一般在淮河流域一带，黄河以北及长江以南大部地区少雨。图 3 1951 2005 年夏季三类雨型年合成图

4、从图 3 可以看出三类雨型主要多雨带位于长江流域或江南一带， 淮河以北大部及 东南沿海地区少雨。文档大全实用标准文案图 4 第一类雨型前期冬季 500hPa 高度场的距平图图 4 显示为中太平洋和北太平洋有显著的反相关性， 或者说是副高的高压和阿留 申低压之间关系， 副高增强的同时，低压也同时增强。 出现这种前期冬季 500hPa 高度场和今年的一类雨型北方性有显著关系。这种相关可能是 WP 型要相关。图 5 第二类雨型前期冬季 500hPa 高度场的距平图图 5 显示为北太平洋和北美洲出现的以北美州西海岸为界限的类轴对称图形。 当 北太平洋的高压脊增强时， 北美大槽也同样增强。 而与两个系统

5、对应的低纬度地 区呈现与之相反的变化关系。我认为这是一种 PNA 和 WP 的杂交型。文档大全实用标准文案图 6 第三类雨型和前期冬季 500hPa 高度场的距平图图 6 显示为副热带太平洋高压加强， 位于阿留申地区的气压加深， 当北太平洋的 高压脊增强时， 北美大槽也同样增强， 这造成了美国的严寒天气， 也造成了我国 的降水主要集中在江淮地区的类型。这是典型的 PNA 型，但低纬地区没有很显 著。附件：附件 1program EPinteger ind(160)real a(160,65),r(144,73),aa(65),bb(65),lat(160),lon(160)real hgt(1

6、44,73,12,65),h1(144,73,65),eu(65)文档大全实用标准文案character*8 id(160)open(1,file=i:duanqi3hgt500.grd,form=binary)open(3,file=i:duanqi3ind.grd,form=binary)open(4,file=i:duanqi3lat_lon.txt)open(5,file=i:duanqi3r.grd,form=binary)do kk=1,65do n=1,12read(1)(hgt(i,j,n,kk),i=1,144),j=1,73)enddoenddodo kk=1,65do j

7、=1,73do i=1,144h1(i,j,kk)=hgt(i,j,1,kk)enddoenddoeu(kk)=-h1(9,59,kk)/4.0+h1(31,59,kk)/2.0-h1(59,53,kk)/4.0enddowrite(3) (eu(i),i=1,65)close(3)do j=1,73do i=1,144文档大全实用标准文案r(i,j)=0.0 do kk=1,65 aa(kk)=eu(kk) bb(kk)=h1(i,j,kk) enddo call Normalized_1D(aa,65) call Normalized_1D(bb,65) do kk=1,65 r(i,j)

8、=r(i,j)+aa(kk)*bb(kk)/65.0 enddo enddo enddo write(*,*)r do j=1,73 do i=1,144 write(5) r(i,j) enddo enddo endSUBROUTINE Normalized_1D(BZH,NUM) REAL BZH(NUM),AVE,VARAVE=0.0文档大全实用标准文案VAR=0.0DO K=1,NUMAVE=AVE+BZH(K)/NUMENDDODO K=1,NUMVAR=VAR+(BZH(K)-AVE)*2ENDDOVAR=SQRT(VAR/NUM)DO K=1,NUMBZH(K)=(BZH(K)-

9、AVE)/VARENDDOEND附件 2program EPinteger ind(160)real a(160,62),r(160),aa(62),bb(62),lat(160),lon(160)real hgt(144,73,12,65),h1(144,73,62),eu(62)character*8 id(160)open(1,file=i:duanqi3hgt500.grd,form=binary)open(2,file=i:duanqi3t1601.txt)open(4,file=i:duanqi3lat_lon.txt)文档大全实用标准文案open(5,file=i:duanqi3

10、r160.grd,form=binary)do kk=1,65do n=1,12read(1)(hgt(i,j,n,kk),i=1,144),j=1,73)enddoenddoclose(1)do kk=1,62do j=1,73do i=1,144h1(i,j,kk)=hgt(i,j,1,kk+3)enddoenddoeu(kk)=-h1(9,59,kk)/4.0+h1(31,59,kk)/2.0-h1(59,53,kk)/4.0 enddoread(2,*)(a(i,j),i=1,160),j=1,62)close(2)do i=1,160read(4,*)id(i),lat(i),lon

11、(i)enddoclose(4)do i=1,160文档大全实用标准文案r(i)=0.0 do kk=1,62 aa(kk)=eu(kk) bb(kk)=a(i,kk) enddocall Normalized_1D(aa,62) call Normalized_1D(bb,62) do kk=1,62 r(i)=r(i)+aa(kk)*bb(kk)/62.0 enddo enddo write(*,*)r do j=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

12、) enddo tim=0.0 nlev=0nflag=1文档大全实用标准文案write(5)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagclose(5)endSUBROUTINE Normalized_1D(BZH,NUM)REAL BZH(NUM),AVE,VARAVE=0.0VAR=0.0DO K=1,NUMAVE=AVE+BZH(K)/NUMENDDODO K=1,NUMVAR=VAR+(BZH(K)-AVE)*2ENDDOVAR=SQRT(VAR/NUM)DO K=1,NUMBZH(K)=(BZH(K)-AVE)/VARENDDOEND附件 3pr

13、ogram EPreal文档大全实用标准文案a6(160,60),a7(160,60),a8(160,60),lat(160),lon(160),rap1(160),rap2(160), rap3(160),ddi(60,3 ),rsum(160),r(160,60),rp(160,60)character*8 id(160)open(3,file=i:duanqi4r1606.txt)open(4,file=i:duanqi4r1607.txt)open(5,file=i:duanqi4r1608.txt) open(6,file=i:duanqi4lat_lon.txt) open(7,f

14、ile=i:duanqi4ddi) open(8,file=i:duanqi4rap1.grd,form=binary) open(9,file=i:duanqi4rap2.grd,form=binary) open(10,file=i:duanqi4rap3.grd,form=binary) read(3,*)(a6(i,j),i=1,160),j=1,60)read(4,*)(a7(i,j),i=1,160),j=1,60)read(5,*)(a8(i,j),i=1,160),j=1,60)do i=1,160read(6,*)lat(i),lon(i)enddodo it=1,60rea

15、d(7,*)(ddi(it,j),j=1,3)enddodo i=1,160rsum(i)=0.0文档大全实用标准文案do it=1,60 rsum(i)=rsum(i)+(a6(i,it)+a7(i,it)+a8(i,it)/3 enddorsum(i)=rsum(i)/60.0do j=1,60rp(i,j)=(a6(i,j)+a7(i,j)+a8(i,j)/3-rsum(i)enddodo j=1,60rp(i,j)=rp(i,j)/rsum(i)enddoenddodo i=1,160rap1(i)=0rap2(i)=0rap3(i)=0do j=1,60rap1(i)=rap1(i)

16、+rp(i,j)*ddi(j,1)/22rap2(i)=rap2(i)+rp(i,j)*ddi(j,2)/19rap3(i)=rap3(i)+rp(i,j)*ddi(j,3)/19enddorap1(i)=rap1(i)*100rap2(i)=rap2(i)*100文档大全实用标准文案rap3(i)=rap3(i)*100enddodo j=1,160id(j)=char(j)tim=0.0nlev=1nflag=1write(8)id(j),lat(j),lon(j),tim,nlev,nflag,rap1(j) enddotim=0.0nlev=0nflag=1write(8)id(j-1

17、),lat(j-1),lon(j-1),tim,nlev,nflag do j=1,160id(j)=char(j)tim=0.0nlev=1nflag=1write(9)id(j),lat(j),lon(j),tim,nlev,nflag,rap2(j) enddotim=0.0nlev=0文档大全实用标准文案nflag=1write(9)id(j-1),lat( j-1),lon(j-1),tim,nlev,nflagdo j=1,160id(j)=char(j)tim=0.0nlev=1nflag=1write(10)id(j),lat(j),lon(j),tim,nlev,nflag,

18、rap3(j)enddotim=0.0nlev=0nflag=1write(10)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagend附件 4program heightreal h(144,73,12,65),ha(144,73,60),ddi(60,3),ave(144,73),ave1(144,73),ave 2(144,73),ave3(144,73),i,j,mo,itopen(1,file=i:duanqi4gdhgt500.grd,form=binary) open(2,file=i:duanqi4gdddi)文档大全实用标准文案open(3,file=i:duanqi4gd1h1.grd,form=binary)open(4,file=i:duanqi4gd2h2.grd,form=binary)open(5,file=i:duanqi4gd3h3.grd,form=binary)read(1) (h(i,j,mo,it),i=1,144),j=1,73),mo=1,12),it=1,65)do it=1