短期气候预测实习报告四.docVIP

.. .. 南京信息工程大学实验（实习）报告 实验（实习）名称 短期气候预测实习 日期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年雨型分类表 1 0 0 一类雨型 0 1 0 二类雨型 0 0 1 三类雨型 实习方法及步骤： （1）降水距平百分率： 为某年夏季降水量， 为1971-2000年夏季降水多年平均值 （2）总体均值的t统计量： ，分别代表样本均值和标准差，为总体均值，n为样本量 根据公式及资料， 一、二、三类雨型的样本数分别为20、17和18。我们认为简单认为它们通过0.01显著性水平的t值分别为2.86，2.92，2.90，通过0.05显著性水平的t值都为2.09，2.12，2.11。 （3）编写程序 program EP real a6(160,60),a7(160,60),a8(160,60),lat(160),h_y(144,73,12,61) real lon(160),rap(160,3),ddi(55,3),rsum(160),r(160,55) real h_m(144,73,12*61),h_p(144,73,3),aa,have(144,73) real ave_p(144,73,3),std(144,73,3),h_p2(144,73,3),sig(144,73,3) integer i,j,k,it character*8 id(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，三类18CCCCCCCCCCCCCCCCCCc do it=1,61 do k=1,12 do j=1,73 do i=1,144 read(2) h_y(i,j,k,it) h_m(i,j,12*(it-1)+k)=h_y(i,j,k,it) enddo enddo enddo enddo read(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) do i=1,160 read(6,*)lat(i),lon(i) enddo !CCCCCCCC t=1951-2005 do it=1,55 read(7,*)(ddi(it,j),j=1,