GrADS实用绘图--站点图函数使用

合集下载

气象绘图软件GrADS使用要点

lat 10( 或为其他固定值 ) t 1 12( 或为其他变化区间 )
d aave( V ar, lon= 0, lo n= 10, lat = 10, lat = 20) 另外 , 在 G rA DS 中定义变量 ( define 命令 ) 时 , 该变量的维数就是当定义该变量时 Gr AD S 的当前维数。所以 , 定义了的变量 ( define 命令 ) 都是有维数的。如果超出了该变量的维数对该变量进行操作 , 将得出意外的结果 ( 或缺测值 ) 。 4 坐标转换在 G rA DS 中多处涉及到坐标的转换问题。在
d ave( V ar, lon= 100, lo n= 120) 又如 , 我们在实际中经常遇到要将一个小区域的面积平均的时间变化曲线画出。一般容易错误地将 X、 Y 方向的维数设为 : set lon 0 10( 或为其他变化区间 ) set lat 10 20( 或为其他变化区间 ) 正确的办法应将 X 、 Y 方向的维数设为固定 , 只有 T 方向的维数是变化的 , 如下 :
50
气象部门使用的绘图分析软件。它不仅能够实现图形的绘制 , 还能对资料进行分析和计算。由于 G rA DS 的后续版本不断推出 , 且有许多气象工作者参与编程和开发 , 所以在实际工作当中我们应当大力应用推广它。参考文献
[ 1] 郭品文 , 刘宣飞 G rAD S 气象绘图系统用户使用手册国家自然科学基金委员会地球科学部南京气象学院大气资料服务中心 [ 2] Brian Dogy , M ike Fi ori no, G rAD S user s manual, Cen t er f or 0cean -Land -A t mosphere St udies, Inst it ue for Gl obal En vironment an d S ociet y

GrADS实用绘图--基本操作

• 关闭数据文件
ga-> close 3 ↵ ga-> close 2 ↵ ga-> close 1 ↵
文件打开从1 开始依次编号, 文件关闭时先关闭序号最大的文件
26
LASG动力论坛版
三. 运行GrADS绘图软件包
• 设置维数环境
ga-> set x 1 73 ↵ ga-> set lat -20 60 ↵ ga-> set lev 850 ↵ ga-> set time 00z01Apr2007 ↵
♥♥♥
18
LASG动力论坛版
二. GrADS绘图软件包的安装
• Windows下的安装(2)
– 设置环境变量 ♥♥♥ “我的电脑”右键菜单属性高级环境变量
• 新建变量 GASDIR c:\PcGrADS\dat GASCRP c:\PcGrADS\lib GAUDFT c:\PcGrADS\udft\udft.list • 修改PATH变量, 增加路径 c:\PcGrADS\win32;c:\PcGrADS\lib;c:\PcGrADS\dat;
8
LASG动力论坛版
一. 引
• GrADS实用文档
言
– LASG编《GrADS实用手册》 ♥♥♥ 中文, 较全面, 可用来自学 – GrADS reference card 英文, 绘图命令和函数速查 – GrADS scripting language reference card 英文, 脚本语言语法速查 – GrADS Documentation Index 联机英文索引, 最全面
12
LASG动力论坛版
二. GrADS绘图软件包的安装
• 操作系统

GrADS_画图入门PPT课件

Descriptor: strike024_041714.ctl Binary: strike024_041714.dat File 2 : TC0801-StrkPro20 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat
命令：query file n
当指定颜色少于等值线条数时，最后一种颜色适用于超出的所有等值线。
等值线标注设置
set clab %.1f 等值线标记方式
set clab on（等值线标数值）/off（不标记）/forded（强制标记）/%gK（在标记后加上字符K）/%g%%（加上字符%）/%.2f(保留2位小数)/ %03.1f（整数部分保留3为，若不足3位，前部用0补齐）
GrADS 画图入门
Grid Analysis and Display System
国家气象中心 2008年12月7日
Outline
• 1、grads功能概述(不讲) • 2、Grads数据准备：
（1）ncep资料的直接读取（2）mm5out_to_Grads,wrfoutput_to_grads （3）二进制格点数据dat的生成和数据描述文件ctl作用（4）站点数据的准备以及站点图形的绘制
strp 0 0 strike probability strn 0 0 strike number
查询命令
命令：query ctlinfo
ga-> q ctlinfo dset strike024_041714.dat title TC0801-StrkPro14 undef 9999 xdef 101 linear 100 0.5 ydef 81 linear 0 0.5 zdef 2 levels 1000 500 tdef 1 linear 00Z19JUL2005 60mn vars 2 strp 0 0 strike probability strn 0 0 strike number

