grads画图坐标设置

合集下载

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画图坐标设置

grads画图坐标设置⼀直听到有⼈抱怨，GrADS的坐标轴怎么那么固定，要设置个起始间隔还不⽀持时间轴，要在经纬度后⾯加个⼩圆圈的度数怎么就那么难，为啥不能四周都弄上坐标轴？好了，你的这些抱怨将会随着你看到这个帖⼦⽽消失，下⾯给出这个脚本的⽤法介绍：注意：在你display变量前需要先设置如下两个命令'set xlpos -20''set ylpos -20'这样能够屏蔽GrADS⾃带的坐标轴，否则会出现坐标轴重叠[code=gs]请先display变量，再运⾏该脚本该脚本主要实现了⾃定义的坐标轴显⽰，您可以免费使⽤该脚本该脚本的⽤法规则如下：#axis -param1 val1 -param2 val2...其中，-param表⽰需要⾃定义的参数类型，val表⽰该参数的具体设定值参数和参数值必须成对出现，例如#axis -type b -position o -sinterval 2上⾯的语句表⽰绘制类型是底部的x轴，刻度位置是朝外，每两个⼤的刻度之间显⽰两个⼩刻度线#注意：-type是必选参数所有参数如下所⽰：-type ：表⽰绘制类型，参数值有：b/l/t/r 分别表⽰底部，左侧，上部，右侧-label：表⽰是否绘制数值标签，参数值有：on/off 分别表⽰绘制和不绘制-position：表⽰刻度的位置，参数值有：o/i/m 分别表⽰外侧，内侧，中间-start：表⽰刻度从该处开始绘制，参数值为该轴范围内的数字当该轴为时间轴时，表⽰开始绘制的时次（1,2,3...）-end：表⽰刻度的结束位置，参数值范围同start，时间轴时表⽰结束的时刻(1,2,3...) -interval：表⽰刻度的间隔，隔多少值绘制⼀个刻度和标签例如，当x轴表⽰经度110-150E，如果interval为10，默认情况下会标出110,120,130,140,150不给定该参数时，默认情况下，会绘制5个刻度，以此确定默认间隔-sinterval：表⽰没两个⼤刻度之间要绘制⼏个⼩的刻度，不给定该参数则不绘制-size：表⽰刻度的长度，单位是英⼨，默认为0.1-color：表⽰坐标轴和刻度标签的颜⾊，默认为1-lfont：表⽰坐标轴标签的字体，默认为当前环境下设置的字体-lsize：表⽰坐标轴标签的字体⼤⼩，默认为0.12-lthick：表⽰坐标轴标签字体的粗细，默认为0.3-langle：表⽰坐标轴标签旋转的⾓度，默认为不旋转-suffix：表⽰需要在每个坐标轴标签后⾯添加的⽂字或者符号后缀，⽐如⽤来添加度数等-asuffix：表⽰当坐标轴为经纬度轴时，是否⾃动添加E/W/N/S/EQ这种标记，默认为添加，参数值为 on/off ，分别表⽰⾃动添加和关闭⾃动添加-hoffset：表⽰坐标轴标签在⽔平⽅向上的偏移量，单位为英⼨，正负均可-voffset：表⽰坐标轴标签在垂直⽅向上的偏移量，单位为英⼨，正负均可-tformat：当所绘制坐标轴为时间轴时产⽣作⽤，表⽰要显⽰的时间格式参数值为y m d h 这四个字母的任意组合，分别表⽰年⽉⽇时如：-tformat my 则会将时间轴标签显⽰为JAN1951这种格式-v：当绘制图形是1-D图（set gxout line等）的时候是必选参数，表⽰当前display的是哪个变量[/code]来具体的看两个例⼦，这两个例⼦只是为了描述脚本的⽤法，所以不⼀定美观。

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快速参考手册

设置y轴标注的颜色<粗细<大小>>
x轴标注位置，offset：英寸；side=b,下|t,上
y轴标注位置，offset：英寸；side=r,右|l,左
数据中z维用对数坐标
设置坐标轴框的颜色<粗细>
y轴数据变化范围从vlo到vhi
x轴数据变化范围从vlo到vhi
把断点数据连上/不连上
在图上方居中写标题string
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息
fwrite
返回当前fwrite命领输出文件的名称
string s
给出经度－纬度坐标值 lon lat, 返回xy坐标值
pp2xy ppx ppy
给出页面坐标 ppx ppy，返回xy坐标值
set imprun script-name
在每个display命令前自动执行一遍script-name文件中的命令组
run file-name <params> file-name <params>
set zlog on|off
设置z坐标取对数坐标,或关闭此设制
set strmden value
设置流线密度[1-10],缺省为5
set stnopts <dig3> <nodig3>
以后3位数标数值
set mdlopts noblank|blank|dig3|nodig3

