GrADS函数大全
- 格式:doc
- 大小:34.50 KB
- 文档页数:2
目录前言 (5)第一章启动GrADS 及打开文件 (5)一.启动GrADS (5)二.打开文件 (6)第二章维数环境 (6)set lat|lon|lev|time val1<val2> (6)set x|y|z|t vall <val2 > (6)第三章显示命令及动画 (7)一、显示命令 (7)二、动画 (7)set loopdim x︱y|z|t (7)set looping on|off (7)第四章图形输出类型 (8)对格点数据: (8)bar: 直方图 (8)barb:风向杆形式绘二维风场 (8)contour:二维等值现图 (8)errbar:单线图及误差分布 (8)frgid:用指定颜色填充二维格点场 (8)fwrite:图形不在屏幕上显示,而是将输出结果存入一个由“set fwrite 文件名”所指定的文件中。
(8)gid:以网格形式在各网点中央标出该点数值 (8)line:单线图 (8)linefill:两单曲线之间填色 (8)scatter:散点图 (8)shaded:二维填色图 (8)stat:输出有关的统计量 (8)stream:流线形式绘二维风场 (8)vector:向量间箭头形式绘二维风场 (8)对站点数据: (8)barb:在各站点绘风向标 (8)findstn:搜索最近的站点(详见描述语言部分) (8)model:以天气图形式将天气观测各分量填放在站点四周 (8)value:在各站点标值 (8)wxsym:绘wx天气符号 (8)Draw wxsym symbol x y size<color<thicks>> (9)对于model,命令为set gxout model ,然后再添图,形式为displayu;v;t;d;slp;delta;cld;wx;vis (9)第五章 GrADS资料格式 (9)格点资料资料描述文件 (9)以下是一个资料描述文件的例子: (9)1、DSET data-set-name (11)2、TITLE string (11)3、UNDEF value (11)4、OPTIONS <keywords> (12)5、XDEF number<LINEAR start increment>或XDEF number<LEVELS value-list> (12)6、YDEF number mapping start<increment>或YDEF number<LEVELS value-list> (13)7、ZDEF number mapping <start increment>或ZDEF number<value-list> (13)8、TDEF number LINEAR start-time increment (14)9、VARS number (14)二. 台站资料资料描述文件 (15)第六章变量名和表达式 (15)一、变量名和表达式 (15)二、表达式 (16)第七章函数 (16)一. aave函数 (16)二.abs函数 (16)三.acos函数 (17)四、asin函数 (17)五、atan2函数 (17)六、ave函数 (17)七、cdiff函数 (17)八、const函数 (18)九、cos函数 (18)十、exp函数 (18)十一、gr2stn函数 (19)十二、hcurl函数 (19)十三、hdivg函数 (19)十四、log函数 (19)十五、log10函数 (19)十六、mag函数 (19)十七、maskout函数 (20)十八、oacres函数 (20)十九、pow函数 (20)二十、sin函数 (20)二十一、skip函数 (20)二十二、smth9函数 (21)二十三、sqrt函数 (21)二十四、stnave函数 (21)二十五、stnmin函数 (21)二十六、stnmax函数 (21)二十七、tan函数 (22)二十八、tloop函数 (22)二十九、tvrh2q函数 (22)三十、tvrh2t函数 (22)三十一、vint函数 (23)第八章图形要素设置 (23)对于图形类型为contour起作用的设置 (23)set ccolor color——设置等值线颜色 (23)set ccolor rainbow——设定等值线颜色用彩虹表示 (23)set cstyle style——设定等值线线型 (23)Set cthick thckns——设定等值线线宽hckns(取值范围1~20) (23)Set cterp on|off——设置样条差值光滑开关 (23)Set clab on|off|forced|string|auto——控制等值线的标记方式。
第九章Grads函数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 1Set t 1 48Set lev 200set lat 0 40Uave=ave(u,lon=0,lon=360,-b)。
求时间平均:set t 1set lon 60 150set lat 0 40set z 1a=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 x 1;set y 1a=aave(u,lon=110,lon=150,lat=0,lat=40)d a2)amean()函数用法用上,不加权。
【转】
GrADS的函数分两类,
一类是对格点/站点数据执行运算的,这一类我们姑且称之为分析函数;
另一类是脚本编程(gs)时使用的, 这后一类我们姑且称之为脚本函数.
第一类分析函数又分为格点分析和站点分析两个子类:
1) 格点分析
绝对值abs(expr)
余弦/反余弦cos(expr) acos(expr)
正弦/反正弦sin(expr) asin(expr)
正切/反正切tan(expr) atan2(expr1, expr2) [expr1/expr2]
指数exp(expr)
对数log(expr) log10(expr)
幂函数pow(expr1, expr2) expr1expr2
直角三角形的弦长/全风速mag(u_expr, v_expr)sqrt(u*u+v*v)
平方根sqrt(expr)
加权/不加权平均ave (expr, dim_expr1,dim_expr2, tincr, { |-b } ) mean()
加权/不加权区域平均aave (expr, { xdim1,xdim2,ydim1,ydim2 | global | g } ) amean()
剔除值后时间求平均tmave (undefine_or_weight_expr, expr, texpr1, texpr2)
加权/不加权区域和asum (expr, { xdim1,xdim2, ydim1, ydim2| global | g } ) asumg()
质量加权垂直积分vint (ps_expr, expr, top)
普通积分gint (expr)
最大/小值Max (expr, dim1, dim2, tincr) min()
最大/小值所在格点Maxloc (expr, dim1,dim2, tincr) minloc()
空间相关/回归scorr (expr1,expr2, {xdim1, xdim2, ydim1, ydim2| global | g}) sregr() 时间相关/回归tcorr (expr1, expr2, tdim1, tdim2) tregr()
常数const (expr, const, { | -u | -a } )
剔除maskout (expr, mask_expr )
跳点skip (expr, skip_x, skip_y )
加权9点平滑smth9 ( expr )1-D 3点平滑
中央差分cdiff ( expr, dim )
时间循环tloop (expr)
极值位置maxloc (expr, dim1, dim2, tincr)minloc()
特征值位置findlvl (expr, expr_to_find, lev1, lev2)
气象上物理量计算的函数:
虚温,相对湿度求比湿tvrh2q(tv_expr, rh_expr)g/g
虚温,相对湿度求温度tvrh2t(tv_expr, rh_expr)degK
垂直涡度hcurl(u_expr, v_expr)
中央差求水平散度hdivg(expr1, expr2)
全风速mag(u_expr, v_expr)
2) 站点分析
从格点插到站点gr2stn(grid_expr, { stn_expr | lon, lat }, -a )
收集生成格点场coll2gr(cnum, { num | -u } )
Cressman插值oacres(grid_expr, stn_expr, radii, 1st_guess)
网格平均oabin(grid_expr, stn_expr, { -f | -c } )
时间平均stnave(expr, dexpr1, dexpr2, -m count)
时间序列最小值stnmin(expr, dexpr1, dexpr2, -m count)
时间序列最大值stnmax(expr, dexpr1, dexpr2, -m count)
转成1-D格点序列s2g1d(expr)
所有数学函数abs,cos,acos,sin,asin,tan,atan2,exp,log,log10,pow,mag,sqrt
第二类, 脚本函数包括有:
abs, log, log10 : math_abs, math_log, math_log10
sin, cos, tan : math_sin, math_cos, math_tan
asin, acos, atan2: math_asin, math_acos, math_atan
math_sinh, math_cosh, math_tanh
math_asinh, math_acosh, math_atanh
pow, sqrt : math_pow, math_sqrt, math_exp
取整: math_nint, math_int
求余: math_fmod, math_mod
字符串处理函数: math_strlen, valnum, wrdpos, strlen
sublin, subwrd, substr
数值格式转换函数: math_format
文件操作函数: read, write, close。