短期实习报告五.docVIP

实习五：预测因子的选择——（合成分析方2） 实习目的 掌握短期气候预测因子的分析和选择，加深对外强迫因子（海温场）在短期气候预测中物理机制的认识 。 实习要求 熟悉资料和方法程序（提供部分子程序）；使用GRADS对结果进行绘图输出和文字分析，完成实习报告。 实习内容 （1）计算前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值及T检验，确定关键区。 北太平洋范围（120°E~60°W，10°S~60°N） （2）根据1951年以来的中国夏季(6－8月)降水距平百分率分布图，着眼于105oE以东的东部地区，把多雨区内降水距平百分率最大的区域作为主要雨带，划分出历年的主要雨型，作为夏季降水的主要预报对象。 四、实习步骤 （1） 资料：NOAA 延伸再建的海表面温度资料 时间范围： 1947年1月－2010年2月 范围（88°S-88°N,0-358°E) 网格距2°×2°,纬向格点数180，经向格点数89 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 注意：陆地上为缺测值32767 资料： 1951~2005年雨型分类表 1 0 0 一类雨型 0 1 0 二类雨型 0 0 1 三类雨型 （2）实习方法： ①不同样本之间的均值t统计量： 、分别代表第i类样本均值和标准差 为第i类样本的样本量 ②总体均值的t统计量： 一、二、三类雨型的样本数分别为20、17和18。我们简单认为任意两类之间通过0.01显著性水平的t值分别为2.73，通过0.05显著性水平的t值都为2.01。 （3）编写程序 i.所用的Fortran程序 program main implicit none real a(180,89,12,64),sst12(180,89,55),ssta(180,89,55)!ssta为海表温距平 real sstave,ddi(55,3),sstaa(180,89,3),num real sstat(180,89,55,3),t(180,89,3) real sas(180,89,3) integer k1,k2,kk1,kk2 integer i,j,k,it,imo open(6,file=e:\duanqi\shixi5\data\sst.grb,form=binary) open(7,file=e:\duanqi\shixi5\data\ddi) open(8,file=e:\duanqi\shixi5\result\sst12.grd,form=binary) !ccccccccccccccc 读取资料（海表面温度、雨型） do it=1,64 do imo=1,12 do j=1,89 do i=1,180 read(6,end=10)a(i,j,imo,it) enddo;enddo;enddo;enddo 10 continue do it=1,55 read(7,*)(ddi(it,j),j=1,3) enddo !ccccccccccccccc 编程求合成 !计算前期12月距平合成图 !先提取60年所有前期12月距平 sst12:年份往后移 do j=1,89 do i=1,180 do it=1,55 sst12(i,j,it)=a(i,j,12,it+4-1) enddo enddo enddo !算距平ssta:sst abnormal do j=1,89 do i=1,180 sstave=0 call cal_ave(sst12(i,j,:),60,sstave) do it=1,55 ssta(i,j,it)=sst12(i,j,it)-sstave enddo enddo enddo !挑年份,求合成sstaa: sst abnormal asemble do j=1,89 do i=1,180 if(ssta(i,j,1)==0)then sstaa(i,j,:)=32767 else do k=1,3 num=0.0 do it=1,55 num=num+ddi(it,k) sstaa(i,j,k)=sstaa(i,j,k)+ssta(i,j,it)*ddi(it,k) enddo sstaa(i,j,k)=sstaa(i,j,k)/num enddo endif enddo enddo !做