短期气候预测基础实习三

实习三：大气遥相关

一、实习目的及要求

掌握大气环流中遥相关型指数的计算及其与大气环流和我国气候关系的分析。

要求运用资料，计算北半球1月遥相关指数，并分析它与环流和我国气候变化的关系；用图形输出指数年际变化曲线、遥相关的空间分布以及与我国气温的相关系数分布，正确分析结果数据，完成实习报告。

二、实习内容

1、计算EU（欧亚）遥相关指数，输出1月份该指数年际变化的时间序列；

2、计算EU遥相关指数与同期环流场（500hPa高度场或海平面气压场）的相关系数；

3、计算EU遥相关指数与同期我国气温的相关系数。

三、实习资料

大气环流场资料：

1、NCEP/NCAR 1948－2008年（61年）的500百帕月平均高度场资料

2、资料范围为（900S-900N,00-3600E)

3、网格距为2.50×2.50,纬向格点数为144，经向格点数为73

4、资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。

站点资料：

1、 我国气候中心整编的160站月平均气温资料；

2、 全国160个台站；

3、 所给的资料是1月份的；

4、 资料为txt格式，参加资料说明。

四、实习步骤

1、编程计算1月份EU遥相关指数（标准化），画出指数图

编写fortran程序：

program main

parameter(nx=144,ny=73,n=nx*ny,nt=61)

real

ind(nt),z(n,nt),eu(nt),hh(144,73,12,61),avey(n),sy(n),r1(n),r(n)

real x,avex,sx

open(10,file='d:\\nyclimate\\sh3\\data\\hgt500.grd',form='binary')

open(20,file='d:\\nyclimate\\sh3\\ind.grd',form='binary')

open(30,file='d:\\nyclimate\\sh3\\eu.grd',form='binary')

open(40,file='d:\\nyclimate\\sh3\\rxy1.grd',form='binary')

do it=1,61

do k=1,12

do j=1,73

do i=1,144

read(10)hh(i,j,k,it)

enddo；enddo；enddo；enddo

do it=1,61

do j=1,73

do i=1,144

z((j-1)*144+i,it)=hh(i,j,1,it)

enddo；enddo；enddo

do it=1,61

eu(it)=-z(9+58*144,it)/4+z(31+58*144,it)/2-z(59+52*144,it)/4

enddo

write(20)(eu(it),it=1,61)

do it=1,61

x=x+eu(it)