GrADS第4章 各类参数设置

第4章各类参数设置

set是一项非常重要的命令，用于设置各种运行环境的参数，包括维数环境、图形类型、图形要素、屏幕显示等等。4.1维数环境的设置

虽然在数据描述文件中给出了各物理变量数组的时空维数范围，但在GrADS

运行环境中还需设定全数据集中参与操作的部分或全部数据集的维数情况，以供以后的表达式、显示命令等使用。这就是维数环境的设置。1.维数环境的概念

GrADS系统视每一个物理变量（VAR）场为一个四维数据集（4Ddataset），即包括三维空间（x，y，z）和一维时间（t）。也可固定其中的一维或几维以获得实际的低于四维的数据子集。2.设置的作用

GrADS

中设置维数用以说明或指定随后的分析或图形操作时参加操作的原始数据集的维数范围。该工作数据集可以是整个原始数据场，也可以是原始数据场的一部分。3.维数环境的定义（1）一种是地球坐标（worldcoordinate）以经纬度为度量单位。形式如下：

setlon|lat|lev|timeval1<val2>

（2）一种是格点坐标（gridcoordinate）以网格点数为度量单位。形式如下：

setx|y|z|tval1<val2>setlon-90

或setx1setlat40

或sety22setlev500

或setz6settime01jan1985或sett5setx1144或setlon0360sety173或setlat-9090setz13或setlev1000850sett13或settime01jan198503jan1985

注意:grads只能显示两维或一维不固定的图！

当所有维数都固定时，得到的是一个单值数据点维数环境的设置一直保持到再次设定时都有效。要查看当前维数情况，可以在运行环境中输入命令：

ga->qdimssetlon180360不固定

setlat090不固定

setlev500sett2dhgtsetlon115setlat15setlev850sett112不固定

du二维平面图一维曲线图4.2图形类型的设置

当维数环境确定后，缺省情况下，一维变量输出的图形为单线图，二维变量为等值线图，改变缺省图形输出类型的命令为：

setgxout

graphics_type

1.格点数据contour:二维等值线图（缺省设置）

setgxoutcontourshaded:二维填色图(通常与contour并用)

先画填色图，再画等值线图；否则等值线会被覆盖。

setgxoutshadeddhgtsetgxoutcontourdhgtsetlon0210setlat070setlev500sett2setgxoutcontourdhgt/10sett1

setgxoutshaded

setclevs-11

setccols1004

d(hgt(t=1)-ave(hgt,t=1,t=50))/10

setgxoutcontour

d(hgt(t=1)-ave(hgt,t=1,t=50))/10grid:二维场不绘图，以网格形式在各网格点中央标出该点数值fgrid:用指定颜色填充二维格点场，对二维场不绘制等值线图，只是将特定值的格点用指定的颜色填充该网格。与命令setfgvalsvalcol合用用法：setgxoutfgridsetfgvalsval1col1setfgvalsval2col2setlon70130setlat1030setlev1000sett1setgxoutgrid

dusetlon70130setlat1030setlev1000sett1setgxoutfgridsetfgvals-14Setfgvals-35duvector:矢量箭头形式绘二维风场（缺省设置）

setgxoutvectordu;v（u,v在同一个数据文件中可这样）或du.1;v.2（u,v在两个数据文件中）stream:流线形式绘制二维风场

setgxoutstreamdu;vbarb:风向杆形式绘制二维风场

setgxoutbarbdu;vline:对一维场绘制单线图setlon120setlat20setlev200sett112setgxoutlinedubar:对一维场不绘单线图，而绘制直方图setlon120setlat20setlev200sett112setgxoutbarsetbargap40Setbarbase0dvlinefill:两单曲线之间填色与命令setlfcolscol1col2合用。用法：setgxoutlinefill

setlfcolscol1col2dline1;line2当line1<line2时，用颜色1（col1）当line1>line2时，用颜色2（col2）'setgxoutline''setcstyle1''setcmark0''setccolor2''dv.4''setgxoutline''setcstyle1''setcmark0''setccolor1''dv.5''setgxoutlinefill''setlfcols410''dv.4;v.5'errbar:单线图误差分布

