第6章 变量和函数.ppt

1、第6章变量和函数是在GrADS执行环境中，操作所涉及的数据对象可以是变量、函数、式形式。 6.1变量名，完全的变量名形式abbrev.file#(dimexpr，dimexpr，)，其中，abbrev :是赋予数据描述文件的变量名的省略形式。 file# :包含此变量的打开文档的序号。 例如，在ga-opend 3360 model.ctl opend 3360 model.le.ctl上的示例中，model.CTL文件的序列号是1，model.le.CTL文件的序列号是2，序列号1是缺省值括号中的dimexpr :当前维环境的局部维设置公式，该公式仅对该变量的维环境进行局部修改，因此仅对固定

2、维有效。 设定形式为x|y|z|t|lat|lev|time=value (value为绝对维度数值)； 相对维数表达式是x|y|z|t|lon|lev|time/- offset (offset是针对当前维数环境设定的偏差维数值)。 例如，对于z.3(lev=500) tv.1(time-12hr) rh默认的当前文件的变量rh q.2(t-1、lev=850 )、第2文件当前时刻的12小时前的时刻的第1文件中的变量tv， dset model.le.datoptionslittle _ endian Cray _ 32bit _ IEEE undef-2.56 e 33title 5大型机

3、模型72.0 zdef 7级别10008500 an 19871 dy var s8PS 099 surfacepressureu 799 uwindsv 799 vwindsz 799 geo端口初始化799临时q 599规范099 surfacete dset model.grbtitlesamplemodeldataforlats4dtutorialundef1e 20 dtypegribindexmodel.gmpxdef 72线性电脑0.0000005.00000 ydef 46线性电脑000 50700500300200 TDE f 5线性0 z1Jan 19871 dy var

4、s8PS 0，1，0，0 surfacepressurehpaua 733，100 eastwardwindm/sva 734，100 northwardwindm/s 100 geopotentialheightmta 711，100 airtemperaturekhus 751，100 specifichumiditykg/kgts 011，105，2曲面(2m ) airtemperaturekpr 059 将0 totalprecipitationrate ga-opend : model.CTL opend 3360 model.le.CTL set 1时间维数固定为第一个时间维数d

5、 z.2(lev=500) set z 3的垂直高度电平设定为第三个电平， 即，设为7000的两个文件： model.ctl文件的序列号为1，model.le.ctl文件的序列号为2 .显示第二个文件的500hPa等压面上的变量z (绝对维数)，第一个文件的第三层次2 .在对相同网格点上的不同变量进行运算时，如果一个变量在某个网格点处的值是缺测，则将该网格点的运算结果除以成为缺测值的运算规则0时的结果也视为缺测。 在计算多个数据文件中的数据时，两个数据网格必须匹配。 也就是说，运算对象的维数的变化范围必须一致。 当一个网格点的变动维数大于其他的网格点时，对于变动维数少的网格点，扩展维数环境以便

