GrADS教学中的几点体会.docVIP

第31卷　第2期气象教育与科技2008年　总第83期 GrAD S 教学中的几点体会 李丽平 (南京信息工程大学大气科学院, 江苏南京　210044 摘要:总结了Gr ADS 软件绘图基本步骤, 提炼出软件的关键命令, 给出实现绘制几种常用重要图形(包括写数据 的控制文件, , 提出了如何应对入门困惑的办法, 旨在帮助Gr ADS 初学者尽快入门。 关键词:Gr ADS; 绘图步骤; 关键命令; 目前有很多流行的绘图软件, 如M I 、NCL 、MAT LAB 等, 其中Gr ADS , 、物理量计算、数据格式转换。和EC WMF 再分析资料 一般为net CDF 格式, 。目前, 该软件已是气象、海洋专业科研工作者和高校学生做宝贵经验[223][1], 一些学者就Gr ADS 处理某些具体问题提出了许多, 成为本科生和函授研究生必须选修的一门课程。 Gr ADS 是一门应用性很强的课程, 内容包含了若干命令, 在学时有限的情况下, 如何让学生在短时间内掌握最核心的精华内容, 做到尽快入门, 已是教学中值得研究的一个重要问题。笔者根据多年教授Gr ADS 课程的教学经验, 从宏观的角度提炼Gr ADS 的几个精华知识点, 试图为学生和其他初学Gr ADS 的科研工作者开辟一条清晰的学习思路。 1　Gr ADS 绘图需要的3类文件 1. 1　Gr ADS 作图的一般步骤 Gr ADS 作图步骤:打开数据描述文件、输入Gr ADS 命令(如维数环境的设置、图形类型的设置、图形要素的设置 、保存图形、退出Gr ADS 。 绘图之前需准备的类型文件:即数据文件(. grd 、数据描述文件(. ctl 、控制文件(. gs 。 1. 2　数据文件 1. 2. 1　数据格式 Gr ADS 能够处理的数据文件常见格式有net CDF 、HDF 2S DS (卫星资料 、GR I B 、B inary (用户自己生成 、站点资料等。美国国家环境预测中心(NCEP /NCAR 和欧洲中期天气预报中心(EC WMF 资料大多为net CDF 格式, 这种格式的数据资料自带数据描述信息, 在Gr ADS 中用sdfopen 命令打开。B inary 是Gr ADS 最基本的数据格式, 常用文件后缀为. dat 、. grd 、. bin, 需要　　作者简介:李丽平(1971— , 女, 山西寿阳人, 副教授, 博士, 研究方向为短期气候预测, lili p ing@nuist. edu . cn . ?1? 。存取方式有直接(ac 2和数据描述(. ctl 一起使用, 是一种无格式文件(f or m =′unfor matted ′ 和间接(access =′ 两种, 间接存取要在. ctl 文件中给以特别说明。cess =′driect ′sequential ′ GR I B 格式资料扩展名一般为. grb, 需和. ctl 文件和. idx 文件一起使用。 1. 2. 2　两个重要概念 关于Gr ADS 数据有两个重要概念, 即数据不同维数之间的关系、同一维数据的存放顺序, 正确理解这两个概念是Gr ADS 绘图或计算准确的保证。 (1 数据不同维之间的关系 Gr ADS 将数据看作是一个五维数据集(图1 , 这是Gr ADS 中的一个非常重要的概念。图1中箭头表示编制程序时五个维数的循环嵌套关系, 越向右, 表示该维越在外循环。无论通过编制Fortran 程序为Gr ADS 绘图和计算准备资料、还是Gr ADS 本身生成数据, 都要遵循这一编程思路, 因为这是Gr ADS 所要求的数据维数关系。 (2 同一维数据的存放顺序 Gr ADS 要求一个x 、y 数据场构成一个记录, 其顺序必须是从西向东, , 从下到上, 即实际大数组以二维数据片存放。如果资料是从北到南、中的op ti ons 选项作特殊说明。一般而言, 描述文件的内容即可知 ; 。 图1　数据资料不同维数的循环嵌套关系 以下给出将ASC II 资料写为无格式直接存取文件(用Fortran 语言实现 的一段程序。重点介绍无格式直接存取文件的生成。 c 打开目标文件, recl 为一个记录的长度, rec 为记录号. open (12, file =′d:\ll p \sy . grd ′, for m =′unf or matted ′, access =′direct ′, recl =nx 3ny 34 cccccccccc 写数据资料ccccccccccc irec =0 do 100it =1, nt do 50iz =1, nz irec =irec +1 write (12, rec =irec ((u (i, j , iz, it , i =1, nx