setgxouterrbardline1;line2'reinit''openD:\grads\data\u.ctl''enableprintD:\grads\data\2.gmf''setgridoff''setparea11018''setgradsoff''setxlopts120.15''setylopts120.15''sett148''setvrange-4060''setgxoutline''setcmark2''du''setgxouterrbar''setccolor2''du-12.8;u+12.8''print''c''disableprint''reinit';scatter散点图setgxoutscatterdv.1;v.2

'reinit'

'opend:\lon.ctl'

'opend:\lat.ctl‘'enableprintD:\grads\data\1.gmf'

'setx1';'set

y1';'setz1';'sett124''setgxoutscatter''dv.1;v.2''print''c''disableprint''reinit';fwrite:将输出结果(数据)存入一个由“setfwrite

文件名”所指定的文件中‘openc:\pcgrads\model\model.le.ctl’

‘setgxoutfwrite’（设置输出类型）

‘setfwrited:\result\u.grd’

it=1while(it<=5)‘setlev500’‘setlon120270’‘setlat–1010’（设置维数环境）‘sett‘it’’‘du’it=it+1endwhile‘disablefwrite’（关闭文件u.grd）注意：

setgxoutfwrite是一种特殊的输出格式，它将向图形输出的数据转向一个二进制文件输出。

由于地图坐标转换向格点坐标转换过程中的误差，需要仔细的确定输出数据的维数。确定输出数据的维数后，可利用Fortran程序进行读取。2.站点数据value:在各站点标值（缺省设置）barb:在各站点绘风向标（缺省设置）findstn:搜索最近的站点stnmark:站点位置分布model:以天气图形式将天气观测各分量填放在站点四周wxsym:绘wx天气符号setgxoutvalue各站点标值'reinit''opene:\bs\rain\rain.qht7.ctl''enableprintD:\grads\data\rain.41.gmf''setgridoff''setxlint5'‘setmpdsetmres''setlat1555''setlon100135''sett1''setgxoutvalue''drain''rune:\bs\rain\test.gs''rune:\bs\rain\southsea.gs''print''c''disable''reinit‘;setgxoutfindstn搜索最近的站点drain;x;y得到的结果：窗口中坐标为（x,y）最近站点的站点号；经纬度；该站点物理量值例如：'setstidon''setgxoutfindstn''drain;7.5;3.3'站点数据变量或者表达式x,y是窗口中的坐标位置

选择显示站点代码'reinit''opene:\bs\rain\rain.qht7.ctl''setlon73136''setlat1754''setstidon''setgxoutfindstn''drain;7.5;3.3'id=subwrd(result,1)x=subwrd(result,2)y=subwrd(result,3)value=subwrd(result,4)say'========================'say'stid:'idsay'value:'valuesay'lon:'xsay'lat:'y'setstatoff''reinit'写出结果将findstn的结果赋值到id，x,y,value'setgxoutstnmark站点位置分布'setstidon''setgxoutstnmark''drain'对于wxsym类型的图形,站点坐标的每一个数值都假定为一种wx天气符号的代码，台风和热带风暴的符号包含在wx天气符号集中。drawwxsymsymbolxysize<color<thickness>>（在一指定的坐标绘一指定的天气符号）其中：symbol—欲绘天气符号的代码x—x坐标（英寸）y—y坐标size—天气符号的大小color—天气符号的颜色，-1为标准色(红色代表热带风暴,蓝色代表降雪等等)thickness—天气符号的线宽当图形输出类型为vector、stream和barb时,在显示命令后需给出两种数据，两者间用分号“;”隔开，前者理解为U分量，后者理解为V分量当图形输出类型为errbar、scatter、linefill时，也须用分号将两种数据分开例如：

displayu；v

显示u、v风场的合成矢量图说明：对于vector和stream类型的图形,可以给出第三种数据，用于指定矢量或流线图的色彩分布，彩色值代表第三分量值。例如：

displayu；v；hcurl(u,v)

（计算垂直涡度）

displayu；v；w

（第三个变量给定矢量或流线图的色彩分布，彩色值代表第三分量值）4.3图形要素的设置

