资源描述：
《短期气候预测实习报告材料4—预测因子选择》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用文档短期气候预测实习报告实习4预测因子选择—合成分析法标准文案实用文档南京信息工程大学实验（实习）报告实验(实习)名称短期气候预测基础课程设计日期2016.10指导教师李忠贤系大气科学专业大气科学班次长望班姓名夏新露一、实习目的掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流异常在短期气候预测中物理机制的认识。二、实习要求熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。三、实习内容计算1951-2005年夏季三类年雨型合成图；计算各类雨型的前期冬季高度场距平合成图，指
2、出可能出现的遥相关型。资料：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）降水距平百分
3、率：某年夏季降水量为1971-2000年夏季降水多年平均值（2）总体均值的t统计量：、分别代表样本均值和标准差为总体均值为样本量标准文案实用文档一、二、三类雨型的样本数分别为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,1
4、2,65)realhj1(144,73,20),hj2(144,73,17),hj3(144,73,18),hjave(144,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,7
5、3doi=1,144read(16)hh(i,j,k,it)enddo;enddo;enddo;enddoread(15,*)((ddi(i,j),j=1,3),i=1,55)print*,((ddi(i,j),j=1,3),i=1,55)a=0b=0c=0标准文案实用文档doi=1,55if(ddi(i,1)==1)thena=a+1标准文案实用文档elseif(ddi(i,2)==1)thenb=b+1elsec=c+1endifenddoprint*,a,b,ce=0f=0g=0doi=1,55if(ddi(i,1)==1
6、)thene=e+1doix=1,144doiy=1,73hj1(ix,iy,e)=hh(ix,iy,1,i)enddoenddoelseif(ddi(i,2)==1)thenf=f+1doix=1,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
7、,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)=hj1(i,j,it)-hjave(i,j)enddoenddo标准文案实用文档enddodoi=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,
8、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)=1
