短期气候预测实习报告四

....南京信息工程大学实验（实习）报告实验（实习）名称短期气候预测实习日期11.15得分指导教师系~专业~年级~班次~姓名Trichtu学号~实习目的：掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流异常在短期气候预测中物理机制的认识。实习要求;熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。实习内容：计算1951-2005年夏季三类年雨型合成图；计算各类雨型的前期冬季高度场距平合成图，指出可能出现的遥相关型。资料：NCEP/NCAR再分析资料1948－2008年1～12月的500百帕月平均高度场资料范围（90°S-90°N,0-360°E)网格距2.5°×2.5°,纬向格点数144，经向格点数73资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。、实习资料：国家气候中心整编的6、7、8月降水量资料（时间段：1951~2010年，资料的格式参见readme.txt文件）1951~2005年雨型分类表100一类雨型010二类雨型001三类雨型实习方法及步骤：（1）降水距平百分率：iR为某年夏季降水量，R为1971-2000年夏季降水多年平均值（2）总体均值的t统计量：X，s分别代表样本均值和标准差，0为总体均值，n为样本量根据公式及资料，一、二、三类雨型的样本数分别为20、17和18。我们认为简单认为它们通过0.01显著性水平的t值分别为2.86，2.92，2.90，通过0.05显著性水平的t值都为2.09，2.12，2.11。%100RRRRPiinsxt0....（3）编写程序programEPreala6(160,60),a7(160,60),a8(160,60),lat(160),h_y(144,73,12,61)reallon(160),rap(160,3),ddi(55,3),rsum(160),r(160,55)realh_m(144,73,12*61),h_p(144,73,3),aa,have(144,73)realave_p(144,73,3),std(144,73,3),h_p2(144,73,3),sig(144,73,3)integeri,j,k,itcharacter*8id(160)open(2,file='e:\copy4\data\hgt500.grd',form='binary')open(3,file='e:\copy4\data\r1606.txt')open(4,file='e:\copy4\data\r1607.txt')open(5,file='e:\copy4\data\r1608.txt')open(6,file='e:\copy4\data\lat_lon.txt')open(7,file='e:\copy4\data\ddi')open(8,file='e:\copy4\rap.grd',form='binary')open(9,file='e:\copy4\hp.grd',form='binary')open(10,file='e:\copy4\sig.grd',form='binary')!ccccccccccccccc读数据(经纬度、160站降水1948-2008、雨型、height(1948-2008))!CCCCCCCCCCCCCCCCCCCC一类20，二类17，三类18CCCCCCCCCCCCCCCCCCcdoit=1,61dok=1,12doj=1,73doi=1,144read(2)h_y(i,j,k,it)h_m(i,j,12*(it-1)+k)=h_y(i,j,k,it)enddoenddoenddoenddoread(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,60)doi=1,160read(6,*)lat(i),lon(i)enddo!CCCCCCCCt=1951-2005doit=1,55read(7,*)(ddi(it,j),j=1,3)enddowrite(*,*)'readdataok'!ccccccccccccccc编程求合成doi=1,160....rsum(i)=0.0doit=21,50rsum(i)=rsum(i)+a6(i,it)+a7(i,it)+a8(i,it)enddorsum(i)=rsum(i)/30.0doit=1,55r(i,it)=(a6(i,it)+a7(i,it)+a8(i,it)-rsum(i))*100/rsum(i)enddoenddodoi=1,160doj=1,3num=0doit=1,55num=num+ddi(it,j)rap(i,j)=rap(i,j)+r(i,it)*ddi(it,j)enddorap(i,j)=rap(i,j)/real(num)enddo;enddo!ccccccccccccc前期高度场dok=1,3doi=1,144doj=1,73aa=.0bb=.0num=0doit=1,55num=num+ddi(it,k)aa=aa+(h_m(i,j,12*(it+1)+12)+h_m(i,j,12*(it+2)+1)+h_m(i,j,12*(it+2)+2))*ddi(it,k)h_p(i,j,k)=aa/3.0bb=bb+h_m(i,j,12*(it+1)+12)+h_m(i,j,12*(it+2)+1)+h_m(i,j,12*(it+2)+2)have(i,j)=bb/3.0enddoave_p(i,j,k)=h_p(i,j,k)/numnum=0doit=1,55num=num+ddi(it,k)aa=(h_m(i,j,12*(it+1)+12)+h_m(i,j,12*(it+2)+1)+h_m(i,j,12*(it+2)+2))*ddi(it,k)/3.0-ave_p(i,j,k)*ddi(it,k)std(i,j,k)=std(i,j,k)+aa**2enddohave(i,j)=have(i,j)/55.0....h_p(i,j,k)=h_p(i,j,k)/num-have(i,j)h_p2(i,j,k)=h_p(i,j,k)/sqrt(std(i,j,k)/num)*sqrt(num/1.0)!write(*,*)h_p2(i,j,k)enddoenddoenddo!CCCCCCCCCCC样本检验dok=1,3doi=1,144doj=1,73selectcase(k)case(1)if(abs(h_p2(i,j,k))=2.52)thensig(i,j,k)=2elseif(abs(h_p2(i,j,k))=1.72)thensig(i,j,k)=1endifcase(2)if(abs(h_p2(i,j,k))=2.57)thensig(i,j,k)=2elseif(abs(h_p2(i,j,k))=1.74)thensig(i,j,k)=1endifcase(3)if(abs(h_p2(i,j,k))=2.55)thensig(i,j,k)=2elseif(abs(h_p2(i,j,k))=1.73)thensig(i,j,k)=1endifendselect!write(*,*)sig(i,j,k)enddoenddoenddo!ccccccccccccccccccc写站点数据doj=1,160id(j)=char(j)tim=0.0nlev=1nflag=1write(8)id(j),lat(j),lon(j),tim,nlev,nflag,(rap(j,i),i=1,3)enddotim=0.0....nlev=0nflag=1write(8)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflagdok=1,3write(9)((h_p(i,j,k),i=1,144),j=1,73)write(10)((sig(i,j,k),i=1,144),j=1,73)enddoEnd对于CTL文件（列举一个）：dsete:\copy4\rap.grddtypestationstnmape:\copy4\rap.mapundef99999.9titlethe160stationwinterranomalytdef1linearjan19841yrvars3r10991rainanomalyr20992rainanomalyr30993rainanomalyendvarsGS文件（列举一个）：'reinit''opene:\copy4\hp.ctl''opene:\copy4\sig.ctl''enableprinte:\copy4\sig.gmf''setgxoutshaded''setclevs0.51.5''setccols02427''setcmin0''dt1.2''setgxoutcontour''dh1''drawtitleFirstkindofraincorralationwith0.05/0.01significance''print''c'....（4）结果输出②绘制1951－2005年夏季三类雨型年合成图：一类雨型：主要多雨带位于黄河流域及其以北地区，江淮流域大范围少雨，梅雨偏弱并常有较明显的伏旱，江南南部至华南一般为次要多雨区。二类雨型：主要多雨带位于黄河至长江之间，雨带中心一般在淮河流域一带，黄河以北及长江以南大部地区少雨。....三类雨型：主要多雨带位于长江流域或江南一带，淮河以北大部及东南沿海地区少雨。③绘制前期冬季500hPa高度场距平合成图，指出可能存在的遥相关型环流特征。第一类雨型和前期冬季500hPa高度场的距平图及相关检验（橙色通过0.05，红色通过0.01）....显示为中太平洋和北太平洋有显著的反相关性，或者说是副高的高压和阿留申低压之间关系，副高增强的同时，低压也同时增强。出现这种前期冬季500hPa高度场和今年的一类雨型北方性有显著关系。这种相关可能是WP型要相关。第二类雨型和前期冬季500hPa高度场的距平图及相关检验（橙色通过0.05，红色通过0.01）显示为北太平洋和北美洲出现的以北美州西海岸为界限的类轴对称图形。当北太平洋的高压脊增强时，北美大槽也同样增强。而与两个系统对应的低纬度地区呈现与之相反的变化关系。我认为这是一种PNA和WP的杂交型。第三类雨型和前期冬季500hPa高度场的距平图及相关检验（橙色通过0.05，红色通过0.01）显示为副热带太平洋高压加强，位于阿留申地区的气压加深，当北太平洋的高压脊增强时，北美大槽也同样增强。这造成了美国的严寒天气，竟然也造成了我国的降水主要集中在江淮....地区的类型，好神奇。这是典型的PNA型，但低纬地区没有很显著，且显著区域的显著性有点低，范围较小，预报性不会很好（Isupposed）。欢迎您的光临，Word文档下载后可修改编辑.双击可删除页眉页脚.谢谢！你的意见是我进步的动力，希望您提出您宝贵的意见！让我们共同学习共同进步！学无止境.更上一层楼。