GrADS实用绘图(基2)
- 格式：pdf
- 大小：2.67 MB
- 文档页数：42
下载文档原格式
合集下载
GrADS绘图与编程(值得收藏)
16
二、文本格式数据
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格式
现有如下资料： 名称： 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学习资料：第2章 数据处理
<litlee-endian>在PC, DEC等机器上生成 的数据 个人电脑
big-endian、little-endian 用于自动改变二进制位存放顺序
<template> 多个时间序列原始数据文件用一个数据描 述文件统一地描述描述这些原始数据时采用的选项， 这些数据文件的文件名形式由 dset 定义的形式命名 文件名，提示所含数据的时次。
例如：
正确的替换为： %y2 两位数年 %y4 4 位数年 %m1 1 或2 位数月 %m2 2 位数月（用 0补齐 1 位数） %mc 3 字符月份缩写 %d1 1 或2 位数天 %d2 2 位数天 %h1 1 或2 位数小时 %h2 2 位小时
设置X方向格点与经度的对应关系
XDEF number <LINEAR Start increment> 或 XDEF number <LEVELS value-list>
number（>=1）给定格点数，整形； LINEAR指明是等间隔分布格点， Start起点坐标，负数表示西经； increment网格间距。 LEVELS参数指明是不等间隔分布格点，列 出具体每个格点的坐标值（以空格分开）。
open(1,file='u.dat')
open(2,file='v.dat')
open(3,file='sst.dat')
! 打开目标文件
open(12,file='mhy.grd',form='binary')
! 把数据文件读入
do 100 it=1,nt do 50 iz=1,nz
read(1,*) ((u(i,j,iz,it),i=1,nx),j=1,ny) 50 continue
big-endian、little-endian 用于自动改变二进制位存放顺序
<template> 多个时间序列原始数据文件用一个数据描 述文件统一地描述描述这些原始数据时采用的选项， 这些数据文件的文件名形式由 dset 定义的形式命名 文件名，提示所含数据的时次。
例如：
正确的替换为： %y2 两位数年 %y4 4 位数年 %m1 1 或2 位数月 %m2 2 位数月（用 0补齐 1 位数） %mc 3 字符月份缩写 %d1 1 或2 位数天 %d2 2 位数天 %h1 1 或2 位数小时 %h2 2 位小时
设置X方向格点与经度的对应关系
XDEF number <LINEAR Start increment> 或 XDEF number <LEVELS value-list>
number（>=1）给定格点数，整形； LINEAR指明是等间隔分布格点， Start起点坐标，负数表示西经； increment网格间距。 LEVELS参数指明是不等间隔分布格点，列 出具体每个格点的坐标值（以空格分开）。
open(1,file='u.dat')
open(2,file='v.dat')
open(3,file='sst.dat')
! 打开目标文件
open(12,file='mhy.grd',form='binary')
! 把数据文件读入
do 100 it=1,nt do 50 iz=1,nz
read(1,*) ((u(i,j,iz,it),i=1,nx),j=1,ny) 50 continue
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使用手册
1
令。详见§6。 (4) .exe 系统命令文件，是 GrADS 系统在 DOS 环境下的各执行文件。如 grads.exe 为 GrADS
图形分析和显示命令；dos4gw.exe 为 DOS 的扩展环境；gxps.exe、gxpsc.exe、gxpscw.exe 都是图元文件转换为 postscript 文件的执行文件；gxtran.exe 是图元文件转换到显示器上 显示的执行文件；gx.exe 是将图元文件转换为各种不带 ps 解释器的打印机输出的执行文 件。详见§1.3。 (5) .gmf（.met） GrADS 系统图元输出文件，格式由 GrADS 内定，文件名随用户自定，
中国科学院大气物理研究所 大气科学和地球流体力学数值模拟国家重点实验室 National Key Laboratory of Numerical Modelling for Atmospheric Sciences and Geophysical Fluid Dynamics ( LASG ) Institute of Atmospheric Physics, Chinese Academy of Sciences
gradsnc
GrADS with readline, printim, LATS, read/write NetCDF
gradshdf GrADS with readline, printim, LATS, read/write HDF, read NetCDF
gribmap
"Maps" a GRIB data set to a GrADS descriptor file
§1.3 基本操作
新版本的 GrADS 交互环境中可以使用上下左右箭头键对曾经使用的命令进行调用和编 辑，但仅限于本次启动 GrADS 交互环境后所使用的命令。
(完整版)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绘图学习技巧与实例
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！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绘图的各种设置
符串中反斜杠表示起新行。
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>
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实用绘图--图形显示
x, y 包含所有的点 z, t 取第一个点 1.GrADS最擅长画2-D图 2. 全球数据的缺省维数♥♥♥
• 改变维数环境
‘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
– 图形类型
• 等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
• • • • • 标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
• 改变维数环境
‘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基础绘图
命令顺序 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中填色标尺cbarn用法
traj.gs 在水平面上绘制当前时刻向前/向后的轨迹线
格式转换
lats4d.gs 输出NetCDF, HDF-SDS 或GRIB 格式数据
动画/缩放
xanim.gs 动画控制
zoom.gs 图形放大
=========================细说cbarn的分割线==================
pinterp.gs 由eta/sigma/等压面等插值到指定等压面
zinterp.gs 由eta/sigma/等压面等插值到指定等高面
物理量计算等
mconv.gs 计算水汽辐合
plotskew.gs 绘制T-logP图
sweat_index.gs 用相对湿度, 温度和风计算凝结指数
使用cbarn.gs只能画出默认的填色标尺，如果想自定义图标效果则使用下面语句：
cbarn sf vert xmid ymid
sf：色标尺寸，1 为全尺寸，0.5为半尺寸；
vert：0为水平，1为垂直；
xmid,ymid：色标中心点的位置。
具体，简单的可省略写为 cbarn 1 1 或 cbarn 0.5 0
cbar_l.gs 折线图图例
cbar_line.填色
makebg.gs 生成具有地形纹理的底图
map.gs 预定义多个可选的地图范围
颜色定义
cmap.gs
define_colors.gs
插值
isen.gs 由等压面插值到特定等熵面
GrADS中填色标尺cbarn用法(2011-09-22 15:56:01)转载标签： 绘图笔记杂谈
===================GrADS自带脚本库==================
气象常用绘图软件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画图入门
ga-> q files File 1 : TC0801-StrkPro024 2008041714 Descriptor: strike024_041714.ctl Binary: strike024_041714.dat File 2 : TC0801-StrkPro024 2008041720 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat
命令：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 当指定颜色少于等值线条数时，最后一种颜色适用于超出的所有等值线。
命令：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)。
line1<line2: #c1 line1>line2: #c2
误差线图 set gxout errbar d line1;line2 颜 色
set ccolor #color
标记颜色
set ccolor #color
标记大小
set digsize #inch
线 宽
set cthick [1~10]
3
LASG动力论坛版
四. 图形显示(D)1
本节内容提要
基本概念
– 数据文件 – .ctl文件 – 维数环境
— 格点绘图
内建图形样式及其定制
– 1-D序列图 – 2-D平面图 – 定制
自由图形 动画*
4
LASG动力论坛版
四. 图形显示(D)1
基本概念
– 数据文件 – 数据描述文件*.ctl
8
LASG动力论坛版
四. 图形显示(D)1
"画布" — 横放 / 纵放
$ grads -l 或 $ grads –p
— 设置布局
虚拟页面尺寸(单位: 英寸!
ga-> set vpage 0 11 0 8.5 ga-> set vpage 0 8.5 0 11
1inch=2.54cm) grads –l grads –p 画布为Letter纸 大小, 不可超出! 请预留坐标和标 题的位置!
scatter
set xlint int set ylint int
18
LASG动力论坛版
四. 图形显示(D)1
等值线平滑
set cterp on/ off
on/ off/ linear
~ csmooth
— 2-D单变量图1
等值线图 set gxout contour ~ ~ shaded
等值线颜色
标记类型
set cmark #marker
线 型
set cstyle #style
17
LASG动力论坛版
四. 图形显示(D)1
errbar
— 1-D双变量图示例
set xaxis strt end int set yaxis strt end int set vrange vlo vhi set vrange2 vlo vhi
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件 在实际空间上的结构
高度
经度
UNIT=20 的数据文 件在实际 空间上的 结构
高度
纬 度
纬 度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
网格坐标系 x, y, z, t ~ 如同数组下标 实际坐标系 lon, lat, lev, time 缺省维数环境
绘图区范围设置
ga-> set parea 1 10.5 1 8
9
LASG动力论坛版
四. 图形显示(D)1
— 实现一页多图
将画布分成几个小的虚页
panels_demo.gs ~ 各虚页相对独立, 不能跨虚页作图写字符
在一个虚页上设置多个绘图区
pages.gs
设置绘图窗口尺寸
'set xsize x y'
四. 图形显示(D)1
一页多图范例
pages.gs
panels_demo.gs
(注意区别)
12
LASG动力论坛版
四. 图形显示(D)1
设置图形类型 1-D图 2-D图
— 设置图形类型
set gxout …
– 维数环境 x | y | z | t 只有1维变化 – 单站的探空曲线, 台风年频数序列… – 维数环境 x | y | z | t 中有2维变化 – 水平分布 / 垂直剖面 / 经向剖面 / 纬向剖面 / 任意剖面 时间演变 /…
'set csmooth on' 'set gxout shaded' 'set black 760 780' 'd var' 'set gxout contour' 'set cint 5' 'd var' 'set clevs 0' 'set ccolor 1' 'd lat'
— 2-D单变量图范例
set ccolor rainbow ~ ~ revrain ~ ~ #color
等值线标记方式
set clab on/ off/ forced/ string/ auto
等值线间隔
set cint int
特征线值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
线 宽
set cthick [1~10]
直方图 set gxout bar
外观属性
set baropts filled ~ ~ outline
颜 色
set ccolor #color
缺测点
set missconn on ~ ~ off
线 宽
set cthick [1~10]
折线图 set gxout line
标记类型
set cmark #mark
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
14
LASG动力论坛版
四. 图形显示(D)1
柱间距
set bargap percent
— 1-D单变量序列图
颜 色
set ccolor #color
基 线
set barbase top ~ ~ bottom ~ ~ value
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
画布/ 虚页/ 绘ຫໍສະໝຸດ 区大小— 格点绘图命令顺序 set display draw print
– 图形类型
等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
set mpdset lowres set mproj latlon set map auto set map 1 1 1 draw map
run cbarn.gs
set grads on set grads off
Set timelab on Set timelab off
11
LASG动力论坛版
errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
display u;v;mag(u,v)
风羽图 set gxout barb 风羽绘制方式
set hempref auto ~ ~ nhem ~ ~ shem
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
箭头大小比例
set arrowhead ±inch
流线图 set gxout stream 流线密度
— 格点绘图
气象要素场按照 一定格式存储 说明数据文件的内容
沿纬线方向的点数 (x方向维数) 沿经线方向的点数 (y方向维数) 沿垂直方向的点数 (z方向维数) 有多少个时间点 (t方向维数) 有多少气象要素/ 每个气象要素的层次
构成 "维数环境"
5
LASG动力论坛版
四. 图形显示(D)1
.ctl文件相当于翻译官——
— 格点绘图
i 经度 j 纬度 k 高度
Write(10,rec=1) ( ( ( temp (i, j, k), i=1,18), j=1,9), k=1,5) Write(20,rec=1) ( ( ( temp (i, j, k), i=1,5), j=1,9), k=1,18)
线 型
set cstyle #style
标记大小
set digsiz #inch
set xlopts #color thick size set ylopts #color thick size
LASG动力论坛版
set xlpos offset side set ylpos offset side
set xlabs b1|b2|b3|… set ylabs b1|b2|b3|…
16
LASG动力论坛版
四. 图形显示(D)1
— 1-D双变量序列图
散点图 set gxout scatter d var1;var2 双线填色图 set gxout linefill d line1;line2 颜 色 set lfcols #c1 #c2
15
四. 图形显示(D)
filled bar linefill
— 1-D序列图示例
line
line
outline bar linefill
'set t 3 9' 'define red=lat-lat' 'set t 1 12' … 'd red'
set xlevs x1 x2 x3… set ylevs y1 y2 y3…
21
LASG动力论坛版
四. 图形显示(D)1
误差线图 set gxout errbar d line1;line2 颜 色
set ccolor #color
标记颜色
set ccolor #color
标记大小
set digsize #inch
线 宽
set cthick [1~10]
3
LASG动力论坛版
四. 图形显示(D)1
本节内容提要
基本概念
– 数据文件 – .ctl文件 – 维数环境
— 格点绘图
内建图形样式及其定制
– 1-D序列图 – 2-D平面图 – 定制
自由图形 动画*
4
LASG动力论坛版
四. 图形显示(D)1
基本概念
– 数据文件 – 数据描述文件*.ctl
8
LASG动力论坛版
四. 图形显示(D)1
"画布" — 横放 / 纵放
$ grads -l 或 $ grads –p
— 设置布局
虚拟页面尺寸(单位: 英寸!
ga-> set vpage 0 11 0 8.5 ga-> set vpage 0 8.5 0 11
1inch=2.54cm) grads –l grads –p 画布为Letter纸 大小, 不可超出! 请预留坐标和标 题的位置!
scatter
set xlint int set ylint int
18
LASG动力论坛版
四. 图形显示(D)1
等值线平滑
set cterp on/ off
on/ off/ linear
~ csmooth
— 2-D单变量图1
等值线图 set gxout contour ~ ~ shaded
等值线颜色
标记类型
set cmark #marker
线 型
set cstyle #style
17
LASG动力论坛版
四. 图形显示(D)1
errbar
— 1-D双变量图示例
set xaxis strt end int set yaxis strt end int set vrange vlo vhi set vrange2 vlo vhi
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件 在实际空间上的结构
高度
经度
UNIT=20 的数据文 件在实际 空间上的 结构
高度
纬 度
纬 度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
网格坐标系 x, y, z, t ~ 如同数组下标 实际坐标系 lon, lat, lev, time 缺省维数环境
绘图区范围设置
ga-> set parea 1 10.5 1 8
9
LASG动力论坛版
四. 图形显示(D)1
— 实现一页多图
将画布分成几个小的虚页
panels_demo.gs ~ 各虚页相对独立, 不能跨虚页作图写字符
在一个虚页上设置多个绘图区
pages.gs
设置绘图窗口尺寸
'set xsize x y'
四. 图形显示(D)1
一页多图范例
pages.gs
panels_demo.gs
(注意区别)
12
LASG动力论坛版
四. 图形显示(D)1
设置图形类型 1-D图 2-D图
— 设置图形类型
set gxout …
– 维数环境 x | y | z | t 只有1维变化 – 单站的探空曲线, 台风年频数序列… – 维数环境 x | y | z | t 中有2维变化 – 水平分布 / 垂直剖面 / 经向剖面 / 纬向剖面 / 任意剖面 时间演变 /…
'set csmooth on' 'set gxout shaded' 'set black 760 780' 'd var' 'set gxout contour' 'set cint 5' 'd var' 'set clevs 0' 'set ccolor 1' 'd lat'
— 2-D单变量图范例
set ccolor rainbow ~ ~ revrain ~ ~ #color
等值线标记方式
set clab on/ off/ forced/ string/ auto
等值线间隔
set cint int
特征线值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
线 宽
set cthick [1~10]
直方图 set gxout bar
外观属性
set baropts filled ~ ~ outline
颜 色
set ccolor #color
缺测点
set missconn on ~ ~ off
线 宽
set cthick [1~10]
折线图 set gxout line
标记类型
set cmark #mark
tserbarb 单站的风羽时间序列 tserwx 单站天气现象时间序列 wxsym 在站点上标记天气符号(set wxopt)
14
LASG动力论坛版
四. 图形显示(D)1
柱间距
set bargap percent
— 1-D单变量序列图
颜 色
set ccolor #color
基 线
set barbase top ~ ~ bottom ~ ~ value
7
LASG动力论坛版
四. 图形显示(D)1
内建图形及其定制
– 布局/ 一页多图
画布/ 虚页/ 绘ຫໍສະໝຸດ 区大小— 格点绘图命令顺序 set display draw print
– 图形类型
等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
标题/ 标值 坐标和坐标轴 线型/颜色/标记 字符/天气符号 地图底图 / 投影方式 / 分辨率
set mpdset lowres set mproj latlon set map auto set map 1 1 1 draw map
run cbarn.gs
set grads on set grads off
Set timelab on Set timelab off
11
LASG动力论坛版
errbar 标误差线 linefill 两线之间填色 scatter 散点图 grfill fgrid grid
按网格填色(马赛克图) stat 按set fgvals指定值填网格颜色 print 网格填值 fwrite
findstn 寻找最近的站点 model 站点填图 stnmark 标记站点位置 value 标站点值
display u;v;mag(u,v)
风羽图 set gxout barb 风羽绘制方式
set hempref auto ~ ~ nhem ~ ~ shem
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
箭头大小比例
set arrowhead ±inch
流线图 set gxout stream 流线密度
— 格点绘图
气象要素场按照 一定格式存储 说明数据文件的内容
沿纬线方向的点数 (x方向维数) 沿经线方向的点数 (y方向维数) 沿垂直方向的点数 (z方向维数) 有多少个时间点 (t方向维数) 有多少气象要素/ 每个气象要素的层次
构成 "维数环境"
5
LASG动力论坛版
四. 图形显示(D)1
.ctl文件相当于翻译官——
— 格点绘图
i 经度 j 纬度 k 高度
Write(10,rec=1) ( ( ( temp (i, j, k), i=1,18), j=1,9), k=1,5) Write(20,rec=1) ( ( ( temp (i, j, k), i=1,5), j=1,9), k=1,18)
线 型
set cstyle #style
标记大小
set digsiz #inch
set xlopts #color thick size set ylopts #color thick size
LASG动力论坛版
set xlpos offset side set ylpos offset side
set xlabs b1|b2|b3|… set ylabs b1|b2|b3|…
16
LASG动力论坛版
四. 图形显示(D)1
— 1-D双变量序列图
散点图 set gxout scatter d var1;var2 双线填色图 set gxout linefill d line1;line2 颜 色 set lfcols #c1 #c2
15
四. 图形显示(D)
filled bar linefill
— 1-D序列图示例
line
line
outline bar linefill
'set t 3 9' 'define red=lat-lat' 'set t 1 12' … 'd red'
set xlevs x1 x2 x3… set ylevs y1 y2 y3…
21
LASG动力论坛版
四. 图形显示(D)1