短期气候预测实习报告五

南京信息工程大学实验（实习）报告实验（实习）名称短期气候预测实习日期11.26得分指导教师系~专业~年级~班次~姓名Trichtu学号~实习目的：掌握短期气候预测因子的分析和选择，加深对外强迫因子（海温场）在短期气候预测中物理机制的认识。实习要求：熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。实习内容：计算前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值及T检验，确定关键区，北太平洋范围（120°E~60°W，10°S~60°N）计算前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值及T检验，确定关键区。实习资料：（1）NOAA延伸再建的海表面温度资料时间范围：1947年1月－2010年2月范围（88°S-88°N,0-358°E)网格距2°×2°,纬向格点数180，经向格点数89资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。注意：陆地上为缺测值32767（2）1951~2005年雨型分类表100一类雨型010二类雨型001三类雨型实习方法：(1)不同样本之间的均值t统计量：iXiS分别代表第i类样本均值和标准差，iN为第i类样本的样本量212122221121112)1()1(NNNNSNSNxxt(2)总体均值的t统计量：一、二、三类雨型的样本数分别为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读数据（经纬度、160站降水1948-2008、sst（1947,1-2010,2））doit=1,61dok=1,12doj=1,89doi=1,180read(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-2005doit=1,55read(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写数据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)enddoend（4）结果输出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差值合成图及其t检验（单位：℃，橙、红阴影表示通过0.05和0.01显著性水平t检验值)前期12月北太平洋SST，二类较三类，其差值场有从东太平洋至白令海峡延美洲东海岸的波裂。€夏季降水第二、三类雨型前期12月北太平洋SST差值合成图及其t检验（单位：℃，橙、红阴影表示通过0.05和0.01显著性水平t检验值)前期12月北太平洋SST，二类较三类，东太平洋的低纬度的海温偏低，而北太平洋在阿留申群岛附近的海温偏高。