短期氣候預(yù)測(cè)實(shí)習(xí)報(bào)告三

南京信息工程大學(xué)實(shí)驗(yàn)（實(shí)習(xí)）報(bào)告實(shí)驗(yàn)（實(shí)習(xí)）名稱短期氣候預(yù)測(cè)實(shí)習(xí)日期11.12得分指導(dǎo)教師系 ~專業(yè)~年級(jí)~班次~姓名Trichtu學(xué)號(hào)~實(shí)習(xí)目的：掌握大氣環(huán)流中遙相關(guān)型指數(shù)的計(jì)算及其與大氣環(huán)流和我國氣候關(guān)系的分析。實(shí)習(xí)要求：要求運(yùn)用資料，計(jì)算北半球1月遙相關(guān)指數(shù)，并分析它與環(huán)流和我國氣候變化的關(guān)系；用圖形輸出指數(shù)年際變化曲線、遙相關(guān)的空間分布以及與我國氣溫的相關(guān)系數(shù)分布，正確分析結(jié)果數(shù)據(jù)，完成實(shí)習(xí)報(bào)告。實(shí)習(xí)內(nèi)容：計(jì)算EU（歐亞）遙相關(guān)指數(shù)，輸出1月份該指數(shù)年際變化的時(shí)間序列；計(jì)算EU遙相關(guān)指數(shù)與同期環(huán)流場(chǎng)（500hPa高度場(chǎng)或海平面氣壓場(chǎng)）的相關(guān)系數(shù)；計(jì)算EU遙相關(guān)指數(shù)與同期我國氣溫的相關(guān)系數(shù)。實(shí)習(xí)資料：（大氣環(huán)流場(chǎng)資料）NCEP/NCAR1948－2008年（61年）的500百帕月平均高度場(chǎng)資料資料范圍為（900S-900N,00-3600E)網(wǎng)格距為2.50×2.50,緯向格點(diǎn)數(shù)為144，經(jīng)向格點(diǎn)數(shù)為73資料為GRD格式，資料從南到北、自西向東排列，每月為一個(gè)記錄，按年逐月排放。(站點(diǎn)資料)：我國氣候中心整編的160站月平均氣溫資料；全國160個(gè)臺(tái)站，所給的資料是1月份的，資料為txt格式，參加資料說明。實(shí)習(xí)方法簡介：EU遙相關(guān)指數(shù)的定義:相關(guān)系數(shù)的計(jì)算：實(shí)習(xí)步驟：編程計(jì)算1月份EU遙相關(guān)指數(shù)（標(biāo)準(zhǔn)化），畫出指數(shù)圖；編程計(jì)算1月份EU遙相關(guān)指數(shù)與500hPa高度場(chǎng)的相關(guān)系數(shù)分布圖，并繪制圖形；programcalculationrealind(61),h500(144,73,12,61),r(144,73)realsumx,sumy,xx,xy,yyopen(2,file='e:\copy3\data\hgt500.grd',form='binary')open(3,file='e:\copy3\ind.grd',form='binary')open(4,file='e:\copy3\corr.grd',form='binary')!readdata1948-2008doit=1,61dok=1,12doj=1,73doi=1,144read(2)h500(i,j,k,it)enddoenddoenddoenddowrite(*,*)'readingisok'!calculationdoit=1,61ind(it)=-0.25*h500(9,59,1,it)+0.5*h500(31,59,1,it)-0.25*h500(59,53,1,it)enddowrite(*,*)'calculationisok'!correlationsumy=.0yy=.0doit=1,61sumy=sumy+ind(it)yy=yy+ind(it)**2enddodoi=1,144doj=1,73sumx=.0xx=.0xy=.0doit=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))enddoenddowrite(*,*)'corralationisok'!writingdoit=1,61write(3)ind(it)enddo!write(*,*)r(2,3)write(4)((r(i,j),i=1,144),j=1,73)write(*,*)'writingisok'End編程計(jì)算EU遙相關(guān)指數(shù)與1月份我國氣溫的相關(guān)系數(shù)分布圖，并繪制圖形。programEP reala(160,60),ind(61),r(160),aa(61),bb(61),lat(160),lon(160) character*8id(160)realavex,avey,xx,xy,yy open(2,file='e:\copy3\data\t1601.txt') open(3,file='e:\copy3\ind.grd',form='binary') open(4,file='e:\copy3\data\lat_lon.txt')open(5,file='e:\copy3\r160.grd',form='binary')ccccccccccccccc讀數(shù)據(jù)（指數(shù)、經(jīng)緯度、160站溫度）ccca:160站氣溫(1951~2010年)ind：指數(shù)序列(1948~2008年)read(2,*)((a(i,j),i=1,160),j=1,60) close(2) doit=1,60 read(3)ind(it) enddodois=1,160 read(4,*)lat(is),lon(is) enddocccccccccccccc編程求相關(guān)年份不同sumy=.0yy=.0doit=1,58sumy=sumy+ind(it+3)yy=yy+ind(it+3)**2enddodois=1,160sumx=.0xx=.0xy=.0doit=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))enddowrite(*,*)'corralationisok'ccccccccccccccccccc寫站點(diǎn)數(shù)據(jù) doj=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) enddo tim=0.0 nlev=0 nflag=1 write(5)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflag end實(shí)習(xí)結(jié)果:1月份EU遙相關(guān)指數(shù)，畫出指數(shù)圖；可以看出EU指數(shù)隨時(shí)間的波動(dòng)，，在1972、1973、1984、1997、2006出現(xiàn)低值，1963、1976、1981、1984、1986、1995、2004出現(xiàn)高值。1月份EU遙相關(guān)指數(shù)與500hPa高度場(chǎng)的相關(guān)系數(shù)分布圖，并繪制圖形；從圖中可以看出EU指數(shù)和亞種中部、北美洲北部及大西洋北部有顯著正相關(guān)，與歐洲、亞洲東部、太平洋北部及北美洲南部有負(fù)相關(guān)。這些系統(tǒng)基本為500hPa的永久性或半永久性環(huán)流中心。而與南半球的相關(guān)性不顯著也反映了這個(gè)指數(shù)是來自于北半球的資料。該指數(shù)反映了北半球500hPa的環(huán)流中心的強(qiáng)弱。EU指數(shù)越強(qiáng)，則位于亞洲中部的位勢(shì)高度越高，即亞洲低壓系統(tǒng)較弱，而北太平洋上的位勢(shì)越低，則反映副高的減弱。因此此圖很好的反映出了兩槽一脊的強(qiáng)度變化，即東亞大槽和歐洲東部的槽加深時(shí)，西伯利亞附近的脊發(fā)展。計(jì)算EU遙相關(guān)指數(shù)與1月份我國氣溫的相關(guān)系數(shù)分布圖，并繪制圖形。