GrADS第4章 各類參數(shù)設(shè)置

第4章各類參數(shù)設(shè)置

set是一項非常重要的命令，用于設(shè)置各種運行環(huán)境的參數(shù)，包括維數(shù)環(huán)境、圖形類型、圖形要素、屏幕顯示等等。4.1維數(shù)環(huán)境的設(shè)置

雖然在數(shù)據(jù)描述文件中給出了各物理變量數(shù)組的時空維數(shù)范圍，但在GrADS

運行環(huán)境中還需設(shè)定全數(shù)據(jù)集中參與操作的部分或全部數(shù)據(jù)集的維數(shù)情況，以供以后的表達式、顯示命令等使用。這就是維數(shù)環(huán)境的設(shè)置。1.維數(shù)環(huán)境的概念

GrADS系統(tǒng)視每一個物理變量（VAR）場為一個四維數(shù)據(jù)集（4Ddataset），即包括三維空間（x，y，z）和一維時間（t）。也可固定其中的一維或幾維以獲得實際的低于四維的數(shù)據(jù)子集。2.設(shè)置的作用

GrADS

中設(shè)置維數(shù)用以說明或指定隨后的分析或圖形操作時參加操作的原始數(shù)據(jù)集的維數(shù)范圍。該工作數(shù)據(jù)集可以是整個原始數(shù)據(jù)場，也可以是原始數(shù)據(jù)場的一部分。3.維數(shù)環(huán)境的定義（1）一種是地球坐標（worldcoordinate）以經(jīng)緯度為度量單位。形式如下：

setlon|lat|lev|timeval1<val2>

（2）一種是格點坐標（gridcoordinate）以網(wǎng)格點數(shù)為度量單位。形式如下：

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

或setx1setlat40

或sety22setlev500

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

注意:grads只能顯示兩維或一維不固定的圖！

當所有維數(shù)都固定時，得到的是一個單值數(shù)據(jù)點維數(shù)環(huán)境的設(shè)置一直保持到再次設(shè)定時都有效。要查看當前維數(shù)情況，可以在運行環(huán)境中輸入命令：

ga->qdimssetlon180360不固定

setlat090不固定

setlev500sett2dhgtsetlon115setlat15setlev850sett112不固定

du二維平面圖一維曲線圖4.2圖形類型的設(shè)置

當維數(shù)環(huán)境確定后，缺省情況下，一維變量輸出的圖形為單線圖，二維變量為等值線圖，改變?nèi)笔D形輸出類型的命令為：

setgxout

graphics_type

1.格點數(shù)據(jù)contour:二維等值線圖（缺省設(shè)置）

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:二維場不繪圖，以網(wǎng)格形式在各網(wǎng)格點中央標出該點數(shù)值fgrid:用指定顏色填充二維格點場，對二維場不繪制等值線圖，只是將特定值的格點用指定的顏色填充該網(wǎng)格。與命令setfgvalsvalcol合用用法：setgxoutfgridsetfgvalsval1col1setfgvalsval2col2setlon70130setlat1030setlev1000sett1setgxoutgrid

dusetlon70130setlat1030setlev1000sett1setgxoutfgridsetfgvals-14Setfgvals-35duvector:矢量箭頭形式繪二維風場（缺省設(shè)置）

setgxoutvectordu;v（u,v在同一個數(shù)據(jù)文件中可這樣）或du.1;v.2（u,v在兩個數(shù)據(jù)文件中）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:將輸出結(jié)果(數(shù)據(jù))存入一個由“setfwrite

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

‘setgxoutfwrite’（設(shè)置輸出類型）

‘setfwrited:\result\u.grd’

it=1while(it<=5)‘setlev500’‘setlon120270’‘setlat–1010’（設(shè)置維數(shù)環(huán)境）‘sett‘it’’‘du’it=it+1endwhile‘disablefwrite’（關(guān)閉文件u.grd）注意：

setgxoutfwrite是一種特殊的輸出格式，它將向圖形輸出的數(shù)據(jù)轉(zhuǎn)向一個二進制文件輸出。

