短期气候预测实习报告
实习4预测因子选择—合成分析法
南京信息工程大学 实验（实习）报告
实验(实习)名称短期气候预测基础课程设计日期2016.10 指导教师李忠贤 系大气科学专业大气科学班次长望班姓名夏新露 一、 实习目的
掌握短期气候预测因子的分析和选择，加深对夏季降水分布、环流异常在短期气候预测中物理机制的认识 。
二、 实习要求
熟悉资料和方法程序（提供部分子程序）；使用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）降水距平百分率：
Ri?RRP?100%i?R
Ri某年夏季降水量
R为1971-2000年夏季降水多年平均值
（2）总体均值的t统计量：
x??0t?ns
x、s分别代表样本均值和标准差
?0 为总体均值 n 为样本量
一、二、三类雨型的样本数分别为22、19和19。我们认为简单认为它们通过0.01显著性水平的t值分别为2.85，2.90，2.90，通过0.05显著性水平的t值都为2.09，2.11，2.11。
六、 实习步骤
程序编写：
（1） 三类雨型合成图程序： real ns
parameter(ns=55) integer ddi(55,3) integer e,f,g real a,b,c
real hh(144,73,12,65)
real hj1(144,73,20),hj2(144,73,17),hj3(144,73,18),hjave(144,73) real hjjp1(144,73,20),hjjp2(144,73,17),hjjp3(144,73,18) real jppj1(144,73),jppj2(144,73),jppj3(144,73) open(15,file='ddi')
open(16,file='hgt48-12500.grd', form='binary',STATUS='OLD') do it=1,65 do k=1,12 do j=1,73 do i=1,144
read(16)hh(i,j,k,it) enddo;enddo;enddo;enddo
read(15,*)((ddi(i,j),j=1,3),i=1,55) print*,((ddi(i,j),j=1,3),i=1,55) a=0 b=0 c=0
do i=1,55
if(ddi(i,1)==1)then a=a+1
else if(ddi(i,2)==1)then b=b+1 else c=c+1 endif enddo
print*,a,b,c e=0 f=0 g=0
do i=1,55
if(ddi(i,1)==1)then e=e+1
do ix=1,144 do iy=1,73
hj1(ix,iy,e)=hh(ix,iy,1,i) end do end do
elseif(ddi(i,2)==1)then f=f+1
do ix=1,144 do iy=1,73
hj2(ix,iy,f)=hh(ix,iy,1,i) enddo enddo else g=g+1
do ix=1,144 do iy=1,73
hj3(ix,iy,g)=hh(ix,iy,1,i) enddo enddo endif enddo
do i=1,144 do j=1,73 hjave(i,j)=0 do it=1,55
hjave(i,j)=hjave(i,j)+hh(i,j,1,it) enddo
hjave(i,j)=1./ns*hjave(i,j) enddo enddo do i=1,144 do j=1,73 do it=1,20
hjjp1(i,j,it)=hj1(i,j,it)-hjave(i,j) enddo enddo enddo
do i=1,144 do j=1,73 do it=1,17
hjjp2(i,j,it)=hj2(i,j,it)-hjave(i,j) enddo enddo enddo
do i=1,144 do j=1,73 do it=1,18
hjjp3(i,j,it)=hj3(i,j,it)-hjave(i,j) enddo enddo enddo
do i=1,144 do j=1,73 jppj1(i,j)=0 do it=1,20
jppj1(i,j)=jppj1(i,j)+hjjp1(i,j,it) enddo
jppj1(i,j)=1./a*jppj1(i,j) enddo enddo
do i=1,144 do j=1,73 jppj2(i,j)=0 do it=1,17
jppj2(i,j)=jppj2(i,j)+hjjp2(i,j,it) enddo
jppj2(i,j)=1./b*jppj2(i,j) enddo enddo
do i=1,144 do j=1,73 jppj3(i,j)=0 do it=1,18
jppj3(i,j)=1./c*jppj3(i,j) enddo enddo
open(21,file='g:\\shixisi\\gedian.grd',form='binary') do j=1,73 do i=1,144
write(21)jppj1(i,j) enddo enddo do j=1,73
do i=1,144
write(21)jppj2(i,j) end
jppj3(i,j)=jppj3(i,j)+hjjp3(i,j,it) enddo
enddo enddo do j=1,73 do i=1,144
write(21)jppj3(i,j) enddo enddo
write(*,*)stid(i),lat(i),lon(i),tim,nlev,nflag,cor(i)
write(*,*)stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
七、 结果输出：
1.1951－2010年夏季1类雨型年合成图：
分析：主要多雨带位于黄河流域及其以北地区，江淮流域大范围少雨，梅雨