第九章Grads函数Word下载.docx

1、2功能：通用的求平均函数。3说明：expr是由dexpr1和dexpr2定义的维数围求平均的量。缺测值不参与求平均运算。4假设求平均的维数是时间，可设定时间增量tincr，其单位与dexpr1和dexpr2一样。5flags为选项-b,表示对每个网格点取同样的权重求平均，包括终端点。例如：求纬向平均Set x 1Set t 1 48Set lev 200set lat 0 40Uave=ave（u,lon=0,lon=360,-b）。求时间平均：set t 1set lon 60 150set z 1a=ave（u,t=1,t=48）d amean函数用法同上,但不加权2、aave函数aave

2、（expr,xdim1,xdim2,ydim1,ydim2）求面积平均。expr: 任何表达式xdim1:X 维数方向的起始维数表达式xdim2: X 维数方向的终止维数表达式ydim1: Y 维数方向的起始维数表达式ydim2: Y 维数方向的终止维数表达式1aave函数采用地球坐标进展运算，运算时考虑了不同纬度的权重。Set z 1 Set x 1;set y 1a=aave（u,lon=110,lon=150,lat=0,lat=40） 2amean函数用法用上，不加权。3、tmave（）高级篇4、asum（）和asumg对指定区域求和二维求和5、sum（）和sumg（） 一维求和函数6

3、、max（）/min（）（高版本）返回指定表达式的最大/最小值maxloc/minloc返回指定表达式的最大/最小值的维数值格点系列使用这个命令后,括号的维数自动固定. 设定t/z/固定，看每个纬带上沿x方向的最大值d max（u,x=1,x=37）d maxloc（u,x=1,x=37）7、vint函数vint（psexpr,expr,top）质量加权垂直积分函数，计算。psexpr表示地面气压变量，expr表示被积变量，top表示积分上限（须为常数，单位hPa），本函数要求定义水平维数环境。如：vint（ps,u,100）表示对u变量从ps到100hPa进展质量加权垂直积分。4gint（e

4、xpr,dim1,dim2）普通积分函数慎用8、scorr（） 求空间相关系数两个场的空间相关系数是一个值，scorr函数之前的x/y维的设置没多大作用，函数会自动固定x/y维，但如果设时间维变动，得到的是相关系数时间序列。*时间序列set x 1;set y 1set z 1set t 1 48d tloop（scoor（u,v,x=1,x=37,y=1,y=17） 剖面z-t：set z 1 2a=tloop（scoor（u,v,x=1,x=37,y=1,y=17）d a9、tcorr（） 求时间相关函数10、空间回归/时间回归reinitopen e:liugddataset z 1se

5、t t 1 set x 1set y 1corr=sregr（u,v,x=1,x=37,y=1,y=17）d corr-求回归系数，一个数值;cuave=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-回归场值-*set z 1 2set t 1 48set x 1set y 1corr=tloop（sregr（u,v,x=1,x=37,y=1,y=17）set x 1 37set y 1 17b=corr*（u-uave）+vaved ba=aave（u,x=1,x=37,y=1,y=17

6、）set t 1*d tregr（a,u,t=1,t=48）-回归系数corr=tregr（a,u,t=1,t=48）aave1=ave（a,t=1,t=48）uave=ave（u,t=1,t=48）d corr*（a-aave1）+uave-回归值四、网格处理/平滑/差分函数1、const函数const（expr,constant设置局部网格点的值取为常数constant。所有非缺测格点处的expr值取为常数，3flag说明：-a: 所有网格点值均设定为指定的常数-u:把缺测格点处的expr值设定为常数。该函数对格点和台站资料均适用。display const（p,0,-u）将变量p在缺测格

7、点处的值设为常数0。又如：在linefill图形类型的输出中，假设须将其中的一条线设为水平直线： set lon -90 set lat -90 90 set gxout linefill set lev 500display const（t,-20）;t-273.162、maskout函数1）格式：maskout（expr,mask）2）功能：标记函数。当mask所在的网格点值取为负值时，对应格点上的expr的值设定为缺测值；当mask所在的网格点值取为正值时,对应格点上的expr值不变。3）说明：该函数对格点和台站资料均适合。file 1中的tsfc表示全球地表温度（在海洋上为SST）资料

