GrADS基础绘图
绘图速度较慢,当数据文件较大时,读入很容易内存溢出 远程处理比较麻烦。
grads的优点; (1),最多能读入20个数据文件,数据读取迅速,出图速度快,但是要写描述文件和gs文 件比较麻烦。 (2),在画二维等值线图时,图形非常平滑,这点是我到现在还没大规模用matlab的主要 原因。 (3),它包含了很多气象上所需要的函数,在gs文件里可以直接调用,因此在用它时也省 了不少事。
Grads在win/Linux下的下载、安装、运行
程序:
Windows:grads-1.8sl11-win32e.exe Linux:grads-bin-1.9b4-linuxRH9.tar.gz grads-1.8sl11-linux.tar.gz data.tar.Z
环境变量:
Windows:
Panoply OpenDX (formerly IBM Data Explorer) NCL (NCAR Command Language) ncBrowse (netCDF File Browser) HDF (Hierarchical Data Format) interface GrADS (Grid Analysis and Display System) FERRET NetCDF Toolbox for MATLAB-5 IDV (Integrated Data Viewer) Origin surfer
命令顺序 grads –b|l|p|c open/run set display draw print quit/reinit
显示类型
bar 直方图 line 折线图 errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill 按网格填色(马赛克图) fgrid 按set fgvals指定值填网格颜色 grid 网格填值 findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值 contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 stat 输出统计数据 print 十进制屏幕输出(set prnopts) fwrite 输出二进制文件(set fwrite) tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
参考资料
1、GrADS Tutorial 2、GrADS Reference Card 3、LASG: GrADS Mannual 4、LASG: GrADS 实用绘图讲义,funny@动 力论坛 5、动力论坛GrADS中文讨论区
常见软件
Grads网络资源
GrADS主页和中文讨论区 – /graபைடு நூலகம்s/ 或ftp:// /grads (软件下载/ 联机文档/ email讨论区) – /cgi-bin/forum/forums.cgi?forum=3 (中文手册下载/ 中文讨论区/ 脚本和其它资源) GrADS实用文档 – LASG编《GrADS实用手册》
显示图形 ga-> display u ↵ ga-> d u; v.2 ↵ ga-> d u; v.2; mag(u,v.2) ↵ • 保存图形 ga-> enable print a.gmf ↵ ga-> print ↵ ga-> disable print ↵ ������ 清屏 ga-> clear ↵ ga-> c ↵ 重置缺省环境, 不关闭 数据文件 ga-> reset ↵ 重新初始化, 关闭 reinit ↵
GrADS基础绘图
何晏春 2008-2-23
提纲
1、grads功能概述 (1)、软件性能介绍:主要是与matlab对比 (2)、 grads在win/Linux下的下载、安装、运行: linux/windows环境变量的设置,安装后各个目录文件的内容 (3)、紧接着进行grads的简单操作:对nc资料的操作,画两张简单的图:contour与shaded叠加 (2)、Grads网络资源 2、Grads的图形显示类型及绘图设置 (1)、一维:直方、折现、散点、误差图.... (2)、二维:等值线、阴影图、风矢量、流线图、站点天气图 (3)、图形设置:坐标、线型、页面、地图投影、地表状况、图形信息开关 3, Grads 函数 4、grads 批处理脚本 3、Grads数据准备 (1)、ncep资料的直接读取 (2)、二进制格点数据的生成、数据描述文件ctl作用 介绍grads格点数据的存储结构,即维数环境,生成二进制的fortran程序段; ctl文件的作用,ctl文件的结构。 (3)、mm5out_to_Grads,wrfoutput_to_grads, (4)、站点数据的准备以及站点图形的绘制 4、grads函数的使用: (1)、数学函数 (2)、统计函数 (3)、网格处理函数 ....... 5、实例示范 (1) precluster应用
自由图形、符号
设置线属性-> 画线-> 画符号-> 画空心|实心方框-> 画多边形-> 画多边形-> 画天气符号-> 设置字体-> 写标题-> 设置字符串属性-> 设置字符串大小-> 写字符串->
‘set line #color #style #thick’ ‘draw line x1 y1 x2 y2’ ‘draw mark marktype x y size’ ‘draw { rec | recf } xlo ylo xhi yhi’ ‘draw ployf x1 y1 x2 y2 ……xn yn’ ‘draw mappoly lon1 lat1 lon2 lat2 … lon1 lat1’ ‘draw wxsym #symbol x y size #color #thick’ ‘set font [0~5]0’ ‘draw title title_string’ ‘set string #color { tl | tc| tr| l | c| r| bl | bc| br} ‘set strsiz width height’ ‘draw string x y string’
set annot #color(1) thick(6) set frame on | off | circle set xyrev on set zlog { on | off | swap | undefine } set coslat {on | off } set tlsupp { year | month } set grid {on| off| horizontal| vertical } linestyle #color set gridln {auto | off | #color } set xlint int set ylint set xaxis starrt end int set yaxis set vrange vlo vhi set vrange2 set xlpos offset(0) { b | t } set ylpos…{l|r} set xlopts #color(1) thick(4) size(0.12) set clopts #col(-1) thick(-1) size(0.09) set ylopts set ylevs set xlevs x1 x2 x3… set ylabs set xlabs b1|b2|b3|… set ylab set xflip { on | off } set yflip set xlab { on | off | aut draw ylab draw title title_string
运行GrADS
打开数据: ga->open *.ctl ga->sdfopen *.nc ga->xdfopen test.ddf 查询信息:query ga->query ctlinfo ga->q dims ga->q file 关闭数据文件 ga-> close 3 ga-> close 2 Ga-> close 1 设置维数环境 ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵ 26 LASG动力论坛版 三. 运行GrADS绘图软件包 • 设置维数环境 ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵
新建变量 GASDIR c:\PcGrADS\dat GASCRP c:\PcGrADS\lib 修改PATH变量, 增加路径 c:\PcGrADS\win32;c:\PcGrADS\lib;c:\PcGrADS\dat;
Linux: 在.cshrc 文件中添加: # 执行文件路径:setenv PATH /usr/local/grads/bin:${PATH} # 字体和底图文件路径:setenv GADDIR /usr/local/grads/dat # 脚本库的路径:setenv GASCRP /usr/local/grads/lib