GrADS绘图软件实用手3

《GrADS绘图软件实用手3》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三章GrADS数据格式每一组GrADS数据应至少包括两组数据文件，数据描述文件—ASCII码和数据文件—二进制，数据的真正存放地。数据文件中只是用户数据的有序排放，而关于数据种类、排放次序等是单独放在一个文件中的称—数据描述文件。而象GRIB和NETCDF等通用数据格式，以上两者是存于同一个文件的—或称为自定义/自解释格式数据。但考虑到GrADS传统，对这类自定义格式数据仍将生成相应的数据描述文件。GrADS并不直接使用“数据文件”，而是通过“描述文件”间接使用“数据文件”。上一章中我们已使用过了

2、这样的一组数据。以此为例，介绍用户如何按GrADS的格式，将自己的数据生成相应的数据文件和数据描述文件。1.格点数据描述文件model.le.ctl文件清单：以*开始的行为注解行。1.1数据描述文件各项解释1.DSET<路经/>数据文件名定义与此数据描述文件相对应的数据文件名。若两者位于同一目录，前面的路经可以省略或以“^”开始，代表两者位于同一目录。若不在同一目录下，应给出路经参数。如：c:/pcgrads/sample/model.le.dat注意路经的给法与DOS不同，而与UNIX环境一致，便

3、于移植！或-113-c:pcgradssamplemodel.le.dat,两种都行。1.TITLE数据文件说明文字串。2.UNDEFvaule定义缺测值。一般给一很大的正/负值，表示，当取值超过这一正值/低于定义的负值，认为该值无效。（GrADS采用跳过或用周围有效点的值处理。）3.OPTIONS这里定义了与二进制存储有关的选项，二进制存储的一大特点是可移植性差，因此通过keywords项来增加可移植性。若keywords省略，则OPTIONS也可省略。

4、>可取：sequential:顺序无格式方式。yrev:Y维与YDEF定义相反方式存放。zrev:Z维与ZDEF定义相反方式存放。big_endian:如数据是在sun,sgi,hpcray机器上生成的，而目前不在此类机器上使用。little_endian:如数据是在iX86,dec机器上生成的，而目前不在此类机器上使用。byteswapped：反序位存放。4.XDEFnumberLINEARX_Startincrementnumber（>=1）给定X方向格点数，其后LINEAR参数指明X方向是等间

5、隔分布格点，X_Start起点坐标，increment网格间距。XDEFLEVELSvalue_listLEVELS参数指明X方向是不等间隔分布格点，因此，其后要给出具体每个格点的坐标值（以空格分开）。5.YDEFnumberLINEARY_StartincrementYDEFLEVELSvalue_list6.ZDEFnumberLINEARZ_StartincrementZDEFLEVELSvalue_list7.TDEFnumberLINEART_StartincrementTDEFLEVEL

6、Svalue_listT_Start和value_list给出时间格式如下：hh:mmZddmmmyyyyhh要以两位数代表小时（缺省00）；“：mm”以两位数代表分钟（缺省：00）；dd以两位或一位数代表日期（缺省1）；mmm为英文月分三个字符的缩写；yyyy以两位或四位数代表年。如：12Z1JAN90（省：mm）14：45Z22JAN1987JUN1960时间increment格式：vvkkvv以两位数字代表；kk取：mn-增量以分为单位；hr-以小时为单位；dy-以天为单位；mo以月为单位；y

7、r以年为单位。8.VARSnumberV_abrevlevsunitsdescription……ENDVARS以VARS指示开始定义变量，number代表变量总的个数；以ENDVARS指示结束变量定义。从VARS以下，每行定义一个变量。V_abrev变量名称，GrADS将用到；levs，-113-数字，代表变量层数，0表示只有一层；units单位，为GRIB预留，给99；description，对变量的文字描述。1.2生成model.le.dat和model.le.ctl文件的程序代码片段progr

8、ammainRealps(72,46,7,5),ts(72,46,7,5),z(72,46,7,5),u(72,46,7,5),v(72,46,7,5)Realt(72,46,7,5),q(72,46,7,5),p(72,46,7,5)C第一部分，读入model.le.dat数据。Open(7,file=’model.le.dat’,form=’unformatted’,access=’direct’,recl=72*46*4)C在微机上，要乘4；用VisualFor