GrADS画图入门

设置 set gxout shaded
（缺省）
命令：
display express
ga-> d strp Contouring: 0 to 100 interval 10
设置：
set *** ###
维数参数设置等值线设置
Set set Set set lon 100 cint 15 120 指定等值线间隔 lat 103 25 指定等值线颜色 ccolor
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取（2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用（4）站点数据的准备以及站点图形的绘制
• 3、图形显示类型及绘图设置：
•
2、Grads数据准备：
（1）ncep资料的直接读取（2）mm5_to_Grads,wrf_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用（4）站点数据的准备以及站点图形的绘制
data open(3,file=filename.ctl,form='formatted',status='unknown') zlev/1000,925,850,700,500,400,300,250,200,150,100/ write(3,'(a)')'dset ^'//filename.dat tim=0.0 ；ilev=11；nflag=0 open(2,file=filename.dat,form='unformatted',access='direct',recl=73*4) write(3,'(a)')'dtype station' irecord=0 write(3,'(a)')'stnmap ^sounding.map' do j=1,irec write(3,'(a)')'undef -999.0' irecord=irecord+1 write(3,'(a)')'title sounding' write(2,rec=irecord)'abc'//sta(j)(1:5),lat(j),lon(j),tim,ilev,nflag, write(3,'(a)')'tdef 1 linear 02JAN1987 1DY ' + (zlev(i),z(j,i),t(j,i),td(j,i),d(j,i),f(j,i),i=1,11) write(3,'(a)')'vars 5' enddo write(3,'(a)')'zz 11 99 '//staa//' height1' ilev=0 write(3,'(a)')'tt 11 99 '//staa//' temper1' irecord=irecord+1 write(3,'(a)')'t_td 11 99 '//staa//' t_td1' write(2,rec=irecord)'abcdefgh',lat(1),lon(1),tim,ilev,nflag, write(3,'(a)')'uu 11 99 '//staa//' u wind1' + write(3,'(a)')'vv (zlev(i),i=1,11),(z(1,i),i=1,11),(t(1,i),i=1,11), 11 99 '//staa//' v wind1' + write(3,'(a)')'endvars‘ (td(1,i),i=1,11),(d(1,i),i=1,11),(f(1,i),i=1,11) close(2) close(3)

GrADS绘图基础

