第四章蒙古高压特征分析.docVIP

第四章:蒙古高压特征分析 4、1 实验实习目得 (1)掌握气象要素气候及异常特征得基本研究方法。 (2)进一步巩固FORTRAN顺序结构与循环结构得程序设计方法与变量、函数、数组得使用方法,重点掌握外部文件对数据进行输入输出操作得方法。 (3)了解并掌握GrADS得使用流程。 (4)学习并掌握数据描述文件(*、ctl文件)得书写方法。 (5)学习书写简单得“*、gs”文件。 (6)掌握GrADS中open、reinit、d、c、quit、q、run等基本命令得使用方法。 (7)掌握GrADS维数环境设置方法。 (8)掌握GrADS绘制“line”、“bar”图形类型得方法,绘制1951-2010年蒙古高压强度异常得时间序列图。 (9)掌握保存GrADS绘制图形得一种方法。 4、2 实验实习内容 4、2、1 问题描述 已知1951-2010年一月蒙古高压强度、面积、位置指数序列,计算蒙古高压各指数得气候值、变率与距平,绘制蒙古高压得强度、面积、位置指数序列得时间序列图,分析冬季蒙古高压得异常规律。 4、2、2 问题分析 已知:1951-2010年一月蒙古高压强度、面积、经度、纬度指数序列资料p、dat、s、dat、lon、dat、lat、dat。 计算:蒙古高压各指数得气候、变率与距平值。 绘制:1951-2010年一月蒙古高压强度、面积、位置指数序列图。 通过分析,根据公式(4、1)、公式(4、2),求得蒙古高压一月环流指数气候及异常值。根据GrADS中line与bar两种图形格式绘制方法,绘制蒙古高压环流指数距平得时间序列图。 4、3 实验实习步骤 4、3、1 蒙古高压环流指数得气候与异常值计算 (1)分析问题,理清算法与程序,设计程序流程图并编写程序。 (2)启动软件开发环境Microsoft Developer Studio。 (3)在I盘上创建新工作区 shixi04。 (4)在工作区shixi04内创建新项目shixi04。 (5)在项目shixi04内创建源程序文件“mh、f90”,编辑输入源程序文本。 (6)在源程序文本中打开数据文件“p、dat”、“s、dat”、“lon、dat”、“lat、dat”,并将其值读入到相应得数组中。 (7)编写计算均值、变率与距平得子程序。 (8)调用子程序分别计算强度、面积、经度、纬度环流指数得均值、变率与距平。 (9)将蒙古高压环流指数得均值与变率写入到“mh1、dat”与“mh1、grd”两个文件中。将蒙古高压环流指数得距平值写入到“mh2、dat”与“mh2、grd”两个文件中。 (10)编译、构建、运行、调试FORTRAN程序。 4、3、2 蒙古高压环流指数距平时间序列图绘制 GrADS就是对数据进行分析、处理与显示得软件,所以在使用该软件进行绘图时,必须有相应得数据文件,并且这些数据必须满足GrADS所要求得数据格式,否则就要进行数据处理;另外,GrADS并不就是直接使用这些数据文件进行操作,而就是通过一个对应得“数据描述文件”间接使用数据文件;对数据得处理与显示就是通过GrADS系统提供得命令完成得;操作完成后,可以在图形输出窗口观察所绘图形,并将图形保存与输出。具体使用流程如1、3、3、3所述,此处不再赘述。 按1、3、3、3 所述得使用流程,完成蒙古高压环流指数距平时间序列图绘制,需要按以下步骤进行: (1)为蒙古高压环流指数距平数据文件“mh2、grd”书写数据描述文件“mh2、ctl”,在此文件中定义四个变量pa,sa,lona,lata、 (2)编写“mh2、gs”可执行文件,利用GrADS基本操作命令与line、bar两种绘图类型得绘图要素设置,以不同颜色与线性显示蒙古高压强度与面积距平时间序列图(曲线),分别以不同颜色显示蒙古高压经度与纬度距平时间序列图(柱状)。 (3)将蒙古高压强度与面积距平时间序列图(曲线)保存到“mhline、gmf”,将蒙古高压经度与纬度距平时间序列图(柱状)分别保存到“mhlonbar、gmf”与“mhlatbar、gmf”中。 (4)启动GrADS,调试、执行“mh2、gs”。 (5)分析蒙古高压气候及其异常特征。 4、4 实验实习关键技术及方法 某数据资料时间序列得距平为数据资料 与其平均值 之差 (4、1) 某数据资料得变率 为其均方差,反应变量围绕平均值得平均变化程度,其计算公式为 (4、2) 4、5 实验实习程序编写 4、5、1 FORTRAN程序编写 以下程序用于蒙古高压环流指数气候及异常值计算。 PROGRAM mh IMPLICIT none inte