由于地圖坐標轉(zhuǎn)換向格點坐標轉(zhuǎn)換過程中的誤差，需要仔細的確定輸出數(shù)據(jù)的維數(shù)。確定輸出數(shù)據(jù)的維數(shù)后，可利用Fortran程序進行讀取。2.站點數(shù)據(jù)value:在各站點標值（缺省設(shè)置）barb:在各站點繪風向標（缺省設(shè)置）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得到的結(jié)果：窗口中坐標為（x,y）最近站點的站點號；經(jīng)緯度；該站點物理量值例如：'setstidon''setgxoutfindstn''drain;7.5;3.3'站點數(shù)據(jù)變量或者表達式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'寫出結(jié)果將findstn的結(jié)果賦值到id，x,y,value'setgxoutstnmark站點位置分布'setstidon''setgxoutstnmark''drain'對于wxsym類型的圖形,站點坐標的每一個數(shù)值都假定為一種wx天氣符號的代碼，臺風和熱帶風暴的符號包含在wx天氣符號集中。drawwxsymsymbolxysize<color<thickness>>（在一指定的坐標繪一指定的天氣符號）其中：symbol—欲繪天氣符號的代碼x—x坐標（英寸）y—y坐標size—天氣符號的大小color—天氣符號的顏色，-1為標準色(紅色代表熱帶風暴,藍色代表降雪等等)thickness—天氣符號的線寬當圖形輸出類型為vector、stream和barb時,在顯示命令后需給出兩種數(shù)據(jù)，兩者間用分號“;”隔開，前者理解為U分量，后者理解為V分量當圖形輸出類型為errbar、scatter、linefill時，也須用分號將兩種數(shù)據(jù)分開例如：

displayu；v

顯示u、v風場的合成矢量圖說明：對于vector和stream類型的圖形,可以給出第三種數(shù)據(jù)，用于指定矢量或流線圖的色彩分布，彩色值代表第三分量值。例如：

displayu；v；hcurl(u,v)

（計算垂直渦度）

displayu；v；w

（第三個變量給定矢量或流線圖的色彩分布，彩色值代表第三分量值）4.3圖形要素的設(shè)置

有些設(shè)置對多數(shù)圖形輸出類型有效，有些設(shè)置只對某一種圖形輸出類型有效。有些設(shè)置一旦設(shè)定，會一直保持不變，有些設(shè)置在鍵入“clear”或“display”命令后會回到缺省設(shè)置狀態(tài)。1.對于圖形類型為contour起作用的設(shè)置（1）setccolorcolor設(shè)置等值線顏色。其中，color為顏色號：0-黑，1-白，2-紅，3-綠，4-蘭，5-青，6-洋紅，7-黃，8-橘黃，15-灰。該設(shè)置在“clear”或“display”命令后即重新設(shè)定顏色。（2）setccolorrainbow

設(shè)定等值線顏色用七彩序列表示。（3）setcstylestyle

設(shè)定等值線線型其中，style為線型號：1-實線，2-長虛線，3-短虛線，5-點線。該設(shè)置在“clear”或“display”命令后即重新設(shè)定。（4）setcthick

thckns

設(shè)定等值線線寬。其中，thckns為線寬值，取值范圍:1～10之間的整數(shù)，線寬大于等于6的線條在屏幕上用粗線顯示，主要用于控制硬拷貝輸出。（5）setcterp

on|off

設(shè)置樣條插值光滑開關(guān)再定義后才重新設(shè)置，填色圖沒有樣條光滑，設(shè)定cterp為off可使填色圖與等值線圖的邊緣重合，也可用csmooth選項來達到上述目的。（6）setclabon|off|forced|string|auto

控制等值線的標記方式。

再定義后才重新設(shè)置。on-快速等值線標記，標記在等值線的水平處(缺省設(shè)置)off-不標記forced-強迫標記所有等值線string-用字符串string替換等值線標記數(shù)值auto–缺省方式例如：setclab%gK

將在等值線標記數(shù)值的最后標記上‘K’setclab%g%%將在等值線標記數(shù)值的最后標記上‘％’setclab%.2f以小數(shù)點后保留2位的格式標記等值線數(shù)值setclab%03.0f以3位整數(shù)(首位0保留)形式標記等值線數(shù)值setclab