显示图形 ga-> display u gaga-> d u; v.2 gaga-> d u; v.2; mag(u,v.2) ga 保存图形 ga-> enable print gaa.gmf ga-> print gaga-> disable print ga 清屏 ga-> clear gaga-> c ga重置缺省环境, 重置缺省环境, 不关闭数据文件 ga-> reset ga重新初始化, 重新初始化, 关闭 reinit
Linux： Linux：在.cshrc 文件中添加: 文件中添加: # 执行文件路径：setenv PATH /usr/local/grads/bin:${PATH} 执行文件路径： # 字体和底图文件路径：setenv GADDIR /usr/local/grads/dat 字体和底图文件路径： # 脚本库的路径：setenv GASCRP /usr/local/grads/lib 脚本库的路径：
GrADS基础绘图 GrADS基础绘图
何晏春 20082008-2-23
提纲
1、grads功能概述 grads功能概述（1）、软件性能介绍:主要是与matlab对比）、软件性能介绍:主要是与matlab对比（2）、 grads在win/Linux下的下载、安装、运行： grads在win/Linux下的下载、安装、运行： linux/windows环境变量的设置，安装后各个目录文件的内容 linux/windows环境变量的设置，安装后各个目录文件的内容（3）、紧接着进行grads的简单操作：对nc资料的操作，画两张简单的图：contour与shaded叠加）、紧接着进行grads的简单操作：对nc资料的操作，画两张简单的图：contour与shaded叠加（2）、Grads网络资源）、Grads网络资源 2、Grads的图形显示类型及绘图设置 Grads的图形显示类型及绘图设置（1）、一维：直方、折现、散点、误差图.... ）、一维：直方、折现、散点、误差图.... （2）、二维：等值线、阴影图、风矢量、流线图、站点天气图（3）、图形设置：坐标、线型、页面、地图投影、地表状况、图形信息开关 3, Grads 函数 4、grads 批处理脚本 3、Grads数据准备 Grads数据准备（1）、ncep资料的直接读取）、ncep资料的直接读取（2）、二进制格点数据的生成、数据描述文件ctl作用）、二进制格点数据的生成、数据描述文件ctl作用介绍grads格点数据的存储结构，即维数环境，生成二进制的fortran程序段；介绍grads格点数据的存储结构，即维数环境，生成二进制的fortran程序段； ctl文件的作用，ctl文件的结构。 ctl文件的作用，ctl文件的结构。（3）、mm5out_to_Grads,wrfoutput_to_grads, ）、mm5out_to_Grads,wrfoutput_to_grads, （4）、站点数据的准备以及站点图形的绘制 4、grads函数的使用： grads函数的使用：（1）、数学函数 (2）、统计函数 (2）、统计函数（3）、网格处理函数 ....... 5、实例示范（1） precluster应用 precluster应用（2） pc演示，用grads的totourial。 pc演示，用grads的totourial。

grads常用命令

第一个命令：这样将得到一个动画显示的序列，现在输入：claerset lon -90set lat -90 90set lev 1000 100set t 1d td u该情况下我们设置y（纬度）和z（高度）变化，所以我们得到一个竖直切面图。

我们显示了两个变量，他们重叠在一起。

你可以让任意多的变量重叠在一起，只要你不使用clear命令就行了。

另一个例子，这种情况X和T是变化的（哈莫图）：cset lon -180 0set lat 40set lev 500set t 1 5d z现在已经知道怎样使用选定的部分数据绘图了，下面我们学习怎样操作数据。

设置维数环境为Z,Y变化：clearset lon -180 0set lat 0 90set lev 500set t 1现在假设我们想显示华氏温度代替开尔文温度。

可以这样转换：display (t-273.16)*9/5 32任何由标准的运算符,-,* 和/ 组成的表达式都是允许的，并且操作数可以包含常量，变量，或者函数。

一个包含函数的例子：d sqrt(u*u v*v)有一个函数用来计算风的级数。

d mag(u,v)另一个内建函数计算平均值：clear d ave(a,t=1,t=5)这种情况我们可以计算5天的平均。

我们也可以从数据中移除平均值（距平值）：d z-ave(z,t=1,t=5)也可以在x方向作平均并求距平：cleard z-ave(z,x=1,x=72)也可以做时间差分：cleard z(t=2)-z(t=1)完整规范的变量名是：name.file(dim |-|=va lue,…) 如果我们打开了两个文件，也许一个是模式输出，另一个是分析，我们应该区分用如下方法二者：display z.2-z.1另一个内置的函数通过有线差分计算水平涡度相关cleard hcurl(u,v)还有另外一个计算数值方向的质量积分：cleard vint(ps,q,275)这儿我们计算了可降水量（单位mm）现在来讨论控制图形输出的话题。

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绘图软件实用手册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绘图的各种设置

符串中反斜杠表示起新行。
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

也可以是gmf格式的

默认输出的图形格式为GrADS meta file,可用附带程序Gv.exe 或者Gv32.exe查看并拷贝到Word文档中。
30
Gv图形处理程序
31
麦卡托投影
Mercator‘s projection
等值线
32
麦卡托投影
Mercator‘s projection等值线+色斑
45
地理网格数据说明
从控制文件可以看到，地理网格资料其实是一个 71×41的数组资料。在GrADS中打开此控制文件（Grid.ctl）, 设定图形输出类型 set gxout grid 可看到所有的数据均由0和1组成。其中需要绘制等值线的区域由1填充，而其他的区域由0填充。

46
47
地理数据的进一步说明
28
常用的绘图控制语句
等值线线形、间隔…… set cstyle style set cint value set cmin(cmax) set black value1 value2 设置绘图区大小 set parea 1 10.5 1 8

29
图形的输出方式

图形输出语句:
printim d:\300.gif gif white x1024 y768
41
对以上流程的检验
需要初步检验一下成果，以确保下一步绘制等值线的资料成功。在GrADS命令窗口中打开相应站点资料控制文件名并显示变量，可看到数据标记在相应的站点上。
42
带有国界的地图：
Set mpdset mres
Set mpdset hires
43
特定的地理网格数据

此文件为创建插值网格而建立。一般而言，我们不需自己创建。对于特定的问题，需要特定的地理网格数据。

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
– 图形类型
• 等值线/ 填色/ 流线/ 风矢/ 站点/…
– 图形要素
• • • • • 标题/ 标值坐标和坐标轴线型/颜色/标记字符/天气符号地图底图 / 投影方式 / 分辨率

第4章 GrADS各类参数设置

1. 对于图形类型为contour起作用的设置
（1）set ccolor color - 设置等值线颜色。
其中，color为颜色号：0 - 黑， 1 - 白， 2 - 红，3 - 绿，4 - 兰， 5 - 青，6 - 洋红，7 - 黄，8 - 橘黄，15 - 灰。该设置在“clear”或“display”命令后即重新设定颜色。
（8）set clskip number – 设置间隔几根(number-1)等值线标示数值。
2. 对于图形类型为contour或 shaded起作用的设置
（1）set cint value - 设置等值线间隔，该设置在“clear”或“display”命令后即重新设置。（2）set clevs lev1 lev2 ...- 设置特定的等值线值，只画lev1，lev2，.......值所在的等值线，用于不等间隔绘图，该设置在“clear”或“display”命令后即重新设置。
（6）set cmax value –设置不画高于此value值的等值线。该设置在 “clear”或“display”命令后即重新设置。（7）set cblack val1 val2–设置不画介于val1和val2之间的等值线。该设置在“clear”或“display”命令后即重新设置。
（4）set csmooth on|off - 设置是否将网格值重新插值。如取“on”，则在绘等值线图前用三次插值将现网格值插到更精细网格上，重新设置才改变本次设置。
（5）set cmin value –设置不画低于此value值的等值线。该设置在 “clear”或“display”命令后即重新设置。
数据描述文件
dset d:\model.le.dat UNDEF -2.56E33 TITLE 5 Days of Sample Model Output XDEF 72 LINEAR 0.0 5.0 YDEF 46 LINEAR -90.0 4.0 ZDEF 7 LEVELS 1000 850 700 500 300 200 100 TDEF 5 LINEAR 02JAN1987 1DY vars 8 ps 0 99 Surface Pressure u 7 99 U Winds v 7 99 V Winds z 7 99 Geopotential Heights t 7 99 Temperature q 5 99 Specific Humidity ts 0 99 Surface Temperature p 0 99 Precipitation ENDVARS

GrADS快速参考手册

set gxout graphic-type
graphic-type可取以下参数：
bar
柱状图形模式
model
不规则分布测站模式
barb
WMO风标表示的风矢量图
scatter
不规则分布数据模式
contour
等值线图模式
shaded
分色图
errbar
误差图模式
stat
按WMO规范显示测站观测数据
fgrid
define
列出当前所有定义的变量名称
dims
列出当前打开的描述文件维数信息
gxinfo
列出图形系统的窗口信息
shades
填色图中每种颜色的序号和它所对应的数值范围
pos
等待用户在图形区单击鼠标，反回一对位置坐标值
time
返回时间维信息
fwrite
返回当前fwrite命领输出文件的名称
string s
set cstyle style
设置线型,style=0,1,2…
set cmark mark_index
设置标注型式
set cterp on|off
打开/关闭样条函数平滑开关
设置标注/不标注/强制标注/字串形式标注/自动标注等直线的值
set clopts thick<size>>
设置等直线标注的颜色<粗细<字大小>>
close file#
关闭数据描述文件,file#=1,2,3…
set dfile number
改变当前ctl文件为number指定描述文件,number=1,2,…
define var=expression
var=expression

Grads使用手册

gxtran
Displays metafiles
stnmap
Maps station data
wgrib
See /wgrib.html
GrADS 英文主页：/grads/
§1.2 基本概念
按磁盘文件记录格式，GrADS 软件包的文件分为如下几类：
gxps，gxpsc，gxpscw
将 GrADS 内 print 命令产生的图元文件分别转换为单色（灰度）、彩色以及反白（即白底）
彩色的 postscript 格式文件，用于硬拷贝输出。语法为：
Technical Report
GrADS究所
大气科学和地球流体力学数值模拟国家重点实验室
（ LASG ）

关于本手册
本手册是在 1995 年张洪编写的《GrADS 气象图形系统使用方法》及其它相关文档的基础上汇编而成。主要介绍了 GrADS 的基本用法、使用技巧、及与 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 内定，文件名随用户自定，
(2) .ctl 原始数据描述文件。该文件为纯 ASCII 文件，用以描述原始数据集的基本信息，包括数据集文件名、数据类型、数据结构、变量描述等等，详见§4。在 GrADS 环境中至少得首先打开（open）一个数据描述文件，以便后续的操作有数据对象。

GrADS实用绘图--图形显示

字符大小
set digsize #inch
特征值马赛克图 set gxout fgrid 指定特征值/颜色
set fgvals v1 #c1 v2 #c2 v3 #c3 …
格线
set gridln auto ~ ~ off ~ ~ #color
特征值/颜色
set clevs v1 v2 … set ccols #c1 #c2…
文件的字节数相同, GrADS通过描述文件(.ctl)正确识别.
UNIT=10的数据文件在实际空间上的结构
高度
经度
UNIT=20 的数据文件在实际空间上的结构
高度
纬度
纬度
经度
6
LASG动力论坛版
四. 图形显示(D)1
– 维数环境
— 格点绘图
• 网格坐标系 x, y, z, t ~ 如同数组下标 • 实际坐标系 lon, lat, lev, time • 缺省维数环境
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…
(李晓峰, 2005)
20
LASG动力论坛版
四. 图形显示(D)1
网格填值图 set gxout grid 小数位数
set dignum number
— 2-D单变量图2
马赛克图 set gxout grfill 颜色

GrADS绘图与编程(值得收藏)

系统的核心文件,由此进入GrADS环境
由GrADS.exe文件调用产生格点资料映射文件(*.idx)
看Grib码资料，转为TXT格式文件
显示*.gmf格式的图形文件产生台站资料的映射文件(*.map)
Grib码资料的解读程序
将*.gmf格式的图形文件转换为*.wmf 格式
10
启动和退出GrADS
100 continue 200 continue 1000 format(2i7) 2000 format(37f6.2) 3000 format(37f8.1) 4000 format(37f7.2) end
三、数据描述文件
Data.ctl 文件内容：
Dset data.grd Undef –9.99e+33 Title Ncep/Ncar reanalysis project Xdef 37 linear 60 2.5 Ydef 17 linear 0 2.5 Zdef 2 levels 850 200 Tdef 48 linear jan1982 1mo Vars 4 U 2 99 u wind (m/s) V 2 99 v wind(m/s) H 1 99 H500 T 1 99 Tsfc data endvars
11 8.5
11
所需预备知识
DOS和Windows基础 FORTRAN编程 WORD基础
学习要求
强调上机编程实践
12
第二章 GrADS 数据格式
• 目标 (1)能将文本格式数据转为GrADS格式
(2)能写出数据描述文件(*.ctl)
重要概念
一、维数环境的概念
1、含义：GrADS视每一个变量(VAR)场为一个四维数据集，即包括三维空间(x,y,z)和一维时间(t)。 2、作用：说明和指定随后的分析或图形操作时参加的原始数据集的维数范围。 3、定义方法： Set lat|lon|lev|time val1 <val2> 地球坐标 Set x|y|z|t val1 <val2> 格点坐标 13

气象绘图软件GrADS使用要点

●计算机与气象仪器　●收稿日期:2002-02-20气象绘图软件GrADS 使用要点高文良　陈忠明　闵文彬(四川省气象科学研究所　成都　610071)摘要:通过介绍GrADS 绘图软件,分析了其优点及难点,指出了GrADS 的使用要点。

关键词:GrADS ;绘图软件;要点中图分类号:TP30214文献标识码:B文章编号:1003-7187(2002)02-0049-021　引言在众多的绘图软件中,针对气象部门的绘图软件GrADS (Grid Analysis and Display System )是一款非常不错的免费共享软件。

它整个系统由C 语言编成,运行速度快,图形显示快速、美观大方,支持多种图形类型的输出。

GrADS 具有丰富的内部函数,可以对资料进行初步的计算和分析处理,而且操作十分简便,易于理解掌握。

它既可以处理格点资料,又可以处理站点资料。

它支持对GriB 码文件、NetCDF 文件的直接读取,并且对气象界最流行的各种数据格式都能正确分析和处理。

用户还可以利用GrADS 的描述语言(GrADS Script )很方便地进行编程、开发。

GrADS 推出后在气象界得到了广泛的使用。

2　基本概念在GrADS 中,每个原始资料数据集都需要一个数据描述文件文件与它相对应。

在数据描述文件中,包括了对原始资料数据文件的名字、数据格式、纬向格点数及间隔、经向格点数及间隔、空中层次数及间隔、时次及间隔、变量数等的说明。

数据描述文件从时间和空间上对数据集的维数作出了明确的规定。

这样,在GrADS 启动后将数据描述文件打开,就可以对整个资料数据集进行处理和分析。

所以,数据描述文件一定要准确,避免任何差错。

数据文件中原始数据的排放次序也有一定的规则。

其排放次序为首先是经度、纬度,然后是高度和物理变量,最后是时次变化。

[1]。

如果改变了数据的排放规则,必须在数据描述文件中进行特别说明。

因此可以将原始数据序列看成是一个五维(或少于五维)的数据场,将每个经向和纬向的数据面看成是一个记录,是整个数据集的一个二维切片。