短期氣候預(yù)測實習(xí)報告五

南京信息工程大學(xué)實驗（實習(xí)）報告實驗（實習(xí)）名稱短期氣候預(yù)測實習(xí)日期11.26得分指導(dǎo)教師系 ~專業(yè)~ 年級~班次~姓名Trichtu學(xué)號~實習(xí)目的：掌握短期氣候預(yù)測因子的分析和選擇，加深對外強迫因子（海溫場）在短期氣候預(yù)測中物理機制的認識。實習(xí)要求：熟悉資料和方法程序（提供部分子程序）；使用GRADS對結(jié)果進行繪圖輸出和文字分析，完成實習(xí)報告。實習(xí)內(nèi)容：計算前期12月北太平洋海溫I和II類雨型合成差值、I和III類合成差值、II和III類合成差值及T檢驗，確定關(guān)鍵區(qū)，北太平洋范圍（120°E~60°W，10°S~60°N）計算前期12月北太平洋海溫I和II類雨型合成差值、I和III類合成差值、II和III類合成差值及T檢驗，確定關(guān)鍵區(qū)。實習(xí)資料：（1）NOAA延伸再建的海表面溫度資料時間范圍：1947年1月－2010年2月范圍（88°S-88°N,0-358°E)網(wǎng)格距2°×2°,緯向格點數(shù)180，經(jīng)向格點數(shù)89資料為GRD格式，資料從南到北、自西向東排列，每月為一個記錄，按年逐月排放。注意：陸地上為缺測值32767（2）1951~2005年雨型分類表100一類雨型010二類雨型001三類雨型實習(xí)方法：(1)不同樣本之間的均值t統(tǒng)計量：分別代表第i類樣本均值和標(biāo)準(zhǔn)差，為第i類樣本的樣本量(2)總體均值的t統(tǒng)計量：一、二、三類雨型的樣本數(shù)分別為20、17和18。我們簡單認為任意兩類之間通過0.01顯著性水平的t值分別為2.73，通過0.05顯著性水平的t值都為2.01(3)編寫程序programssttrealsst(180,89,12,61),ddi(55,3),sst_m(180,89,12*61),sstave(180,89),sst_p(180,89,3)realcha(180,89,3),sig(180,89,3),ave_p(180,89,3),std(180,89,3),sst_p2(180,89,3)integeri,j,k,itcharacter*8id(160)open(2,file='e:\copy5\data\sst.grb',form='binary')open(7,file='e:\copy5\data\ddi')open(9,file='e:\copy5\sst_p.grd',form='binary')open(10,file='e:\copy5\cha.grd',form='binary')open(11,file='e:\copy5\sig.grd',form='binary')!ccccccccccccccc讀數(shù)據(jù)（經(jīng)緯度、160站降水1948-2008、sst（1947,1-2010,2））doit=1,61 dok=1,12 doj=1,89 doi=1,180 read(2)sst(i,j,k,it)if(sst(i,j,k,it)==32767)thensst(i,j,k,it)=0endifsst_m(i,j,12*(it-1)+k)=sst(i,j,k,it) enddoenddoenddoenddo!CCCCCCCCt=1951-2005 doit=1,55 read(7,*)(ddi(it,j),j=1,3) enddowrite(*,*)'readdataok'!ccccccccccccccc編程求合成!ccccccccccccc前期海溫場dok=1,3doi=1,180doj=1,89aa=.0bb=.0num=0doit=1,55num=num+ddi(it,k)aa=aa+sst_m(i,j,12*(it+2)+12)*ddi(it,k)sst_p(i,j,k)=aabb=bb+sst_m(i,j,12*(it+2)+12)sstave(i,j)=bbenddoave_p(i,j,k)=sst_p(i,j,k)/numnum=0doit=1,55num=num+ddi(it,k)aa=(sst_m(i,j,12*(it+2)+12)-ave_p(i,j,k))*ddi(it,k)std(i,j,k)=std(i,j,k)+aa**2enddosstave(i,j)=sstave(i,j)/55.0sst_p(i,j,k)=sst_p(i,j,k)/num-sstave(i,j)enddoenddoenddodoi=1,180doj=1,89cha(i,j,1)=sst_p(i,j,1)-sst_p(i,j,2)cha(i,j,2)=sst_p(i,j,1)-sst_p(i,j,3)cha(i,j,3)=sst_p(i,j,2)-sst_p(i,j,3)enddoenddo!CCCCCCCCCCCCCCCCCCCC雙變量t檢驗doi=1,180doj=1,89aa=sqrt((19*std(i,j,1)/20+16*std(i,j,2)/17)/35.0)sst_p2(i,j,1)=(ave_p(i,j,1)-ave_p(i,j,2))/sqrt(1/20+1/17.0)/aawrite(*,*)sst_p2(i,j,1)aa=sqrt((19*std(i,j,1)/20+17*std(i,j,2)/18)/36.0)sst_p2(i,j,2)=(ave_p(i,j,1)-ave_p(i,j,3))/sqrt((1/20+1/18.0))/aaaa=sqrt((16*std(i,j,2)/17+17*std(i,j,3)/18)/33.0)sst_p2(i,j,3)=(ave_p(i,j,2)-ave_p(i,j,3))/sqrt((1/17+1/18.0))/aaenddoenddodok=1,3doi=1,180doj=1,89selectcase(k)case(1)if(abs(sst_p2(i,j,k))>=2.44)thensig(i,j,k)=2elseif(abs(sst_p2(i,j,k))>=1.69)thensig(i,j,k)=1endifcase(2)if(abs(sst_p2(i,j,k))>=2.44)thensig(i,j,k)=2elseif(abs(sst_p2(i,j,k))>=1.69)thensig(i,j,k)=1endifcase(3)if(abs(sst_p2(i,j,k))>=2.44)thensig(i,j,k)=2elseif(abs(sst_p2(i,j,k))>=1.69)thensig(i,j,k)=1endifendselectenddoenddoenddo!其實都一樣，就可以不用select了!ccccccccccccccccccc寫數(shù)據(jù)dok=1,3write(9)((sst_p(i,j,k),i=1,180),j=1,89)write(10)((cha(i,j,k),i=1,180),j=1,89)write(11)((sig(i,j,k),i=1,180),j=1,89)enddo end（4）結(jié)果輸出1繪制繪制I、II和III類雨型的前期12月北太平洋海溫距平合成圖。夏季降水第一類雨型前期12月北太平洋SSTA合成圖（單位：℃）夏季降水第二類雨型前期12月北太平洋SSTA合成圖（單位：℃）夏季降水第三類雨型前期12月北太平洋SSTA合成圖（單位：℃）2繪制繪制前期12月北太平洋海溫I和II類雨型合成差值、I和III類合成差值、II和III類合成差值圖（包括它們的t檢驗）夏季降水第一、二類雨型前期12月北太平洋SST差值合成圖及其t檢驗（單位：℃，橙、紅陰影表示通過0.05和0.01顯著性水平t檢驗值)前期12月北太平洋SST，一類較二類，從西北太平洋到北太平洋海溫更高，而東太平洋的海溫偏低。夏季降水第一、三類雨型前期12月北太平洋SST差值合成圖及其