Foo

用Foo標記所有等值線（7）setcloptscolor<thickness<size>>設(shè)置等值線標記的顏色。其中：

color是顏色號，-1為省缺，

表示采用等值線的顏色進行標記；

thickness為標記的線寬，-1為省缺；

size為標記的大小，0.09為省缺。該設(shè)置在下一個setclopts命令前一直有效。（8）setclskipnumber

設(shè)置間隔幾根等值線標示數(shù)值。2.對于圖形類型為contour或shaded起作用的設(shè)置

（1）setcintvalue

設(shè)置等值線間隔，該設(shè)置在“clear”或“display”命令后即重新設(shè)置。例如：間隔3m/s畫風速線：setcint3du（2）setclevslev1lev2...-設(shè)置特定的等值線值，只畫lev1，lev2，值所在的等值線，用于不等間隔繪圖，該設(shè)置在“clear”或“display”命令后即重新設(shè)置。例如：畫零風速線：setclevs0du

畫赤道:setclevs0dlat

（3）setccolscol1col2...

設(shè)置對應于“setclev”命令設(shè)定的特定等值線的顏色，該設(shè)置在“clear”或“display”命令后即重新設(shè)置。缺省的彩虹顏色號序列為：9,14,4,11,5,13,3,10,7,12,8,2,6。（4）setcsmooth

on|off

設(shè)置是否將網(wǎng)格值重新插值。如取“on”，則在繪等值線圖前用三次插值將現(xiàn)網(wǎng)格值插到更精細網(wǎng)格上，重新設(shè)置才改變本次設(shè)置。（5）setcminvalue設(shè)置不畫低于此value值的等值線。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（6）setcmaxvalue設(shè)置不畫高于此value值的等值線。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（7）setcblackval1val2設(shè)置不畫介于val1和val2之間的等值線。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。3.對于圖形類型為contour,shaded,vector,stream起作用的設(shè)置

（1）setstrmdenvalue設(shè)置流線密度其中，value的取值范圍為1—10的整數(shù)，5為缺省設(shè)置值。（2）setrbcolscolor1color2<color3>...

設(shè)置新的彩虹顏色序列。其中，color1，color2，可以用“setrgb”命令定義新的顏色號，該新的彩虹顏色序列在隨后的彩虹顏色調(diào)用中取代原缺省的彩虹顏色序列，重新設(shè)置后才改變原設(shè)置。（3）setrbcolsauto

起用內(nèi)定的彩虹顏色。重新設(shè)置后才改變原設(shè)置。（4）setrbrangelowhigh

設(shè)置彩虹顏色序列對應的等值線的取值范圍，缺省時，最低值和最高值對應取為變量場的最小值和最大值，“clear”命令后即重新設(shè)置。4.對于圖形類型為line起作用的設(shè)置（1）setccolorcolor--設(shè)置線條的顏色號。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（2）setcstylestyle設(shè)置線條類型。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（3）setcmarkmarker

設(shè)置線條上的標記符號。其中：

marker為標記符號值：0-無標記，1-叉號，2-空心圓，3-實心圓，4-空心方框，5-實心方框。該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（4）setmissconn

on|off

缺省設(shè)置時，線條在缺測資料點斷開，設(shè)置“setmissconnon”將連接缺測資料點。5.對于圖形類型為bar起作用的設(shè)置

（1）setbargap

val

以百分比值設(shè)定直方條之間的間隔。其中，val取值為0～100，省缺值為0，即無間隔，當val取100時直方圖退化為垂直線條直方圖。（2）setbarbase

val|bottom|top

如給出val值，則各直方條從該值處起畫(向上和向下)，所畫直方條取值于y軸坐標尺度之內(nèi)；如給bottom，各直方條從圖框的底邊向上繪出；如給top，直方條從圖框頂邊向下繪出。從坐標最底值開'setbarbase32''setbargap20'setlon120setlat20setlev200sett112setgxoutbarsetbargap40Setbarbase0dv6.對于圖形類型為grid起作用的設(shè)置（1）setdignumnumber

