短期气候预测基础实习四.docVIP

实习四：预测因子的选择——合成分析方法（1） 一、实习目的 掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流异常在短期气候预测中物理机制的认识 。 二、实习要求 熟悉资料和方法程序（提供部分子程序）； 使用GRADS对结果进行绘图输出和文字分析，完成实习报告。 三、实习内容 计算1951-2001年夏季三类年雨型合成图； 计算各类雨型的前期冬季高度场距平合成图，指出可能出现的遥相关型。 四、资料介绍 A、大气环流场资料： NCEP/NCAR再分析资料1948－2008年（61年）的500百帕月平均高度场资料；资料范围为（90°S-90°N,0°-360°E)；网格距为2.50×2.50；纬向格点数为144，经向格点数为73。 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 B、站点资料： 国家气候中心整编的6、7、8月降水量资料（时间段：1951~2010年，资料的格式参见readme.txt文件） C. 1951~2005年雨型分类表:ddi 1 0 0 一类雨型 0 1 0 二类雨型 0 0 1 三类雨型 五、方法介绍 1、降水距平百分率： Ri某年夏季降水量； 为1971-2000年夏季降水多年平均值 2、总体均值的t统计量： 分别代表样本均值和标准差 为总体均值 n为样本量 一、二、三类雨型的样本数分别为20、17和18。我们认为简单认为它们通过0.01显著性水平的t值都为2.88，通过0.05显著性水平的t值都为2.10。 六、实习步骤： 编程计算并绘制1951－2005年夏季三类雨型年合成图； 编程计算并绘制前期冬季500百帕高度场距平合成图，指出可能存在的遥相关型环流特征。 A.源程序: 1.三类雨型及前期高度场合成; program main implicit none real,dimension(160)::sumr,aver,sumrp1,sumrp2,sumrp3,rp1,rp2,rp3,sumrp,averp,sumrp11,sumrp22,sumrp33,jfangcha1,jfangcha2,jfangcha3,t1,t2,t3 integer::rsty(55,3),i,j,it,imo real,dimension(160,55)::r6,r7,r8,summer,rp real,dimension(144,73,55)::winter,jph real::a(144,73,12,61) real,dimension(144,73)::aveh,sumh,sumjph1,sumjph2,sumjph3,jph1,jph2,jph3,sumjph11,sumjph22,sumjph33,jfangcha11,jfangcha22,jfangcha33,th1,th2,th3,sumjph,avejph open(1,file=g:\r1606) open(2,file=g:\r1607) open(3,file=g:\r1608) open(4,file=g:\ddi) open(5,file=g:\rp1.txt) open(6,file=g:\rp2.txt) open(7,file=g:\rp3.txt) open(8,file=g:\t1.txt) open(9,file=g:\t2.txt) open(10,file=g:\t3.txt) open(11,file=g:\h5004808.grd,form=binary) open(12,file=g:\th1.grd,form=binary) open(13,file=g:\th2.grd,form=binary) open(14,file=g:\th3.grd,form=binary) open(15,file=g:\jph1.grd,form=binary) open(16,file=g:\jph2.grd,form=binary) open(17,file=g:\jph3.grd,form=binary) read(4,*) ((rsty(i,j),j=1,3),i=1,55) read(1,*) ((r6(i,j),i=1,160),j=1,55) read(2,*) ((r7(i,j),i=1,160),j=1,55) read(3,*) ((r8(i,j),i=1,160),j=1,55) do it=1,61