20131344075王彩虹实验四.docVIP

南京信息工程大学 实验（实习）报告 实验（实习）名称 蒙古高压与我国气温异常的关系分析 实验日期 2016.5.27 指导教师 马利 系 计软院 专业 软件工程 年级 13 班次 3 姓名 王彩虹 学号 20131344075 一、实验目的 （1）掌握气象要素相互影响研究的基本方法。 （2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。 （3）了解站点资料的数据结构。 二、实验内容 1 问题描述 利用1951-2010年1月蒙古高压标准化强度指数（data/h-p.dat）和1951-2010年1月中国160站气温资料(data/t1601.dat)，根据3.5节计算1月蒙古高压强度与中国160站气温的相关关系。要求以“*.dat”和“*.grd”两种格式保存1951-2010年1月蒙古高压强度与我国气温的同期相关数据。 2 问题分析 已知：1951-2010年1月蒙古高压标准化强度指数序列和1951-2010年中国160站逐月气温站点资料。计算：1951-2010年1月蒙古高压强度与中国160站气温同期相关系数。 三、实验要求 （1）分析问题，理解所给出的方法与技术，理清问题思路，分解任务，设计算法。 （2）输出数据采用有格式输入输出，使输出数据规范、醒目、简洁。 （3）用注释的方法指出程序中函数调用的起始和结束位置，并为函数进行注释。 （4）巩固顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。 （5）学习并掌握同期相关系数的求解方法，利用函数求1951-2010年1月蒙古高压强度与我国气温同期相关的160个相关系数，并以*.dat和*.grd两种格式保存。 四、实验步骤 蒙古高压环流指数的气候和异常值计算 1.分析问题，理清算法和程序，设计程序流程图并编写程序。 2.启动软件开发环境Microsoft Developer Studio。 3.在D盘上创建新工作区shixi04。 4.在工作区shixi3内创建新项目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程序。 蒙古高压环流指数距平时间序列图绘制 1.为蒙古高压环流指数距平数据文件“mh2.grd”书写数据描述文件“mh2.ctl”，在此文件中定义四个变量p、s、lo、la。 2.编写“mh2.gs”可执行文件，利用GrADS基本操作命令和line、bar两种绘图类型的绘图要素设置，以不同颜色和线形显示蒙古高压强度和面积时间序列图（曲线），分别以不同颜色显示蒙古高压经度和纬度时间序列图（柱状）。 3.将蒙古高压强度和面积时间序列图（曲线）保存到“mhline.gmf”，将蒙古高压经度和纬度时间序列图（柱状）分别保存到“mhlonbar.gmf”和“mhlatbar.gmf”中。 4.启动GrADS，调试、执行“mh2.gs”。 5.分析蒙古高压气候及其异常特征。 实习关键技术及方法 的n对观测资料和，则样本的相关系数可这样计算 其中x 表示压强 ，y表示 气温 ，r 表示相关系数 。 程序设计 以下程序用于蒙古高压环流指数气候及异常值计算。（Fortran程序编写） PROGRAM EX034 IMPLICIT none integer,parameter:: ny=60 real p(ny),s(ny),lon(ny),lat(ny),pa(ny),sa(ny),lona(ny),lata(ny),pav,sav,lonav,latav,pd,sd,lond,latd integer i,j,k open(1,file=D:\shixi4\p.dat,form=formatted) open(2,file=D:\shixi4\s.dat,form=formatted) open(3,file=D:\shixi