設(shè)置小數(shù)點后的位數(shù)為number。

（2）setdigsizesize

設(shè)置數(shù)字的字符大小，size單位為英寸，通常取為0.1-0.15。7.對于圖形類型為fgrid起作用的設(shè)置

setfgvalsvaluecolor<valuecolor><valuecolor>...

對取值為value的網(wǎng)格點用顏色為color的色塊標記該網(wǎng)格，每個格點的值取法是四舍五入，要繪出的值點須逐個舉出，未列出的值不繪圖。見第12頁ppt8.對于圖形類型為vector起作用的設(shè)置

（1）setarrsclsize<magnitude>

設(shè)置矢量箭頭的長度為size英寸(虛頁英寸)。通常size取為0.5～1.0；選項magnitude為設(shè)定矢量箭頭的大小。缺省時所有矢量同長，該設(shè)置在“clear”或“display”命令后即重新設(shè)置。（2）setarrowheadsize

設(shè)置箭頭大小。

size值通常取為0.05，如取為0，則不畫箭頭的頭；如取為負值，箭頭大小與矢量值成比例（張角的大?。?。4.4坐標要素的設(shè)置

（1）setvrangey1y2

setvrangex1x2

分別設(shè)置y坐標軸和x坐標軸標尺的取值范圍，“clear”命令后即重新設(shè)置。'setvrange3246'不設(shè)置setvrange（2）setzlog

on|off

對z維數(shù)方向取對數(shù)尺度的開關(guān)。

on表示z維數(shù)方向取為對數(shù)尺度。重新設(shè)置后才改變原設(shè)置。setzlogonsetzlogoff（3）setxaxis|yaxisstartend<incr>

設(shè)置坐標軸（x軸或y軸）的坐標從起始值start到結(jié)束值end，并用incr作為刻度的增量，標尺可與所繪的數(shù)據(jù)和維數(shù)無關(guān)。'setxaxis1204''setyaxis1102'（4）setgridon|off|horizontal|vertical|linestyle|color

控制是否繪網(wǎng)格線。

on繪網(wǎng)格(缺省)；

off不繪網(wǎng)格（一般不繪制網(wǎng)格）；

color和linestyle為網(wǎng)格線的顏色和線型，缺省時，color為15(灰)，linestyle為5(點線)；

horizontal表示只畫水平網(wǎng)格線；

vertical表示只畫垂直網(wǎng)格線。缺省時或setgridon'setgridhorizontal22'線型顏色（５）setxloptscolor<thickness<size>>

setyloptscolor<thickness<size>>設(shè)置x，y軸的顏色、線寬和字符大小。其中：

xlopts控制x坐標軸；

ylopts控制y坐標軸；

color為坐標軸標尺的顏色號（缺省為1）；

thickness為坐標軸標尺的線寬（缺省為4）；

size為坐標軸刻度的大?。ㄈ笔?.12）。'setxlopts120.2''setylopts120.2''setxlopts240.3''setylopts240.3'（６）setxlevslab1lab2...setylevslab1lab2...

設(shè)置x，y坐標軸標尺上要標記的值，本設(shè)置不適用于時間坐標軸，“clear”命令后即重新設(shè)置。'setxlevs607085105120135''setylevs0152040'（7）setxlintintervalsetylintinterval

設(shè)置坐標軸的標記間隔interval。

setxlevs/ylevs可再控制標記的分布，clear命令后即重新設(shè)置。setxlint10setylint10setxlint10setylint5（8）

setgradson|off

開關(guān)選擇是否打印出GrADS的標注。一般不打印，要不然插入word時還要裁剪圖片。4.5圖注設(shè)置

（１）drawtitlestring

使用“draw”命令在圖形頂部寫一字符串string作為圖的標題，字符串中如有反斜桿“/”，表示另起新行。'drawtitle1pentad'（２）setxlabon|off|auto|string

setylabon|off|auto|string同setclab命令原理。例如：setylab%.2f將以小數(shù)點后兩位的格式標記y軸坐標。“clear”命令后即重新設(shè)置。drawxlabLongitudedrawylabpentad（３）setannotcolor<thickness>

