短期气候预测实习报告五

实验（实习）名称 短期气候预测实习 日期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统计量： x1x2 t2(N11)S12(N21)S211 N1N22N1N2

Xi Si 分别代表第i类样本均值和标准差， Ni为第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))*ddi(it,k) std(i,j,k)=std(i,j,k)+aa**2 end do

sstave(i,j)=sstave(i,j)/55.0

sst_p(i,j,k)=sst_p(i,j,k)/num-sstave(i,j) end do end do end do

do i=1,180 do j=1,89

cha(i,j,1)=sst_p(i,j,1)-sst_p(i,j,2) cha(i,j,2)=sst_p(i,j,1)-sst_p(i,j,3) cha(i,j,3)=sst_p(i,j,2)-sst_p(i,j,3) end do end do

!CCCCCCCCCCCCCCCCCCCC双变量t检验

do i=1,180 do j=1,89

aa=sqrt((19*std(i,j,1)/20+16*std(i,j,2)/17)/35.0)

sst_p2(i,j,1)=(ave_p(i,j,1)-ave_p(i,j,2))/sqrt(1/20+1/17.0)/aa write(*,*) sst_p2(i,j,1)

aa=sqrt((19*std(i,j,1)/20+17*std(i,j,2)/18)/36.0)

sst_p2(i,j,2)=(ave_p(i,j,1)-ave_p(i,j,3))/sqrt((1/20+1/18.0))/aa aa=sqrt((16*std(i,j,2)/17+17*std(i,j,3)/18)/33.0)

sst_p2(i,j,3)=(ave_p(i,j,2)-ave_p(i,j,3))/sqrt((1/17+1/18.0))/aa end do

end do

do k=1,3 do i=1,180 do j=1,89

select case(k) case(1)

if(abs(sst_p2(i,j,k))>=2.44)then sig(i,j,k)=2

elseif(abs(sst_p2(i,j,k))>=1.69)then sig(i,j,k)=1 endif case(2)

if(abs(sst_p2(i,j,k))>=2.44)then sig(i,j,k)=2

elseif(abs(sst_p2(i,j,k))>=1.69)then sig(i,j,k)=1 endif case(3)

if(abs(sst_p2(i,j,k))>=2.44)then sig(i,j,k)=2

elseif(abs(sst_p2(i,j,k))>=1.69)then sig(i,j,k)=1 endif end select end do end do end do

!其实都一样，就可以不用select了 !ccccccccccccccccccc写数据 do k=1,3

write(9) ((sst_p(i,j,k),i=1,180),j=1,89) write(10) ((cha(i,j,k),i=1,180),j=1,89) write(11) ((sig(i,j,k),i=1,180),j=1,89) end do end

（4）结果输出

1绘制绘制I、II和III类雨型的前期12月北太平洋海温距平合成图。

夏季降水第一类雨型前期12月北太平洋SSTA合成图（单位：℃）

夏季降水第二类雨型前期12月北太平洋SSTA合成图（单位：℃）

夏季降水第三类雨型前期12月北太平洋SSTA合成图（单位：℃）

2绘制绘制前期12月北太平洋海温I和II类雨型合成差值、I和III类合成差值、II和III类合成差值图（包括它们的t检验）

€ 夏季降水第一、二类雨型前期12月北太平洋SST差值合成图及其t检验（单位： ℃ ，

橙、红阴影表示通过0.05和0.01显著性水平t检验值)

前期12月北太平洋SST，一类较二类，从西北太平洋到北太平洋海温更高，而东太平洋的海温偏低。

€ 夏季降水第一、三类雨型前期12月北太平洋SST差值合成图及其t检验（单位： ℃ ，

橙、红阴影表示通过0.05和0.01显著性水平t检验值)

前期12月北太平洋SST，二类较三类，其差值场有从东太平洋至白令海峡延美洲东海岸的波裂。

€ 夏季降水第二、三类雨型前期12月北太平洋SST差值合成图及其t检验（单位： ℃ ，

橙、红阴影表示通过0.05和0.01显著性水平t检验值)

前期12月北太平洋SST，二类较三类，东太平洋的低纬度的海温偏低，而北太平洋在阿留申群岛附近的海温偏高。