GrADS绘图系统在站点等值线中的应用.docVIP

GrADS绘图系统在站点等值线中的应用 第 2001姜期JOURNAlO广FG西UAN气GX象IMETEOROLOGY年6月V01.22No.2Jun.2001 文章编号:1001—5191(2001)02-0053-04 GRADS绘图系统在站点等值线中的应用 况雪源 (广西气候中心,广西南宁530021) 摘要;姒一个实例详细介绍了利用GRADS绘图系统绘制具有不规则边界的站点资料等值线的方法,使站点等值 线的绘制变得简单易行,图形输出更为美观.并对GRADS系统存在的一些问豚进行了讨论. 荚镑调;GRADS系统;站点等值线绘图 中囝分类号:TP317,4文献标识码B ApplicationofGRADSDrawingSysteminStationIsoline KUNXL!e—yuan (GuangxiClimateCenter,Nanning530022,China) Abstract:Bytakinganactualexample,thispaperintroducesindetailtheapplicationofGRADS drawingsystemtomappinginformationisolineforthestationswithirregularboundary.This methodmakesthemappingofstationisolineeasyandsimple,andmorehandsomeimageoutput.It alsodiscussestheexistingproblemsofGRADSdrawingsystem. Keywards:GRADSsystem;stationisoline;mapping 引言1基本思路 GrADS(GridAnalysisandDisplaySystem)绘 图系统是一组多平台应用及气象数据图形分析和显 示系统,该系统具有气象数据分析功能强,地图投 影丰富,高级编程语言容易,图形显示快速,并具 有彩色动画功能等特点,并已成为国际气象数据显 示的标准平台之一. GRADS系统提供了一个交互操作的气象格点 数据和站点数据的分析和显示环境,但相对于格点 资料,GRADS处理站点资料的能力较弱,且系统无 法直接采用站点资料来绘制等值线.对于一般用户 来说,如何利用GRADS绘制格点等值线功能及输 出图形美观等优点,将一定区域的站点资料以等值 线方式显示,并使等值线不越过不规则边界呢?刘 宣飞Lj对此作了一些介绍,但较为简略.本文将针对 这些问题作进一步探讨,并以绘制广西88个站点 2000年6月平均气温分布图为例作详细介绍.文中 程序部分用Basic语言编写,对GRADS系统的基本 使用可参看文献]. GRADS系统无法采用站点资料来绘制等值线, 而需将站点资料内插成格点资料后才能绘制等值 线,GRADS提供的oacres内部函数可实现此内插 功能,同时,内部函数maskout可将符合一定条件 的格点值设为缺测,为此,我们采取以下方式进行 站点等值线的绘制. 首先,形成一份格点资料文件,此文件提供插 值格点区域,并将边界外资料设为负值,边界内资 料设为正值,然后,采用oacres函数将站点资料内 插成格点资料,并利用maskout函数将边界外的插 值场设为缺测后绘制等值线,最后,利用广西区界, 地区界经纬度资料在图上定位以绘制边界. 2GRADS绘制站点等值线步骤 2.1所需基本资料及格式 为绘制广西88个站点2000年6月平均气温分 布图形,需以下资料(表1); 收稿日期:2000—07—31 作者简升;况雪源(1971--)t女.贵州赫章人,硬士.广西气候中心预测工程师,现从事短期气候珂工作 54广酉气象22卷 表l绘制等值线圈资料 2.2站点数据文件及其描述文件 2.2.1站点文件格式 站点资料采用二进制文件进行存储,按每个时 次一个报告的顺序直接记录各站资料,每个时次的 站点报告组叉分成两部分:地面报和高空报.站点 报数据的排放顺序为: ①提供站点经纬度的头记录 ②地面报变量 ③高空报变量 ④一个时次完成后加上一个特殊头记录(没有 数据组),表示本时次数据报结束. 头记录共占28个字节,包括: STAT:站点标识,字符型(string或char),占8个 字节;其赋值可以是任意数字,也可以是字 母. RLON:经度,单精度型(single),占4个字节. RLAT:纬度,单精度型(single),占4个字节. TIM:本报的时次,单精度型(single),占4个字节, NLEV:头记录后的数据组数,单精度型(single), 占4个字节,包括一个地面报组加上各高 度层数,取0时标志一个时次报文的结束. NFLAG:单精度型(single),占4个字节