短期气候实习报告六.doc

南京信息工程大学短期气候实验（实习）报告 实验名称 夏季区域降水的定量预测 日期 2017.6.6得分 指导教师 系 大气科学 专业 大气科学 年级 2014级 班次 2姓名 车楚玉 学号 20141301043 目的要求： 目的 ：掌握短期气候预测中物理统计预测的基本步骤。 要求 ：能运用提供的资料和方法子程序，编写或补充完成程序当中的部分片断，了解区域降水的预测方法及其建立过程，输出实验要求的相应结果，并就方法对区域降水的拟合及试验预测效果进行分析 。 资料和方法： 资料： 1、前期1月的Nino3.4指数（来自CPC）； 2、西太平洋副高脊线、西太平洋副高西伸脊点、亚洲极涡面积、南方涛动指数（来自中国气象局整编的74个环流指数）； 3、夏季华北区域10站的降水量，距平百分率。 方法： 回归分析(mregrssion.for)是用来寻找若干变量之间统计关系的一种方法，利用所找到的统计关系对某一变量作出未来时刻的估计，称为回归预报值。 实习步骤： ① 说明所用资料方法； ② 计算方法的简单介绍； ③ 输出反映回归效果的参数及回归系数，并就相关参数分析回归效果； ④预测量与回归方程计算的估计值和观测值的历年曲线变化图（1952~2001年），并附简单的说明； ⑤输出独立预测试验的观测与预测值。 结果（图表并解释说明）： 程序如下： PROGRAM MAIN INTEGER,PARAMETER::N=50 INTEGER,PARAMETER::K=5 REAL,DIMENSION(K,N)::X REAL,DIMENSION(N)::Y REAL,DIMENSION(K+1)::A REAL,DIMENSION(K+1,K+1)::B REAL,DIMENSION(K)::V REAL Q,S,R,U,ind(6,60),year(50),y1(50),pre(7) integer i,j C x(5,50) y(50) a(6) b(6,6),v(5) C OPEN THE INPUT DATA FILE OPEN(10,FILE=D:\duanqi\shixi6\shixi.txt) C READ THE DATA and give data to X and Y do i=1,50 read(10,*) year(i),x(1:5,i),y(i) enddo MM=K+1 call DYHG(X,Y,K,MM,N,A,Q,S,R,V,U,B,DYY) write(*,88) A(1) 88 format(/1x,b 0=,f19.5) do 89 j=2,MM 89 write(*,100) j-1,A(j) 100 format(1x,b,i2,=,f9.5) open(11,FILE=D:\duanqi\shixi6\guji.txt) open(12,FILE=D:\duanqi\shixi6\guji.grd,form=binary) open(13,FILE=D:\duanqi\shixi6\prediction.txt) do j=1,50 y1(j)=a(1)+a(2)*x(1,j)+a(3)*x(2,j)+a(4)*x(3,j)+a(5)*x(4,j)+a(6) * *x(5,j) enddo pre(1)=a(1)+a(2)*26.50+a(3)*13.00+a(4)*100.00+a(5)*196.00+a(6)* * 2.00 pre(2)=a(1)+a(2)*27.76+a(3)*12.00+a(4)*90.0+a(5)*197.00+a(6)*(-1.00) pre(3)=a(1)+a(2)*26.74+a(3)*12.0+a(4)*120.00+a(5)*199.00+a(6)* *(-11.00) pre(4)=a(1)+a(2)*27.10+a(3)*13.00+a(4)*110.00+a(5)*220.00+a(6)*3.00 pre(5)=a(1)+a(2)*25.64+a(3)*14.00+a(4)*105.00+a(5)*215.00+a(6)*13.00 pre(6)=a(1)+a(2)*27.26+a(3)*15.00+a(4)*90.00+a(5)*185.00+a(6)* *(-8.00) pre(7)=a(