有些设置对多数图形输出类型有效，有些设置只对某一种图形输出类型有效。有些设置一旦设定，会一直保持不变，有些设置在键入“clear”或“display”命令后会回到缺省设置状态。1.对于图形类型为contour起作用的设置（1）setccolorcolor设置等值线颜色。其中，color为颜色号：0-黑，1-白，2-红，3-绿，4-兰，5-青，6-洋红，7-黄，8-橘黄，15-灰。该设置在“clear”或“display”命令后即重新设定颜色。（2）setccolorrainbow

设定等值线颜色用七彩序列表示。（3）setcstylestyle

设定等值线线型其中，style为线型号：1-实线，2-长虚线，3-短虚线，5-点线。该设置在“clear”或“display”命令后即重新设定。（4）setcthick

thckns

设定等值线线宽。其中，thckns为线宽值，取值范围:1～10之间的整数，线宽大于等于6的线条在屏幕上用粗线显示，主要用于控制硬拷贝输出。（5）setcterp

on|off

设置样条插值光滑开关再定义后才重新设置，填色图没有样条光滑，设定cterp为off可使填色图与等值线图的边缘重合，也可用csmooth选项来达到上述目的。（6）setclabon|off|forced|string|auto

控制等值线的标记方式。

再定义后才重新设置。on-快速等值线标记，标记在等值线的水平处(缺省设置)off-不标记forced-强迫标记所有等值线string-用字符串string替换等值线标记数值auto–缺省方式例如：setclab%gK

将在等值线标记数值的最后标记上‘K’setclab%g%%将在等值线标记数值的最后标记上‘％’setclab%.2f以小数点后保留2位的格式标记等值线数值setclab%03.0f以3位整数(首位0保留)形式标记等值线数值setclab

Foo

用Foo标记所有等值线（7）setcloptscolor<thickness<size>>设置等值线标记的颜色。其中：

color是颜色号，-1为省缺，

表示采用等值线的颜色进行标记；

thickness为标记的线宽，-1为省缺；

size为标记的大小，0.09为省缺。该设置在下一个setclopts命令前一直有效。（8）setclskipnumber

设置间隔几根等值线标示数值。2.对于图形类型为contour或shaded起作用的设置

（1）setcintvalue

设置等值线间隔，该设置在“clear”或“display”命令后即重新设置。例如：间隔3m/s画风速线：setcint3du（2）setclevslev1lev2...-设置特定的等值线值，只画lev1，lev2，值所在的等值线，用于不等间隔绘图，该设置在“clear”或“display”命令后即重新设置。例如：画零风速线：setclevs0du

画赤道:setclevs0dlat

（3）setccolscol1col2...

设置对应于“setclev”命令设定的特定等值线的颜色，该设置在“clear”或“display”命令后即重新设置。缺省的彩虹颜色号序列为：9,14,4,11,5,13,3,10,7,12,8,2,6。（4）setcsmooth

on|off

设置是否将网格值重新插值。如取“on”，则在绘等值线图前用三次插值将现网格值插到更精细网格上，重新设置才改变本次设置。（5）setcminvalue设置不画低于此value值的等值线。该设置在“clear”或“display”命令后即重新设置。（6）setcmaxvalue设置不画高于此value值的等值线。该设置在“clear”或“display”命令后即重新设置。（7）setcblackval1val2设置不画介于val1和val2之间的等值线。该设置在“clear”或“display”命令后即重新设置。3.对于图形类型为contour,shaded,vector,stream起作用的设置

（1）setstrmdenvalue设置流线密度其中，value的取值范围为1—10的整数，5为缺省设置值。（2）setrbcolscolor1color2<color3>...

设置新的彩虹颜色序列。其中，color1，color2，可以用“setrgb”命令定义新的颜色号，该新的彩虹颜色序列在随后的彩虹颜色调用中取代原缺省的彩虹颜色序列，重新设置后才改变原设置。（3）setrbcolsauto

起用内定的彩虹颜色。重新设置后才改变原设置。（4）setrbrangelowhigh

设置彩虹颜色序列对应的等值线的取值范围，缺省时，最低值和最高值对应取为变量场的最小值和最大值，“clear”命令后即重新设置。4.对于图形类型为line起作用的设置（1）setccolorcolor--设置线条的颜色号。该设置在“clear”或“display”命令后即重新设置。（2）setcstylestyle设置线条类型。该设置在“clear”或“display”命令后即重新设置。（3）setcmarkmarker

