GrADS实用绘图(基4)
GrADS绘图与编程(值得收藏)
二、文本格式数据
GrADS格式
现有如下资料： 名称： U850，V850，U200，V200，H500 和TSFC 范围： 60~150E，0~40N 分辨率：2.5×2.5 时间： 1982年1月~1985年12月的逐月资料。
如何将上述资料写成GrADS下的数据格式？
lon) y (或lat) z (或lev) t (或time)
从西向东的水平坐标 从南到北的水平坐标 从地面到高空的垂直坐标 时序坐标。
如 ： set lon 40 160 ；set lat 0 60 定义了水平变化范围 set lat 30；set lon 0 180 定义了沿30N的纬向变化范围 set time jan81 dec94 定义了从81年1月到94年12月的时段
z
(X,Y) (X,Y)
v
t
• T=1 • T=2 • T=3
U V H T
U200
分(z)
时(v)
5月20日(t)
秒(x,y)
18
FORTRAN源程序
Parameter(ii=37,jj=17) Real var(ii,jj) Open(1,file=„u850.dat‟) Open(2,file=„v850.dat‟) Open(3,file=„u200.dat‟) Open(4,file=„v850.dat‟) Open(5,file=„h500.dat‟) Open(6,file=„tsfc.dat‟)
7
简
介
• GrADS 是Grid Analysis and Display System的缩写,它由美国马里兰大学气象 系Brian E. Doty开发。利用该系统可实现 包括格点数据和站点数据的彩色气象图 形分析和显示。具有操作简便、功能丰 富、图形美观、显示快速的特点。在国 内外气象界得到广泛应用。 • 版本：V1.5, V1.6.8, V1.7.7，V1.7.9等
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
显示类型
(完整版)GrADS绘图软件使用手册
GrADS绘图软件实用手册2002年1月目录第一章GrADS绘图软件概述1.GrADS绘图软件简介2.Internet上的GrADS资源2.1GrADS在Internet上的主页2.2windows环境下GrADS资源3.GrADS绘图软件的安装（windows环境）3.1在windows环境下安装GrADS软件包X server 的安装第二章GrADS绘图模板1.GrADS示例演示1.1 启动GrADS1.2 退出GrADS1.3 示例演示GrADS命令的使用2.GrADS绘图模板3.GrADS模板的高级应用GrADS描述语言GrADS高级模板的应用第三章GrADS数据格式1.格点数据描述文件1.1 数据描述文件各项解释1.2 生成model.le.dat和model.le.ctl文件的程序代码片段2.站点数据的格式附录1.如何精确控制图形输出的尺寸—Landscape纸型2.台站资料的显示3.Linux环境下的安装第二章GrADS绘图软件概述1GrADS绘图软件简介The Grid Analysis and Display System(GrADS) 是一套应用广泛、使用方便的科学数据绘图软件包。
其主要特点：●GrADS属于自由软件，可以从Internet上免费获得。
●可运行于各种Windows 和Unix工作平台。
●GrADS可用于4D数据的分析。
既经度、纬度、层（气压层、高度层等）和时间/xyzt 4维。
数据可以是格点化的数据或离散点数据。
GrADS特别适用于气象类数据的分析。
但也完全可以用于更广泛类型的数据分析。
●GrADS有多种显示方式：等值线、流线、矢量图、风矢量图、站点填图、折线图、直方图等多种两维图形。
●可处理多种数据格式的数据。
GRIB、NetCDF、HDF-SDS等通用数据格式和系统自定义的一种二进制数据格式。
●采用命令行输入的方式交互式地显示图形。
并有多种命令对数据进行再加工。
GrADS绘图学习技巧与实例
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
GrADS实用绘图--站点图 函数使用
存放在 系统变量 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
• 注意点
♥♥♥
– 为避免插值出现虚假结果, 参考网格的间隔与站点 间距相近为好; – 具体情况具体分析, 并不一定Cressman函数更好!
12
LASG动力论坛版
六. 图形显示(D)2
• 变量定义技巧
— 站点绘图扩展
♥♥♥
– define命令只适用于格点数据; – 在当前缺省文件为格点数据时, 可以用define 命令将站点数填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
天气符号颜色
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)
第二讲GrADS绘图的各种设置
draw xlab string draw ylab string 分别在水平坐标轴或垂直坐标轴左侧写字符串
string，作为x轴或y轴的说明。
set annot color <thickness> 设置上述图注所用的颜色和线宽，缺省为白色，
线宽为6，该命令同时设置了坐标轴线以及刻 度的颜色和线宽，坐标刻度和标尺的线宽为图 注设置的线宽再减1。
例：draw title Geopotential Hight (October 2010)
set xyrev on 交换水平和垂oj 设置当前地图投影方式。
常用proj取值包括： latlon：缺省设置，用固定的投影角进行Lat/lon投影
； scaled：用不固定的投影角进行latlon投影，地图比例
失效； nps | sps：北半球|南半球极地投影； robinson：robinson投影，x：-180 180；y：-90 90 off：不画出地图，坐标轴也不代表lat/lon。
（缺省时即按数据描述文件的维数进行绘图。）
特别地， 当数据文件为全球多层次多时次资料时，例如全球
144*73*12的多时次资料，缺省时维数为x、y方向为 （0-360E，90S-90N），z方向和t为第一个点。
GrADS最擅长画2-D图，即全球数据的缺省维数。
维数说明分为两类：
set lon | lat | lev | time val1 <val2> set x | y | z | t val1 <val2> 两者对应于同一组数据，前者是地球坐标
set grid on|off|value|horizontal|vertical <style> <color>
GrADS第4章 各类参数设置
'set gxout line' 'set cmark 2' 'd u' 'set gxout errbar'
'set ccolor 2'
'd u-12.8;u+12.8'
'print' 'c' 'disable print' 'reinit' ;
scatter 散点图 'reinit'
set gxout scatter 'open d:\lon.ctl' 'open d:\lat.ctl‘ d v.1;v.2 'enable print D:\grads\data\1.gmf' 'set x 1 ';'set y 1';'set z 1 '; 'set t 1 24' 'set gxout scatter ' 'd v.1;v.2 ' 'print' 'c' 'disable print' 'reinit' ;
vector: 矢量箭头形式绘二维风场（缺省设置） set gxout vector d u;v （u,v在同一个数据文件中可这样） 或 d u.1;v.2 （u,v在两个数据文件中） stream: 流线形式绘制二维风场 set gxout stream d u;v barb: 风向杆形式绘制二维风场 set gxout barb d u;v
站点数据变量或者表达式
GrADS实用绘图--图形显示
• 改变维数环境
‘set ‘set ‘set ‘set x 1 10’ y 1 10’ z 6’ t 15’ ‘set lon 120 180’ ‘set lat -20 20’ ‘set lev 500’ ‘set time 00z01jan2007’
3
LASG动力论坛版
四. 图形显示(D)1
本节内容提要
• 基本概念
– 数据文件 – .ctl文件 – 维数环境
— 格点绘图
• 内建图形样式及其定制
– 1-D序列图 – 2-D平面图 – 定制
• 自由图形 • 动画*
4
LASG动力论坛版
四. 图形显示(D)1
基本概念
– 数据文件 – 数据描述文件*.ctl
• 绘图区范围设置
ga-> set parea 1 10.5 1 8
9
LASG动力论坛版
四. 图形显示(D)1
— 实现一页多图
• 将画布分成几个小的虚页
panels_demo.gs ~ 各虚页相对独立, 不能跨虚页作图写字符
• 在一个虚页上设置多个绘图区
pages.gs
• 设置绘图窗口尺寸
‘set xsize x y’
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
• 画布/ 虚页/ 绘图区大小
— 格点绘图
命令顺序 set display draw print
– 图形类型
• 等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
• • • • • 标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
GrADS绘图学习技巧与实例(阿木)
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！（笔者：阿木）1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。
软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推荐学习者安装2.0版本，选择默认安装路径就可以。
2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。
3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20) sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,nz,nt)real sstopen(15,file='sst.grd',form='binary') ！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt') ！新建txt文件do it=1,ntdo iz=1,nzread(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！read后只有文件号，数据是无格式的enddoenddodo it=1,ntdo iz=1,nzwrite(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddoclose(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2 ，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。
GrADS第4章各类参数设置.
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout grid du
set lon 70 130 set lat 10 30 set lev 1000 set t 1 set gxout fgrid set fgvals -1 4 Set fgvals -3 5 du
3. 维数环境的定义
（1）一种是地球坐标（world coordinate） 以经纬度为度量单位。
形式如下： set lon|lat|lev|time val1 <val2>
（2）一种是格点坐标（gridcoordinate） 以网格点数为度量单位。 形式如下：
set x | y | z | t val1 <val2>
1. 维数环境的概念
GrADS系统视每一个物理变量（VAR） 场为一个四维数据集（4D data set）， 即包括三维空间（x，y，z）和一维时间 （t）。也可固定其中的一维或几维以获 得实际的低于四维的数据子集。
2.设置的作用
GrADS 中设置维数用以说明或指定随 后的分析或图形操作时参加操作的原 始数据集的维数范围。该工作数据集 可以是整个原始数据场，也可以是原 始数据场的一部分。
或 set lev 1000 850
set t 1 3
或 set time 01jan1985 03jan1985
注意: grads只能显示两维或一维不固定的图！
当所有维数都固定时，得到的是一个单 值数据点
维数环境的设置一直保持到再次设定时 都有效。要查看当前维数情况，可以在 运行环境中输入命令： ga-> q dims
第4章 各类参数设置
GrADS精致绘图说讲解
GrADS 精致绘图说【转】GrADS 精致绘图说(一)分类:技术文摘2007.2.13 21:37 作者：步行者| 评论：0 | 阅读：2872[程序见GrADS Script Library]1. basemap.gs:basemap L | O | U <fill_color> <out_color> <hi/lo>在低分辨率海岸廓线范围内用颜色覆盖陆地／海洋。
适用于各种投影方式，需lpoly.asc, lpoly_hires.asc, lpoly_US.asc, opoly.asc, opoly_hires.asc文件。
其中：L(l)：覆盖陆地，O(o)：覆盖海洋，U(u)：覆盖20N-50N的墨西哥和加拿大领土（低分辨率，适用美国），fill_color：填充色号，缺省为15，out_color：廓线的颜色号，缺省为15，hi/lo：高分辨率('set mpdset hires'，仅对15N-53N, 130W-60W区域)/低分辨率。
2. cbar.gs、cbarn.gs、cbarc.gs、cbar_l.gs和cbar_line.gs:cbarn sf vert xmid ymidcbarc center_x center_y back_color绘制'set gxout shaded'图形的填色标尺。
sf：色标尺寸，1为全尺寸，0.5为半尺寸；vert：0为水平，1为垂直；xmid,ymid：色标中心点的位置。
cbar_l -x X -y Y -n number -t text -pcbar_line -x X -y Y -c color -m mark -l linestyle -t text -p加'set gxout line'的图例说明。
其中：-x,-y：图中x和y的位置，-n：线条的数目（最多可为10条），-t：文字说明的内容（最多10条，需双引号括起），-c：线和标记的颜色，-m：定义标记；-l：定义线型，-p：用户可在图中点击给定图例的放置位置。
GrADS基础绘图
命令顺序 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)
自由图形、符号
设置线属性-> 画线-> 画符号-> 画空心|实心方框-> 画多边形-> 画多边形-> 画天气符号-> 设置字体-> 写标题-> 设置字符串属性-> 设置字符串大小-> 写字符串->
‘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’
气象常用绘图软件grads
27
常用的绘图控制语句
设置数据维数 set lat 15 55 set lev 850 设置坐标 set xlopts color thickness size set xlopts 2 4 0.18 set zlog on set xlabs x1|x2|x3|x4……
X坐标设置
Z坐标取对数坐标
45
地理网格数据说明
从控制文件可以看到，地理网格资料其实是一个 71×41的数组资料。 在GrADS中打开此控制文件（Grid.ctl）, 设定图形输出类型 set gxout grid 可看到所有的数据均由0和1组成。其中需要绘制等 值线的区域由1填充，而其他的区域由0填充。
46
47
地理数据的进一步说明
8
Set mpdset yunn
9
中国政区图
Set mpdset cn
10
主要内容
气象中常用的绘图方法
Grads, surfer, matlab, mapinfo……
11
简介
GrADS 是Grid Analysis and Display System的缩写, 它由美国马里兰大学气象系Brian E. Doty开发。利用 该系统可实现包括格点数据和站点数据的彩色气象 图形分析和显示。具有操作简便、功能丰富、图形 美观、显示快速的特点。在国内外气象界得到广泛 应用。 版本：V1.5, V1.6.8, V1.7.7，V1.7.9等 目前使用得比较多的V1.8sl11
33
麦卡托投影
Mercator‘s projection风矢量
34
极地投影
等值线、色斑
35
地图投影设置
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 } )
范例: 输出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的格式
在%后和格式字符d、o、x 等之间可插入以下附加符号:
l(L) 用于长整型整数,可加在d, o, x, u 前面, m (正整数m)数据最小的宽度, .n (正整数n)对实数输出n 位小数, 对字符串截取n 个字符, － 输出的数字或字符串在域内向左靠.
18
LASG动力论坛版
六. 图形显示(D)2
天气符号大小
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 ’
GrADS画图入门
命令：query dims
ga-> q dims Default file number is: 1 X is varying Lon = 100 to 150 X = 1 to 101 Y is varying Lat = 0 to 40 Y = 1 to 81 Z is fixed Lev = 1000 Z = 1 T is fixed Time = 00Z19JUL2005 T = 1
命令：open filename2.ctl
ga-> open strike024_041720 Scanning description file: strike024_041720 Data file strike024_041720.dat is open as file 2
查询命令
命令：query files
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取 （2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用 （4）站点数据的准备以及站点图形的绘制
• 3、图形显示类型及绘图设置：
set cmax 90 不画超过90（不含）以上的等值线 set clevs 495 523 534 556 560 564 572 584 588 set ccols 1 2 3 4 5 当指定颜色少于等值线条数时，最后一种颜色适用于超出的所有等值线。
第一讲-GrADS绘图软件基础知识
GrADS软件的优点
1）软件的专业性强，它是专门为从事气象研究 的科研和业务人员开发的，因此特别适用于气象 数据资料的分析和显示；
2）软件适用于多种操作系统，如Windows7， Linux和工作站下的Unix等； 3）软件采用命令输入和图形显示交互式的操作 方式
GrADS操作界面由命令窗口和图形显示窗口组成， 在命令窗输入命令，则在图形窗出现相应图形，这种 操作方式简单直观，容易理解和使用。
GrADS主要功能
1）可以处理四维数据，包含经度、纬度、层次和 时间； 2）可以处理多种数据格式，包括二进制数据格式 ，GRIB格式，NetCDF格式等通用数据格式； 3）可以对数据进行再加工，如fwrite语句； 4）可以显示多种类型图形； 5）可以精细设置各种绘图要素； 6）可以按多种方式存储图形文件； 7）提供多种后期处理的软件工具。
授课方式
本课程一半学时为老师讲授，另一半学 时为学生上机实习，老师在一旁指导。
每周上课边讲授边上机实习，有利于学 生及时消化掌握老师讲授的知识，使教 学效果达到最好。
总评成绩
平时成绩：40%
包括考勤、课堂作业、课后作业等
期末成绩：60%
机考
课程内容
第一章 气象绘图软件GrADS基础知识 第二章 GrADS绘图的各种设置 第三章 变量和函数的定义及功能 第四章 站点数据格式及绘制站点图 第五章 气象数据的处理 第六章 其他气象常用绘图软件简介（选讲）
多时刻多个层次多个变量的数据存放次序
数据文件的读取和写入
Readuwnd.f90
数据描述文件（即.ctl文件）（例）：
DSET ua.dat
TITLE Upper Air Data
DTYPE grid
GrADS实用绘图--脚本编程
LASG动力论坛版
八. 批处理和脚本语言(S)
— 脚本编程关键运算范例 范例2－相应的.ctl文件
15
范例2－格式转换
‘sdfopen uwnd.2000.nc’ ‘set gxout fwrite’ ‘set fwrite uwnd.500.feb1.bin’ ‘set lev 500’ ‘set time 1feb2000’ ‘set x 1 144’ * 缺省维数是 1~145 ! ‘set y 37 73’ * 取北半球 ‘d uwnd’ ‘disable fwrite’
LASG动力论坛版
1
GrADS实用绘图
第五讲 脚本编程/图形保存
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
LASG动力论坛版
2
内容提要
请注意标记♥♥♥的地方
第 五 讲
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1 — 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2 — 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
下拉菜单
‘set dropmenu fc bc oc1 oc2 tfc tbc toc1 toc2 bfc bbc boc1 boc2 soc1 soc2 thick ’ ‘draw dropmenu #widget[0~64] xc yc width height text ’
对话框
‘set dialog #strcol #fcol #bkcol #linecol #thk { | numeric | n }’ ‘q dialog { | xc yc width height } { prompt_string | prompt_and_initial_string }’ say result
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
5
LASG动力论坛版
六. 图形显示(D)2
三类维数变化
x, y 变化 z变化 t 变化
— 缺省维数
命令顺序
set display 变量引用方式 draw var ( stid = 'station_id' ) print
6
LASG动力论坛版 SLP显示方式
set stnopts dig3| nodig3
本节内容回顾
站点数据处理
– 图样式及其定制 – 技巧 defiine 站点客观分析函数
统计和输出(set gxout其它选项)
– 统计信息 – ASCII和Binary数据输出(文件格式转换)
19
LASG动力论坛版
内容提要
请注意标记的地方
一. 引言 二. GrADS绘图软件包安装 三. 运行GrADS绘图软件包 四. 图形显示(D)1— 格点绘图 五. 数据准备(Gr) 六. 图形显示(D)2— 站点绘图 七. 数据分析(A) — 函数的使用 八. 批处理和脚本语言(S) 九. 图形保存和转换
天气符号颜色
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)
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
天气现象序列图 set gxout tserwx 注意:
1. 只能时间维变化 2. 变量注明站号stid
天气符号大小
set digsize #inch
返回信息
stid lat lon value
天气符号颜色
set wxcols #c1 #c2 …
10
LASG动力论坛版
六. 图形显示(D)2
— 站点绘图范例
注意点
– 为避免插值出现虚假结果, 参考网格的间隔与站点 间距相近为好; – 具体情况具体分析, 并不一定Cressman函数更好!
12
LASG动力论坛版
六. 图形显示(D)2
变量定义技巧
— 站点绘图扩展
– define命令只适用于格点数据; – 在当前缺省文件为格点数据时, 可以用define 命令将站点数据的客观分析结果保存到变量 中.
范例: 输出Binary格式数据 file='output.grd' 'set gxout fwrite' 'set fwrite 'file 'd tem' 'disen数据文件 write数据 close数据文件
17
LASG动力论坛版
六. set gxout其它选项(D)2—print的格式
20
LASG动力论坛版
七. 数据分析(A)
本节内容提要
数学函数 统计函数 网格处理函数 滤波平滑函数 差分函数 物理量计算函数 站点数据专用函数
— 函数的使用
21
LASG动力论坛版
七. 数据分析(A)
功 能
abs(expr) cos(expr) sin(expr) tan(expr) exp(expr) log(expr) pow(expr1, expr2)
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
数据格式
%d %o %x %u %c %s %f %e %g 以带符号的10进制形式输出整数(正整数无符号), 以8 进制无符号形式输出整数(不输出前导符0), 以16 进制无符号形式输出整数(不输出前导符0x), 以无符号10进制形式输出整数, 以字符形式输出,只输出一个字符, 输出字符串, 以小数形式输出单, 双精度数,隐含输出6 位小数, 以标准指数形式输出单, 双精度数,数字部分小数位数为6 位, 选用%f 或%e 格式中输出宽度较短的一种格式,不输出无意义的0.
11
LASG动力论坛版
六. 图形显示(D)2
客观分析函数的使用
– Cressman分析
— 站点绘图扩展
oacres(grid_expr, stn_expr, { radii | 10,7,4,2,1 } )
– 栅格平均
oabin(grid_expr, stn_expr, { | -f | -c } )
— 数学函数
函
数
备 注
绝对值 余弦/反余弦 正弦/反正弦 正切/反正切 指数 对数 幂函数
acos(expr) asin(expr) atan2(expr1, expr2) [expr1/expr2] log10(expr) expr1expr2 sqrt(u*u+v*v)
存放在 系统变量 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
天气符号大小
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 '
1
LASG动力论坛版
GrADS实用绘图
第四讲 站点图/函数使用
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
=Funny @ LASG动力论坛=
2
LASG动力论坛版
内容提要
请注意标记的地方
第 四 讲
一. 引言 二. GrADS绘图软件包安装 三. 运行GrADS绘图软件包 四. 图形显示(D)1— 格点绘图 五. 数据准备(Gr) 六. 图形显示(D)2— 站点绘图 七. 数据分析(A) — 函数的使用 八. 批处理和脚本语言(S) 九. 图形保存和转换
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
输出格式
set fwrite { | -be | -le } { | -sq | -st } filename
范例: 输出ASCII格式数据 file='output.txt' 'set gxout print' 'd tem' rc=write(file, result) rc=close(file)
改变当前缺省文件的命令: 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
bar line 直方图 折线图 errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
六. 图形显示(D)2
— 站点图类型
contour 等值线图 shaded 填色等值线图 vector 风矢量图 barb 风羽 stream 流线图 输出统计数据 二进制屏幕输出(set prnopts) 输出二进制文件(set fwrite)