资源描述：
《短期气候预测实习报告4预测因子选择》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、短期气候预测实习报告实习4预测因子选择—合成分析法南京信息工程大学实验（实习）报告实验(实习)名称短期气候预测基础课程设计日期2016.10指导教师李忠贤系大气科学专业大气科学班次长望班姓名夏新露一、实习目的掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流异常在短期气候预测中物理机制的认识。二、实习要求熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。三、实习内容计算1951-2005年夏季三类年雨型合成图；计算各类雨型的前期冬季高度场距平合成图，指出可能出现的遥相关型。资料：NCEP/NCAR
2、再分析资料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）降水距平百分率：某年夏季降水量为1971-2000年夏季降水多年平均值（2）总体均值的t统计量：
3、、分别代表样本均值和标准差为总体均值为样本量一、二、三类雨型的样本数分别为22、19和19。我们认为简单认为它们通过0.01显著性水平的t值分别为2.85，2.90，2.90，通过0.05显著性水平的t值都为2.09，2.11，2.11。一、实习步骤程序编写：（1）三类雨型合成图程序：realnsparameter(ns=55)integerddi(55,3)integere,f,greala,b,crealhh(144,73,12,65)realhj1(144,73,20),hj2(144,73,17),hj3(144,73,18),hjave(1
4、44,73)realhjjp1(144,73,20),hjjp2(144,73,17),hjjp3(144,73,18)realjppj1(144,73),jppj2(144,73),jppj3(144,73)open(15,file='ddi')open(16,file='hgt48-12500.grd',form='binary',STATUS='OLD')doit=1,65dok=1,12doj=1,73doi=1,144read(16)hh(i,j,k,it)enddo;enddo;enddo;enddoread(15,*)((ddi(i,j
5、),j=1,3),i=1,55)print*,((ddi(i,j),j=1,3),i=1,55)a=0b=0c=0doi=1,55if(ddi(i,1)==1)thena=a+1elseif(ddi(i,2)==1)thenb=b+1elsec=c+1endifenddoprint*,a,b,ce=0f=0g=0doi=1,55if(ddi(i,1)==1)thene=e+1doix=1,144doiy=1,73hj1(ix,iy,e)=hh(ix,iy,1,i)enddoenddoelseif(ddi(i,2)==1)thenf=f+1doix=1
6、,144doiy=1,73hj2(ix,iy,f)=hh(ix,iy,1,i)enddoenddoelseg=g+1doix=1,144doiy=1,73hj3(ix,iy,g)=hh(ix,iy,1,i)enddoenddoendifenddodoi=1,144doj=1,73hjave(i,j)=0doit=1,55hjave(i,j)=hjave(i,j)+hh(i,j,1,it)enddohjave(i,j)=1./ns*hjave(i,j)enddoenddodoi=1,144doj=1,73doit=1,20hjjp1(i,j,it)=h
7、j1(i,j,it)-hjave(i,j)enddoenddoenddodoi=1,144doj=1,73doit=1,17hjjp2(i,j,it)=hj2(i,j,it)-hjave(i,j)enddoenddoenddodoi=1,144doj=1,73doit=1,18hjjp3(i,j,it)=hj3(i,j,it)-hjave(i,j)enddoenddoenddodoi=1,144doj=1,73jppj1(i,j)=0doit=1,20jppj1(i,j)=jppj1(i,j)+hjjp1(i,j,it)enddojppj1(i,j)
8、=1./a*jppj1(i,j)enddoenddodoi=1,144doj=1,73jppj2(i,j)=
