您好,欢迎访问三七文档
当前位置:首页 > 商业/管理/HR > 销售管理 > 短期气候预测实习报告五
南京信息工程大学实验(实习)报告实验(实习)名称短期气候预测实习日期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,二类较三类,东太平洋的低纬度的海温偏低,而北太平洋在阿留申群岛附近的海温偏高。
本文标题:短期气候预测实习报告五
链接地址:https://www.777doc.com/doc-5583040 .html