报告彩虹实验报告.docx

报告彩虹实验报告.docx

《报告彩虹实验报告.docx》由会员分享，可在线阅读，更多相关《报告彩虹实验报告.docx（14页珍藏版）》请在冰豆网上搜索。

报告彩虹实验报告

【关键字】报告

彩虹实验报告

篇一：

XX1344075王彩虹实验四

　　南京信息工程大学实验（实习）报告实验（实习）名称蒙古高压与我国气温异常的关系分析实验日期指导教师马利系计软院专业软件工程年级13班次3姓名王彩虹学号XX1344075

　　一、实验目的

（1）掌握气象要素相互影响研究的基本方法。

（2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。

　　（3）了解站点资料的数据结构。

　　二、实验内容

　　1问题描述

　　利用1951-XX年1月蒙古高压标准化强度指数（data/h-p.dat）和1951-XX年1月中国160站气温资料（data/t1601.dat），根据3.5节计算1月蒙古高压强度与中国160站气温的相关关系。

要求以“*.dat”和“*.grd”两种格式保存1951-XX年1月蒙古高压强度与我国气温的同期相关数据。

　　2问题分析

　　已知：

1951-XX年1月蒙古高压标准化强度指数序列和1951-XX年中国160站逐月气温站点资料。

计算：

1951-XX年1月蒙古高压强度与中国160站气温同期相关系数。

　　三、实验要求

（1）分析问题，理解所给出的方法与技术，理清问题思路，分解任务，设计算法。

（2）输出数据采用有格式输入输出，使输出数据规范、醒目、简洁。

　　（3）用注释的方法指出程序中函数调用的起始和结束位置，并为函数进行注释。

　　（4）巩固顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。

　　（5）学习并掌握同期相关系数的求解方法，利用函数求1951-XX年1月蒙古高压强度与我国气温同期相关的160个相关系数，并以*.dat和*.grd两种格式保存。

　　四、实验步骤

　　蒙古高压环流指数的气候和异常值计算

　　1.分析问题，理清算法和程序，设计程序流程图并编写程序。

　　2.启动软件开发环境MicrosoftDeveloperStudio。

　　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.分析蒙古高压气候及其异常特征。

　　五、实习关键技术及方法

　　x、y的n对观测资料x1,x2,...,xn和y1,y2,...,yn，则样本的相关系数rxy可这样计算

　　rxy?

　　1n（xt?

）（yt?

）?

其中x表示压强，y表示气温，r表示相关系数。

　　程序设计

　　以下程序用于蒙古高压环流指数气候及异常值计算。

（Fortran程序编写）

　　PROGRAMEX034

　　IMPLICITnone

　　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

　　integeri,j,k

　　open（1,file='D:

\shixi4\p.dat',form='formatted'）open（2,file='D:

\shixi4\s.dat',form='formatted'）

　　open（3,file='D:

\shixi4\lon.dat',form='formatted'）open（4,file='D:

\shixi4\lat.dat',form='formatted'）doi=1,ny

　　read（1,*）,p（i）

　　read（2,*）,s（i）

　　read（3,*）,lon（i）

　　read（4,*）,lat（i）

　　enddo

　　close

（1）

　　close

（2）

　　close（3）

　　close（4）

　　callcha（ny,p,pa,pav,pd）

　　callcha（ny,s,sa,sav,sd）

　　callcha（ny,lon,lona,lonav,lond）

　　callcha（ny,lat,lata,latav,latd）

　　open（5,file='D:

\shixi4\mh1.dat',form='formatted'）write（5,*）,pav,pd

　　write（5,*）,sav,sd

　　write（5,*）,lonav,lond

　　write（5,*）,latav,latd

　　close（5）

　　open（6,file='D:

\shixi4\mh1.grd',form='binary'）

　　write（6）,pav,pd

　　write（6）,sav,sd

　　write（6）,lonav,lond

　　write（6）,latav,latd

　　close（6）

　　open（7,file='D:

\shixi4\mh2.dat',form='formatted'）write（7,*）,pa（60）

　　write（7,*）,sa（60）

　　write（7,*）,lona（60）

　　write（7,*）,lata（60）

　　close（7）

　　open（8,file='D:

\shixi4\mh2.grd',form='binary'）Do

　　write（8）,pa（i）write（8）,sa（i）

　　write（8）,lona（i）write（8）,lata（i）

　　enddo

　　close（8）end

　　subroutinecha（ny,x,xa,xav,xd）

　　integer:

:

ny

　　integeri

　　real:

:

x（ny）,xa（ny）,xav,xd,sum

　　sum=0

　　doi=1,ny

　　sum=sum+x（i）

　　Enddo

　　xav=sum/ny

　　xd=0

　　doi=1,ny

　　xa（i）=x（i）-xav

　　xd=xa（i）*xa（i）+xd

　　enddo

　　xd=sqrt（xd/ny）

　　Returnend

　　mh2.ctl”

　　dsetD:

\shixi4\mh2.grd

　　undef-9.99E+33

　　titleTheMongolianhighpressure

　　xdef1levels1

　　ydef1levels1

　　zdef1levels1

　　tdef60linearJAN19511yrGrADS程序编写）以下程序用于蒙古高压环流指数距平值绘图。

（“

　　vars4

　　p00qiangdu

　　s00mianji

　　lo00jingdu

　　la00weidu

　　endvars

　　“mh2.gs”

　　'reinit''openD:

\shixi4\mh2.ctl'

　　'enableprintD:

\shixi4\mhline.gmf''setlat1'

　　'setlon1'

　　'setlev1'

　　'sett160'

　　'setgxoutline'

　　'setccolor1'

　　'setcstyle1'

　　'setcthick6'

　　'setcmark1''dp'

　　'setccolor2'

　　'setcstyle2'

　　'setcthick6'

　　'setcmark2''ds'

　　'print'

　　'disableprint''c'

　　'enableprintD:

\shixi4\mhlonbar.gmf'

　　'setgxoutbar'

　　'setbarbase0'

　　'setbargap30'

　　'setccolor3'

　　'print'

　　'disableprint''c'

　　'enableprintD:

\shixi4\mhlatbar.gmf''setgxoutbar'

　　'setbarbase0'

　　'setbargap30'

篇二：

XX1344075王彩虹实验五

　　南京信息工程大学实验（实习）报告实验（实习）名称与Grads相关的应用

（2）实验（实习）日期马利系计软院专业软件工程年级13班次3姓名王彩虹学号XX1344075

　　一、实验目的

（1）掌握气象要素相互影响研究的基本方法。

（2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。

　　（3）了解站点资料的数据结构。

　　二、实验内容

　　问题描述

　　利用1948-XX年NCAR/NCEP月平均气温和降水再分析资料（nc格式文件），分析1月气温和降水气候特征。

要求利用Fortran提取1948-XX年1月数据并进行计算，以“*.grd”格式保存1948-XX年1月气温、降水气候场数据，用GrADS生成1948-XX年1月气温、降水气候二维等值线和二维填色图两种格式叠加的图像，书写标题，并在高温、低温中心标记“H”和“L”。

保存并分析之。

　　1）问题分析

　　已知：

1948-XX年NCAR/NCEP月降水和气温再分析资料。

　　绘制：

1948-XX年1月气温、降水气候图。

　　通过分析，首先用Grads将从NCAR/NCEP月气温和降水再分析资料（nc格式

　　文件）中提取气温和降水1月逐年数据，并保存为为grd格式文件，再用Fortran程序利用grd格式资料按照公式求得气温和降水的气候值，用GrADS绘制二者的气候图，分析1948-XX年1月蒙古高压气候特征。

　　三、实验步骤

　　分析问题，理清算法和程序，设计和编写程序。

　　①书写“air1948-XX-1.gs”和“pre1948-XX-1.gs”从NCAR/NCEP再分析资料（nc格式），利用fwrite提取1月逐年数