設(shè)置上述圖注所用的顏色color和線寬thicknesses。省缺時，顏色為白色，線寬為6。

該命令同時也設(shè)置了坐標軸的邊框、坐標軸標記和刻度的顏色及線寬值，坐標刻度和標尺的實際繪圖線寬為thickness值減1。

（４）setxyrevon

交換水平和垂直坐標軸所代表的維數(shù)方向。（５）setxfliponsetyflipon

設(shè)置水平坐標軸或垂直坐標軸的維數(shù)方向取反向。以上坐標軸的設(shè)置當“clear”或“setvpage”命令后即重新設(shè)置。4.6地圖投影設(shè)置（１）setmproj

proj

設(shè)置當前地圖投影方式。proj取值包括：latlon-缺省設(shè)置，用固定的投影角進行Lat/lon投影；scaled-用不固定的投影角進行l(wèi)atlon投影，地圖比例失效；nps-北半球極地投影；sps-南半球極地投影；off-同“scaled”設(shè)置，但不畫出地圖，坐標軸也不代表lat/lon；robinson

–robinson偽圓柱投影（x:-180180；y:-9090；）；lambert-lambert投影（等角圓錐投影）；orthogr：正投影；mollweide：等面積偽圓柱投影。（2）setmpvals

lonmin

lonmax

latmin

latmax

設(shè)置極地投影時經(jīng)度和緯度值的取值范圍，缺省時取為當前維數(shù)環(huán)境。

setmpvals-140-401580

setmprojnps

setmpdsetmres

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

設(shè)置地圖數(shù)據(jù)集。

lowres（缺?。榈头直媛实貓D集，

mres（hires）為中（高）分辨率地圖集，

nmap為北美地圖集。（4）setpoli

on|off

在mres或hires地圖集中開關(guān)選擇是否使用行政邊界，省缺設(shè)置為on。（5）setmapcolorstylethickness

設(shè)置地圖背景的顏色color、線型style和線寬thickness。（6）setmpdraw

on|off

選off則不繪地圖背景，但地圖標尺仍起作用，非經(jīng)緯度數(shù)據(jù)需此項。4.7繪圖區(qū)域設(shè)置

GrADS

的繪圖工作區(qū)分三個層次，一層是實際頁（realpage）；一層是虛擬頁（virtualpage）缺省時虛頁等同于實頁；第三層是在虛頁中指定繪圖區(qū)域（parea），即缺省時等同于實際英寸，當設(shè)置虛頁后按比例度量。（1）setvpage

xmin

xmax

ymin

ymax

通過定義在實頁上一個或多個虛頁來控制繪圖的數(shù)目和大小。本命令在實頁上用xmin，xmax，ymin，ymax（英寸）設(shè)置了一個虛頁，隨后的所有圖形都輸出到這張?zhí)擁撋希ㄌ擁搯挝粸橛⒋纾?，直到下個“setvpage”命令出現(xiàn)。

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

set

vpage01108.5set

vpage08.5011set

vpage1816（２）setvpageoff回到缺省的實虛頁相同的狀態(tài)。（３）setparea

xmin

xmax

ymin

ymax

在虛頁中定義了一塊區(qū)域parea

用于GrADS

的繪圖，但該區(qū)域不包括標題、坐標軸標記等。要預留出寫標題和坐標軸標記的位置。設(shè)置的區(qū)域用于等值線繪圖、地圖繪圖、單線繪圖，該區(qū)域內(nèi)以虛頁英寸為單位。缺省時，自動按圖形類型設(shè)置繪圖區(qū)域.

可實現(xiàn)一頁多圖

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

（４）setpareaoff回到缺省狀態(tài)。4.8屏幕顯示設(shè)置

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

設(shè)置顯示模式，省缺為彩色，填色圖和等直線圖用彩虹的顏色填繪。grey–灰色顯示；greyscale-等直線用單一的灰色顯示，填色圖用greyscale序列填色。color–選項black或white用于設(shè)定硬拷貝的背景色，省缺為black。（２）setframe