預(yù)測因子的選擇(短期氣候預(yù)測實(shí)習(xí)報(bào)告)

1、短期氣候預(yù)測實(shí)習(xí)報(bào)告實(shí)習(xí)4預(yù)測因子的選擇南京信息工程大學(xué)短期氣候預(yù)測實(shí)習(xí)報(bào)告實(shí)驗(yàn)名稱預(yù)測因子的選擇日期2016年11月22日得分指導(dǎo)教師系大氣科學(xué)專業(yè)大氣科學(xué)年級(jí)班次學(xué)號(hào)實(shí)習(xí)目的：掌握短期氣候預(yù)測因子的分析和選擇，加深對(duì)夏季降水分布、環(huán)流異常在短期氣候預(yù)測中物理機(jī)制的認(rèn)識(shí)。實(shí)習(xí)要求：熟悉資料和方法程序提供部分子程序；使用GRADS寸結(jié)果進(jìn)行繪圖輸出和文字分析，完成實(shí)習(xí)報(bào)告。實(shí)習(xí)資料：NCEP/NCAR1分析資料19482012年112月的500百帕月平均高度場資料范圍90°S-900N,0-360°E)網(wǎng)格距2.5°X2.5°,緯向格點(diǎn)數(shù)144,經(jīng)向格點(diǎn)

2、數(shù)73資料為GRD&式，資料從南到北、自西向東排列，每月為一個(gè)記錄，按年逐月排放。國家氣候中心整編的6、7、8月降水量資料1951-2010年雨型分類表100一類雨型010二類雨型001三類雨型實(shí)習(xí)方法：降水距平百分率：RRRPi100%RR某年夏季降水量R為19712000年夏季降水多年平均值？實(shí)習(xí)步驟：用FORTRAN寫相關(guān)程序程序見附件3和附件4利用繪圖軟件繪制19512010年夏季三類雨型年合成圖；繪制前期冬季500hPa高度場距平合成圖，指出可能存在的遙相關(guān)型環(huán)流特征。實(shí)習(xí)結(jié)果及討論:圖119512005年夏季一類雨型年合成圖從圖1可以看出一類雨型主要多雨帶位于黃河流域及其以

3、北地區(qū)，江淮流域大范圍少雨，梅雨偏弱并常有較明顯的伏旱，江南南部至華南一般為次要多雨區(qū)。20141.543022*hNMH.'hiMJUM圖219512005年夏季二類雨型年合成圖從圖2可以看出二類雨型主要多雨帶位于黃河至長江之間，雨帶中心一般在淮河流域一帶，黃河以北及長江以南大部地區(qū)少雨。圖319512005年夏季三類雨型年合成圖從圖3可以看出三類雨型主要多雨帶位于長江流域或江南一帶，淮河以北大部及東南沿海地區(qū)少雨。20141343022圖4第一類雨型前期冬季500hPa高度場的距平圖圖4顯示為中太平洋和北太平洋有顯著的反相關(guān)性，或者說是副高的高壓和阿留中低壓之間關(guān)系，副高增強(qiáng)的同時(shí)

4、，低壓也同時(shí)增強(qiáng)。出現(xiàn)這種前期冬季500hPa高度場和今年的一類雨型北方性有顯著關(guān)系。這種相關(guān)可能是WPffl要相關(guān)。2014134.3022圖5第二類雨型前期冬季500hPa高度場的距平圖圖5顯示為北太平洋和北美洲出現(xiàn)的以北美州西海岸為界限的類軸對(duì)稱圖形。當(dāng)北太平洋的高壓脊增強(qiáng)時(shí)，北美大槽也同樣增強(qiáng)。而與兩個(gè)系統(tǒng)對(duì)應(yīng)的低緯度地區(qū)呈現(xiàn)與之相反的變化關(guān)系。我認(rèn)為這是一種PN用口WP勺雜交型。20141343022圖6第三類雨型和前期冬季500hPa高度場的距平圖圖6顯示為副熱帶太平洋高壓加強(qiáng)，位于阿留申地區(qū)的氣壓加深，當(dāng)北太平洋的高壓脊增強(qiáng)時(shí)，北美大槽也同樣增強(qiáng)，這造成了美國的嚴(yán)寒天氣，也造成了

