《气象图形图象软件2.ppt》由会员分享，可在线阅读，更多相关《气象图形图象软件2.ppt（58页珍藏版）》请在三一文库上搜索。
1、气象图形图像软件,课程安排,GraDS绘图软件 Origin绘图软件 Compaq Array Viewer介绍,GrADS,Grid Analysis and Display System,GrADS主页 http:/www.iges.org/grads/ LASG中文讨论区: http:/ Linux Macintosh OSX Sun OS MS Windows SGI / IRIX DEC IBM / AIX 版本： GrADS 2.0（a0、a1、a2、a3）、GrADS 1.8s11、 GrADS 1.9b4、 GrADS 1.7. ftp:/iges.org/grads/ 安装,
2、1运行环境和基本操作,操作界面,提示：选择风景画（横放，11*8.5）还是肖像画（竖放， 8.5*11 ）(硬拷贝输出) 直接回车横放 no 竖放,1运行环境和基本操作,以versions 1.9为例，各程序功能：,1运行环境和基本操作,基本概念,三种文件： .dat: 原始数据文件 .ctl ：原始数据描述文件 .gs：GrADS控制文件,1运行环境和基本操作,(1).dat 二进制无格式记录的原始数据文件，其数据既可以是格点数据，也可以是站点数据。它们是从其它气象数据(如站点气象报、格点气象报、模式格点输出结果)转换生成的，格式是二进制无格式数据。 整个数据集是一个大的五维数据场，包括三维
3、物理空间、一维物理变量、一维时间变量，存放时以二维数组片的形式按照水平、垂直、物理变量、时间序列的顺序排放。 第一维：从西到东 第二维；从南到北 第三维：从低层到高层 第四维： 变量 第五维： 时间,GrADS数据数据类型与格式 支持的数据格式： netCDF数据格式 HDF 数据格式 Binary数据格式 站点数据格式 GRIB数据格式,数据存放：,1运行环境和基本操作,DSET model.dat OPTIONS little_endian UNDEF -2.56E33 TITLE 5 Days of Sample Model Output XDEF 72 LINEAR 0.0 5.0 Y
4、DEF 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,(2).ctl 原始数据描述
5、文件。是一个文本文件，用于描述原始数据文件的基本信息，包括数据集文件名、数据类型、数据结构、变量描述等。在GrADS中至少得首先打开一个数据描述文件。,格点数据的ctl文件,1运行环境和基本操作,DSET data-set-name（给定二进制原始数据文件的文件名（包含路径） DTYPE （数据类型）(station / grib) （可省略） OPTIONS （可省略） 选项： yrev: y维数方向反向 zrev: z维数方向反向 template：用于多个时间序列原始文件仅用一个描述 UNDEF value（定义缺测值） TITLE string（用字符串简略描述数据文件的内容） （可省
6、略）,Model.ctl,1运行环境和基本操作,XDEF number LINEAR start increment XDEF number LEVELS value-list,时空维数环境设置,number：x方向格点数，用整数； LINEAR：格距均匀；start:起始经度，负数表示西经； increment：格距； LEVELS：格距不均匀，用枚举法列出各网格点对应的经度值； value-list：顺序列出各格点的经度值，可在下一行续行。,例：XDEF 72 LINEAR 0.0 5.0,1运行环境和基本操作,YDEF number mapping start YDEF number m
7、apping Mapping: 映射方式：,1运行环境和基本操作,ZDEF number mapping ZDEF number mapping ,例如： ZDEF 7 LEVELS 1000 850 700 500 300 200 100,TDEF number LINEAR start-time increment number:时次数 start-time格式：hh:mmZddmmmyyyy,例： 14:20Z22JAN1987 12Z1JAN1990 JUN1960,hh:两位数的小时,mm：两位数的分钟,dd：为一或二位数的日期,mmm：三个字符的月份缩写,yyyy：为两位或四位数年
8、份（两位代表年份：1950-2049）,注：不给出时，hh缺省为00时，mm缺省为00分，dd缺省为1号 ，月年值不能缺省。整个时间串不能有空格。,increment为时间增量，格式为vvkk,其中vv为增量值，用或位整型数表示， kk为增量类型，有如下几种：,1运行环境和基本操作,VARS number: 表示变量描述开始，并给出变量个数 每个变量描述记录格式如下：,abrev:由个字符组成的变量名，要求以字母开头 levs：表示该变量在本数据集中含有的垂直层次数，整型数，若 该量为地表变量，仅一层，则设为 units：预留量，暂设为常数99 description: 为一段说明该变量的字符
9、串，最多40个字符,ENDVARS:表示数据描述文件结束,PS 0 99 Surface Pressure,Model.ctl,1运行环境和基本操作,(3).gs GrADS控制文件，用GrADS命令 run执行。这也是一个纯文本文件， 内含用描述语言 script language写成的批处理 GrADS系统设置和命令， 可集成处理 GrADS命令。,运行*.gs文件：直接在ga- 输入run .gs文件(含路径)，run 可省略,reinit open h:2model.ctl set lat 20 40 set lon 90 120 set z 5 set t 1 d t,*.gs文件简
10、单范例：,reinit *同时关闭所有打开的文件，并释放所有定义的对象 enable print h:21.gmf *与print合用，输出图形 open h:2model.ctl set lat 20 40 *设置纬度 set lon 90 120 *设置经度 set z 5 *设置高度 set t 1 *设置时间 d z *显示图形命令：d（或display） c *清屏（或clear） d u;v *显示图形命令：d（或display） print *将当前显示的图形保存到输出文件中） disable print *关闭该输出文件,(1)二维等值线绘图 (2)二维填色等值线绘图 (3)标
11、示网格点数值 (4)矢量箭头形式绘二维风场 (5)天气填图 (6)直方图 (7)单线图 (8)用指定颜色填充特定点 (9)风向杆形式绘二维风场,根据需要，GrADS可画以下不同类型的图,练习写*.ctl文件,exer.dat 数据，保存在d:model路径下，数据类型为格点资料,缺测值为99999， x方向 80个格点，从西经140度开始，从西到东排列，格距为1度， y方向 50个格点，从北纬20度开始，从南到北排列，格距为1度， 垂直方向7层，线性排列1000hpa到400hpa，间隔100hpa： 时间从2008年4月10日00时开始到4月12日00，每12小时一次 变量 u(风向的东西分
12、量，7层) v(风向的南北分量，7层) rh(相对湿度，5层) slp（海平面气压，0层）,简单演示,DSET model.dat OPTIONS little_endian 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 W
13、inds 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,model.dat 文件描述文件如下：,2基本绘图指令,一 图形类型设置 命令： set gxout graphics-type graphics-type包括以下图形类型：,coutour 二维等值线绘图 shaded 二维填色等值线绘图 grid 二维场不画图，在网格中央标示网格点数值 vecto
14、r 矢量箭头形式绘二维风场 stream 流线形式绘二维风场 barb 风向杆形式绘二维风场 bar 直方图 line 单线图 fgrid 二维场不画等值线，仅用指定颜色填充特定点 model 天气填图,二 图形要素设置,对图形类型为等值线（gxout=contour）起作用的设置： set ccolor color 设置等值线颜色，颜色号color为0黑、1白、2红、3绿、4兰、5青、 6紫、7黄、8桔黄、9灰。clear或display即重新设定颜色 set ccolor rainbow 设置等值线用七色序列表示 set ccolor revrain 设置等值线用反序的七色序列表示 set
15、 cstyle style 设置等值线线形， style：1实线、2长虚线、3短虚线、4长短虚线、5点线，clear或display即重新设定,set cthick thckns 设置等值线线宽， 取值1-10之间的整数，缺省为4 set clab on/off/forced/string/auto 设置等值线标记方式， 一直持续到重新设置,2. 对图形类型为等值线（gxout=contour或shaded）起作用的设置： set cint value 设置等值线间隔。 clear或display即重新设定或回到缺省设置 set clevs lev1 lev2 设置特定等值线值，只画lev1
16、lev2值所在的等值线， clear或d即重新设定 set cmin value 不画低于此value的等值线， c或d即重新设定 set cmax value 不画高于此value的等值线， c或d即重新设定 set black val1 val2 不画值介于val1和val2之间的等值线， c或d即重新设定,*仅对gxout=shaded起作用的设置： cbarn sf vert xmid ymid 绘制shaded 图形的填色标尺 sf：色标尺寸，1为全尺寸，0.5为半英寸 vert：0为水平，1为竖直 xmid ymid:色标中心点的位置 例： cbarn 1 1 8.8 4.3,3.
17、,4.,5.,什么是虚页？,6 如何画矩形、直线等命令,6 如何画矩形、直线等命令,三 坐标要素控制,实习要求： 一、熟悉基本环境、打开、运行相关文件，熟练简单绘图命令 二、学会画等值线图、阴影图、流场图、单线图、熟练选择各种线形、颜色、标记等,3变量和函数的定义、功能,完整的变量名形式为： abbreb.file#(dimexpr,dimexpr,) 其中 abbreb:是数据描述文件给出的变量名缩写 file#为包含此变量的已打开的文件序号 dimexpr为当前维数环境进行的局地维数设置表达式，其中绝对维数表达式为： x|y|z|t|lon|lat|lev|time= value 例： z
18、.2(lev=500): 表示文件2中高度为500hPa等压面上的变量z u.2(t=2,lev=500): 表示文件2中时次为2高度为 500hPa等压面上的变量u,相对维数表达式（相对于当前维数环境设置）： x|y|z|t|lon|lat|lev|time+/-offset offset：为相对偏差维数值 例： ps.1(time-12hr): 相对于当前时间之前12小时时刻的1号文件中的变量ps q.2(t-1,lev=850): 2号文件中相对于当前时刻前一时刻，高度为850hpa面上的变量q,表达式 GraDS的变量可进行运算：+ - * /,即可以在*.gs中输入命令：d z-z(
19、t-1),定义临时变量,例：define zave=(z(lev=850)+z(lev=500)/2 d zave,可定义一些新的临时变量供以后的操作使用，定义方式如下： define varname=expression,GraDS可调用的内部函数,4 如何把文本数据写成二进制数据,文本数据,程序： real p400(150,91) open(1,file=f:208092420.000) open(2,file=f:208092420.dat,form=binary) read(1,*) read(1,*) read(1,*) do j=1,91 read(1,*)(p400(i,j),
20、i=1,150) enddo do j=1,91 write(2)(p400(i,j),i=1,150) enddo end,相应的描述文件： DSET n:208092420.dat TITLE 400 height undef 999999 XDEF 150 LINEAR 30 1.0 YDEF 91 LINEAR 0.0 1.0 ZDEF 1 LEVELS 400 1 TDEF 1 LINEAR 20z24sep2008 1DY VARS 1 z 1 99 Surface Pressure ENDVARS,5台站资料的使用,台站资料,站点数据排放次序: micaps第三类数据通用格式（共
21、计14行） diamond 3 05年05月18日20时6小时降水量 05 05 18 20 -1 6 1 5 10 25 50 100 1 25 39 122.7 37.3 122.5 36.7 120.8 36.2 119.3 35.0 120.4 34.3 121.3 38.5 121.2 40.3 118.5 38.6 120.8 38.1 1 582,54115 118.07 43.60 800 0 54134 121.28 43.60 242 .4 58634 118.25 28.68 114 7.0 58718 116.38 27.20 95 .4 59102 115.65 24
22、.95 299 .6 57894 114.16 26.58 838 .01 58519 116.68 29.00 41 6.0 57793 114.38 27.80 129 8.0 57598 114.58 29.03 147 12.0 58626 117.21 28.30 49 22.0 58506 115.98 29.58 1165 46.0 58813 116.33 26.85 142 .01 58715 116.65 27.58 82 3.0 57799 114.97 27.12 78 7.0 58527 117.20 29.30 60 5.0 54857 120.33 36.07 7
23、7 0 53446 109.85 40.66 1069 0,台站号 经度 纬度 海拔高度 降水,地面报变量,台站型离散资料画图的一般原则 1、把台站资料r.dat写成二进制文件r.grd 2、编写r.grd数据的数据说明文件r.ctl 3、用r.ctl生成r.map文件 4、生成插值用的背景格点文件 5、编写grid.grd的.ctl文件：grid.ctl 6、把r.grd插值到grid.grd的网格点上，并画图,1、把台站资料*.000写成二进制文件*.dat 运行fortran程序shikuang1.f,parameter(n=582) character*8 stid(n) real r
24、lat(n),rlon(n),rain(n),nh(n) open(1,file=H:wwwstation05051820.000) open(10,file=H:wwwstation05051820.dat,form=bin -* ary) do i=1,14 read(1,*) enddo do i=1,n read(1,*,end=1111) stid(i),rlon(i),rlat(i),nh(i),rain(i) 1111 continue enddo close(1),tim=0.0 *本次报告的相对时刻 nlev=1 *一组地面报 nflag=1 *=1：表示有地面报；=0则无
25、do i=1,n write(10)stid(i),rlat(i),rlon(i),tim,nlev,nflag,rain(i) enddo nlev=0 *表示该时刻记录结束 write(10)stid(n-1),rlat(n-1),rlon(n-1),tim,nlev,nflag close(10) end,2、编写*.dat数据的数据说明文件*.ctl,dset H:wwwstation05051820.dat dtype station stnmap H:wwwstation05051820.map undef -999.0 tdef 1 linear 00z16aug2002 1dy
26、 vars 1 rain 0 99 rain endvars,rain.ctl,3、用*.ctl生成*.map文件,“开始”/“所有程序”/“附件”/打开“命令提示符”，进入dos 系统，进入“C:Program FilesPCGrADSwin32e”， 运行“stnmap”提示输入“stn ctl filename”如“H:wwwstationrain.ctl”生成相应的*.map文件,4、生成插值用的网格点文件(已有),grid.grd,5、编写grid.grd的.ctl文件：grid.ctl（已有）,grid.ctl,dset H:wwwstationgrid.grd undef -99
27、9.0 xdef 40 linear 95 1 ydef 40 linear 5 1 zdef 1 linear 1 1 tdef 1 linear 00z16aug2002 1dy vars 1 g 0 99 grid data is prepared for oacres function endvars,6、把*.dat插值到grid.grd的网格点上，并画图(编写*.gs),reinit open H:wwwstationgrid.ctl open H:wwwstationrain.ctl enable print H:wwwstation05051820.gmf set lon 105 126 set lat 23 38 define ro=oacres(g,rain.2,2.5,0.8) set grid off set xlint 1 set ylint 1 set xlopts 1 4 0.1 set ylopts 1 4 0.1 set cint 10 set grads off set gxout contour set clab forced d ro print disable print,oacres函数： 功能：Cressman内插函数， 将台站资料插至格点上。,6小时降水分布图,