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 页