短期气候预测实习报告五.doc

南京信息工程大学实验（实习）报告 实验（实习）名称 短期气候预测实习 日期11.26得分 指导教师 系 ~ 专业 ~ 年级 ~ 班次~姓名 Trichtu 学号~ 实习目的： 掌握短期气候预测因子的分析和选择，加深对外强迫因子（海温场）在短期气候预测中物理机制的认识 。 实习要求： 熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。 实习内容： 计算前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值及T检验，确定关键区，北太平洋范围（120°E~60°W，10°S~60°N） 计算前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值及T检验，确定关键区。 实习资料： （1）NOAA 延伸再建的海表面温度资料时间范围： 1947年1月－2010年2月 范围（88°S-88°N,0-358°E) 网格距2°×2°,纬向格点数180，经向格点数89 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 注意：陆地上为缺测值32767 （2）1951~2005年雨型分类表 1 0 0 一类雨型 0 1 0 二类雨型 0 0 1 三类雨型 实习方法： (1)不同样本之间的均值t统计量： 分别代表第i类样本均值和标准差， 为第i类样本的样本量 (2) 总体均值的t统计量： 一、二、三类雨型的样本数分别为20、17和18。我们简单认为任意两类之间通过0.01显著性水平的t值分别为2.73，通过0.05显著性水平的t值都为2.01 (3) 编写程序 program sstt real sst(180,89,12,61),ddi(55,3),sst_m(180,89,12*61),sstave(180,89),sst_p(180,89,3) real cha(180,89,3),sig(180,89,3),ave_p(180,89,3),std(180,89,3),sst_p2(180,89,3) integer i,j,k,it character*8 id(160) open(2,file=e:\copy5\data\sst.grb,form=binary) open(7,file=e:\copy5\data\ddi) open(9,file=e:\copy5\sst_p.grd,form=binary) open(10,file=e:\copy5\cha.grd,form=binary) open(11,file=e:\copy5\sig.grd,form=binary) !ccccccccccccccc 读数据（经纬度、160站降水1948-2008、sst（1947,1-2010,2 ）） do it=1,61 do k=1,12 do j=1,89 do i=1,180 read(2) sst(i,j,k,it) if(sst(i,j,k,it)==32767)then sst(i,j,k,it)=0 endif sst_m(i,j,12*(it-1)+k)=sst(i,j,k,it) enddo enddo enddo enddo !CCCCCCCC t=1951-2005 do it=1,55 read(7,*)(ddi(it,j),j=1,3) enddo write(*,*)read data ok !ccccccccccccccc 编程求合成 !ccccccccccccc前期海温场 do k=1,3 do i=1,180 do j=1,89 aa=.0 bb=.0 num=0 do it=1,55 num=num+ddi(it,k) aa=aa+sst_m(i,j,12*(it+2)+12)*ddi(it,k) sst_p(i,j,k)=aa bb=bb+sst_m(i,j,12*(it+2)+12) sstave(i,j)=bb end do ave_p(i,j,k)=sst_p(i,j,k)/num num=0 do it=1,55 num=num+ddi(it,k) aa=(sst_m(i,j,12*(it+2)+12)-ave_p(i,j,k))*