设置线条上的标记符号。其中：

marker为标记符号值：0-无标记，1-叉号，2-空心圆，3-实心圆，4-空心方框，5-实心方框。该设置在“clear”或“display”命令后即重新设置。（4）setmissconn

on|off

缺省设置时，线条在缺测资料点断开，设置“setmissconnon”将连接缺测资料点。5.对于图形类型为bar起作用的设置

（1）setbargap

val

以百分比值设定直方条之间的间隔。其中，val取值为0～100，省缺值为0，即无间隔，当val取100时直方图退化为垂直线条直方图。（2）setbarbase

val|bottom|top

如给出val值，则各直方条从该值处起画(向上和向下)，所画直方条取值于y轴坐标尺度之内；如给bottom，各直方条从图框的底边向上绘出；如给top，直方条从图框顶边向下绘出。从坐标最底值开'setbarbase32''setbargap20'setlon120setlat20setlev200sett112setgxoutbarsetbargap40Setbarbase0dv6.对于图形类型为grid起作用的设置（1）setdignumnumber

设置小数点后的位数为number。

（2）setdigsizesize

设置数字的字符大小，size单位为英寸，通常取为0.1-0.15。7.对于图形类型为fgrid起作用的设置

setfgvalsvaluecolor<valuecolor><valuecolor>...

对取值为value的网格点用颜色为color的色块标记该网格，每个格点的值取法是四舍五入，要绘出的值点须逐个举出，未列出的值不绘图。见第12页ppt8.对于图形类型为vector起作用的设置

（1）setarrsclsize<magnitude>

设置矢量箭头的长度为size英寸(虚页英寸)。通常size取为0.5～1.0；选项magnitude为设定矢量箭头的大小。缺省时所有矢量同长，该设置在“clear”或“display”命令后即重新设置。（2）setarrowheadsize

设置箭头大小。

size值通常取为0.05，如取为0，则不画箭头的头；如取为负值，箭头大小与矢量值成比例（张角的大小）。4.4坐标要素的设置

（1）setvrangey1y2

setvrangex1x2

分别设置y坐标轴和x坐标轴标尺的取值范围，“clear”命令后即重新设置。'setvrange3246'不设置setvrange（2）setzlog

on|off

对z维数方向取对数尺度的开关。

on表示z维数方向取为对数尺度。重新设置后才改变原设置。setzlogonsetzlogoff（3）setxaxis|yaxisstartend<incr>

设置坐标轴（x轴或y轴）的坐标从起始值start到结束值end，并用incr作为刻度的增量，标尺可与所绘的数据和维数无关。'setxaxis1204''setyaxis1102'（4）setgridon|off|horizontal|vertical|linestyle|color

控制是否绘网格线。

on绘网格(缺省)；

off不绘网格（一般不绘制网格）；

color和linestyle为网格线的颜色和线型，缺省时，color为15(灰)，linestyle为5(点线)；

horizontal表示只画水平网格线；

vertical表示只画垂直网格线。缺省时或setgridon'setgridhorizontal22'线型颜色（５）setxloptscolor<thickness<size>>

setyloptscolor<thickness<size>>设置x，y轴的颜色、线宽和字符大小。其中：

xlopts控制x坐标轴；

ylopts控制y坐标轴；

color为坐标轴标尺的颜色号（缺省为1）；

thickness为坐标轴标尺的线宽（缺省为4）；

size为坐标轴刻度的大小（缺省为0.12）。'setxlopts120.2''setylopts120.2''setxlopts240.3''setylopts240.3'（６）setxlevslab1lab2...setylevslab1lab2...

设置x，y坐标轴标尺上要标记的值，本设置不适用于时间坐标轴，“clear”命令后即重新设置。'setxlevs607085105120135''setylevs0152040'（7）setxlintintervalsetylintinterval

设置坐标轴的标记间隔interval。

setxlevs/ylevs可再控制标记的分布，clear命令后即重新设置。setxlint10setylint10setxlint10setylint5（8）

setgradson|off

