气象统计方法实习报告 下载本文

1 气象统计方法课程实践内容2013

目录

实习一 求500hPa高度场气候场、距平场和均方差场 ------------------------------------------------- 3

1、资料介绍 ------------------------------------------------------------------------------------------------- 3 2.要求 -------------------------------------------------------------------------------------------------------- 3 3、实习结果 ------------------------------------------------------------------------------------------------- 3 1)、FORTRAN源程序 ---------------------------------------------------------------------------------- 3

(2)、grads文件 --------------------------------------------------------------- 5 (3)、实习结果------------------------------------------------------------------------------------- 6

实习二 计算给定数据资料的简单相关系数和自相关系数 ------------------------------------------- 11

1、资料介绍 ------------------------------------------------------------------------------------------------ 11 2、要求 ------------------------------------------------------------------------------------------------------- 11 3、实习结果 ------------------------------------------------------------------------------------------------ 12

(1)、Fortran源程序 ----------------------------------------------------------------------------- 12 (2)、程序运行结果: -------------------------------------------------------------------------- 14

实习三 分析中国夏季降水线性趋势的分布特征 -------------------------------------------------------- 15

1.资料介绍及要求: -------------------------------------------------------------------------------------- 15 2.实习结果 ------------------------------------------------------------------------------------------------- 15

(1).matlab程序 ----------------------------------------------------------------------------------- 15 (2).程序运行结果 -------------------------------------------------------------------------------- 16

实习四 求给定数据的一元线性回归方程 ----------------------------------------------------------------- 17

1、资料介绍及要求 --------------------------------------------------------------------------------------- 17 2、实习结果 ------------------------------------------------------------------------------------------------ 18

第 1 页 共 36 页

2 气象统计方法课程实践内容2013

(1)、MATLAB程序 ---------------------------------------------------------------------------- 18 (2)、程序运行结果 ------------------------------------------------------------------------------ 18 (3)、结果分析 ------------------------------------------------------------------------------------ 19

实习五 对给定的海温数据进行EOF分析 ---------------------------------------------------------------- 21

1、资料介绍 ------------------------------------------------------------------------------------------------ 21 2、要求 ------------------------------------------------------------------------------------------------------- 21 3、实习结果: --------------------------------------------------------------------------------------------- 21

(1)、FORTRAN源程序 ------------------------------------------------------------------------ 21 (2)空间场和时间序列的ctl文件 ------------------------------------------------------------ 23 (3)运行结果 --------------------------------------------------------------------------------------- 24 (4)分析 --------------------------------------------------------------------------------------------- 26

实习三(附加) 计算给定数据的11年滑动平均和累积距平 -------------------------------------- 28

1、资料介绍 ------------------------------------------------------------------------------------------------ 28 2、要求 ------------------------------------------------------------------------------------------------------- 28 3、实习结果 ------------------------------------------------------------------------------------------------ 28 实习四(附加) 求给定数据的多元线性回归方程 ----------------------------------------------------- 31

1、说明 ------------------------------------------------------------------------------------------------------- 31 2、要求 ------------------------------------------------------------------------------------------------------- 31 3、实习结果: --------------------------------------------------------------------------------------------- 31

(1)Matlab源程序 -------------------------------------------------------------------------------- 31 (2)运行结果 --------------------------------------------------------------------------------------- 35

第 2 页 共 36 页

3 气象统计方法课程实践内容2013

实习一 求500hPa高度场气候场、距平场和均方差场

1、资料介绍

有一500hPa高度场资料,文件名h500.dat,范围:60~150E,0~40N. 时段:1982.1~1985.12共48个月。水平分辨率:2.5*2.5,格点数:37*17。

2.要求

编fortran程序,求500hPa高度场的 (1)气候场; (2)距平场; (3)均方差场。

并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图,图注要清楚,即要注明是哪个时间的图形,并做简单分析。

注:h500.For给出了如何用fortran读取ASCII码资料h500.dat.

3、实习结果

1)、FORTRAN源程序

program ex_grads

implicit none

第 3 页 共 36 页

4 气象统计方法课程实践内容2013

integer,parameter::nx=37,ny=17,nz=4,nt=12 integer i,j,iz,it

real var(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt)

! Opening file

open(10,file='g:\\gradsdata\\h500.dat') do iz=1,nz do it=1,nt read(10,1000)

read(10,3000) ((var(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo

1000 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) close(10) !Output

open(16,file='g:\\gradsdata\\h500.grd',form='binary') do iz=1,nz do it=1,nt

write(16) ((var(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo

!Calculating the Climatological Field do it=1,nt do i=1,nx do j=1,ny sum=0

do iz=1,nz

sum=sum+var(i,j,iz,it) enddo

cl(i,j,it)=sum/4 enddo enddo

enddo

! Output climate-file

open(12,file='g:\\gradsdata\\climate.grd',form='binary') do it=1,nt

write(12) ((cl(i,j,it),i=1,nx),j=1,ny) enddo

第 4 页 共 36 页

5 气象统计方法课程实践内容2013

!Calculating the Anomaly

do iz=1,nz do it=1,nt do i=1,nx do j=1,ny

jp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it) enddo enddo enddo enddo

open(13,file='g:\\gradsdata\\anomaly.grd',form='binary')

!Output anomaly-file do iz=1,nz do it=1,nt

write(13) ((jp(i,j,iz,it),i=1,nx),j=1,ny) enddo enddo

!Calculating the Mean-square Deviation do it=1,nt do i=1,nx do j=1,ny sum=0

do iz=1,nz

sum=sum+(jp(i,j,iz,it))**2 enddo

jfc(i,j,it)=sqrt(sum/4) enddo enddo

enddo

!Output mean-square deviation-file

open(14,file='g:\\gradsdata\\deviation.grd',form='binary') do it=1,nt

write(14) ((jfc(i,j,it),i=1,nx),j=1,ny) enddo end

(2)、grads文件

'open g:\\gradsdata\\***.ctl' (***为所求场对应的ctl文件名) 'set lat 0 40'

第 5 页 共 36 页