6、于运算操作。 例如，z-z (t-1 ) t (lev=200 )-t (lev=850 ) ave (z，t=1，t=5)z- .调用方式采用直接引用方式，形式例如，函数名(残奥仪表1、残奥仪表2、)、2 .常用函数、(1)ave函数例如，ga-open model.ctlsetlev 500 Dave (z，t=1，t=5)表示500hPa等压面上的从第一次到第五次的变量z的平均量，setp area 1816 setgridoffsetgradsoffs e a=ave(v，t=7，t=660，12 ) daprintc，例如求出纬向的平均： ave(expr，lon=00，lon=。

7、如果对终点取一半的权重，则可以使用ave(expr、lon=0、lon=360、-b )。setxyrevonsetlev 500 set t1set lat 090 setx1definea=ave (v、lon=0、lon=360、-b)dasetcmark0setcolor (2) const函数格式： const (退出恒定(constant )功能：将部分网格点的值设为常数的恒定说明：将所有非缺陷网格点的expr值设为常数，将flag设为可选项，加上选项-a，此函数同时适用于网格点和站数据。 (3)aave函数格式：求出aave(expr，xdim1，xdim2，ydim1，ydim

8、2)功能：面积平均。 说明： expr:任意式xdim1: X维方向的开始维数学式xdim 23360 x维方向的结束维数学式ydim 13360 y维方向的开始维数学式ydim 23360 y维方向的结束维数学式set lev 850 sett 148 setx1set y1definea=aavv lon=100、lon=120、lat=30、lat=40) d a、设定lon 60 200设定lat-2080 d mag (u.1、v.2 )、(5)hcurl函数格式： hcurl(uexpr、vexpr )功能边界上的涡度值设定为缺失测量。 6)hdivg函数格式： hdivg(uex

9、pr，vexpr )功能：修正水平分散度。 说明： uexpr是u风速成分，vexpr是v风速成分，风速单位是m/s。 set lon 60 200 set lat -20 80 d hdivg (u，v )、(7)skip函数格式： skip(expr，skipx，skipy )功能：设定样本的采样密度。 说明： skipx、skipy的数值决定x方向和y方向的采样密度(取值1时可以省略) *此函数主要用于向量场的稀疏化显示。 d u.1； v.2、德斯基普(u.1、2、2 )； skip (v.2，2，2 )，(8)cdiff函数格式： cdiff(expr，dim )功能：中央差分函数说

10、明： expr表示进行差分运算的量，dim表示进行差分运算的次元的方向。 边界网格的运算值被设定为预设值。 例如，通过校正温度平流，设定定义dtx=c diff (t，x )定义DTY=c diff (t，y )定义dx=c diff (lon，x )。 y ) *3. 1416/180 d-1 * (u * dtx )/(cos (lat *3. 1416/180 ) * dx ) v * DTY/dy )此函数主要用于改善图形的输出品质，使其更美观。 格点数据为一维时，9点平滑化劣化为3点平滑化。 d smth9(b1)、(11)vint函数格式： vint(psexpr，expr，top

11、 )功能：质量加权垂直积分函数说明： psexpr是地面气压变量，expr是vint(psexpr，expr，top )=f/g * sopene : bsuq.ctlopene : bsp RES.ctlsetgxoutfwritesetfwritee : bsmonzuq.grdi=1while (I=660 )设置1定义常见问题=vint SETI setx 1144 sety 173 setz1d fuqi=I 1生成endwhiledisaad (2) map文件(3)创建网格文档(创建后也可以继续使用) (4)通过插值函数将站点数据插值到网格文件中，并绘制文件的时间必须与站点的数据

12、时间一致。 范围： lat:15 55 N lon:70 140 E格式：1*1 grid.grd文件中的每个点一个，格式文件grid.grd : 分配在第2章数据处理中写入70-73页的网格数据grid.grrd。与该网格文件对应的数据描述文件grid.ctl， 作为dsetd 3360 datagrid.grd undef-999.0 titlesamplegriddataxdef 71线性化70 ydef 41线性化15 z de f 1线性化的f 1线性化jul 19511移动1 g 099 gr 如果有该网格文件和与其对应的数据描述文件，则可通过GrADS插值函数将站点数据插值到该网

13、格文件中来使用。 (12) oacres函数格式： oacres(gexpr，sexpr )功能： cressman插值函数说明：站点数据插值中经常使用的网格数据。 为了避免gexpr大于站点数据范围的网格数据变量名称sexpr站点数据变量名称radii影响半径，缺省值为：10、7、4、2、1插值造成的虚假结果，参照网格的间隔最好接近站点间隔。 将站点数据插值到作为插值函数的网格文件后使用。 中的组合图层性质变更选项。 与此网格文件相对应的数据描述文件grid.ctl是dsetd 3360 datagrid.grd undef-999.0 titlesamplegriddataxdef 71线性文件70 ydef