20131344075王彩虹实验五.docxVIP

南京信息工程大学 实验（实习）报告实验（实习）名称与Grads相关的应用(2)实验（实习）日期 2016.6.3得分指导教师 马利系 计软院专业 软件工程年级 13班次 3 姓名 王彩虹 学号 20131344075一、实验目的（1）掌握气象要素相互影响研究的基本方法。（2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。（3）了解站点资料的数据结构。二、实验内容问题描述利用1948-2010年NCAR/NCEP月平均气温和降水再分析资料（nc格式文件），分析1月气温和降水气候特征。要求利用Fortran提取1948-2010年1月数据并进行计算，以“*.grd”格式保存1948-2010年1月气温、降水气候场数据，用GrADS生成1948-2010年1月气温、降水气候二维等值线和二维填色图两种格式叠加的图像，书写标题，并在高温、低温中心标记“H”和“L”。保存并分析之。1）问题分析已知：1948-2010年NCAR/NCEP月降水和气温再分析资料。绘制：1948-2010年1月气温、降水气候图。通过分析，首先用Grads将从NCAR/NCEP月气温和降水再分析资料（nc格式文件）中提取气温和降水1月逐年数据，并保存为为grd格式文件，再用Fortran程序利用grd格式资料按照公式求得气温和降水的气候值，用GrADS绘制二者的气候图，分析1948-2010年1月蒙古高压气候特征。三、实验步骤分析问题，理清算法和程序，设计和编写程序。书写“air1948-2010-1.gs”和“pre1948-2010-1.gs”从NCAR/NCEP再分析资料（nc格式），利用fwrite提取1月逐年数据，生成气温、降水1月的二进制数据“air1.grd”、“pre1.grd”，编写1月气温、降水资料“air1.grd”、“pre11.grd”的数据描述文件“air1.ctl”、“pre1.ctl”。启动软件开发环境Microsoft Developer Studio。在D盘上创建新工作区shixi05。在工作区shixi2内创建新项目shixi05。在项目shixi21内创建源程序文件“shixi05.f90”，编辑输入源程序文本，打开1948-2010年1月气温、降水资料“air1.grd”、“pre1.grd”，根据求均值公式求得1948-2010年1月气温、降水气候场值。编译、构建、运行、调试Fortran程序，生成1948-2010年1月气温、降水气候场数据“tc1.grd”、“pc1.grd”。打开记事本，编写“tc1.grd”和“pc1.grd”的数据描述文件“tc1.ctl”和“pc1.ctl”。利用GrADS的基本绘图命令和contour、shaded两种图形类型的图形设置方法和图形要素设置方法以及基本绘图命令，利用printim和enable print两种输出图像方法及输出图像格式设置，编写“tc1.gs”和 “pc1.gs”。启动GrADS，运行、调试“tc1.gs”和 “pc1.gs”，保存气温、降水气候图“tc1.gmf”、“pc1. gmf”或者“tc1.gif”、“pc1. gif”。分析1948-2010年1月蒙古高压气候特征。程序设计提取NCAR/NCEP再分析资料中1月气温、降水二进制数据“air1948-2010-1.gs”reinitsdfopen D:\shixi5\air.mon.mean.ncset gxout fwriteset fwrite D:\shixi5\air1.grdset x 1 144set y 1 73i=1while(i=756)set t id airi=i+12endwhiledisable fwrite; “pre1948-2010-1.gs”j=1reinitsdfopen D:\shixi5\pr_wtr.eatm.mon.mean.ncset gxout fwriteset fwrite D:\shixi5\pre1.grdset x 1 144set y 1 73i=1while(i=756)set t id pr_wtri=i+12endwhiledisable fwrite 编写“air1.grd”、“pre1.grd”数据描述文件“air1.ctl”和“pre1.ctl”。“air1.ctl”dset D:\shixi5\air1.grdtitle air temperature of NCEP Reanalysi s in Jan xdef 144 linear 0 2.5ydef 73 linear -90 2.5zdef 1 linear 0 1tdef 63 li