GrADS实用绘图(基4)

11
LASG动力论坛版
六. 图形显示(D)2
客观分析函数的使用
– Cressman分析
— 站点绘图扩展
oacres(grid_expr, stn_expr, { radii | 10,7,4,2,1 } )
– 栅格平均
oabin(grid_expr, stn_expr, { | -f | -c } )
5
LASG动力论坛版
六. 图形显示(D)2
三类维数变化
x, y 变化 z变化 t 变化
— 缺省维数
命令顺序
set display 变量引用方式 draw var ( stid = 'station_id' ) print
6
LASG动力论坛版 SLP显示方式
set stnopts dig3| nodig3
20
LASG动力论坛版
七. 数据分析(A)
本节内容提要
数学函数 统计函数 网格处理函数 滤波平滑函数 差分函数 物理量计算函数 站点数据专用函数
— 函数的使用
21
LASG动力论坛版
七. 数据分析(A)
功 能
abs(expr) cos(expr) sin(expr) tan(expr) exp(expr) log(expr) pow(expr1, expr2)
存放在 系统变量 result中
16
六. set gxout其它选项(D)2—print / fwrite
输出格式
set prnopts format(%g) line_n(8) blank_n(1) { |u}
LASG动力论坛版
数据写屏输出 set gxout print binary数据输出 set gxout fwrite
本节内容回顾
站点数据处理
– 图样式及其定制 – 技巧 defiine 站点客观分析函数
统计和输出(set gxout其它选项)
– 统计信息 – ASCII和Binary数据输出(文件格式转换)
19
LASG动力论坛版
内容提要
请注意标记的地方
一. 引言 二. GrADS绘图软件包安装 三. 运行GrADS绘图软件包 四. 图形显示(D)1— 格点绘图 五. 数据准备(Gr) 六. 图形显示(D)2— 站点绘图 七. 数据分析(A) — 函数的使用 八. 批处理和脚本语言(S) 九. 图形保存和转换
15
LASG动力论坛版
六. set gxout其它选项(D)2
站点数据值输出
set stnprint on | off
ga-> set stat on ga-> d tem
—stat
统计信息 set gxout stat
ga-> set gxout stat ga-> d tem Data Type = grid Dimensions = 0 1 I Dimension = 36 to 109 Linear -2.5 5 J Dimension = 1 to 36 Linear -87.5 5 Sizes = 74 36 2664 Undef value = -32767 Undef count = 0 Valid count = 2664 Min, Max = -48.8 32 Cmin, cmax, cint = -40 30 10 Stats[sum,sumsqr,root(sumsqr),n]: 8574.9 1.14741e+06 1071.17 2664 Stats[(sum,sumsqr,root(sumsqr))/n]: 3.21881 430.708 20.7535 Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 3.22002 430.87 20.7574 Stats[(sigma,var)(n)]: 20.5024 420.348 Stats[(sigma,var)(n-1)]: 20.5062 420.506
改变当前缺省文件的命令: set dfile n
13
LASG动力论坛版
六. 图形显示(D)2
内容回顾
站点数据处理
–图样式及其定制 –技巧 站点客观分析函数 defiine 的使用
14
LASG动力论坛版
六. set gxout其它选项(D)2
set gxout 图形类型
bar line 直方图 折线图 contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 输出统计数据 二进制屏幕输出(set prnopts) 输出二进制文件(set fwrite) errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
7
LASG动力论坛版
六. 图形显示(D)2
天气符号填图 set gxout wxsym 绘图选项
set wxopt wxsym|mark|char
— 站点绘图2
站点位置分布图 set gxout stnmark 站号显示
set stid on | off
单要素填图 set gxout value 站号显示
9
LASG动力论坛版
六. 图形显示(D)2
(x,y)的最近站点 set gxout findstn d stnexpr; x; y 注意:
1. 维数只能X,Y变化 2. x, y 为虚页坐标
— 站点绘图3
风羽序列图 set gxout tserbarb d uexpr; vexpr 注意:
1. 只能时间维变化 2. 变量注明站号stid
六. 图形显示(D)2
站点填图 set gxout model
— 站点绘图1
数值显示方式
set mdlopts noblank|blank ~ ~ dig3| nodig3 set dignum number
'd u;v;t;d;slp;Δp;cloud;wx;vis' t vis wx d slp Δp
天气符号颜色
set wxcols #c1 #c2 …
字符大小
set digsize #inch
粗 细
set cthick #thick
颜 色
set ccolor #color
u,v,t,d,slp,Δp: 风,温度,露点,地面气压,变压 cloud: [1~9] 标记的类型, [20~25]云量: 20晴, 21少云, 22中云, 23多云, 24阴; 25缺测(M)
注意点
– 为避免插值出现虚假结果, 参考网格的间隔与站点 间距相近为好; – 具体情况具体分析, 并不一定Cressman函数更好!
12
LASG动力论坛版
六. 图形显示(D)2
变量定义技巧
— 站点绘图扩展
– define命令只适用于格点数据; – 在当前缺省文件为格点数据时, 可以用define 命令将站点数据的客观分析结果保存到变量 中.
在%后和格式字符d,o,x 等之间可插入以下附加符号:
l(L) 用于长整型整数,可加在d, o, x, u 前面, m (正整数m)数据最小的宽度, .n (正整数n)对实数输出n 位小数, 对字符串截取n 个字符, - 输出的数字或字符串在域内向左靠.
18
LASG动力论坛版
六. 图形显示(D)2
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
天气符号大小
set digsize #inch
标记颜色
set rbcols #col | auto
天气符号颜色
set wxcols #c1 #c2 …
标记类型
set cmark #mark
8
LASG动力论坛版
六. 图形显示(D)2
—天气符号
'set wxcols #color1 #color2 …… ' 'draw wxsym #symbol x y size #color #thick '
— 数学函数
函
数
备 注
绝对值 余弦/反余弦 正弦/反正弦 正切/反正切 指数 对数 幂函数
acos(expr) asin(expr) atan2(expr1, expr2) [expr1/expr2] log10(expr) expr1expr2 sqrt(u*u+v*v)
范例: 输出Binary格式数据 file='output.grd' 'set gxout fwrite' 'set fwrite 'file 'd tem' 'disable fwrite'
输出到binary文件 open数据文件 write数据 close数据文件
17
LASG动力论坛版
六. set gxout其它选项(D)2—print的格式
输出格式
set fwrite { | -be | -le } { | -sq | -st } filename
范例: 输出ASCII格式数据 file='output.txt' 'set gxout print' 'd tem' rc=write(file, result) rc=close(file)
数据格式
%d %o %x %u %c %s %f %e %g 以带符号的10进制形式输出整数(正整数无符号), 以8 进制无符号形式输出整数(不输出前导符0), 以16 进制无符号形式输出整数(不输出前导符0x), 以无符号10进制形式输出整数, 以字符形式输出,只输出一个字符, 输出字符串, 以小数形式输出单, 双精度数,隐含输出6 位小数, 以标准指数形式输出单, 双精度数,数字部分小数位数为6 位, 选用%f 或%e 格式中输出宽度较短的一种格式,不输出无意义的0.
3
LASG动力论坛版
六. 图形显示(D)2
本节内容提要
— 站点绘图
适用站点资料的内建图形样式及其定制 站点处理的局限性和技巧
变量的引用格式 – defiine – 函数的使用
–
set gxout其它选项*
– stat/ print/ fwrite
4
LASG动力论坛版
set gxout 图形类型
1
LASG动力论坛版
GrADS实用绘图
第四讲 站点图/函数使用
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
=Funny @ LASG动力论坛=
2
LASG动力论坛版
内容提要
请注意标记的地方
第 四 讲
一. 引言 二. GrADS绘图软件包安装 三. 运行GrADS绘图软件包 四. 图形显示(D)1— 格点绘图 五. 数据准备(Gr) 六. 图形显示(D)2— 站点绘图 七. 数据分析(A) — 函数的使用 八. 批处理和脚本语言(S) 九. 图形保存和转换
set stid on | off
标值颜色
set ccolor #color
站号显示
set stid on | off
站号颜色
set ccolor #col
标值大小
set digsize #inch
站号颜色
set ccolor #color
标记/站号大小
set digsize #inch
小数位数
set dignum #number
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
bar line 直方图 折线图 errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
六. 图形显示(D)2
— 站点图类型
contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 输出统计数据 二进制屏幕输出(set prnopts) 输出二进制文件(set fwrite)
天气现象序列图 set gxout tserwx 注意:
1. 只能时间维变化 2. 变量注明站号stid
天气பைடு நூலகம்号大小
set digsize #inch
返回信息
stid lat lon value
天气符号颜色
set wxcols #c1 #c2 …
10
LASG动力论坛版
六. 图形显示(D)2
— 站点绘图范例