8、，file 2中的mask表示海陆标记（海洋上为0，陆地上为1）资料，假设需显示SST的分布，那么可：maskout（tsfc,0.5 - mask.2）3、skip函数skip（expr,skipx,skipy）设定样本的取样密度。display skip（u,2）表示在x,y方向上对变量u隔一个网格点取一次值构成要显示的场，display skip（u,1,3）;v表示在y方向隔二个点取一次值，而x方向取所有的格点进展显示操作。该函数主要用于对矢量场的稀疏化显示。5、cdiff函数cdiff（expr,dim）中央差分函数expr表示进展差分运算的量，dim 表示进展差分运算的维数方向，为

9、X,Y,Z中的任一个字符。边界格点的运算值设为缺测。例如，以下运算是计算水平涡度，其效果与hcurl函数完全一样： define dv=cdiff（v,x） define dx=cdiff（lon,x）*3.1416/180 define du=cdiff（u*cos（lat*3.1416/180）,y）define dy=cdiff（lat,y）*3.1416/180display （dv/dx - du/dy）/（6.37e6*cos（lat*3.1416/180）以下运算是计算温度平流： define dtx=cdiff（t,x） define dty=cdiff（t,y） define

10、 dy=cdiff（lat,y）*3.1416/180 display -1*（u*dtx）/（cos（lat*3.1416/180）*dx）+v*dty/dy）/6.37e66、smth9函数smth9（expr）9点平滑函数。expr表示须平滑的表达式。该函数主要用于改善图形的输出质量，使其更加美观。当格点资料为1维时，9点平滑退化为3点平滑。7、tloop函数tloop（expr）通过循环时间维数逐个逐个局地固定时间维来求expr的值，即对expr固定其时间维数求值，并在求值过程这对时间维进展逐个循环。典型用法是求面积平均的时间序列，由于aave函数要求时间维固定，这时必须通过tloop

11、函数才能实现这一目的，见下例： set x 1 set y 1 set z 3 set t 1 31 d tloop（aave（ts,lon=0,lon=360,lat=-90,lat=90）为求ts变量的全球面积平均的时间序列，须用前四行语句来设定维数环境，即固定x,y,z,而让时间t变化，再用tloop函数，让aave在求面积平均时固定时间维，并逐时刻循环求值。8、寻找特征值五、物理量计算/站点函数1、gr2stn函数gr2stn（grid_expr,stn_expr）将格点值回插至台站值。grid_expr为代表格点资料的表达式，插将采用该格点资料进展，stn_expr为代表台站资料的表

12、达式，插结果将放至其中。为检验客观分析资料的质量，可将客观分析后得到的格点资料回插到原来台站，看其与原始台站资料之间的差值大小，即：d t.3-gr2stn（t.1,t.3）其中file 1 为格点资料，file 3 为台站资料。2、Coll2gr（）3、hcurl函数hcurl（uexpr,vexpr）计算垂直涡度。uexpr表示U分量，vexpr表示V分量，风速单位用m/s。边界上的涡度值设定为缺测。d hcurl（u,v）4、hdivg函数hdivg（uexpr,vexpr）计算水平散度。5、tvrh2q函数tvrh2q（tvexpr,rhexpr）由虚温tvexpr和相对湿度rhexp

13、r计算比湿。tvexpr的单位为0K，rhexpr的单位为百分比（取值围0100），运算结果的单位为g/g。该函数仅适合用格点资料。6、tvrh2t函数tvrh2t（tvexpr,rhexpr）由虚温tvexpr和相对湿度rhexpr计算温度。运算结果的单位为0K，其余说明参见tvrh2q函数。7、oacres函数oacres（grid_expr,stn_exprCressman插函数，将台站资料插至格点上。grid_expr表示格点资料表达式，stn_expr表示台站资料表达式，选项radii表示影响半径，一般采用多个影响半径，它们之间用逗号分开，如：10,7,4,2,1（省缺值）。oacr

14、es（ts,ts.2） , oacres（ts,ts.2,12,8,5,4,3,2,1）8、stnave函数stnave（expr,dexpr1,dexpr22） 功能：台站资料时间平均函数。3） 说明：expr表示欲平均的台站资料表达式，dexpr1表示起始时间，dexpr2表示完毕时间，选项-mt 为质量控制项，表示参与平均运算的台站资料样本数最小不得少于cnt，否那么运算结果设定为缺测。cnt的省缺值为1。stnave（ts,t=1,t=20,-m 10）表示对头20个时刻的ts进展平均，假设样本数少于10，那么平均的结果为缺测。9、stnmin函数stnmin（expr,dexpr1,dexpr2时段台站资料最小值。参看stnave函数的说明。10、stnmax函数1） 格式：stnmax（expr,dexpr1,dexpr2时段台站资料最大值。11、oabin