GrADS实用绘图--数据准备
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绘图实习一数据文件的转换及数据描述文件的建立
实习一数据文件的转换及数据描述文件的建立一、实习资料：现在“data”文件夹下有ASCII码数据资料文件：月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)、月平均高度场 (hgt500.dat)，2002.1-2005.12共48个月，资料的水平网格范围是M×N个网格点（M=37，N=17）,分辨率为2.5°×2.5°，自西向东经度60-150°E，由南至北纬度0-40°N。
二、实习要求：编写出将ASCII码数据资料文件：月平均风场和月平均高度场分别转换成二进制无格式直接存取（GrADS数据格式）文件（uv.grd和hgt.grd文件）的Fortran程序，给出相应的数据描述文件（uv.ctl和hgt.ctl文件），并利用GrADS基本操作命令显示该地区850hPa和200hpa uv风场，最后保存图形。
三、实习目的：通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。
四、实习步骤：4.1 利用“写字板”程序打开*.dat数据文件，熟悉该数据资料；4.2 按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为mhy.grd；4.3 通过“写字板”或“记事本”程序编写相应的数据描述文件，保存为mhy.ctl；4.4 画出2002年1月850hpa风场图，与图1比较验证数据正确性；4.5 利用GrADS基本操作命令（open,display,print等）显示2002年7月850hpa和200hpa uv 风场图形；4.5 完成实习报告1）说明所用资料资料名称：月平均风场(u200.dat、u850.dat、v200.dat、v850.dat)；月平均高度场(hgt500.dat)范围：60~150E，0~40N.层次：风场，850 hPa、200hPa；高度场，500hPa；时段：2002.1~2005.12共48个月。
《GrADS基础绘》课件
GrADS的主要应用领 域
GrADS被广泛应用于气象 学、气候学、海洋学、环 境科学等领域的数据处理 和可视化。
安装及环境配置
1
GrADS的安装步骤
详细介绍如何下载、安装和配置GrADS软件。
2
GrADS环境的配置
探索如何设置和调整GrADS的环境变量以适应个人需求。
3
GrADS的基本设置
了解GrADS的基本设置和首选项，以优化您的工作流程。
GrADS与其他软件集成
GrADS和Python的 数据交互
探索如何结合GrADS和 Python的强大功能，进行数 据交互和自动化处理。
GrADS与Matlab的 数据联动
了解如何在GrADS和Matlab 之间无缝传递数据，实现更 多分析和可视化功能。
GrADS在GIS领域的 应用
发现如何将GrADS与地理信 息系统（GIS）工具结合，实 现空间数据的分析和可视化。
总结
GrADS的优缺点评价
综合评价GrADS的优点和不足，帮助您了解 该工具是否适合自己的数据可视化需求。
GrADS使用心得交流
分享使用GrADS的一些心得体会和实用技巧， 促进经验交流和学习。
GrADS高阶绘图技能
GrADS高级绘图基础
学习如何创建精美且高度可定 制的数据可视化图表，彰显数 据的特点和趋势。
GrADS图像调整及美化
发掘GrADS图像调整和美化工 具，使您的图表更具吸引力和 专业性。
GrADS数据可视化技巧
了解如何创建交互式、动态和 多维数据可视化，提升数据沟 通的效果。
GrADS案例分析
1
GrADS数据处理应用案例
2
分享一些使用GrADS进行数据处理和
GrADS画图入门
（缺省）
命令：
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实用绘图--站点图 函数使用
存放在 系统变量 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绘图软件使用手册
GrADS绘图软件实用手册2002年1月目录第一章 GrADS绘图软件概述1.GrADS绘图软件简介2.Internet上的GrADS资源2.1GrADS在Internet上的主页2.2 windows环境下GrADS资源3.GrADS绘图软件的安装（windows环境）3.1在windows环境下安装GrADS软件包3.2X server 的安装第二章 GrADS绘图模板1.GrADS示例演示1.1 启动GrADS1.2 退出GrADS1.3 示例演示GrADS命令的使用2.GrADS绘图模板3.GrADS模板的高级应用3.1GrADS描述语言3.2GrADS高级模板的应用第三章 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) 是一套应用广泛、使用方便的科学数据绘图软件包。
其主要特点：z GrADS属于自由软件，可以从Internet上免费获得。
z可运行于各种Windows 和Unix工作平台。
z GrADS可用于4D数据的分析。
既经度、纬度、层（气压层、高度层等）和时间/xyzt 4维。
数据可以是格点化的数据或离散点数据。
GrADS特别适用于气象类数据的分析。
但也完全可以用于更广泛类型的数据分析。
z GrADS有多种显示方式：等值线、流线、矢量图、风矢量图、站点填图、折线图、直方图等多种两维图形。
z可处理多种数据格式的数据。
GRIB、NetCDF、HDF-SDS等通用数据格式和系统自定义的一种二进制数据格式。
z采用命令行输入的方式交互式地显示图形。
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 (Grid Analysis and Display System)概要内容提要一、GrADS绘图软件包的安装和运行二、GrADS绘图基本步骤三、GrADS绘图需要的三类型文件四、GrADS基本命令五、Fortran中如何读写二进制文件六、数据描述文件和gs文件简介一、GrADS绘图软件包的安装和运行1、演示安装2、启动GrADS1）双击grads图标2）绘画形式的选择风景画(横放)：11⨯8.5 英寸(缺省)肖像画形式(竖放)：肖像画形式是8.5 ⨯11 英寸3）GrADS的工作窗口分为两类：文本窗口：键入GrADS命令图形窗口。
显示所产生的图形两窗口之间可以切换。
二、GrADS绘图的基本步骤1、准备好了绘图资料（包括数据资料*.grd和数据描述文件*.ctl）2、启动GrADS3、打开数据描述文件（*.ctl）4、设置维数环境(set)5、生成图形文件6、查看图形文件(gv32.exe)例如：打开描述文件1)Open filename（盘符：路径/数据描述文件名）注意：启动GrADS后首先需打开至少一个数据描述文件2)显示变量图像D varname三、GrADS绘图需要的三类型文件.grd文件：这种文件格式为二进制无格式文件，可以是格点或站点资料，是从其他气象数据转化而来。
格点资料可以是直接存取或者顺序存取。
对格点资料而言，是一个五纬数据场，即时空、变量维（x/y/z/t/var）.ctl文件：数据描述文件，纯ASCII码文件。
是对数据信息的描述，包括数据名称、时空范围、变量说明等。
.gs文件：GrADS控制文件，也是一个纯文本文件。
由GrADS描述（脚本）语言（Scripts Language）写成的批处理GrADS系统设置和命令。
.gmf文件：GrADS系统图元输出文件，由gv.exe或者gv32.exe查看。
*.grd:数据文件，无格式直接存取文件。
*.ctl:数据描述文件*.gs:GrADS命令批处理文件四、基本操作命令命令浏览1. open 打开数据描述文件2. d 显示一个GrADS "表达式“3. set 设置在绘图中所需要的环境变量1）地球坐标(world coordinates)：set lat val1 <val2>set lon val1 <val2>set lev val1 <val2>set time val1 <val2>2）格点坐标(grid coordinates)：set x val1 <val2>set y val1 <val2>set z val1 <val2>set t val1 <val2>当键入一个值时，该维--"固定"；当键入两个值时，该维--"变化"。
第二讲GrADS绘图的各种设置
地图投影设置：
set mproj proj 设置当前地图投影方式。
常用proj取值包括： latlon：缺省设置，用固定的投影角进行Lat/lon投影
； scaled：用不固定的投影角进行latlon投影，地图比例
失效； nps | sps：北半球|南半球极地投影； robinson：robinson投影，x：-180 180；y：-90 90 off：不画出地图，坐标轴也不代表lat/lon。
对图形为矢量图起作用的 设置：
set arrscl size <magnitude> 设置矢量箭头的长度为size（
英寸）， 选项magnitude的值为矢量的
大小（例10m/s）。
set arrowhead size 设置箭头大小，缺省为0.05，
取0不画头，取负与矢量值 成比例。
例：set arrowhead 0.06
，后者是网格坐标。 例：set lon -180 0
set lat 0 90 set lev 500
set t 1 换成set time 00:00z01Jan2001
图形类型设置：
当维数环境确定后，缺省情况下，一维变量输出的图形为单线 图，二维变量为等值线图，若想进行特殊指定，键入命令：
set gxout graphics-type（图形类型） 对格点数据： contour: 二维数据等值线图（缺省） shaded: 二维填色等值线图 grid: 二维场不绘图，以网格形式在各网格点中央标出该场点
数值 vector: 矢量箭头形式绘二维风场 stream: 流线形式绘二维风场 barb: 风向杆形式绘二维风场 bar: 对一维场不绘单线图，而绘直方图 line: 对一维场绘单线图（缺省） 例：set gxout shaded
GrADS实用绘图(基5)
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实用绘图--图形显示
• 改变维数环境
‘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基础绘图
命令顺序 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实用绘图(基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绘图软件使用手3
GrADS绘图软件使用手3第三章GrADS数据格式每一组GrADS数据应至少包括两组数据文件，数据描述文件—ASCII 码和数据文件—二进制，数据的真正存放地。
数据文件中只是用户数据的有序排放，而关于数据种类、排放次序等是单独放在一个文件中的称—数据描述文件。
而象GRIB和NETCDF等通用数据格式，以上两者是存于同一个文件的—或称为自定义/自解释格式数据。
但考虑到GrADS传统，对这类自定义格式数据仍将生成相应的数据描述文件。
上一章中我们已使用过了这样的一组数据。
以此为例，介绍用户如何按GrADS的格式，将自己的数据生成相应的数据文件和数据描述文件。
1.格点数据描述文件model.le.ctl文件清单：以某开始的行为注解行。
1.1数据描述文件各项解释1.DSET数据文件名定义与此数据描述文件相对应的数据文件名。
若两者位于同一目录，前面的路经可以省略或以“^”开始，代表两者位于同一目录。
若不在同一目录下，应给出路经参数。
如：c:/pcgrad/ample/model.le.dat56注意路经的给法与DOS不同，而与UNI某环境一致，便于移植！或c:\\pcgrad\\ample\\model.le.dat,两种都行。
2.TITLE数据文件说明文字串。
3.UNDEFvaule定义缺测值。
一般给一很大的正/负值，表示，当取值超过这一正值/低于定义的负值，认为该值无效。
（GrADS采用跳过或用周围有效点的值处理。
）4.OPTIONS这里定义了与二进制存储有关的选项，二进制存储的一大特点是可移植性差，因此通过keyword项来增加可移植性。
若keyword省略，则OPTIONS也可省略。
可取：equential:顺序无格式方式。
yrev:Y维与YDEF定义相反方式存放。
zrev:Z维与ZDEF定义相反方式存放。
big_endian:如数据是在un,gi,hpcray机器上生成的，而目前不在此类机器上使用。
GrADS实用绘图--基本操作
17
LASG动力论坛版
二. GrADS绘图软件包的安装
• Windows下的安装(1)
– 安装X server (仅X-win32版本需要)
• MicroImages X server (轻巧,但图形窗内容被覆盖后无法刷新恢复) • X manager • Cygwin-X (软件“冲突”的解决) (LASG动力论坛精华版)
21
LASG动力论坛版
二. GrADS绘图软件包的安装
• 本节内容回顾
– 版本选择
• 仿UNIX的图形环境 ♥♥♥
– 安装路径 – 环境变量设置 – 备份和搬家
♥♥♥ ♥♥♥
• 本节推荐的是以方便使用(个人观点)为原则的设置
22
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方 第 一 讲
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1 — 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2 — 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
20
LASG动力论坛版
二. GrADS绘图软件包的安装
• 关于地图底图
– GrADS自带的底图mres和hires在我国西南国界与 实际不符!! (请以实际行动爱国哦)
– 新增的底图文件下载
/cgi-bin/forum/forums.cgi?forum=3
19
LASG动力论坛版
二. GrADS绘图软件包的安装
• 备份
– dat\ 目录的内容 – lib\ 目录的内容 – udft\ 目录的内容
最具扩展性的内容都保存在这几个目录下, 这也是分目 录存放不同文件的好处.
气象图形软件应用-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绘图与编程(值得收藏)
由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绘图学习技巧与实例
以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。
不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！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，大气的数据要根据数据的层次确定几层。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
zdef 7 levels 1000 850 700 500 300 200 100
• 文件路径推荐用UNIX惯例写
– d:/mydata/test.grd
10
LASG动力论坛版
五. 数据准备(Gr)
– – – – – – – – – yrev zrev sequential 365_day_calendar template byteswapped big_endian little_endian cray_32bit_ieee
17
LASG动力论坛版
五. 数据准备(Gr)
%ix1 %ix3 %iy2 %iy4 %im1 %im2 %imc %id1 %d1 %id2 %ih1 %ih2 %ih3 %in2
— .ctl文件语法9-4
tdef 中起始时间的替代方式
起始的1 位数年代(v1.9) 起始的3 维数年代 (v1.9) 起始的2 位数年 起始的4 位数年 起始的1/2 位数月 起始的2 位数月(必要时前置0) 起始的3 字符月 起始的1/2 位数天(必要时前置0) 1 位数字天 起始的2 位数天 起始的1/2 位数时 起始的2 位数时 起始的3 位数时 起始的2 位数分(必要时前置0)
— .ctl文件语法7-1
pdef关键字
－lcc | lccr －nps | sps －BILIN | file －eta.u | pse －ops
－特殊投影坐标系下的数据 Lambert投影 极射赤面投影 其它直线网格投影 US NMC eta模式 | SSM/I极射赤面投影 CSU RAMS斜极射赤面投影
3
LASG动力论坛版
五. 数据准备(Gr)
本节内容提要
– .ctl文件语法
– – – – – –
GrADS支持的数据格式 Binary数据格式 netCDF数据格式 GRIB数据格式 站点数据格式 进阶向导
4
LASG动力论坛版
五. 数据准备(Gr)
• 3-D空间数据排列方式
– 自西向东 – 从南到北 – 从低层到高层
15
LASG动力论坛版
五. 数据准备(Gr)
文 件 名 中 时 间 格 式 说 明 符
%x1 %x3 %y2 %y4 %m1 %m2 %mc %d1 %d2 %h1 %h2 %h3 %f2 %f3 %n2
— .ctl文件语法9-2
(如: 1910, 1920, 1930) (如: 1980, 1990, 2000) (如: 1, 2, …, 11,12) (如: 01, 03, …, 12) (如: Jan, Feb, Mar, …) (如: 1, 2, …, 31) (如: 01, 02, 10, 31) (如: 1z, 2z, …, 12z) (如: 01, 02, 10, 31) (如: 120, 012) (如: 12, 24, 48, 72, 120) (如: 012, 024, 072, 120) (如: 01, 35, 60)
HDF4(SDS)
20
LASG动力论坛版
五. 数据准备(Gr)
Binary格式数据
— Binary格式1 ♥♥♥
– GrADS最基本, 我们最容易准备 – 常用后缀 .dat .grd .bin ,…… – 需要 .ctl 文件一起使用 – 是无格式文件 (form=‘unformatted’) – 若用Fortran 读写
– 方向: 自西向东, 自南向北 (永远不变) – 单位: 度 纬度/经度 – 正负值: 以东为正/ 以北为正
7
LASG动力论坛版
五. 数据准备(Gr)
• Z 方向的维数说明
— .ctl文件语法3
zdef 7 levels 1000 850 700 500 300 200 100 列举法 给出 坐标值
– 方向: 自下而上 – 单位: hPa / 其他
可用linear 关键字替换
8
LASG动力论坛版
五. 数据准备(Gr)
• T 方向的维数说明
— .ctl文件语法4
1mo 时间 步长 3mn 6hr 1dy 1mo 1yr 分钟 小时 天 月 年
– tdef 12 linear 00:00z21jan2000
6
LASG动力论坛版
五. 数据准备(Gr)
• X, Y 方向的维数说明
– xdef
144 linear 坐标值 线性 增加
— .ctl文件语法2
0 起始 坐标 值 2.5 坐标 间隔 可用levels 关键字替换
Y: gausT62 gausR15 gausR20 gausR40
X(lon) 共有 方向 144个 维数 格点 定义
12
LASG动力论坛版
五. 数据准备(Gr)
pdef关键字
－使用注意
— .ctl文件语法7-2
－u,v分量需与GRIB数据格式中风场变量33和34的说明格式 相同, 以便执行矢量的旋转运算
u 18 V 18 u10 0 v10 0 33,100 34,100 33,105 34,105 U-Wind Components on Pressure Levels V-Wind Components on Pressure Levels 10 Meter U Wind 10 Meter V Wind
— .ctl文件语法6
• options 是可选的关键字, 其选项有
自北向南的数据 自高层到低层的数据 顺序存贮的数据 任何年份均是365天 多个文件 文件存储的字节次序与本机相反 在SUN, SGI等机器上生成的数据 在PC, DEC等机器上生成的数据 在Cray机上生成的数据
11
LASG 数据准备(Gr) 五. 动力论坛版
— .ctl文件语法1
• 格点数据描述文件 .ctl 的写法
dset d:/mydata/test.grd title This is a GrADS test data options sequential yrev undef -9999.0 xdef 144 linear 0 2.5 ydef 73 linear -90 2.5 zdef 7 levels 1000 850 700 500 300 200 100 tdef 12 linear 00:00z01jan2000 1mo vars 5 slp 0 99 Sea Surface Temperature (hPa) t2m 0 99 2m temperature( degC) temp 7 99 upper level temperature (degC) uwnd 7 99 upper level U component of wind (m/s) vwnd 7 99 upper level V component of wind (m/s) endvars
1位数字表示的年代 3位数字表示的年代 2位数字表示的年 4位数字表示的年 1位或2位数字表示的月 2位数字表示的月 3位字符表示的月 1位或2位数字表示的日 2位数字表示的日 1位或2位数字表示的小时 2位数字表示的小时 3位数字表示的小时 2位或3位数字表示的预报时效 3位数字表示的预报时效 2位数字表示的分钟
dset ^pr.%ch.nc chsub 1 600 1851-1900 chsub 601 1800 1901-2000 options template dtype netcdf tdef 1800 linear jan1851 1mo
若两个文件分别位于不同的磁盘上:
dset %ch chsub 1 600 /disk1/pr.1851-1900.nc chsub 601 1800 /disk2/pr.1901-2000.nc
1
LASG动力论坛版
GrADS实用绘图
第三讲 数据准备
＝Funny @ LASG动力论坛＝
— 2007年8月22日—
版权所有, 欢迎传播, 请保持文档完整性
2
LASG动力论坛版
内容提要
请注意标记♥♥♥的地方
一． 引言 二． GrADS绘图软件包安装 三． 运行GrADS绘图软件包 四． 图形显示(D)1 — 格点绘图 五． 数据准备(Gr) 六． 图形显示(D)2 — 站点绘图 七． 数据分析(A) — 函数的使用 八． 批处理和脚本语言(S) 九． 图形保存和转换
16
LASG动力论坛版
五. 数据准备(Gr)
— .ctl文件语法9-3
dset %ch 和 chsub t1 t2 string 若有如下NetCDF 数据文件, pr.1851-1900.nc包含了50 年 的逐月资料(即第1-600月), pr.1901-2000.nc包含100 年(第 601-1800月), 则:
pdef isize jsize { lcc | lccr } latref lonref iref jref Slat Nlat std_lon dx dy lcc 表示数据中风分量已经旋转成经纬度坐标, lccr 表示风分量仍然是lambert 正型投影,需GrADS 对其进行旋转. isize, jsize X, Y 方向的总点数; latref, lonref 参考点的经, 纬度; iref, jref 参考点的(i,j)位置; Slat, Nlat Lambert 投影m=1 时南, 北两个纬度; std_lon 标准经度; dx, dy 以M 为单位的X 和Y 方向的网格距.
— .ctl文件语法9-1
• 用一个 .ctl 文件描述多个数据文件
dset %y4%m2%d2.grd options template … … tdef 365 linear %id2%im2%iy4 1dy …
.ctl
• %y4%m2%d2 是时间的替代格式符
%y4 | %iy4 4位数的年份 %m2 | %im2 2位数的月份 %d2 | %id2 2位数的日期