短期气候预测实习报告4—预测因子选择.docxVIP

短期气候预测实习报告4—预测因子选择

LT

（1）降水距平百分率：

某年夏季降水量

为1971-2000年夏季降水多年平均值

（2）总体均值的t统计量：

、分别代表样本均值和标准差

为总体均值

为样本量

一、二、三类雨型的样本数分别为22、19和19。我们认为简单认为它们通过0.01显著性水平的t值分别为2.85，2.90，2.90，通过0.05显著性水平的t值都为2.09，2.11，2.11。

实习步骤

程序编写：

三类雨型合成图程序：

realns

parameter(ns=55)

integerddi(55,3)

integere,f,g

reala,b,c

realhh(144,73,12,65)

realhj1(144,73,20),hj2(144,73,17),hj3(144,73,18),hjave(144,73)

realhjjp1(144,73,20),hjjp2(144,73,17),hjjp3(144,73,18)

realjppj1(144,73),jppj2(144,73),jppj3(144,73)

open(15,file=ddi)

open(16,file=hgt48-12500.grd, form=binary,STATUS=OLD)

doit=1,65

dok=1,12

doj=1,73

doi=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

doi=1,55

if(ddi(i,1)==1)then

a=a+1

elseif(ddi(i,2)==1)then

b=b+1

else

c=c+1

endif

enddo

print*,a,b,c

e=0

f=0

g=0

doi=1,55

if(ddi(i,1)==1)then

e=e+1

doix=1,144

doiy=1,73

hj1(ix,iy,e)=hh(ix,iy,1,i)

enddo

enddo

elseif(ddi(i,2)==1)then

f=f+1

doix=1,144

doiy=1,73

hj2(ix,iy,f)=hh(ix,iy,1,i)

enddo

enddo

else

g=g+1

doix=1,144

doiy=1,73

hj3(ix,iy,g)=hh(ix,iy,1,i)

enddo

enddo

endif

enddo

doi=1,144

doj=1,73

hjave(i,j)=0

doit=1,55

hjave(i,j)=hjave(i,j)+hh(i,j,1,it)

enddo

hjave(i,j)=1./ns*hjave(i,j)

enddo

enddo

doi=1,144

doj=1,73

doit=1,20

hjjp1(i,j,it)=hj1(i,j,it)-hjave(i,j)

enddo

enddo

enddo

doi=1,144

doj=1,73

doit=1,17

hjjp2(i,j,it)=hj2(i,j,it)-hjave(i,j)

enddo

enddo

enddo

doi=1,144

doj=1,73

doit=1,18

hjjp3(i,j,it)=hj3(i,j,it)-hjave(i,j)

enddo

enddo

enddo

doi=1,144

doj=1,73

jppj1(i,j)=0

doit=1,20

jppj1(i,j)=jppj1(i,j)+hjjp1(i,j,it)

enddo

jppj1(i,j)=1./a*jppj1(i,j)

enddo

enddo

doi=1,144

doj=1,73

jppj2(i,j)=0

doit=1,17

jppj2(i,j)=jppj2(i,j)+hjjp2(i,j,it)

enddo

jppj2(i,j)=1./b*jppj2(i,j)

enddo

enddo

doi=1,144

doj=1,73

jppj3(i,j)=0

doit=1,18

jppj3(i,j)=jppj3(i,j)+hjjp3(i,j,it)