GrADS绘图软件实用手.docVIP

AUTHOR zhu he - PAGE 102 - - PAGE 101 - GrADS数据格式 每一组GrADS数据应至少包括两组数据文件，数据描述文件—ASCII码和数据文件—二进制，数据的真正存放地。数据文件中只是用户数据的有序排放，而关于数据种类、排放次序等是单独放在一个文件中的称—数据描述文件。而象GRIB和NETCDF等通用数据格式，以上两者是存于同一个文件的—或称为自定义/自解释格式数据。但考虑到GrADS传统，对这类自定义格式数据仍将生成相应的数据描述文件。 GrADS并不直接使用“数据文件”，而是通过“描述文件”间接使用“数据文件”。 上一章中我们已使用过了这样的一组数据。以此为例，介绍用户如何按GrADS的格式，将自己的数据生成相应的数据文件和数据描述文件。 格点数据描述文件 model.le.ctl文件清单： 以*开始的行为注解行。 1.1 数据描述文件各项解释 DSET 路经/数据文件名 定义与此数据描述文件相对应的数据文件名。若两者位于同一目录，前面的路经可以省略或以“^”开始，代表两者位于同一目录。若不在同一目录下，应给出路经参数。如：c:/pcgrads/sample/model.le.dat 注意路经的给法与DOS不同，而与UNIX环境一致，便于移植！或c:\pcgrads\sample\model.le.dat,两种都行。 TITLE 数据文件说明文字串。 UNDEF vaule 定义缺测值。一般给一很大的正/负值，表示，当取值超过这一正值/低于定义的负值，认为该值无效。（GrADS采用跳过或用周围有效点的值处理。） OPTIONS keywords 这里定义了与二进制存储有关的选项，二进制存储的一大特点是可移植性差，因此通过keywords项来增加可移植性。若keywords省略，则OPTIONS也可省略。 keywords可取： sequential: 顺序无格式方式。 yrev: Y维与YDEF定义相反方式存放。 zrev: Z维与ZDEF定义相反方式存放。 big_endian: 如数据是在sun,sgi,hp cray机器上生成的，而目前不在此类机器上使用。 little_endian: 如数据是在iX86,dec机器上生成的，而目前不在此类机器上使用。 byteswapped：反序位存放。 XDEF number LINEAR X_Start increment number（=1）给定X方向格点数，其后LINEAR参数指明X方向是等间隔分布格点，X_Start起点坐标，increment网格间距。 XDEF LEVELS value_list LEVELS参数指明X方向是不等间隔分布格点，因此，其后要给出具体每个格点的坐标值（以空格分开）。 YDEF number LINEAR Y_Start increment YDEF LEVELS value_list ZDEF number LINEAR Z_Start increment ZDEF LEVELS value_list TDEF number LINEAR T_Start increment TDEF LEVELS value_list T_Start和value_list给出时间格式如下：hh:mmZddmmmyyyy hh要以两位数代表小时（缺省00）；“：mm”以两位数代表分钟（缺省：00）；dd以两位或一位数代表日期（缺省1）；mmm为英文月分三个字符的缩写；yyyy以两位或四位数代表年。如： 12Z1JAN90（省：mm） 14：45Z22JAN1987 JUN1960 时间increment格式： vvkk vv以两位数字代表；kk取：mn-增量以分为单位；hr-以小时为单位；dy-以天为单位；mo以月为单位；yr以年为单位。 VARS number V_abrev levs units description …… ENDVARS 以VARS指示开始定义变量，number代表变量总的个数；以ENDVARS 指示结束变量定义。从VARS以下，每行定义一个变量。 V_abrev变量名称，GrADS将用到；levs，数字，代表变量层数，0表示只有一层；units单位，为GRIB预留，给99；description，对变量的文字描述。 1.2 生成model.le.dat和model.le.ctl文件的程序代码片段 program main Real ps(72,46,7,5),ts(72,46,7,5),z(72,46,7,5),u(72,46,7,5),v(72,46,7,5) Real t(72,46,7,5),q(7