GRADS绘图实例教程

500mb高度场等直线图.gs文件** Draw the COR.COEF SST and nhc000*'reinit''enable print h9601.31''clear''open hs.ctl''set dfile 1''set vpage 0.5 10. 0.2 8.5''set lon 0. 360.''set lat -90. 90.''set mproj latlon''set mpdset mres''set poli off''set ylint 10.''set xlint 20.''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set clopts -1 -1 0.06''set clab forced''set cint 40.''set gxout contour''set cthick 4''set grads off''d rsst55''set string 3 c 5 0''set strsiz 0.14''draw string 4.5 7.5 1996.1.31 Global 500mb Geopotential Height Field' 'print'pull dummy.ctl文件DSET h9601.31aTITLE heightFORMAT yrevUNDEF -9999.00XDEF 145 LINEAR 0. 2.5YDEF 73 LINEAR -90. 2.5ZDEF 1 LINEAR 1 1TDEF 1 LINEAR jan1960 1moVARS 1RSST55 0 99 COR.COEFENDVARS图形降水异常场.gs文件** Draw Precipitation Anomalous Field*'reinit''enable print spad6754.met''clear''open prp.ctl''set dfile 1''set vpage 0.5 11.0 0.5 8.5''set lon -180. 180.''set lat -90. 90.''set mproj latlon''set mpdset mres''set poli off''set ylint 10.''set xlint 20.''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set gxout shaded''set cmin 0''set grads off''d ps''set cint 1.''set clopts -1 -1 0.05''set xlint 20.''set ylint 10.''set gxout contour''set ccolor rainbow''set clab forced''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''d ps''set string 3 c 4 0''set strsiz 0.14''draw string 5.2 8. The Summer Precipitation Difference of 54-63 and 67-76 (interval = 1.) ''print'pull dummy.ctl文件DSET spad54.67aTITLE PrecipitationUNDEF -9999.0XDEF 145 LINEAR -180. 2.5YDEF 73 LINEAR -90. 2.5ZDEF 1 LINEAR 1 1TDEF 1 LINEAR jan1960 1moVARS 1ps 0 99 COR.COEFENDVARS图形风速矢量场.gs文件'reinit''enable print uvd56.850''clear''open uw.ctl''open vw.ctl''set dfile 1''set grads off''set vpage 1. 8.5 0.2 5.8 ''set lat -90 90''set lon 0 360''set t 1''set mproj scaled''set xlint 30''set ylint 30''set gxout vector''set arrscl 0.5 8'*'set arrowhead 0.2'*'set clab forced''set grads off''set grid off''set ylopts 1 5 0.12''set xlopts 1 5 0.12''set map 1 1 1''d skip(-u.1,2);(-v.2)'*'d u.3;v.4''set string 4 tl 4''set strsiz 0.12 0.14''draw string 0.5 0.3 850mb ud*vd for 1958-63 and 1967-76' 'print'pull dummyuw.ctl文件DSET ud85.56aTITLE SURFACE PRESSURE AND WIND ANOMALYFORMAT yrevUNDEF -9996.00XDEF 145 LINEAR 0 2.5YDEF 73 LINEAR -90 2.5ZDEF 1 LINEAR 1 1TDEF 1 linear JAN00 1moVARS 1u 0 99 u windENDVARSvw.ctl文件DSET vd85.56aTITLE SURFACE PRESSURE AND WIND ANOMALYFORMAT yrevUNDEF -9996.00XDEF 145 LINEAR 0 2.5YDEF 73 LINEAR -90 2.5ZDEF 1 LINEAR 1 1TDEF 1 linear JAN00 1moVARS 1v 0 99 v windENDVARS图形剖面流线图.gs文件'reinit''enable print uw35.met''clear''open uu.ctl''open wu.ctl''set vpage 0.5 11. 0.5 8.5''set lat 10''set lon 40 240''set vrange 1000 100 ''set ylint -100''set zlog on''set gxout stream''set strmden 6''set z 1 12''set grads off''d uu.1;(-wu.2*30)''set strsiz 0.18 0.15''draw string 1.8 0.3 u*(-w*30) for 35N in July for 1980-1995 ''print'pull dummyuu.ctl文件DSET um.35aTITLE SURFACE PRESSURE AND WIND ANOMALYUNDEF -9996.00YDEF 1 LINEAR 10 1XDEF 81 LINEAR 40 2.5ZDEF 12 Levels 1000 925 850 700 600 500 400 300 250 200 150 100 TDEF 1 linear JAN00 1moVARS 1uu 12 99 v windENDVARSwu.ctl文件DSET wm.35aTITLE SURFACE PRESSURE AND WIND ANOMALYUNDEF -9996.00YDEF 1 LINEAR 10 1XDEF 81 LINEAR 40 2.5ZDEF 12 Levels 1000 925 850 700 600 500 400 300 250 200 150 100 TDEF 1 linear JAN00 1moVARS 1wu 12 99 w windENDVARS图形时间演变图lat.gs文件** Draw the Latitude-year section*'reinit''enable print tlatr20.met''clear''open lat.ctl''set dfile 1''set vpage 0.5 11. 0.5 8.5''set t 1 40''set lat -90 90''set lon 357.5''set ylint 5.''set xlint 5.''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set gxout shaded''set cmax 0.''set grads off''set xyrev on''set cint 0.5''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lon=0,lon=357.5)''set clopts -1 -1 0.05''set xlint 5.''set ylint 5.''set gxout contour''set clab forced''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 4''set cint 0.5''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lon=0,lon=357.5)''set string 3 c 4 0''set strsiz 0.18''set strsiz 0.10''draw string 6. 8.2 1958-96 summer 200mb Temperature Zonal Mean (0-357.5E,90S- 90N),interval=0.5 C''print'pull dummylat.Ctl文件DSET t20ar.40aTITLE COR.COEFFORMAT yrevUNDEF -327.68XDEF 144 LINEAR 0 2.5YDEF 73 LINEAR -90 2.5ZDEF 1 LEVELS 1000TDEF 40 LINEAR jul1958 1yrVARS 1ps 0 99 COR.COEFENDVARS图形** Draw the COR.COEF SST and nhc000 *'reinit''enable print tlonr20.met''clear''open lon.ctl''set dfile 1''set vpage 0.5 11. 0.5 8.5''set t 1 40''set lon 0 357.5''set lat -90''set ylint 10.''set xlint 5.''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set gxout shaded''set cmax 0.''set grads off''set xyrev on''set cint .3''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lat=-90,lat=90)''set clopts -1 -1 0.05''set xlint 5.''set ylint 10.''set gxout contour'*'set clab forced''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 4''set cint .3''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lat=-90,lat=90)''set string 3 c 4 0''set strsiz 0.13''draw string 6. 8.2 1958-96 200mb Temperature Meridional Mean (0-357.5E,90S- 90N),interval=0.3 C''print'pull dummylon.gs文件** Draw the COR.COEF SST and nhc000*'reinit''enable print tlonr20.met''clear''open lon.ctl''set dfile 1''set vpage 0.5 11. 0.5 8.5''set t 1 40''set lon 0 357.5''set lat -90''set ylint 10.''set xlint 5.''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set gxout shaded''set cmax 0.''set grads off''set xyrev on''set cint .3''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lat=-90,lat=90)''set clopts -1 -1 0.05''set xlint 5.''set ylint 10.''set gxout contour'*'set clab forced''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 4''set cint .3''set ccols 0 0 5 5 5 5 5 5 5 5 5''d ave(ps,lat=-90,lat=90)''set string 3 c 4 0''set strsiz 0.13''draw string 6. 8.2 1958-96 200mb Temperature Meridional Mean (0-357.5E,90S- 90N),interval=0.3 C''print'pull dummylon.ctl文件DSET t20ar.40aTITLE COR.COEFFORMAT yrevUNDEF -327.68XDEF 144 LINEAR 0 2.5YDEF 73 LINEAR -90 2.5ZDEF 1 LEVELS 1000TDEF 40 LINEAR jul1958 1yrVARS 1ps 0 99 COR.COEFENDVARS图形曲线图.gs文件'reinit''enable print sahan2.met''clear''open sah1.ctl''set dfile 1''set vpage 0.5 11. 0.5 8.5''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set gxout line''set t 1 14''set vrange 52.0 62.0''set xaxis 1984 1997 5''set cmark 2''set cstyle 2''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''set grads off''d pr4''set string 3 c 4 0''set strsiz 0.20''print'pull dummy.ctl文件DSET sahr2.01aTITLE COR.COEFUNDEF -9999.99XDEF 1 LINEAR 1 1YDEF 1 LINEAR 1 1ZDEF 1 LINEAR 1 1TDEF 14 LINEAR jan1900 1yr VARS 1pr4 0 99 COR.COEF ENDVARS图形一页多图.gs文件* Draw line*'Reinit’'enable print awi2.met''clear''open aam.ctl''open am1.ctl''open am2.ctl''open am3.ctl''set dfile 1''set vpage 0.5 10.5 0.5 8.0''set gxout line''set t 1 90''set vrange -110. 190.''set xaxis 1900 1989 10'*'set yaxis -140. 200. 50.''set cmark 0''set cstyle 1''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''d rsst55''set dfile 2''set vpage 0.5 10.5 0.5 8.0''set gxout line''set t 1 90''set vrange -110. 190.''set xaxis 1900 1989 10'*'set yaxis -140. 200. 50.''set cmark 0''set cstyle 2''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''d rsst55''set dfile 3''set vpage 0.5 10.5 0.5 8.0' 'set gxout line''set t 1 90''set vrange -110. 190.''set xaxis 1900 1989 10'*'set yaxis -130. 120. 100.' 'set cmark 0''set cstyle 3''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''d rsst55''set dfile 4''set vpage 0.5 10.5 0.5 8.0' 'set gxout line''set t 1 90''set vrange -110. 190.''set xaxis 1900 1989 10'*'set yaxis -130. 120. 100.' 'set cmark 0''set cstyle 5''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set csmooth on''set cthick 6''d rsst55''set string 3 c 4 0''set strsiz 0.17''draw string 6.0 8.0 Asia-Africa Summer Precipitation Anomalies Time Series ''set grads off''print'pull dummyaam.ctl文件DSET uam.39aTITLE COR.COEFUNDEF -9999.99XDEF 1 LINEAR 1 1YDEF 1 LINEAR 1 1ZDEF 1 LINEAR 1 1TDEF 39 LINEAR jan1900 1yrVARS 1RSST55 0 99 COR.COEFENDVARSam1.ctl文件DSET ssr.01aTITLE COR.COEFUNDEF -9999.99XDEF 1 LINEAR 1 1YDEF 1 LINEAR 1 1ZDEF 1 LINEAR 1 1TDEF 96 LINEAR jan1900 1yrVARS 1sar 0 99 COR.COEFENDVARSam2.ctl文件DSET csr.01aTITLE COR.COEFUNDEF -9999.99XDEF 1 LINEAR 1 1YDEF 1 LINEAR 1 1ZDEF 1 LINEAR 1 1TDEF 96 LINEAR jan1900 1yrVARS 1chr 0 99 COR.COEFENDVARSam3.ctl文件DSET wwi1.40aTITLE COR.COEFUNDEF -9999.99XDEF 1 LINEAR 1 1YDEF 1 LINEAR 1 1ZDEF 1 LINEAR 1 1TDEF 40 LINEAR jan1900 1yr VARS 1RSST55 0 99 COR.COEF ENDVARS。