5、我國的降水主要集中在江淮地區(qū)的類型。這是典型的PNA®,但低緯地區(qū)沒有很顯著。附件：附件1programEPintegerind(160)reala(160,65),r(144,73),aa(65),bb(65),lat(160),lon(160)realhgt(144,73,12,65),h1(144,73,65),eu(65)character*8id(160)open(1,file='i:duanqi3hgt500.grd',form='binary')open(3,file='i:duanqi3ind.grd',form=

6、9;binary')open(4,file='i:duanqi3lat_lon.txt')open(5,file='i:duanqi3r.grd',form='binary')dokk=1,65don=1,12read(1)(hgt(i,j,n,kk),i=1,144),j=1,73)enddoenddodokk=1,65doj=1,73doi=1,144h1(i,j,kk)=hgt(i,j,1,kk)enddoenddoeu(kkenddowrite(3)(eu(i),i=1,65)close(3)doj=1,73doi=1,144r(

7、i,jdokk=1,65aa(kk)=eu(kk)bb(kk)=h1(i,j,kk)enddocallNormalized_1D(aa,65)callNormalized_1D(bb,65)dokk=1,65r(i,j)=r(i,j)+aa(kk)*bb(kkenddoenddoenddowrite(*,*)rdoj=1,73doi=1,144write(5)r(i,j)enddoenddoendSUBROUTINENormalized_1D(BZH,NUM)REALBZH(NUM),AVE,VARDOK=1,NUMAVE=AVE+BZH(K)/NUMENDDODOK=1,NUMVAR=VAR

8、+(BZH(K)-AVE)*2ENDDOVAR=SQRT(VAR/NUM)DOK=1,NUMBZH(K)=(BZH(K)-AVE)/VARENDDOEND附件2programEPintegerind(160)reala(160,62),r(160),aa(62),bb(62),lat(160),lon(160)realhgt(144,73,12,65),h1(144,73,62),eu(62)character*8id(160)open(1,file='i:duanqi3hgt500.grd',form='binary')open(2,file='i:d

9、uanqi3t1601.txt')open(4,file='i:duanqi3lat_lon.txt')open(5,file='i:duanqi3r160.grd',form='binary')dokk=1,65don=1,12read(1)(hgt(i,j,n,kk),i=1,144),j=1,73)enddoenddoclose(1)dokk=1,62doj=1,73doi=1,144h1(i,j,kk)=hgt(i,j,1,kk+3)enddoenddoeu(kkenddoread(2,*)(a(i,j),i=1,160),j=1

10、,62)close(2)doi=1,160read(4,*)id(i),lat(i),lon(i)enddoclose(4)doi=1,160r(idokk=1,62aa(kk)=eu(kk)bb(kk)=a(i,kk)enddocallNormalized_1D(aa,62)callNormalized_1D(bb,62)dokk=1,62r(i)=r(i)+aa(kk)*bb(kkenddoenddowrite(*,*)rdoj=1,160id(j)=char(j)tim=nlev=1nflag=1write(5)id(j),lat(j),lon(j),tim,nlev,nflag,r(j

11、)enddotim=nlev=0nflag=1write(5)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagclose(5)endSUBROUTINENormalized_1D(BZH,NUM)REALBZH(NUM),AVE,VARDOK=1,NUMAVE=AVE+BZH(K)/NUMENDDODOK=1,NUMVAR=VAR+(BZH(K)-AVE)*2ENDDOVAR=SQRT(VAR/NUM)DOK=1,NUMBZH(K)=(BZH(K)-AVE)/VARENDDOEND附件3programEPreala6(160,60),a7(160,60),a8(

12、160,60),lat(160),lon(160),rap1(160),rap2(160),rap3(160),ddi(60,3),rsum(160),r(160,60),rp(160,60)character*8id(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,file='i:d

13、uanqi4ddi')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

14、,60)doi=1,160read(6,*)lat(i),lon(i)enddodoit=1,60read(7,*)(ddi(it,j),j=1,3)enddodoi=1,160rsum(idoit=1,60rsum(i)=rsum(i)+(a6(i,it)+a7(i,it)+a8(i,it)/3enddorsum(i)=rsum(idoj=1,60rp(i,j)=(a6(i,j)+a7(i,j)+a8(i,j)/3-rsum(i)enddodoj=1,60rp(i,j)=rp(i,j)/rsum(i)enddoenddodoi=1,160rap1(i)=0rap2(i)=0rap3(i)=0

15、doj=1,60rap1(i)=rap1(i)+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)*100rap3(i)=rap3(i)*100enddodoj=1,160id(j)=char(j)tim=nlev=1nflag=1write(8)id(j),lat(j),lon(j),tim,nlev,nflag,rap1(j)enddotim=nlev=0nflag=1write(8)i

16、d(j-1),lat(j-1),lon(j-1),tim,nlev,nflagdoj=1,160id(j)=char(j)tim=nlev=1nflag=1write(9)id(j),lat(j),lon(j),tim,nlev,nflag,rap2(j)enddotim=nlev=0nflag=1write(9)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagdoj=1,160id(j)=char(j)tim=nlev=1nflag=1write(10)id(j),lat(j),lon(j),tim,nlev,nflag,rap3(j)enddotim=nle

17、v=0nflag=1write(10)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagend附件4programheightrealh(144,73,12,65),ha(144,73,60),ddi(60,3),ave(144,73),ave1(144,73),ave2(144,73),ave3(144,73),i,j,mo,itopen(1,file='i:duanqi4gdhgt500.grd',form="binary")open(2,file='i:duanqi4gdddi')open(3,file=&

18、#39;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)doit=1,60read(2,*)(ddi(it,j),j=1,3)enddodoi=1,144doj=1,73doit=1,60ha(i,j,it)=(h(i,j,12,it+2)+h(i,j,1,it+3)+h(i,j,2,it+3)/3enddoenddoenddodoi=1,144doj=1,73doit=1,60ave(i,j)=ave(i,j)+ha(i,j,it)/60enddoenddoenddodoi=1,144doj=1,73doit=1,60if(ddi(it,1)=1)thenave1(i,j)=ave1(i,j)+ha(i,j,it)/22elseif(dd