GrADS实习报告6

《GrADS实习报告6》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、南京信息工程大学实验（实习）报告实验（实习）名称蒙古高压与中国气温关系分析日期2016.10得分指导教师孙晓娟系大气科学班次长望实验班姓名夏新露一、实习目的（1）掌握气象要素互相影响研究的基本方法。（2）进一步巩固FORTRAN顺序结构和循环结构的程序设计方法和变量、函数、数组的使用方法，重点掌握外部文件对数据进行输入输出的方法。（3）掌握FORTRAN中WLF函数库的使用方法。（4）了解站点资料的数据结构。（5）掌握站点资料数据描述文件和站点映射文件的建立方法。（6）掌握站点资料生成个点资料的方法。（7）学习并掌握GRADS函数的使用方法，加强对oacre

2、s（）、maskout（）、smth9（）等函数的理解。（8）掌握GRADS变量的定义和使用方法。（9）巩固GRADS数据处理流程、绘图要素设置、基础绘图命令的使用方法，进一步掌握描述语言的应用方法。二、实习内容l问题描述利用1951-2010年1月蒙古高压标准化强度指数（p.dat）和1951-2010遍1月中国160站气温资料（t160.dat），根据6.5节相关系数的计算方法计算1951-2010年1月蒙古高压与中国气温的相互影响关系。要求以“*.dat”和“*.grd”两种格式保存1951-2010年1月蒙古高压强度与中国气温同期相关数据；用GRAD

3、S绘制1951-2010年1月蒙古高压强度与中国气温同期相关等值线图，写上标题，画出黄河、长江，并用阴影标记出显著相关区，并分析之。l问题分析已知: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月蒙古高压强度

4、指数“p.dat”，根据同期相关系数的计算方法，编写计算相关系数的函数计算1951-2010年1月蒙古高压强度与中国160站气温同期相关系数，并保存为“ptcorl.dat”和“ptcor1.grd”文件。根据教材中站点资料数据描述文件和站点映射文件的建立方法，结合rain.map建立站点的资料“ptcor1.grd”的数据描述文件和站点映射文件；生成中国160站的格点文件“grid.grd”，注意“grid.grd”的数据描述文件中的时间描述说明一定要与“ptcor1.grd”的数据描述文件相一致。用GRADS编写.gs文件，完成将二进制站点资料“ptco

5、r1.grd”插值到格点文件“grid.grd”的网格点上，并按照要求显示和保存图形。一、实习步骤（1）分析问题，绘制程序流程图，设计算法，编写程序。（2）启动软件开发环境MicrosoftDeveloperStudio。（3）在E盘上创建新工作区shixi6（4）在工作区shixi6内创建新项目shixi06（5）在项目shixi06内创建计算1951-2010年1月蒙古高压强度与中国气温同期相关的源程序文件“ptcor1.f90”，编辑输入源程序文本。（6）拷贝WFL到C:，配置FORTRAN使用WFL函数库的环境。（7）编译、构建、运行、调试FORT

6、RAN程序，生成“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”，并书写

7、该文件的数据描述文件“grid.ctl”，注意“grid.ctl”中的时间说明要与“mh-t-gr.ctl”中时间一致。（11）在GRADS运行环境中运行“!stanmap”命令，在提示下输入“mh-t-gr.ctl”的路径和文件名，生成“mh-t-gr.grd”站点映射文件mh-t-gr.map。（12）利用插值函数oacres（）、标记函数maskout（）、平滑函数smth9（）等编写“ptcor1.gs”，完成将二进制站点资料文件“ptcor1.grd”插值到格点文件“grid.grd”的网格点上，按要求绘制1951-2010年1月蒙古高压与我国气温

8、同期相关图。（13）启动GRADS运行、调试“ptc