(完整版)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实用绘图(基5)

绝对值, 对数三角函数反三角函数双曲函数反双曲函数幂函数平方根, e的幂最接近整数, 截去小数余数
LASG动力论坛版
8
res = functionname(args) ;
* res为函数值
math_abs(expr), math_log(expr), math_log10(expr) math_sin(radian), math_cos(radian), math_tan(radian) math_asin(expr), math_acos(expr), math_atan(expr) math_sinh(radian), math_cosh(radian), math_tanh(radian) math_asinh(expr), math_acosh(expr), math_atanh(expr) math_pow(expr, exponent) math_sqrt(expr), math_exp(exexpr) math_fmod(expr1,expr2), math_mod(expr) math_strlen(string) strlen(string) valnum(string) wrdpos(string,int) sublin(string,n) Subwrd(string,n) substr(string,start,length)
res = write(filename, record, { | append } ) res: 状态值 =0 OK =1 打开错误 =8 文件为read打开 res = close(filename) res: 状态值 =0 OK =1 文件未打开
LASG动力论坛版
11
八. 批处理和脚本语言(S) —控制语句

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

利用Grads画站点图（contour、shaded、grfill等）

利用Grads画站点图（contour、shaded、grfill等）1.利用Fortran程序将数据输出为grd格式。

