第九篇Grads函数.doc

GRADS (Grid Analysis and Display System) 讲义（9）函数 内容提要 概况 函数详细介绍 概况 1、函数的调用 1）通过函数名直接引用， 2）参数放在括号中用逗号分开， 3）可以嵌套调用， 4）有些函数在运算时会改变维数环境。 2、函数的类型 二、数学函数 注意：sin/cos/tan等三角函数单位是弧度； asin等反三角函数结果单位是弧度。 三、统计积分函数 1、ave函数 1）格式： ave(expr,dexpr1,dexpr2,tincr,flags) 2）功能： 通用的求平均函数。 3）说明：expr是由dexpr1和dexpr2定义的维数范围内求平均的量。缺测值不参与求平均运算。 4）若求平均的维数是时间，可设定时间增量tincr，其单位与dexpr1和dexpr2一样。 5）flags为选项-b,表示对每个网格点取同样的权重求平均，包括终端点。 例如：求纬向平均 Set x 1 Set t 1 48 Set lev 200 set lat 0 40 Uave=ave(u,lon=0,lon=360,-b)。 求时间平均： set t 1 set lon 60 150 set lat 0 40 set z 1 a=ave(u,t=1,t=48) d a 注意：mean函数用法同上,但不加权 2、aave函数 1）格式： aave(expr,xdim1,xdim2,ydim1,ydim2) 2）功能：求面积平均。 3）说明： expr: 任何表达式 xdim1:X 维数方向的起始维数表达式 xdim2: X 维数方向的终止维数表达式 ydim1: Y 维数方向的起始维数表达式 ydim2: Y 维数方向的终止维数表达式 注意： 1）aave函数采用地球坐标进行运算，运算时考虑了不同纬度的权重。 例如： Set t 1 48 Set z 1 Set x 1;set y 1 a=aave(u,lon=110,lon=150,lat=0,lat=40) d a 2）amean（）函数用法用上，不加权。 3、tmave()高级篇 4、asum()和asumg 对指定区域求和（二维求和） 5、sum()和sumg() 一维求和函数 6、max()/min()(高版本) 返回指定表达式的最大/最小值 maxloc/minloc（） 返回指定表达式的最大/最小值的维数值（格点）系列 使用这个命令后,括号内的维数自动固定. 例如： 设定t/z/固定，看每个纬带上沿x方向的最大值 d max(u,x=1,x=37) d maxloc(u,x=1,x=37) 7、vint函数 1）格式：vint(psexpr,expr,top) 2）功能：质量加权垂直积分函数，计算 。 3）说明：psexpr表示地面气压变量，expr表示被积变量，top表示积分上限(须为常数，单位hPa)，本函数要求定义水平维数环境。如：vint(ps,u,100)表示对u变量从ps到100hPa进行质量加权垂直积分。 4）gint(expr,dim1,dim2)普通积分函数（慎用） 8、scorr() 求空间相关系数 注意：两个场的空间相关系数是一个值，scorr函数之前的x/y维的设置没多大作用，函数会自动固定x/y维，但如果设时间维变动，得到的是相关系数时间序列。 *时间序列 ‘set x 1’;’set y 1’ ‘set z 1’ ‘set t 1 48’ ‘d tloop(scoor(u,v,x=1,x=37,y=1,y=17)) 剖面z-t： ‘set x 1’;’set y 1’ ‘set z 1 2’ ‘set t 1 48’ ‘a=tloop(scoor(u,v,x=1,x=37,y=1,y=17)) ‘d a’ 9、tcorr() 求时间相关函数 10、空间回归/时间回归 reinit open e:\liugd\data set z 1 set t 1 ‘set x 1’ ‘set y 1’ corr=sregr(u,v,x=1,x=37,y=1,y=17) d corr----求回归系数，一个数值 ; c uave=aave(u,x=1,x=37,y=1,y=17) vave=aave(v,x=1,x=37,y=1,y=17) d corr*(u-uave)+vave----回归场值 ; ------------ reinit open e:\liugd\data set z 1 *s