资源描述：
《短期气候预测基础实习二》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用文档实习二：大气环流分型一、实习目的及要求本实习的目的是掌握大气环流分型的基本方法－－EOF（经验正交函数分解）；要求熟悉EOF方法和程序的应用：用相关GRADS气象绘图系统，编写数据描述文件以及GRADS执行程序，用图形输出空间向量场；并能正确分析结果数据，完成实习报告。二、实习内容对1948-2008年1月欧亚（20-700N；40o-140oE）500hPa位势高度场的标准化序列进行自然正交展开，绘图主要的异常环流型环流型，输出讨论EOF分析主要参数标。三、实习步骤1、实习资料NCEP/NCAR1948－2008年1月的500百帕月平
2、均高度场资料，资料范围为（900S-900N,00-3600E)网格距为2.50×2.50,纬向格点数为144，经向格点数为41资料为GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。2、实习方法EOF功能：从一个气象场多次观测资料中识别出典型空间场型及其时间演变规律。EOF展开就是将气象变量场分解为空间函数(V)和时间函数(T)两部分的乘积之和：X=VT。应用步骤：①数据输入（主程序）②资料预处理（距平或标准化处理）callTRANSF(N,M,F,AVF,DF,KS)③计算协方差矩阵callFORMA(N,M,MNH,F
3、,A)④用Jacobi方法计算协方差矩阵的特征值与特征向量JCB(mnh,A,S,EPS）⑤将特征值从大到小排列callARRANG(KV,MNH,A,ER,S)⑥计算特征向量的时间系数callTCOEFF(KVT,KV,N,M,MNH,S,F,V,ER)⑦计算每个特征向量的方差贡献callOUTER(KV,ER)⑧结果输出（主程序）3、编写程序（1)绘制1948-2008年1月欧亚地区500hPa位势平均高度场图实习配置的GRADS数据描述文件：dsetd:yclimatesh2hgt500.grdundef-9.99E+33titl
4、emonthlymeanhgtfromtheNCEPReanalysis标准文案实用文档xdef144linear0.0002.500ydef73linear-90.0002.500zdef1levels500tdef732linearjan19481movars1hgt1-999monthlymeanhgtendvars;gs文件：'reinit''opend:yclimatesh2hgt500.ctl''enableprintd:yclimatesh2have1.gmf''setlev500''setlat-9090''se
5、tlon0360''sett1''definehave=ave(hgt,t+0,t=732,12)''dhave''print''disableprint';（2）编写关于EOF数据处理的FORTRAN程序：doit=1,61doi=1,41doj=1,21h(i,j,it)=hh(i+16,j+44,1,it)enddo;enddo;enddodoit=1,61doj=1,21doi=1,41f((j-1)*41+i,it)=h(i,j,it)enddo;enddoenddo（3）编写GRADS数据描述文件：dsetd:yclimate
6、sh2egvt.grdundef-9.99E+33titlemonthlymeanhgtfromtheNCEPReanalysisxdef41linear40.0002.500ydef21linear20.0002.500标准文案实用文档zdef1levels500tdef8linearjan19481movars1egvt199egvtendvars;（4）编写GRADS的执行（gs）文件：'reinit''opend:yclimatesh2egvt.ctl''enableprintd:yclimatesh2egvt.gmf
7、'i=1while(i<=3)'sett'i'''degvt''drawtitleNO.-'i'ofEOF''print''c'i=i+1endwhile'disableprint';4、结果输出1948-2008年1月500hPa高度平均图（单位：gpm）分析：上图表示的是1948-2008年标准文案实用文档1月份500hPa的高度场，从图中可以看出北半球冬季极地地区为极地低压，而中高纬地区呈三槽三脊的形势，三槽分别位于：亚洲东部、北美东部和欧洲东部，其中欧洲东部的槽很弱；三脊分别位于：美国阿拉斯加地区、西欧沿岸和青藏高原东部，可以看出脊的强
8、度比槽弱得多。在低纬地区，可以看到除了北美和东亚大槽向南伸到较低纬度外，在地中海、东太平洋都有比较明显的槽。而副热带高压强度小，且中心都位于海上。第一