参考程序如下：parameter(num=160) （站点数）character*8 sta(num) (站名名数组，可任意）dimension xlon(num),ylat(num),rc(num) (经、纬度数组及其上对应的数值) open(30,file='cor.grd',form='unformatted') (工作站)open(30,file='cor.grd',form='binary') (微机)do 100 lev=1,20 (共输出20个时间上的观测或结果)tim=0.0ilev=1nflag=1do 10 i=1,numsta(i)=char(i)write(30) sta(i),ylat(i),xlon(i),tim,ilev,nflagwrite(30) rc(i)10 continueilev=0write(30)sta(num),ylat(num),xlon(num),tim,ilev,nflag (每个时次的结束) 100continueclose(30)2.ctl文件(创建与步骤1中输出的站点数据配对的station.ctl文件)：dset cor.grdformat sequentialdtype stationstnmap cor.mapundef -9.99e33tdef 20 linear jun1958 1movars 1r 0 99 correlationendvars3.为插值函数准备格点数据grid.grd及对应的grid.ctldset grid.grdundef -99.0xdef 71 linear 70 1ydef 41 linear 15 1zdef 1 linear 1 1tdef 20 linear jun1958 1movars 1g 0 99 grid data prepared for oacres functionendvars注意：两个ctl（station.ct l和grid.ctl）文件中的时间要严格一致。

[转]grads中站点处理和oacresmaskout函数应用心得

[转]grads中站点处理和oacresmaskout函数应⽤⼼得SJ 发表于: 2011-1-05 13:22 来源:昨⼉刚刚帮同事画了两张降⽔图（不同时间段的），很是费劲，中间碰到若⼲问题，但好在最后基本达到了满意的效果，下⾯把这次画图中遇到的问题及解决办法贴出来，供⼤家参考或指正：1，数据处理的问题由于micaps的24⼩时降⽔资料中，并不是每个时次的⽂件⾥的站点个数都是相同的，就更不⽤提两个⽂件⾥有完全相同的站点了，这样就给各个站点降⽔的累加处理造成了极⼤的困难（如果⼀开始就有⼈或专门的机构来做这件事，将⼤⼤提⾼科研⼈员的⼯作效率），但本帖不打算过多描述数据处理的情形，因此这个问题就到此为⽌。

2，站点数据中打⽹格和maskout的联系降⽔的站点资料处理完毕并转成grads能识别的⼆进制格式后，就需要给站点数据打⽹格了，这个⽹格是⽤于站点资料差值使⽤的（如果⽤grads，就要⽤格点资料），此时就涉及到要打⼀个什么分辨率的⽹格了，在这次画图过程中，我⼀共打了两个全球⽹格，⼀个是0.1X0.1（细⽹格）的另⼀个是0.5X0.5（粗⽹格）的。

然后按部就班的使⽤了如下的语句：'set mpdset china' * 设定⼀个地图，⽤于只画中国范围内的物理量'define a=oacres(g,tt.2,XXXX)' * 插值函数，把变量tt插值到⽹格g上，其中XXXX代表影响半径'define a1=maskout(a,a-8)' * maskout（n1,n2）⾥第⼀个变量n1是我们要画的物理量，第⼆个变量n2的意思是只画当其值⼤于零时的物理量n1，* 当然n1和n2是⼀⼀对应的，⽽这个例⼦中的意思是，只画物理量a的值⼤于4的部分'define aa=smth9(a1)' * 平滑函数'run province-basemap china aa' * 只画中国范围内物理量下⾯通过格点的分辨率-2种，和oacres中的影响半径做组合实验：（1）粗⽹格+XXXX（50,30,20,10,5,1）见图⼀（2）细⽹格+XXXX（20,10,5,1）见图⼆（3）细⽹格+XXXX（50,30,20,10）见图三（4）粗⽹格+XXXX（50,30,20）50是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中的函数

GrADs中的函数GrADS的函数分两类,一类是对格点/站点数据执行运算的,这一类我们姑且称之为分析函数;另一类是脚本编程(gs)时使用的,这后一类我们姑且称之为脚本函数.第一类分析函数又分为格点分析和站点分析两个子类:1)格点分析绝对值abs(expr)余弦/反余弦cos(expr)acos(expr)正弦/反正弦sin(expr)asin(expr)正切/反正切tan(expr)atan2(expr1,expr2)[expr1/expr2]指数exp(expr)对数log(expr)log10(expr)幂函数pow(expr1,expr2)expr1expr2直角三角形的弦长/全风速mag(u_expr,v_expr)sqrt(u*u+v*v)平方根sqrt(expr)加权/不加权平均ave(expr,dim_expr1,dim_expr2,tincr,{|-b})mean()加权/不加权区域平均aave(expr,{xdim1,xdim2,ydim1,ydim2|global|g}) amean() 剔除值后时间求平均tmave(undefine_or_weight_expr,expr,texpr1,texpr2)加权/不加权区域和asum(expr,{xdim1,xdim2,ydim1,ydim2|global|g}) asumg() 质量加权垂直积分vint(ps_expr,expr,top)普通积分gint(expr)最大/小值Max(expr,dim1,dim2,tincr)min()最大/小值所在格点Maxloc(expr,dim1,dim2,tincr)minloc()空间相关/回归scorr(expr1,expr2,{xdim1,xdim2,ydim1,ydim2|global| g})sregr() 时间相关/回归tcorr(expr1,expr2,tdim1,tdim2)tregr()常数const(expr,const,{|-u|-a})剔除maskout(expr,mask_expr)跳点skip(expr,skip_x,skip_y)加权9点平滑smth9(expr)1-D3点平滑中央差分cdiff(expr,dim)时间循环tloop(expr)极值位置maxloc(expr,dim1,dim2,tincr)minloc()特征值位置findlvl(expr,expr_to_find,lev1,lev2)气象上物理量计算的函数:虚温,相对湿度求比湿tvrh2q(tv_expr,rh_expr)g/g虚温,相对湿度求温度tvrh2t(tv_expr,rh_expr)degK垂直涡度hcurl(u_expr,v_expr)中央差求水平散度hdivg(expr1,expr2)全风速mag(u_expr,v_expr)2)站点分析从格点插到站点gr2stn(grid_expr,{stn_expr|lon,lat},-a)收集生成格点场coll2gr(cnum,{num|-u})Cressman插值oacres(grid_expr,stn_expr,radii,1st_guess)网格平均oabin(grid_expr,stn_expr,{-f|-c})时间平均stnave(expr,dexpr1,dexpr2,-m count)时间序列最小值stnmin(expr,dexpr1,dexpr2,-m count)时间序列最大值stnmax(expr,dexpr1,dexpr2,-m count)转成1-D格点序列s2g1d(expr)所有数学函数abs,cos,acos,sin,asin,tan,atan2,exp,log,log10,pow,mag,sqrt 第二类,脚本函数包括有:abs,log,log10:math_abs,math_log,math_log10sin,cos,tan:math_sin,math_cos,math_tanasin,acos,atan2:math_asin,math_acos,math_atan math_sinh,math_cosh,math_tanhmath_asinh,math_acosh,math_atanhpow,sqrt:math_pow,math_sqrt,math_exp取整:math_nint,math_int求余:math_fmod,math_mod字符串处理函数:math_strlen,valnum,wrdpos,strlen sublin,subwrd,substr数值格式转换函数:math_format文件操作函数:read,write,close。

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实用绘图--站点图函数使用

