短期气候预测实习报告四

南京信息工程大学实验（实习）报告南京信息工程大学实验（实习）报告 实验（实习）名称 短期气候预测实习 日期 11.15 得分指导教师 系 专业 年级 班次姓名 Trichtu 学号 实习目的实习目的掌握短期气候预测因子的分析和选择， 加深对夏季降水分布、 环流异常在短期气 候预测中物理机制的认识 。 实习要求实习要求; ; 熟悉资料和方法程序（提供部分子程序） ；使用 GRADS 对结果进行绘图输出和文 字分析，完成实习报告。 实习内容实习内容 计算 1951-2005 年夏季三类年雨型合成图； 计算各类雨型的前期冬季高度场距平合成图，指出可能出现的遥相关型。 资料NCEP/NCAR 再分析资料 1948－2008 年 1～12 月的 500 百帕月平均高度场资料 范围（90S-90N,0-360E 网格距 2.52.5,纬向格点数 144，经向格点数 73 资料为 GRD 格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 、 实习资料实习资料 国家气候中心整编的 6、7、8 月降水量资料（时间段19512010 年，资料的格式参见 readme.txt 文件） 19512005 年雨型分类表 1 0 0一类雨型 0 1 0二类雨型 0 0 1三类雨型 实习方法及步骤实习方法及步骤 （1）降水距平百分率 R i R 100RP i  R R i 为某年夏季降水量， R 为 1971-2000 年夏季降水多年平均值 （2）总体均值的 t 统计量 0 x  t  s n 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。 （3）编写程序 program EP real a6160,60,a7160,60,a8160,60,lat160,h_y144,73,12,61 real lon160,rap160,3,ddi55,3,rsum160,r160,55 real h_m144,73,12*61,h_p144,73,3,aa,have144,73 real ave_p144,73,3,std144,73,3,h_p2144,73,3,sig144,73,3 integer i,j,k,it character*8 id160 open2,filee\copy4\data\hgt500.grd,binary open3,filee\copy4\data\r1606.txt open4,filee\copy4\data\r1607.txt open5,filee\copy4\data\r1608.txt open6,filee\copy4\data\lat_lon.txt open7,filee\copy4\data\ddi open8,filee\copy4\rap.grd,binary open9,filee\copy4\hp.grd,binary open10,filee\copy4\sig.grd,binary ccccccccccccccc 读数据经纬度、160 站降水 1948-2008、雨型、height1948-2008 CCCCCCCCCCCCCCCCCCCC一类 20，二类 17，三类 18CCCCCCCCCCCCCCCCCCc do it1,61 do k1,12 do j1,73 do i1,144 read2 h_yi,j,k,it h_mi,j,12*it-1kh_yi,j,k,it enddo enddo enddo enddo read3,*a6i,j,i1,160,j1,60 read4,*a7i,j,i1,160,j1,60 read5,*a8i,j,i1,160,j1,60 do i1,160 read6,*lati,loni enddo CCCCCCCC t1951-2005 do it1,55 read7,*ddiit,j,j1,3 enddo write*,*read data ok ccccccccccccccc 编程求合成 do i1,160 rsumi0.0 do it21,50 rsumirsumia6i,ita7i,ita8i,it enddo rsumirsumi/30.0 do it1,55 ri,ita6i,ita7i,ita8i,it-rsumi*100/rsumi enddo enddo do i1,160 do j1,3 num0 do it1,55 numnumddiit,j rapi,jrapi,jri,it*ddiit,j enddo rapi,jrapi,j/realnum enddo;enddo ccccccccccccc 前期高度场 do k1,3 do i1,144 do j1,73 aa.0 bb.0 num0 do it1,55 numnumddiit,k aaaah_mi,j,12*it112h_mi,j,12*it21h_mi,j,12*it22*ddiit, k h_pi,j,kaa/3.0 bbbbh_mi,j,12*it112h_mi,j,12*it21h_mi,j,12*it22 havei,jbb/3.0 end do ave_pi,j,kh_pi,j,k/num num0 do it1,55 numnumddiit,k aah_mi,j,12*it112h_mi,j,12*it21h_mi,j,12*it22*ddiit,k/ 3.0-ave_pi,j,k*ddiit,k stdi,j,kstdi,j,kaa**2 end do havei,jhavei,j/55.0 h_pi,j,kh_pi,j,k/num-havei,j h_p2i,j,kh_pi,j,k/sqrtstdi,j,k/num*sqrtnum/1.0 write*,* h_p2i,j,k end do end