短期气候预测基础实习二.docVIP

实习二：大气环流分型 一、实习目的及要求 本实习的目的是掌握大气环流分型的基本方法－－EOF（经验正交函数分解）；要求熟悉EOF方法和程序的应用：用相关GRADS气象绘图系统，编写数据描述文件以及GRADS执行程序，用图形输出空间向量场；并能正确分析结果数据，完成实习报告。 二、实习内容 对1948-2008年1月欧亚（20-700N；40o-140oE）500hPa位势高度场的标准化序列进行自然正交展开 ，绘图主要的异常环流型环流型，输出讨论EOF分析主要参数标。 三、实习步骤 1、实习资料 NCEP/NCAR 1948－2008年1月的500百帕月平均高度场资料，资料范围为（900S-900N,00-3600E) 网格距为2.50×2.50,纬向格点数为144，经向格点数为41 资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。 2、实习方法 EOF功能：从一个气象场多次观测资料中识别出典型空间场型及其时间演变规律。EOF展开就是将气象变量场分解为空间函数(V)和时间函数(T)两部分的乘积之和：X=VT。 应用步骤： 数据输入（主程序） 资料预处理（距平或标准化处理） call TRANSF(N,M,F,AVF,DF,KS) 计算协方差矩阵call FORMA(N,M,MNH,F,A) 用Jacobi方法计算协方差矩阵的特征值与特征向量 JCB(mnh,A,S,EPS） 将特征值从大到小排列call ARRANG(KV,MNH,A,ER,S) 计算特征向量的时间系数call TCOEFF(KVT,KV,N,M,MNH,S,F,V,ER) 计算每个特征向量的方差贡献call OUTER(KV,ER) 结果输出（主程序） 3、编写程序 （1)绘制1948-2008年1月欧亚地区500hPa位势平均高度场图 实习配置的GRADS数据描述文件： dset d:\nyclimate\sh2\hgt500.grd undef -9.99E+33 title monthly mean hgt from the NCEP Reanalysis xdef 144 linear 0.000 2.500 ydef 73 linear -90.000 2.500 zdef 1 levels 500 tdef 732 linear jan1948 1mo vars 1 hgt 1 -999 monthly mean hgt endvars ; gs文件： reinit open d:\nyclimate\sh2\hgt500.ctl enable print d:\nyclimate\sh2\have1.gmf set lev 500 set lat -90 90 set lon 0 360 set t 1 define have=ave(hgt,t+0,t=732,12) d have print disable print ; （2） 编写关于EOF数据处理的FORTRAN程序： do it=1,61 do i=1,41 do j=1,21 h(i,j,it)=hh(i+16,j+44,1,it) enddo;enddo;enddo do it=1,61 do j=1,21 do i=1,41 f((j-1)*41+i,it)=h(i,j,it) enddo;enddo enddo （3）编写GRADS数据描述文件： dset d:\nyclimate\sh2\egvt.grd undef -9.99E+33 title monthly mean hgt from the NCEP Reanalysis xdef 41 linear 40.000 2.500 ydef 21 linear 20.000 2.500 zdef 1 levels 500 tdef 8 linear jan1948 1mo vars 1 egvt 1 99 egvt endvars ; （4）编写GRADS的执行（gs）文件： reinit open d:\nyclimate\sh2\egvt.ctl enable print d:\nyclimate\sh2\egvt.gmf i=1 while(i=3) set t i d egvt draw title NO.-i of EOF print c i=i+1 endwhile disable print ; 4、结果输出 1948-2008年1月500hPa高度平均图（单位：gpm） 分析：上图表示的是1948-2008年1月份500hPa的高度场，从图中可以看出北半球冬季极地地区为极地低压，而中高纬地区呈三槽三脊的形势，三槽分别位于：亚