内容发布更新时间 : 2024/6/8 10:38:18星期一 下面是文章的全部内容请认真阅读。
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 页