grads 相关系数_实习-大气科学与环境气象试验教学中心.DOC

实习-大气科学与环境气象试验教学中心

实 习

实习一 大气环流状况的表征

在短期气候预测中，经常用大气环流场的时间平均图、距平图及纬偏图来表征大气环流的基本状况。下面以500hPa高度场的时间平均图、距平图及纬偏图为例，学习大气环流的基本状态的表征及分析方法。

(1)目的和要求

要求学生回顾FORTRAN语言和GRADS气象绘图系统的相关知识，在了解资料格式(特别要注意的是所读取资料的范围)基础上，应用FORTRAN语言编写程序，正确读写资料，计算500hPa高度场的时间平均、距平及纬偏值，并使用GRADS编写数据描述文件以及执行文件，对计算结果进行绘图并完成实习报告。目的是掌握大气环流基本状况的表征及分析方法。

(2)实习内容

计算1948-2007年(60年)1月的平均高度场，绘制环流平均图；

计算2008年1月的高度距平，绘制环流距平图(相对于1948-2007年共60年的平均)；

计算2008年1月的高度场纬偏值，绘制环流纬偏图。

(3)步骤

1)熟悉资料

资料为NCEP/NCAR 61年(1948-2008年)1～12月的500hPa月平均高度场资料，资料范围为(900S-900N；00E-3600E)，网格距为2.50*2.50，纬向格点数为144，经向格点数为73。资料为GRD格式，按照从南到北、自西向东的顺序排列，每月为一个记录，按年逐月排放，注意读取方式以及记录长度。

本次实习应用NCEP/NCAR 1948-2008年1月500百帕的平均高度场资料。

2)方法

① 求时间平均

，

其中(i，j)为经、纬格点，m为样本长度，k为年份。

② 求距平

距平图就是用给定时段的平均图减去同时段的多年平均图所得的差值图。距平图表示这时段的平均状况对其多年平均状况的偏差。

，

其中(i，j)为经、纬格点，m为样本长度，k为年份。

③ 求纬偏

纬偏图是用给定时段的平均图减去同时段的纬圈平均图所得的差值图。纬偏图上的槽脊位置与该时段的平均槽脊位置一致，且两者沿纬圈的相对强弱也是定性一致的。

其中，(i，j)为经、纬格点，m为纬圈格点数，k为年份。

3)编写程序

编写计算1948-2007年(共60年)1月500hPa高度平均程序；编写计算2008年1月500hPa高度场距平程序(相对于1948-2007年共60年的平均)程序；以及编写2008年1月500hPa高度场纬向平均和偏差程序；编写GRADS数据描述文件及执行文件。

4)结果输出

通过GRADS绘制图形。其中包括1948-2007年(共60年)1月500hPa高度平均图；2008年1月500hPa高度场距平图；2008年1月500hPa高度场纬偏图。

5)完成实习报告

① 说明所用资料方法；

② 绘制环流平均图、距平图及纬偏图；

③ 用所学的知识分析图表。

实习二 大气环流分型

为了更充分地描述环流场，短期气候预测中利用自然函数正交展开(EOF)方法来分析环流的结构特征。下面以欧亚地区500hPa高度场进行EOF方法为例，来学习大气环流的模态法。

(1)目的和要求

目的是要掌握大气环流分型的基本方法。要求学生熟悉EOF方法以及相关程序的应用；使用GRADS完成图形的绘制，并能正确分析结果数据，完成实习报告。

(2)实习内容

对1月份欧亚地区(200-700N；400-1400E)500hPa高度场进行自然正交展开 (EOF)，输出EOF的前三个特征向量，以及主要参数指标；绘制环流型图。

(3)步骤

1)熟悉资料

资料为NCEP/NCAR 61年(1948-2008年)1～12月的500hPa月平均高度场资料，资料范围为(900S-900N；00E-3600E)，网格距为2.50*2.50，纬向格点数为144，经向格点数为73。资料为GRD格式，按照从南到北、自西向东的顺序排列，每月为一个记录，按年逐月排放，注意读取方式以及记录长度。

本次实习应用NCEP/NCAR (200-700N；400-1400E)1948-2008年1月500百帕的平均高度场资料。

2)方法(自然正交函数分解EOF)

具体方法介绍参见第二章。应用步骤主要包括：

资料预处理(距平或标准化处理)

计算协方差矩阵

用Jacobi方法计算协方差矩阵的特征值与特征向量

将特征值从大到小排列

计算特征向量的时间系数

计算每个特征向量的方差贡献

结果输出

3)编写程序

FORTRAN程序参见附录EOF.for

编写GRADS数据描述文件和执行文件

4)结果输出

通过GRADS绘制图形，绘制EOF的前三个特征向量场。并通过表格形式给出前三个特征向量场各自的方差贡献和总的方差贡献。

5)完成实习报告

①说明所用资料方法；

②绘制前三个特征向量场；给出前三个特征向量场各自的方差贡献和总的方差贡献；

③用所学的知识分析图表资料，讨论它们各自