资源描述：
《grads绘图与编程(值得收藏)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、GrADS绘图与编程图例11月份500hPa位势高度场2图例2500hPa高度场（1、4、7、10月）3图例3亚洲季风区7月850hPa流场图4图例4SST和SeaLevel分布图5图例5（130E，25N）处的感热和潜热的时间演变6第一章概述★简介★GrADS的安装★GrADS软件的组成★GrADS的启动和退出★所需预备知识★学习本课程的要求7简介GrADS是GridAnalysisandDisplaySystem的缩写,它由美国马里兰大学气象系BrianE.Doty开发。利用该系统可实现包括格点数据和站点数据的彩色气象图形分析和显示。具有操作简便、功
2、能丰富、图形美观、显示快速的特点。在国内外气象界得到广泛应用。版本：V1.5,V1.6.8,V1.7.7，V1.7.9等8安装1、将系统软件拷入硬盘某一子目录下，如：c:grads2、修改autoexec.bat文件，加上pathc:grads;%path% setgaddir=c:grads setgascrp=c:grads9GrADS软件包的组成Grads.exe系统的核心文件,由此进入GrADS环境Dos4gw.exe由GrADS.exe文件调用Gribmap.exe产生格点资料映射文件(*.idx)Gribscan.exe看Grib码资
3、料，转为TXT格式文件Gxtran.exe显示*.gmf格式的图形文件Stnmap.exe产生台站资料的映射文件(*.map)Wgrib.exeGrib码资料的解读程序Gv.exe将*.gmf格式的图形文件转换为*.wmf格式10启动和退出GrADS在DOS下键入：grads既可进入；键入:quit即可退出grads,返回DOS。GrADS在初始化绘图环境之前，将提示用户选择风景画(Landscape)或肖像画(Portrait)形式。风景画形式的大小为11×8.5英寸，肖像画为8.5×11英寸。启动GrADS系统时，有以下选择：-b以批处理形式运行Gr
4、ADS。-l以风景画形式运行GrADS。-p以肖像画形式运行GrADS。-c进入GrADS后首先执行随后跟在RUN命令后的描述文件，如：grads–c“runprofile.gs”上述选项可结合起来使用，如：grads–blc“runbatch.gs”118.58.51111所需预备知识DOS和Windows基础FORTRAN编程WORD基础学习要求强调上机编程实践12第二章GrADS数据格式目标(1)能将文本格式数据转为GrADS格式(2)能写出数据描述文件(*.ctl)一、维数环境的概念重要概念1、含义：GrADS视每一个变量(VAR)场为一个四维数
5、据集，即包括三维空间(x,y,z)和一维时间(t)。2、作用：说明和指定随后的分析或图形操作时参加的原始数据集的维数范围。3、定义方法：Setlat
6、lon
7、lev
8、timeval1地球坐标Setx
9、y
10、z
11、tval1格点坐标13x(或lon)从西向东的水平坐标y(或lat)从南到北的水平坐标z(或lev)从地面到高空的垂直坐标t(或time)时序坐标。如：setlon40160；setlat060定义了水平变化范围setlat30；setlon0180定义了沿30N的纬向变化范围settimejan81dec94定义了从81年1
12、月到94年12月的时段 所以，固定的维数环境和变化的维数环境相结合，就构成了当前的维数环境。14(1)x,y,z,t均固定时，得到一个单值数据点(2)x,y,z,t中只有一维变化时，得到一维曲线(3)二维发生变化时，对应于二维剖面图，如X-Y平面图，X-Z，Y-Z，X-T，Y-T，Z-T剖面图(4)当三维发生变化时，GrADS以动画方式显示二维切片。15说出下列维数环境的含义Setlon60150Setlat030Sett7Setlev200Setlon120Setlat30Setz112Sett1Setlat30Setlon60130Sett112S
13、etz1Setlon60150Setlat030Sett112Setlev20016二、文本格式数据GrADS格式现有如下资料：名称：U850，V850，U200，V200，H500和TSFC范围：60~150E，0~40N分辨率：2.5×2.5时间：1982年1月~1985年12月的逐月资料。如何将上述资料写成GrADS下的数据格式？17顺序：(x,y)zvtT=1T=2T=3UVHTU850(X,Y)U200(X,Y)时(v)分(z)秒(x,y)5月20日(t)18Parameter(ii=37,jj=17)Realvar(ii,jj)Open(1,
14、file=‘u850.dat’)Open(2,file=‘v850.dat’)O
