GrADS程序设计第6章变量和函数解析.ppt

第6章 变量和函数 在GrADS运行环境中，参与操作的数据对象可以是变量、函数和表达式形式。 6.1 变量名 完全的变量名形式 abbrev.file#（dimexpr， dimexpr，…） 其中，abbrev：是数据描述文件中给出的变量名缩写；file#：为包含此变量的已打开的文件序号。 例如： ga-open model.ctl open model.le.ctl 上例中，model.ctl文件的序列号为1，model.le.ctl文件的序列号为2。序列号1为缺省情况，可以省略。 例如: ga d z.2 (显示已打开的2号文件中的变量Z) 括号内的dimexpr：是对当前维数环境进行的局域维数设置表达式，该表达式仅对该变量的维数环境进行局域修正，因而仅对固定的维有效。 设置形式可以采用以下两种方式： 绝对维数表达式为： X|Y|Z|T|LON|LAT|LEV|TIME = value （value为绝对维数值）； 相对维数表达式为： X|Y|Z|T|LON|LAT|LEV|TIME +/- offset （offset为相对于当前维数环境设置的偏差维数值）。 例如： z.3（lev=500） （表示文件3中高度为 500hPa等压面上的变量Z） tv.1（time-12hr） （相对于当前时刻之前12 小时时刻的1号文件中的变量tv） rh （缺省的当前文件中的变量rh） q.2（t-1,lev=850） （2号文件中相对于当前 时刻的前一时刻，高度为 850hPa等压面上的变量q） 2. 举例说明 ga-open model.ctl open model.le.ctl （同时打开两个文件，model.ctl 文件的序列号为1，model.le.ctl 文件的序列号为2） set t 1 （将时间维数固定在第1个时次） d z.2(lev=500) （显示2号文件中500hPa等压面 上的变量z(绝对维数)） set z 3 （设置垂直高度层次为第3个层次， 即700hPa） d z(z+1) （显示1号文件中相对于第3个层次 的后一个层次，即500hPa面上的 变量z （相对维数）） 6.2 表达式 表达式的构成 运算符： + 加，- 减，* 乘，/ 除； 运算域： 变量, 函数和常数； 括号：用于控制运算的次序。 2. 运算规则 对相同网格点上的不同变量进行运算时，只要有一个变量在某格点的值为缺测，则该网格点的运算结果为缺测值；当被0除时的结果也设为缺测。 对多个数据文件的数据做运算时，两种数据网格要一致，即运算对象的维数变化范围要一致。 如某一格点的变动维数多于其他格点，则具有较少变动维数的格点其维数环境将被拓展，以便于运算操作。 例如： z - z(t-1) (表示变量Z的时间变化) t(lev=500)-t(lev=850) (表