20131344075王彩虹实验四

《20131344075王彩虹实验四》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、南京信息工程大学实验（实习）报告实验（实习）名称蒙古高压与我国气温异常的关系分析实验日期2016.5.27指导教师马利系计软院专业软件工程年级13班次3姓名王彩虹学号20131344075一、实验目的（1）掌握气象要素相互影响研究的基本方法。（2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。（3）了解站点资料的数据结构。二、实验内容1问题描述利用1951-2010年1月蒙古高压标准化强度指数（data/h-p.dat）和1951-2010年1月中国160站气温资料(data/t1601

2、.dat)，根据3.5节计算1月蒙古高压强度与中国160站气温的相关关系。要求以“*.dat”和“*.grd”两种格式保存1951-2010年1月蒙古高压强度与我国气温的同期相关数据。2问题分析已知：1951-2010年1月蒙古高压标准化强度指数序列和1951-2010年中国160站逐月气温站点资料。计算：1951-2010年1月蒙古高压强度与中国160站气温同期相关系数。三、实验要求（1）分析问题，理解所给出的方法与技术，理清问题思路，分解任务，设计算法。（2）输出数据采用有格式输入输出，使输出数据规范、醒目、简洁。（3）用注释的方法指出程序中函数调用的起始和结束

3、位置，并为函数进行注释。（4）巩固顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。（5）学习并掌握同期相关系数的求解方法，利用函数求1951-2010年1月蒙古高压强度与我国气温同期相关的160个相关系数，并以*.dat和*.grd两种格式保存。四、实验步骤蒙古高压环流指数的气候和异常值计算1.分析问题，理清算法和程序，设计程序流程图并编写程序。2.启动软件开发环境MicrosoftDeveloperStudio。3.在D盘上创建新工作区shixi04。4.在工作区shixi3内创建新项目shixi04。5.在项

4、目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.为蒙

5、古高压环流指数距平数据文件“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

6、”。5.分析蒙古高压气候及其异常特征。五、实习关键技术及方法的n对观测资料和，则样本的相关系数可这样计算其中x表示压强，y表示气温，r表示相关系数。程序设计以下程序用于蒙古高压环流指数气候及异常值计算。（Fortran程序编写）PROGRAMEX034IMPLICITnoneinteger,parameter::ny=60realp(ny),s(ny),lon(ny),lat(ny),pa(ny),sa(ny),lona(ny),lata(ny),pav,sav,lonav,latav,pd,sd,lond,latdintegeri,j,kopen(1,file='

7、D:shixi4p.dat',form='formatted')open(2,file='D:shixi4s.dat',form='formatted')open(3,file='D:shixi4lon.dat',form='formatted')open(4,file='D:shixi4lat.dat',form='formatted')doi=1,nyread(1,*),p(i)read(2,*),s(i)read(3,*),lon(i)read(4,*),lat(i)enddoclose(1)close(2)close(3)close(4)