GrADS任意方向垂直剖面.docVIP

GrADS下任意方向垂直剖面的实现 缪明 余清平 （解放军理工大学气象学院） 廖德敏 （93966部队气象室） 引言 GrADS（Grid Analysis and Display System）是当前气象学中应用最为广泛的图形图像处理软件之一，它功能强大，使用方便而深受用户的欢迎。利用它可直接实现图像的动画显示，在三维数值模拟研究中，在固定模式高度的情况下，可较为方便地显示某一模式预报量在该层次上的水平分布；在固定经（纬）度的情况下，还可方便地显示某一模式预报量沿该经（纬）度的经（纬）向垂直剖面分布。但在实际问题的三维数值模拟研究中，用户关心的不仅仅是模式预报量沿某一经纬向的垂直剖面分布，更多的时候用户关心的是沿某一方向的垂直剖面分布，但在GrADS下如何实现这一功能，这可能是困扰一般的GrADS用户的问题。其实，GrADS除了强大的图形图象的分析处理能力外，还有较强的数据处理能力，我们可利用GrADS的交互性和数据处理的特点，利用GrADS脚本语言的编程技术，按GrADS的数据格式，通过提取所要分析的模式预报量沿某一方向上的值（ascII码形式），然后再将提取ascII码值转换成GrADS支持的二进制格式，并为转化后的数据编写相应的数据描述文件（ctl文件），在GrADS下打开它就可实现某一模式预报量沿该方向的垂直剖面分布的分析。 实现的思想和步骤 在介绍实现的思想和gs编程（GrADS Script programing）之前先说明GrADS下某些基本的常识和某些基本命令的含义和功能。 GrADS基本常识和基本命令的简介 （1）进入GrADS后，在“ga -”后每执行一个GrADS命令，在文本窗口都会返回1行或2行的的文本代码，返回的代码即为GrADS的默认变量“result”值，每执行一次GrADS命令，原来的“result”值将被刷新成刚执行的GrADS命令返回码。在gs编程中，可通过代码“say result ”看到默认变量“result”对应的返回文本代码。 （2）命令“q bpos”： 功能：确定鼠标在图形窗口中所点击处的位置（页面坐标） 返回码：执行完该命令后，有一行形如“position = x0 y0 1 0”的文本返回码，其中x0、y0即为用户在图形窗口中点击处对应的页面坐标。 注意：执行该命令后，GrADS将等待用户的响应，直到用户用鼠标点击图形窗口后才可继续执行其它的命令或返回到“ga -”状态下。 （3）在GrADS下成功打开一个.ctl描述文件对应的数据文件后，在某一水平显示环境下（即t、z一定的情况下）一旦执行了“d f”(f为.ctl文件中的某个变量名)，图形窗口任何一点就在页面坐标、经纬坐标和格点坐标三者间建立起一一对应的关系，这种对应关系只有在一定的缩放比例环境下才有意义，.ctl描述文件对应的数据文件后，执行命令“d f ”的一个目的就是要在页面坐标与经纬坐标间建立一定的缩放比例。执行“clear”命令或“reset” 命令后，页面坐标与经纬坐标间建立一定的缩放比例自动终止。 （4）命令“q xy2w x0 y0”： 功能：实现页面坐标（x0，y0）与对应的经纬坐标（lon0，lat0）之间的转换； 返回码：在页面坐标与经纬坐标间建立一定的缩放比例后，执行该命令，在文本窗口中将有一行形如“Lon = lon0 Lat = lat0”文本返回码，其中lon0、lat0为与页面坐标（x0，y0）对应的经纬度。 注意：只有在页面坐标与经纬坐标间建立起一定的缩放比例的前提下，才会有“Lon = lon0 Lat = lat0”的返回码，也才能实现页面坐标（x0，y0）与对应的经纬坐标（lon0，lat0）之间的转换。 （5）几个GrADS脚本编程语言的内部函数（intrinsic function） sublin (string, nth_line): 功能：从字符串“string”中提取第“nth_line”行子字符串。 注意：“nth_line”为正整数。 举例：假设string=“Demonstrate how to use the intrinsic function！ ” 则：sublin(string,2)=“the intrinsic function！” subwrd (string,nth_word)： 功能：从字符串“string”中提取第“nth_word”个“字(word)”。 注意：“nth_word”为正整数；在GrADS 下，用空格符来作为“字(word)”与“字(word)”间的分隔符，因此“字(word)”是指任何两相邻两空格间的非空隔字符或非空隔字符串。 举例：如：string的值同上，则：subwrd