GrADS实习报告6.docVIP

南京信息工程大学 实验（实习）报告 实验（实习）名称蒙古高压与中国气温关系分析日期 2016.10 得分 指导教师孙晓娟系大气科学班次长望实验班姓名夏新露 一、 实习目的 （1） 掌握气象要素互相影响研究的基本方法。 （2） 进一步巩固 FORTRAN 顺序结构和循环结构的程序设计方法和变量、 函 数、数组的使用方法，重点掌握外部文件对数据进行输入输出的方法。 （3） 掌握 FORTRAN 中 WLF 函数库的使用方法。 （4） 了解站点资料的数据结构。 （5） 掌握站点资料数据描述文件和站点映射文件的建立方法。 （6） 掌握站点资料生成个点资料的方法。 （7） 学习并掌握 GRADS 函数的使用方法，加强对 oacres（）、maskout（）、 smth9（）等函数的理解。 （8） 掌握 GRADS 变量的定义和使用方法。 （9） 巩固 GRADS 数据处理流程、 绘图要素设置、 基础绘图命令的使用方法， 进一步掌握描述语言的应用方法。 二、 实习内容 问题描述 利用 1951-2010年 1 月蒙古高压标准化强度指数（ p.dat）和 1951-2010 遍 1 月中国 160站气温资料（t160.dat），根 据 6.5节相关系数的计算方法计算 1951-2010 年 1 月蒙古高压与中国气温的相互影响关系。要求以“ *.dat”和“*.grd”两种格 式保存 1951-2010年 1 月蒙古高压强度与中国气温同期相关数据；用 GRADS 绘 制 1951-2010年 1 月蒙古高压强度与中国气温同期相关等值线图，写上标题，画 出黄河、长江，并用阴影标记出显著相关区，并分析之。 问题分析 已知:1951-2010 年 1 月蒙古高压标准化强度指数序列和 1951-2010 年中国 160 站 1 月气温站点资料。 计算：1951-2010年 1月蒙古高压强度与中国 160站气温同期相关系数。 绘制：1951-2010年 1 月蒙古高压与中国 160站气温同期相关图， 写上标题， 并用阴影标记显著相关区。 通过分析，首先用 FORTRAN 打开从 1951-2010年中国 160 站1 月气温数据 “t160.dat”和 1951-2010年1 月蒙古高压强度指数“ p.dat”，根据同期相关系数 的计算方法， 编写计算相关系数的函数计算 1951-2010年 1 月蒙古高压强度与中 国 160 站气温同期相关系数，并保存为“ ptcorl.dat”和“ptcor1.grd”文件。 根据教材中站点资料数据描述文件和站点映射文件的建立方法，结合 rain.map建立站点的资料“ ptcor1.grd”的数据描述文件和站点映射文件；生成中 国 160 站的格点文件“ grid.grd”，注意“grid.grd”的数据描述文件中的时间描述 说明一定要与“ ptcor1.grd”的数据描述文件相一致。 用 GRADS编写.gs 文件，完成将二进制站点资料“ ptcor1.grd”插值到格点 文件“ grid.grd”的网格点上，并按照要求显示和保存图形。 三、实习步骤 （1） 分析问题，绘制程序流程图，设计算法，编写程序。 （2） 启动软件开发环境 Microsoft Developer Studio。 （3） 在 E盘上创建新工作区 shixi6 （4） 在工作区 shixi6 内创建新项目 shixi06 （5） 在项目 shixi06 内创建计算 1951-2010年 1 月蒙古高压强度与中国气 温同期相关的源程序文件“ ptcor1.f90”，编辑输入源程序文本。 （6） 拷贝WFL 到 ，配置 FORTRAN 使用 WFL 函数库的环境。 （7）编译、构建、运行、调试FORTRAN 程序，生成“ mh-t-1.dat”和 “mh-t-1.grd”。 （8） 用 FORTRAN读取存放全国 160 个站台经纬度记录的文本文件 “china.dat”和“ptcor01.dat”，编写 FORTRAN 程序“cor-stntogrd.f90”， 调试运行， 将站点数据文件 “mh-t-1.dat”转换成按 GRADS 要求排列 格点数据“ mh-t-gr.grd”。 （9）为“ mh-t-gr.grd”编写数据描述文件“ mh-t-gr.ctl” （10）将站点资料通过插值函数插值到格点文件上，首先建立“ grid,f90” 程序，生成 160站的格点文件“ grid.grd”，并书写该文件的数据描述 文件“ grid.ctl”，注意“ grid.ctl”中的时间说明要与“ mh-t-gr.ctl”中 时间一致。 （11）在 GRADS 运行环境中运行“!stanmap”命令，在提示下输入“mh-t-gr.ctl” 的路径