3
LASG动力论坛版
六. 图形显示(D)2
本节内容提要
— 站点绘图
• 适用站点资料的内建图形样式及其定制 • 站点处理的局限性和技巧
变量的引用格式 – defiine – 函数的使用
–
• set gxout其它选项*
– stat/ print/ fwrite
4
LASG动力论坛版
set gxout 图形类型
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

几个有用的GRADS函数

几个有用的GRADS函数monmask多年时间序列中显示每年的相同月的资料<run> monmask.gs 5 define mmask for May (5)d (expr)*mmask plot (expr) only for Maycbar_line: 显示图例exp: cbar_line -x num(x方向位置) -y num(y方向位置) -c num(线条颜色) -l num(线条类型) -n num(线条数，max=10) -m num(标记) -t "(文本)" -p (加上-p，就不要-x, -y, 鼠标点到哪里就画在那)cbar: 对shaded图形有效，做图例cbarn: 同上，可加参数，exp: cbarn 1(图例粗细) 1(图例竖放，横放=0) <5(图例位置，x方向)> <5(图例位置，y方向)>cbarc: 同上，默认图例在右上角，可加参数，exp: cbarc <5(图例位置，x方向)> <5(图例位置，y方向)>connect_the_dots: 用鼠标画线，点击图的下方退出，无参数font: 设置字体，参数0-5，默认font=0, ROMAN字体为font=1loop: 循环作图，exp: loop x(or y，z，t，lat，lon，lev) min(起始值) max(最终值) step(步长)page: 选择画图区域，usage:run page top plot on top half of screenrun page bottom plot on bottom half of screenrun page left plot on left half of screenrun page right plot on right half of screenrun page q1 plot on quadrent 1run page q2 plot on quadrent 2run page q3 plot on quadrent 3run page q4 plot on quadrent 4run page xy XY (x,y,X,Y) are single digitsxy = location, XY = sizeonly 1 blank between xy and XYquadrents:-----------| q1 | q2 |-----------| q3 | q4 |-----------Matrix notation:xy is the location of the pageXY is number of subdivisionsThe idea here is that the page is divided into a grid of X by Yrectangles. Plotting is done in the x-y rectangle. In thefollowing example, the page is divided into two columns (X=2)and 3 rows (Y=2).-----------| 13 | 23 |-----------| 12 | 22 | xy values for XY=23-----------| 11 | 21 |-----------To draw into the top left rectangle, you "run page 13 23".To get the bottom right rectangle, you "run page 21 23".The following commands are equivalennt:"run page q1" "run page 12 22""run page q2" "run page 22 22""run page top" "run page 12 12""run page bottom" "run page 11 12"Limitations: X and Y must have values from 1 to 9.关掉多图形显示：page offsave: 保存图形，exp: save <路径>文件名.gmfwrdlin: 用经纬度画直线，exp: wrdlin minlon(起始经度) minlat(起始纬度) maxlon(最后经度) maxlat(最后纬度)wrdmrk: 在经纬度上标记，exp: wrdmrk num(标记类型) lon(经度) lat(纬度) siz(大小)wrdrec: 在图上画矩形，exp: wrdrec minlon(起始经度) minlat(起始纬度)[矩形左下角坐标] maxlon(最后经度) maxlat(最后纬度)[矩形右上角坐标]wrdrecf: 掩盖矩形内的图形，exp: wrdrec minlon(起始经度) minlat(起始纬度)[矩形左下角坐标] maxlon(最后经度) maxlat(最后纬度)[矩形右上角坐标]wrdstr: 在图上写字符串，exp: wrdstr lon(经度) lat(纬度) "(字符串)"map: 设置地图投影，nps (north-pole stereographic) lat='20 90' lon='-270 90'nps2 (north-pole stereographic) lat='20 90' lon='-180 180'nps3 (north-pole stereographic) lat='20 90' lon='-60 300'sps (sorth-pole stereographic) lat='-90 -20' lon='-270 90'lola (lat-lon) (0..360)'lola2 (lat-lon) (-180..180)'usa 'latlon' lat='24 52' lon='-127 -65'usa2 'nps' lat='15 80' lon='-150 -45'n_amer 'nps' lat='5 90' lon='-270 90's_amer 'latlon' lat='-60 20' lon='-90 -30'africa 'latlon' lat='-40 50' lon='-20 60'europe 'nps' lat='5 90' lon='-180 180'euro2 'nps' lat='5 90' lon='-180 180'asia 'latlon' lat='0 80' lon='40 170'aust 'latlon' lat='-50 0' lon='100 180'robinson 'robinson' lat='-90 90' lon='-180 180'c_pac 'latlon' lat='-45 45' lon='120 290'n_pac 'nps' lat='0 90' lon='90 270'custom [lon0 dlon lat0 dlat] (lat0,lon0=left-bottom corner)。