开关选择是否打印出GrADS的标注。一般不打印，要不然插入word时还要裁剪图片。4.5图注设置

（１）drawtitlestring

使用“draw”命令在图形顶部写一字符串string作为图的标题，字符串中如有反斜杆“/”，表示另起新行。'drawtitle1pentad'（２）setxlabon|off|auto|string

setylabon|off|auto|string同setclab命令原理。例如：setylab%.2f将以小数点后两位的格式标记y轴坐标。“clear”命令后即重新设置。drawxlabLongitudedrawylabpentad（３）setannotcolor<thickness>

设置上述图注所用的颜色color和线宽thicknesses。省缺时，颜色为白色，线宽为6。

该命令同时也设置了坐标轴的边框、坐标轴标记和刻度的颜色及线宽值，坐标刻度和标尺的实际绘图线宽为thickness值减1。

（４）setxyrevon

交换水平和垂直坐标轴所代表的维数方向。（５）setxfliponsetyflipon

设置水平坐标轴或垂直坐标轴的维数方向取反向。以上坐标轴的设置当“clear”或“setvpage”命令后即重新设置。4.6地图投影设置（１）setmproj

proj

设置当前地图投影方式。proj取值包括：latlon-缺省设置，用固定的投影角进行Lat/lon投影；scaled-用不固定的投影角进行latlon投影，地图比例失效；nps-北半球极地投影；sps-南半球极地投影；off-同“scaled”设置，但不画出地图，坐标轴也不代表lat/lon；robinson

–robinson伪圆柱投影（x:-180180；y:-9090；）；lambert-lambert投影（等角圆锥投影）；orthogr：正投影；mollweide：等面积伪圆柱投影。（2）setmpvals

lonmin

lonmax

latmin

latmax

设置极地投影时经度和纬度值的取值范围，缺省时取为当前维数环境。

setmpvals-140-401580

setmprojnps

setmpdsetmres

dz（3）setmpdset<lowres|mres|hires|nmap>

设置地图数据集。

lowres（缺省）为低分辨率地图集，

mres（hires）为中（高）分辨率地图集，

nmap为北美地图集。（4）setpoli

on|off

在mres或hires地图集中开关选择是否使用行政边界，省缺设置为on。（5）setmapcolorstylethickness

设置地图背景的颜色color、线型style和线宽thickness。（6）setmpdraw

on|off

选off则不绘地图背景，但地图标尺仍起作用，非经纬度数据需此项。4.7绘图区域设置

GrADS

的绘图工作区分三个层次，一层是实际页（realpage）；一层是虚拟页（virtualpage）缺省时虚页等同于实页；第三层是在虚页中指定绘图区域（parea），即缺省时等同于实际英寸，当设置虚页后按比例度量。（1）setvpage

xmin

xmax

ymin

ymax

通过定义在实页上一个或多个虚页来控制绘图的数目和大小。本命令在实页上用xmin，xmax，ymin，ymax（英寸）设置了一个虚页，随后的所有图形都输出到这张虚页上（虚页单位为英寸），直到下个“setvpage”命令出现。

1英寸=2.54cm不能超出画布大小8.5*11或11*8.5

set

vpage01108.5set

vpage08.5011set

vpage1816（２）setvpageoff回到缺省的实虚页相同的状态。（３）setparea

xmin

xmax

ymin

ymax

在虚页中定义了一块区域parea

用于GrADS

的绘图，但该区域不包括标题、坐标轴标记等。要预留出写标题和坐标轴标记的位置。设置的区域用于等值线绘图、地图绘图、单线绘图，该区域内以虚页英寸为单位。缺省时，自动按图形类型设置绘图区域.

可实现一页多图

'setvpage05.508''sett1''du''setylaboff''setvpage5.51108''sett2''du'

（４）setpareaoff回到缺省状态。4.8屏幕显示设置

（１）setdisplaygrey|greyscale|color<black|white>

设置显示模式，省缺为彩色，填色图和等直线图用彩虹的颜色填绘。grey–灰色显示；greyscale-等直线用单一的灰色显示，填色图用greyscale序列填色。color–选项black或white用于设定硬拷贝的背景色，省缺为black。（２）setframe