气象图形软件应用-GrADS绘图软件的使用

气象图形软件应用——G r A D S绘图软件的使用成都信息工程学院大气科学系巩远发主要内容一、GrADS的一般介绍二、GrADS的数据结构详解三、GrADS数据数据处理和绘图详解四、GrADS绘图编程（脚本）语言第一讲 GrADS的一般介绍什么是GrADS相关的关于GrADS的网络资源启动和退出GrADS最基本的GrADS操作（命令）GrADS使用文件的分类GrADS的主要功能介绍1. 什么是GrADS？G r A D S——Grid data Analysis and Display System 格点数据分析和显示系统GrADS（Grid data Analysis and Display System）是一个全32 位的交互操作的气象格点数据和站点数据的分析和显示环境。

系统具有气象数据分析功能强、地图投影坐标丰富、高级编程语言使用容易、图形显示快速，并具有彩色动画功能等特点。

目前是国内外气象数据显示的标准平台之一。

2. 相关的关于GrADS的网络资源WWW Sites/grads/head.html/home.htmlGrADS论坛动力论坛专业绘图软件格点资料下载http://www.ecmwf.int//cdc/reanalysis/3. 启动和退出GrADSgrads 启动GrADSga-> GrADS的系统提示符help给出GrADS的主要命令quit 退出GrADS4. 最基本的GrADS操作Open :打开一个网格点或站点数据文件供GrADS使用(sdfopen: 打开netCDF数据)Query (或仅用q) :查看当前GrADS环境的各种参数Display (或仅用d) :显示(绘制)一幅GrADS “表达式”的图形Set : 设置绘图的操作方式(包括绘什么图，在那里绘图，怎样绘图……)Clear (或仅用c) : 清除当前显示的图形5. GrADS使用文件的分类数据文件: 二进制无格式记录的原始数据，可以是格点的，也可以是站点的。

气象绘图软件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]。

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

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

GrADS实用绘图--站点图 函数使用

气象绘图软件GrADS使用要点

GrADS实用绘图--基本操作

GrADS_画图入门PPT课件

GRADS绘图实例教程

(完整版)GrADS绘图软件使用手册

GrADS实用绘图(基5)

GrADS实用绘图--图形显示

利用Grads画站点图（contour、shaded、grfill等）

[转]grads中站点处理和oacresmaskout函数应用心得

气象常用绘图软件grads

GrADs中的函数

Grads使用手册

GrADS实用绘图--站点图 函数使用

几个有用的GRADS函数

气象图形软件应用-GrADS绘图软件的使用

气象绘图软件GrADS使用要点

GrADS实用绘图--站点图函数使用

GrADS实用绘图--站点图函数使用