GrADS繪圖軟件使用手

10其次章 GrADS繪圖模板初學(xué)者可以通過“初學(xué)指南”既參考GrADS軟件的“Tutorial”聯(lián)機(jī)幫助文開頭程序—>“tutorial”Tutorial”聯(lián)機(jī)幫助文件。tutorial的條件需要有一組演示數(shù)據(jù)：Model.le.dat(數(shù)據(jù)文件—二進(jìn)制)，Model.le.ctl(描述文件—ASCII碼)Model.grb，Model.ctl，Model.gmp(GRIB碼數(shù)據(jù))GRIB格式存儲(chǔ)的數(shù)據(jù)，而*.ctl文件〔ASCII碼〕是對(duì)應(yīng)數(shù)據(jù)文件的描述文件。GrADS用“數(shù)據(jù)文件GrADS中“翻開一個(gè)數(shù)據(jù)文件”既是指翻開一個(gè)數(shù)據(jù)描述文件。關(guān)于GrADS數(shù)據(jù)格式和用戶如何生成該格式的數(shù)據(jù)，將有特地章節(jié)表達(dá)；關(guān)于GRIB、NETCDF等數(shù)據(jù)格式的使用也將有特地論述。以下以第一組〔或其次組〕GrADS的使用。在使用數(shù)據(jù)前，先Model.le.ctl和Model.le.dat文件存于硬盤上的某一名目下，如c:\pcgrads\sample。首先來了解一下數(shù)據(jù)的內(nèi)容。Model.le.dat是一組模式輸Model.le.ctl清單：GrADS例如演示目的：GrADS的根本使用方法。GrADS的常用命令。GrADSxmanagerMS-DOS窗口。在MS-DOS窗口輸入命令：cdc:\pcgrads\sample再輸入命令:grads，顯示如下提示：首先顯示GrADS的一般信息，最底行提示是用“Landscape”模式grads-b以批處理方式運(yùn)行，不顯示圖形輸出grads-l模式運(yùn)行。grads-p模式運(yùn)行。grads-c進(jìn)入GrADSgrads-cl以上選項(xiàng)可聯(lián)合使用。〔11X8.5英寸〕或portrait〔8.5X11英寸—GrADS中長(zhǎng)度和大小單位都用英寸〕模式顯示圖形輸出窗口。輸入 L回車或回車是模式；輸入P回車是用portrait模式。之后進(jìn)入GrADS的命令交互模式 等待用戶輸入命令此時(shí)你的桌面應(yīng)如以下圖所示的直接進(jìn)入“Landscape”模式或grads-b以批處理方式運(yùn)行，不顯示圖形輸出grads-l模式運(yùn)行。grads-p模式運(yùn)行。grads-c進(jìn)入GrADSgrads-cl以上選項(xiàng)可聯(lián)合使用。左上角一個(gè)窗口是圖形顯示區(qū)，GrADS的全部圖形輸出結(jié)果在此窗口下顯示；MS-DOS窗口，此時(shí)正在運(yùn)行GrADga->，說明正處于GrADS命令等待狀態(tài)。GrADS的全部命令都只能通過該窗口輸出，全部文字信息也都由此窗口輸出。GrADS命令承受以下格式：ga->命令<參數(shù)<>>式”符號(hào)命令交互模既是在GrADSGrADS命令產(chǎn)生各種圖形。式GrADSga->quitga->openmodel.le.ctl進(jìn)入的第一步，翻開一個(gè)描述文件〔擴(kuò)展名ctlga->openmodel.le.ctl進(jìn)入的第一步，翻開一個(gè)描述文件〔擴(kuò)展名ctl可省略，GrADS據(jù)文件ga->qfilega->dps顯示GrADS翻開數(shù)據(jù)的內(nèi)容，顯示結(jié)果應(yīng)與描述文件說明全都。(ga->qq命令所代的參數(shù)及多項(xiàng)功能解釋。)示”命令ps”是model.le.ctl量名稱，代表地面氣壓。2緊接上例輸入以下命令：ga->c 去除圖面。如不去除圖面，GrADS后續(xù)顯示的圖形將與已存在的圖形產(chǎn)生疊加。ga->setlat40ga->setlon-90ga->setlev500ga->sett1GrADS視全部物理量為四維數(shù)據(jù)，以坐標(biāo)〔lon|lat〕自動(dòng)認(rèn)為是經(jīng)度/緯度坐標(biāo)。以x|y|z|t方式表示維數(shù)時(shí)，每一維都是一組從1開頭的序列數(shù)，對(duì)應(yīng)網(wǎng)格的序號(hào)。本例中：x從1到73〔代表從00到3600共72個(gè)格點(diǎn)第73個(gè)格點(diǎn)即第1個(gè)格點(diǎn)46個(gè)格點(diǎn)〕;z從1到7〔從1000hPa開頭，共7層〕;t從1到5〔5天〕運(yùn)行“open model.le.ctl”命令后，lon|lat或x|y是可變的〔lon:00~360;lat:-900~90。而垂直和時(shí)間維取固定值lev=1000hPa或z=1; 或t=1。維數(shù)參數(shù)數(shù)值1<2>”命令轉(zhuǎn)變當(dāng)前維數(shù)設(shè)定值。當(dāng)取“數(shù)值12”時(shí)，表示該維是變化的；而只取一個(gè)值時(shí)，表示該維取固定值，對(duì)于網(wǎng)格坐標(biāo)可以是分?jǐn)?shù)值，系統(tǒng)會(huì)自動(dòng)作內(nèi)插。維數(shù)參數(shù)始終ga->qdimsga->dz 勢(shì)高度。緊接上例輸入以下命令：ga->c緊接上例輸入以下命令：ga->c去除圖面。ga->setga->dlonz-1800顯示沿400N,，1800W至0度，500hPa層，1987.1.2.0hr位勢(shì)高度剖面圖。4ga->cga->setlat去除圖面。090ga->d zLon、lev、t設(shè)定值承受前例使用值不變。5ga->cga->setgaga->cga->setgaga->d去除圖面。tz15動(dòng)畫顯示西北半球1987.1.2.0hr位勢(shì)高度。1987.1.2.0hr至ga->c去除圖面。ga->setlat-9090ga->setlon-90ga->setlev1000100ga->sett1ga->dga->dtu溫度剖面圖。在上圖根底上再疊加上東西風(fēng)重量。ga->cga->setga->cga->setga->setga->setga->setga->d去除圖面。latlonlevtz40-180500015500hPa400N，高度的時(shí)間剖面。8ga->c去除圖面。ga->dsqrt(u*u+v*v)顯示500hPa全風(fēng)速值?；騡a->dmag(u,v)GrADS內(nèi)部定義了多種函數(shù)。9ga->c去除圖面。ga->dave(z,t=1,t=5)500hPa5天平均。ga->dz-ave(z,t=1,t=5)1987.1.2.0hr高度與平均的偏差。10ga->cga->dz-ave(z,x=1,x=72)去除圖面。高度與緯向平均值的偏差。11ga->c去除圖面。ga->dga->dz(t=2)-z(t=1)z(t+1)-z兩個(gè)時(shí)刻高度的差?；驁?zhí)12ga->c去除圖面。ga->dga->setz(lev=500)-z(lev=700)z 3兩個(gè)高度間的厚度?；驁?zhí)ga->dz(z+1)-zga->openga->openmodel.ctlga->openga->openmodel.ctlmodel.le.ctl同時(shí)翻開兩個(gè)文件。ga->dga->setz.2(lev=500)-z.1(lev=700)z 3兩個(gè)高度間的厚度。或執(zhí)ga->dz.1(z+1)-z.212所示下：變量名·文件序號(hào)〔維數(shù)參數(shù)變量名·文件序號(hào)〔維數(shù)參數(shù)+/-/=某一數(shù)值，…〕個(gè)文件，并為每個(gè)翻開的文件編一個(gè)序列號(hào)〔1開頭。上例中，序列號(hào)1lat/lon/lev/timex/y/z/t。下畫線局部可省或局部省略。對(duì)于有兩個(gè)不同文件中的數(shù)據(jù)作運(yùn)算時(shí)，兩種數(shù)據(jù)網(wǎng)格要全都對(duì)于有兩個(gè)不同文件中的數(shù)據(jù)作運(yùn)算時(shí)，兩種數(shù)據(jù)網(wǎng)格要全都。ga->reinitga->openga->reinitga->openGrADS狀態(tài)。model.ctlga->setga->setga->dlatlon090-1800vint(ps,q,275)vint(ps,q,top)=1topq的垂直質(zhì)量積分，計(jì)算可降水量。qdp, ps和top：hPa。gps15ga->resetga->setga->setga->setga->dopen命令后的全部設(shè)置。latlonlecurl(u,v)由風(fēng)場(chǎng)導(dǎo)出渦度場(chǎng)。ga->cga->setgxoutshadedga->cga->setgxoutshadedga->d hcurl(u,v)ga->’cbarn 1ga->setgxout0’contour去除圖面。以分色圖形方式輸出。由風(fēng)場(chǎng)導(dǎo)出渦度場(chǎng)。畫圖例。1:相對(duì)長(zhǎng)短，>1放大；0:水平；1:垂直。以等值線方式輸出〔缺省方式。ga->dz疊加500hPa高度場(chǎng)。留意疊加次序。ga->draw title 500hPaHeightsandVorticity寫圖標(biāo)題。18ga->18ga->cga->setgxoutstreamga->d u;v;q去除圖面。以流線方式表示矢量場(chǎng)。風(fēng)流線場(chǎng)?！?lt;>”內(nèi)的局部只起標(biāo)顏色的作用。ga->cga->gxoutvector去除圖面。以箭頭方式表示矢量場(chǎng)。ga->du;v;q風(fēng)矢量場(chǎng)顯示矢量時(shí)，d19ga->cga->setgxoutbarbga->d u;v;q去除圖面。以WMO風(fēng)標(biāo)方式表示矢量場(chǎng)。風(fēng)標(biāo)。192019ga->c去除圖面。ga->setgxoutgrid直接輸出網(wǎng)格點(diǎn)數(shù)值。ga->setdignum0設(shè)置保存小數(shù)位數(shù)。ga->setdigsize0.1設(shè)置數(shù)字大小〔英寸。ga->du東西風(fēng)網(wǎng)格點(diǎn)數(shù)值。202222ga->c去除圖面。ga->setgxout contourga->setlat 15 80ga->setlon-140-40ga->setlev500ga->setmpdraw on〔非經(jīng)緯度數(shù)據(jù)需此項(xiàng)〕ga->setpoli onoff不畫國(guó)界省界等。ga->setmap auto由系統(tǒng)自定地圖背景的顏色、線型和粗細(xì)。而set map color <line_style <thickness>>由用戶自定義。ga->set mpvals -140 -40 15 80 極射投影經(jīng)緯度范圍。ga->setmprojnps北半球極射投影(其它取值：latlon:等比經(jīng)緯度投影〔缺省scaled經(jīng)緯度投影并沖滿整個(gè)sps:Lambert:lambert方式投影;robinson:robinson式投影；off:關(guān)閉ga->setmpdset地圖投影，缺省是翻開的。)mres lowres/mres/hires分別代表可使用低分〔缺省ga->dzga->cga->cga->gxoutga->setlatga->setlonga->setlevga->mprojga->mpdset去除圖面。contourcaledmresga->d z23ga->c 去除圖面。ga->gxout contourga->setlat -90 ga->setlon -180 180ga->setlev 500ga->mproj robinsonga->mpdset mresga->d zga->setgxoutga->setgxoutvectorga->du;v;q結(jié)果如封面所示。ga->c去除圖面。ga->setgxoutga->setlat contour90ga->setlon-180180ga->setlev500ga->mpvalsga->mprojga->mpdsetga->d z-180npsmres1800-90GrADS繪圖模板上述使用方法是GrADS的根本使用方法，其特點(diǎn)是在GrADS系統(tǒng)提GrADS命令來完成繪畫。但很簡(jiǎn)潔造成輸入錯(cuò)誤，特別是有些命令設(shè)定后，假設(shè)不再重設(shè)置，是永久有效，而有些命令只是一次有效。假設(shè)用戶對(duì)初次繪畫的效果不滿足，要增加一些GrADS的根本方法使用會(huì)有需多重復(fù)命令，因此效所為“繪圖模板”就GrADS下以批處理方式執(zhí)行。這種文件被稱為“模板”或叫“描述語言”文件—“script比照上例編制一些簡(jiǎn)潔的模板設(shè)計(jì)。1md01.gs文件清單〔1：’open‘d；model.le.ctl’’open‘d；model.le.ctl’ps’﹡本例中，維數(shù)參數(shù)承受翻開數(shù)據(jù)文件后的缺省值。﹡模板文件中以“﹡”開頭的行是注解行。﹡命令必需用單引號(hào)括起來，行尾肯定要有回車。﹡DOS命令提示符下輸入命令：c:\pcgrads\sample\grads –cl sample01.gsc:\pcgrads\sample\grads –cl“sample0112…”(.gs擴(kuò)張名可以省略)c:\pcgrads\sample\grads –bcl sample01.gsGrADSGrADS命令提示符下輸入命令ga->runsample01.gs 12…或ga->sample01 12…〔在在GrADS命令提示符下，可反復(fù)執(zhí)行run命令。但考慮到之間可以先運(yùn)行“ga->reinit/resetrun〕GrADS在退出前，假設(shè)滿足，可把圖形存于文件。有多種方法：ga->printimfile.pngwhite存于名為“file.png”的文件(png格式)。printim命令格式：ga->printimfile選項(xiàng)1，選項(xiàng)2…選項(xiàng)：gif—gif存于“file”文件〔缺省為。png格式〕white—白色背景。black—黑色背景。xnnn如：ga->printimfile.gifgifwhitex800y600另一種方式是：ga->enableprintfile-name.gmfga->print執(zhí)行輸出，結(jié)果存于file-name.gmf—用戶指定的文件，metafile格式,最好代gmf后綴，ga->disableprintdisable命令后，print命令的結(jié)果才真正存metafile文件可用gv32查看或轉(zhuǎn)換成wmfgxps/gxeps轉(zhuǎn)成ps/epsmetafile和ps/eps格式圖象的印刷效果是最正確的。最終退出GrADSga->quit2‘reinit’’openc:/pcgrads/sample/model.le.ctl’DOS不同‘reinit’’openc:/pcgrads/sample/model.le.ctl’DOS不同DOS式的路徑。‘set‘set‘set*‘setlat 40’lon -180lev 500’t 1’0’*‘setccolor 0’*設(shè)定顏色[注1]。缺省取前風(fēng)光。*‘setcmark 3’*設(shè)定折線圖節(jié)點(diǎn)標(biāo)記[注2]。缺省取2。*‘setcstyle 1’*設(shè)定線型[注3]1。*‘setcthick 1’*設(shè)定線粗細(xì)[注4]1。*‘setgridon33’*設(shè)定是否畫網(wǎng)格線[注5]gridon?！甦z’‘printimsp02.png white’*把圖形存于文件[注]0：1：234：567‘printimsp02.png white’*把圖形存于文件2]：標(biāo)記取值：0：不作標(biāo)記；1：＋；2：○；3：●；4：□；5：■；6：X；7：

；8：△；9：不作標(biāo)記；10：①。3]：線型：1：實(shí)線；2：長(zhǎng)斷線；3：短斷線；4：一長(zhǎng)一短斷線；5：點(diǎn)線；6：點(diǎn)劃線；7：點(diǎn)點(diǎn)劃線。[注4]：線粗細(xì)取值：1-10〔個(gè)點(diǎn)陣。取值1-5〔含〕在屏幕上不顯示變化，5以上才顯示。因此，對(duì)于Printim，5變化。on <線型值> <顏色值>;畫網(wǎng)格并可指定線型和顏〔之一或全部或按缺省。gridhorizontal<線型值><顏色值>;只畫水平網(wǎng)格。gridvertical<線型值><顏色值>;只畫垂直網(wǎng)格。gridoff不畫網(wǎng)格。[注6]以“*”開使的命名行表示用戶可選擇是否使用。很多設(shè)置可首先考由GrADS自行打算，如不滿足，再由用戶自定義，以免畫蛇添足。3’open‘set‘set‘setmodel.le.ctl’lat’open‘set‘set‘setmodel.le.ctl’latlonlev090’-180500’0’*‘set*‘set*‘set*‘set*‘set*‘setcsmoothrgbrgb1617ccoloron’15622233’156234133’16’用戶自定義顏色[注1]。clab %.1f’clskip 2等值線標(biāo)記方式[注2]。每隔一條等值線標(biāo)記數(shù)值。*‘setcint8’指定等值線間隔。*‘setcterpon’樣條平滑“ON”或“OF。*‘setclopts30.10.2’等值線標(biāo)記的顏色<粗細(xì)<大小>>‘dz/9.8’‘‘printim；md03.gufgif white’[注1]：‘set rgb 用戶自定義顏色序號(hào)〔16到99—1到15系統(tǒng)定義〕紅綠藍(lán)‘“紅綠藍(lán)”取值范圍0-255[注2]：set clab on〔等值線標(biāo)數(shù)值〕/off〔不標(biāo)記〕/forded〔強(qiáng)制標(biāo)記加上字符保存%03·1f〔33位，前部用0補(bǔ)齊〕模板4’open‘set‘set‘setmodel.le.ctl’lat’open‘set‘set‘setmodel.le.ctl’latlonlev090’-180 0’500’*畫指定值的等值線?！畇et‘set‘d；clevsccolsz/9.8’495523534556560564572584588’123’可以分別指定以上每一條線的顏色。模板5md05.gs文件清單：’open’open‘setmodel.le.ctl’lat090’*‘setcmin0’不畫低于0〔含〕以下的等值線。*’set‘dcmax 0’ 0〔不含〕以上的等值線。hdivg(u,v)*1.e5’ 散度*‘setblackoff’‘setclevs0’‘setccols1’‘setcthick8’‘set‘set‘set‘set‘set‘setlonlevblack-180 0’500’-0.10.1’跳過-0.10.1〔0線〕不畫。black“d“of特別處理“0”線畫法?！甦；hdivg(u,v)*1.e5’’open‘set‘set‘set*‘set’open‘set‘set‘set*‘set‘set*‘setc:/pcgrads/sample/model.le.ctl’latlonlevtccolorgrid40’-180 0’500’1’3’on33’以直方圖方式輸出。直方圖間隔〔20。在直方圖中的矩形中填色outline”只畫矩形框而不填色。直方圖中的矩形從底部向上畫to”從‘setgxoutbar’‘setbargap20’‘setbaroptsfilled’‘setbarbasebottom’頂部向下畫畫。*GrADSGrADS*2D坐標(biāo)，X軸指垂直軸。*與數(shù)據(jù)中定義的4D坐標(biāo)是兩個(gè)概念?！畇etxaxis-180010’‘setyaxis51505850 50’‘setylint200’‘setxlint20’‘setxflipon’‘setyflipoff’‘setxlpos0 b’‘setannot8 8’‘setxlopts‘setxaxis-180010’‘setyaxis51505850 50’‘setylint200’‘setxlint20’‘setxflipon’‘setyflipoff’‘setxlpos0 b’‘setannot8 8’‘setxlopts4 80.2’‘d；z’X軸翻轉(zhuǎn)。[注3]X軸位置。4]坐標(biāo)軸顏色和線粗細(xì)。[注5]標(biāo)記特性。6][注1]:set xaxis 開頭完畢<間隔>/set yaxis 開頭完畢<間隔；設(shè)置X軸/Y軸標(biāo)記范圍，留神“開頭完畢開頭—束范圍要與維數(shù)參數(shù)的范圍一set lon -180 如不一制圖形數(shù)據(jù)將沒有代表，此時(shí)你給什么，它標(biāo)什么，而圖形不變。如不能預(yù)先定出范圍，可由先由GrADS 自己打算或調(diào)整間隔〔xlint/ylin使用xaxis/yaxis的狀況是當(dāng)你處理非經(jīng)緯度數(shù)據(jù)時(shí)，GrADS總是把水平坐標(biāo)處理成經(jīng)緯度來標(biāo)記。[注2]：set xlint/ylint 間隔，如在此設(shè)置間隔，將取代xaxis/yaxis設(shè)置的間隔。[注3]：X/Y軸，同時(shí)圖像也作相應(yīng)的翻轉(zhuǎn)。[注4]：set xlpos offset〔+/-英寸〕b〔底部〕/t〔頂部〕set ylpos offset〔+/-英寸〕l〔左〕/r〔右〕[注5]：set annot 顏色〈粗細(xì)；邊框的顏色和線粗細(xì)。[注6]：set xlopts/ylopts 顏色〈粗細(xì)〈大小—英寸；標(biāo)記數(shù)字的顏色、線粗細(xì)和字大小。7’open‘set’open‘set‘set‘set‘setc:/pcgrads/sample/model.le.ctl’latlonlevt160’42’10003’100’*分左右放置兩幅圖‘set‘set‘set‘set‘set‘set‘dvxyrevannotxloptsyloptszlogt’05.506’on’8413on’8’120.15’0.2’設(shè)置虛頁面大小〔英寸。[注1]坐標(biāo)軸特性將始終保持。標(biāo)記特性將始終保持。Z坐標(biāo)取對(duì)數(shù)坐標(biāo)〔將始終保持?！畇et‘set‘set‘set‘d；vxlintxloptsyloptst’5.510’1341106’120.15’0.2’標(biāo)記特性。[注1]：一個(gè)物理頁是指：水平11英寸X垂直8.5英寸(Landscape)或18.5X11(portrait)，其上可以畫分出多塊區(qū)域，實(shí)現(xiàn)多圖同時(shí)輸出。Set v 水平開頭 完畢 垂直開頭 完畢8’open‘set’open‘set‘set‘set‘setc:/pcgrads/sample/model.le.ctl’latlonlevt0160’42’700’1’*分左右放置兩幅圖。左圖中Y軸尺度由GrADS自定。‘set v 0 5.5 0 6’連接斷點(diǎn)。設(shè)置Y〔vrangeX軸的尺度。*‘setxyrevon’‘setannot8 8’‘setxlopts4 10.2’‘set‘d‘set‘dyloptst’yloptsz/9.8’1 23’0.2’‘setv5.511 0 6’*‘setxyrevon’‘setmissonnon’‘setylopts1’‘d‘dvranget’z/9.8’240 330’‘cbar_l‘cbar_l-x2.5-y1.5–t“temp”“heigh”’標(biāo)圖例。-x2.5-y1.5位置；–t“temp“heigh”標(biāo)名稱。md09.gs文件清單：md09.gs文件清單：‘set*‘set‘set*‘set‘run‘set‘setv 3803.4’arrlabon’onarrowhead-0.2’arrsclsp17’vstrmden0.1’17。3 83’3.46.8’調(diào)整流線密度〔1稀—10疏?！畆un‘setsp18.gs’v 3 8 6.811’18。*‘setdigsize 0.04’調(diào)整風(fēng)向標(biāo)大小?！畆unsp19.gs’19?！畃rintimmd09.png white’sp17.gssp17.gs文件清單”open”set”set”set”set”d*’dc:/pcgrads/sample/model.ctl”lonlatlev-180 0”090”500”gxoutu;v;q”vector”skip(X_skip,Y_skip),通過X_skip,Y_skip”XY方向的取樣密度〔1可省略不給u,vXYvector或sp19.gssp19.gs文件清單”open”set”set”set”setc:/pcgrads/sample/model.ctl”lonlatlev-180 0”090”500”gxoutbarb””d u;v;q”barbskip函數(shù)調(diào)整密度。sp18.gssp18.gs文件清單”open”set”set”set”setc:/pcgrads/sample/model.ctl”lon-180 0”latlev090”500”gxoutstream””d u;v;q”10md10.gsmd10.gs文件清單：‘open‘set‘set‘setmodel.le.ctl’lonlatlev0 180’30’500’”define zv=ave(z/9.8,lat=20,lat=40)”定義zv20~40N緯度帶間的平均位勢(shì)高度?！眘et gxout linefill”**setlfcols23’”d z/9.8;zv”z/9.8zvz/9.8<zv〔z/9.8>zv時(shí)用綠色〔。md11.gs文件清單：md11.gs文件清單：‘open‘set‘set‘set‘setmodel.le.ctl’lonlatlev00180’90’500’rbcols914411513310712826’設(shè)置彩虹色。‘d z/9.8’*以下給出GrADS一組根本繪圖指令的使用方法。‘*第一組：‘set font2’選擇字型庫(kù)〔０－５〕顏色、位置、粗細(xì)、角度。字水平大小、垂直大小?！甦rawstring83 firtgroup’。位置可先運(yùn)行‘qpos’命令打算?！甦raw title first line\second line’“”起分行作用?！痵etstring1c245’‘setstrsiz0.30.6’‘drawline3457’‘drawrec3456’‘drawrecf2133’*其次組：‘set line346’‘draw‘draw‘drawmark 10polyf 5wxsym 26177 0.5’6 2 6設(shè)置線的顏色、線型、粗細(xì)到〔5,7〕畫直線。畫矩形(3,4),(5,6)對(duì)角點(diǎn)。畫填充色矩形〔當(dāng)前線色〕=5。70.54’-1畫多邊形8’*其次組：‘set line346’‘draw‘draw‘drawmark 10polyf 5wxsym 26177 0.5’6 2 6設(shè)置線的顏色、線型、粗細(xì)到〔5,7〕畫直線。畫矩形(3,4),(5,6)對(duì)角點(diǎn)。畫填充色矩形〔當(dāng)前線色〕=5。70.54’-1畫多邊形8’顏色為-18的氣象符號(hào)=2。set strsiz 字水平大小 垂直大小>set line 線的顏色<線型 粗細(xì)>>set string 字串的顏色 <位置 <粗細(xì) <角度>>>位置取以下值：tl/tc/tr/l/c/r/bl/bc/br其中，c代表以起點(diǎn)位置居中；r從起點(diǎn)向右，等等。角度取正值逆時(shí)針，取負(fù)順時(shí)針轉(zhuǎn)。draw wxsym 氣象符〔1-41〕x y 大小<顏色 <粗細(xì)>>顏色取負(fù)值時(shí)表示將由系統(tǒng)自定。預(yù)知每個(gè)符號(hào)的圖形，運(yùn)行runwxsym.gs，結(jié)果如以下圖：GrADS模板的高級(jí)應(yīng)用模板12‘openmodel.le.ctl’21-24，設(shè)計(jì)通用地圖投影圖模板?！畂penmodel.le.ctl’21-24，設(shè)計(jì)通用地圖投影圖模板。lon1=-140;lon2=-40。lat1=15;mp=nps*mp=spslat2=80?。〔恍栌脝我?hào)括起來*mp=lambert*mp=robinson‘set lev 500’‘set‘set‘set‘sett 1’lon ’lon1‘’lon2lat’lat1‘’lat2只有GrADS命令局部需要用單引號(hào)括起。?！甦efinevar=zzzz=z/9.8’ 必需在setlon/lat后定義留意以上兩種“變量定義”方式的不同，不能直接定義*maps(zz,mp,lon1,lon2,lat1,lat2)maps(var,mp,lon1,lon2,lat1,lat2)件中。調(diào)用函數(shù)?。『瘮?shù)必需在同一文functionmaps(vars,mpj,ln1,ln2,lt1,lt2)定義函數(shù)。‘set‘set‘setmpdrawon’polimap*’setmapoff’auto’556’*‘setmpvals‘set mproj‘set mpdset’ln1‘’ln2‘’lt1‘’lt2’mpjmres’‘display;’vars函數(shù)和調(diào)用不在同一文件中‘openmodel.le.ctl’21-24，設(shè)計(jì)通用地圖投影圖模板。lat1=15;mp=npslat2=80！！不需用單引號(hào)括起來*mp=lambert‘set‘set‘set‘setlevt 1’lon500’’lon1‘’lon2lat’lat1‘’lat2只有GrADS命令局部需要用單引號(hào)括起。?！甦efinevar=zzzz=z/9.8’ 必需在setlon/lat后定義*‘md132 ’var‘’mp‘’lon1‘’lon2‘’lat1‘’lat2‘md132 ’var‘’mp??！調(diào)用函數(shù)不再是以函數(shù)名，而是以函數(shù)所在的文件名。參數(shù)可以從后向前省略;md132.gsmd132.gs文件清單—函數(shù)局部：function maps(args)定義函數(shù)。args,代表了全部參數(shù)。vars=subwrd(args,1) 取args的第1個(gè)字。mpj=subwrd(args,2) 取args的第2個(gè)字ln1=subwrd(args,3)ln2=subwrd(args,4)lt1=subwrd(args,5)lt2=subwrd(args,6)‘set‘set‘setmpdrawon’polimap*’setmapoff’auto’556’*‘setmpvals ’ln1‘’ln2‘’lt1‘’lt2‘set mproj ’mpj‘set mpdset mres’‘display;’var設(shè)計(jì)在一頁紙上輸出多幅圖功能的模板md14.gs文件清單：FunctionFunctionall_in_one_(args)*setgradsoff”setimprundefault.gs”在每次”d”命令前運(yùn)行default.gs模板。usage1=“putseveraltempletsinonereal.arrangthemintwocollons/rows.“usage2=“md14<file.gmf>>“templet1templet2...<-printim<file.png>>/<-printgeteachsinglemapfromcommandlineargsi=0pt=0out_file=””wrd=subwrd(args,1)if(wrd=””);sayusage1;say“usage:“usage2;return;endifwhile(wrd!=””)getoutputoptionsif(wrd=”-printim”)pt=”-printim”out_file=subwrd(args,i+2)breakendifif(wrd=”-print”)pt=”-print”out_file=subwrd(args,i+2)breakendifgettempletnamesmd.i=wrdsaymd.i” i=”ii=i+1wrd=subwrd(args,i+1)endwhileputthelasttempletaszeromd.i=0;num=math_nint(i/2)*2say“numberoftemplets“i“in“num“virtualsintwocollons/rows.“say“outputoption-pt=“pt“out_file=“out_filegetinformation”querygxinfo”result代表querygxinfo命令產(chǎn)生的一個(gè)多行輸出結(jié)果。rec2=sublin(result,2) sublin(result,2)取result的第2行=rec2xsiz=subwrd(rec2,4) 取rec2的其次個(gè)字。ysiz=subwrd(rec2,6)*say“width=“xsiz“andhiegh=“ysizif(ysiz<xsiz)=“l(fā)andscape“else=“portrait“endif*setvandploteachsinglemapineachvirtual.if(=“l(fā)andscape“)sayif(num<4)i=0x1=0x2=xsiz/numwhile(i<num)y1=0y2=ysiz”setv”x1””x2””y1””y2”run”md.isay”setv”x1””x2””y1””y2x1=x2x2=x2+xsiz/numi=i+1endwhileif(0) drawaboxoneachvirtul”setvoff”i=1x1=xsiz/numwhile(i<num)”drawline”x1”0”x1””ysizi=i+1x1=x1+xsiz/numendwhileendifelsei=num/2x1=0x2=xsiz/iy1=ysiz/2y2=ysizk=0while(k<i)”setv”x1””x2””y1””y2*say”setv”x1””x2””y1””y2”run”md.k”setv”x1””x2””0””y1*say”setv”x1””x2””0””y1j=k+iif(md.j);”run”md.j;endif* say“ “x1=x2x2=x2+xsiz/ik=k+1endwhileif(0)drawaboxoneachvirtul”setvoff””drawline0”ysiz/2””xsiz””ysiz/2k=1x1=xsiz/iwhile(k<i)”drawline”x1”0”x1””ysizk=k+1x1=x1+xsiz/iendwhileendifendifendifif(=“portrait“)sayif(num<4)i=0y1=0y2=ysiz/numwhile(i<num)x1=0x2=xsiz”setv”x1””x2””y1””y2say”setv”x1””x2””y1””y2j=num-1-i”run”md.jy1=y2y2=y2+ysiz/numi=i+1endwhileif(0) drawaboxoneachvirtul”setvoff”i=1y1=ysiz/numwhile(i<num)”drawline0”y1””xsiz””y1y1=y1+ysiz/numi=i+1endwhileendifelsei=num/2x1=0x2=xsiz/2y1=0y2=ysiz/ik=0j=num-1while(k<i)”setv”x1””x2””y1””y2say”setv”x1””x2””y1””y2l=j-1”run”md.l”setv”x2””xsiz””y1””y2say”setv”x2””xsiz””y1””y2if(md.j);”run”md.j;endifsay“ “y1=y2y2=y2+ysiz/ij=j-2k=k+1endwhileif(0) drawaboxoneachvirtul”setvoff””drawline”xsiz/2”0”xsiz/2””ysizk=1y1=ysiz/iwhile(k<i)”drawline0”y1””xsiz””y1y1=y1+ysiz/ik=k+1endwhileendifendifendif”setvoff””drawrec00.03”xsiz-0.01””ysiz”drawrec00.0”xsiz ””ysizif(pt=”-printim”)if(out_file!=””)”printim”out_file”white”else”printimtmp.pngwhite” ”printimtmp.gifgifwhite”endifreturnendifif(pt=”-print”)if(out_file!=””)”enableprint”out_fileelse”enableprinttmp.gmf”endif”print””disableprint”returnendifsay“entertheprintim<file.png>orprint<file.gmf>forprintouttoafile. “say“ORJUSTPRESSENTERKEYFORNOTHING. “pullargwrd=subwrd(arg,1)if(wrd=”printim”)out_file=subwrd(arg,2)if(out_file!=””)”printim”out_file”white”else”printimtmp.pngwhite” ”printimtmp.gifgifwhite”endifreturnendifif(wrd=”print”)out_file=subwrd(arg,2)if(out_file!=””)”enableprint”out_fileelse”enableprinttmp.gmf”endif”print””disableprint”returnendif*say”out_file=”out_filereturn;‘setgrads‘setgradsoff’;執(zhí)行模板：dos>grads–cl“md14sp21sp22sp23sp24sp16”grads后執(zhí)行：ga->md14 sp12sp22sp23sp24sp16處顯示出由于多個(gè)模板同時(shí)運(yùn)行，某些設(shè)置相互干擾〔單獨(dú)運(yùn)行都是正常的。圖中圖片大大小小，規(guī)格不能統(tǒng)一。如何準(zhǔn)確掌握?qǐng)D片大小是個(gè)難題。setimprundefault.gsdisplay命令前執(zhí)行一個(gè)模板default.gs或其它display命令后，setgradson自動(dòng)打開。GrADS描述語言GrADS描述語言GrADS的一行寫一變量間需要用空格分開〔GrADS單/1個(gè)或多個(gè)空格單/雙引號(hào)的方式來。定義變量變量名〔1~8字母數(shù)字〕=數(shù)值/“字串”/‘字串’名2=變量名1 使用變量。while(表達(dá)式)變量名·i=值 定義數(shù)組“”取一數(shù)值范圍。endwhile以下畫線開頭的變量為全局變量。_variable=123 變量才能超越函數(shù)，超越文件使用，但不能作函數(shù)的參數(shù)。以下變量名用戶不要定義，為系統(tǒng)保存字。lat,lon,lev,result,recfunction語句方式1——函數(shù)過程體與函數(shù)過程的調(diào)用在同一文件中定義函數(shù)過程functionname(variable1,variable2,…)緊接函數(shù)過程體；〔函數(shù)過程體完畢標(biāo)志〕函數(shù)調(diào)用過程語句…name(variable1,variable2…)name 假設(shè)沒有定義變量語句…目前GrADS要求函數(shù)過程的調(diào)用與函數(shù)過程體在同一文件中。方式2——函數(shù)過程體與函數(shù)過程的調(diào)用不在同一文件中定義函數(shù)過程函數(shù)過程文件清單functionname(variable)緊接函數(shù)過程體；〔函數(shù)過程體完畢標(biāo)志〕函數(shù)調(diào)用過程函數(shù)調(diào)用過程文件清單：語句…件名參數(shù)1‘’參數(shù)2‘’…語句留意單引號(hào)和強(qiáng)制參加空格運(yùn)算符和表達(dá)式數(shù)學(xué)運(yùn)算：+，-，*，/規(guī)律運(yùn)算：=,!=,>,>=,<,<=, |,&,!,%(連接符)表達(dá)式即由運(yùn)算符連接起來的變量，其結(jié)果為0—假或非0—即真。掌握語句if(express)語句elseendifif(express);語句;語句;endifwhile(express)語句endwhile內(nèi)部函數(shù)subwrd(string,nwd)從stringnwd〔數(shù)字〕個(gè)字。subline(string,nl) 從stringnl〔數(shù)字〕行一整行。Substr(string,start,nl)從字串中取從start字開頭nl長(zhǎng)的字串。雜項(xiàng)輸出line=”PeterPan,theflyingone”;say“Shesaiditislineprompt“Shesaiditis“l(fā)ine輸入pullvariable遇到pull語句后，程序等待用戶輸入，variable將代表用戶的輸入的值。query/q<參數(shù)>命令是一個(gè)格外重要的命令，用它可以獲得很多圖原信息。以下片段說明其一種用法及與描述語言其它局部信息交換確實(shí)方式?！眖uerygxinfo”result代表querygxinfo命令產(chǎn)生的一個(gè)多行輸出結(jié)果。rec2=sublin(result,2) sublin(result,2)取result的第2行=rec2xsiz=subwrd(rec2,4) 取rec2的其次個(gè)字。ysiz=subwrd(rec2,6) 〔xsiz,ysiz〕是圖形左上角的坐標(biāo)〔英寸。GrADS高級(jí)模板的應(yīng)用pcgrads\lib以作為有用的工具和學(xué)習(xí)參考。由于該名目是一個(gè)特別名目，因此這些工具可ascii碼文件，具體使用方法可直接參考文件說明。basemap.gs—陸地或海洋填色背景圖Usage:basemapL(and)/O(cean)<fill_color><outline_color><hi/lo>調(diào)用數(shù)據(jù)：lpoly.asc/lpoly_hires.asc;opoly.asc/opoly_hires.ascbasemap.gs文件，在數(shù)據(jù)文件前加上路徑?！眔pen”openc:/pcgrads/sample/model.ctl””setlon-180180””setlat-9090””setlev500””sett 1””setmpdrawoff” 不畫海陸線”setcmin9999” 用basemap前要求畫任意圖，”dz/9.8” setcmin9999使等值線并不顯示出來。D命令后自動(dòng)矢效。”basemapl3” 給陸地填綠色。”basemapo4” 海洋填藍(lán)色?！眘etccolorrainbow” 設(shè)為彩虹色，d命令后彩虹色自動(dòng)矢效?！眃z/9.8”;cbar_l.gs(cbar_line.gs)—為線條圖畫圖例。Usage:cbar_l<-xX-yY-n#>-ttext<p>見模板8，cbar_line.gs與其作用一樣。cbarn.gs(cbar.gs,cbarc.gs)—為填色圖畫圖例Usage: cbarn<sf<vert<xmidymid>>>cmap.gs—調(diào)整顏色gs->cmap numberconnect_the_dots.gs—通過點(diǎn)擊鼠標(biāo)畫折線。define_colors.gs—定義顏色值：21~79defval_demo.gs—演示“qdefval”和“setdefval…”的使用?？蓪?duì)壞點(diǎn)數(shù)據(jù)進(jìn)展修改。font.gs—演示字庫(kù)樣本。gs->font <0/1/2/3/4/5>isen.gs—等熵面圖。md16.gs”openc:/pcgrads/sample/model.ctl””setlon50 140””setlat0 90””sett 1””setlev1000100””defineue=”isen(u,t,lev,320)320K等熵面層風(fēng)場(chǎng)ue,ve;u,v,t三維要素場(chǎng).”defineve=”isen(v,t,lev,320) lev氣壓層(hPa)”set z 1””d ”isen(z,t,lev,320) 320等熵面高度.”setgxoutbar””due;ve;q”以下要把isen.gs文件的內(nèi)容包括內(nèi)。functionisen(field,tgrid,pgrid,tlev)．．map.gs—設(shè)置地圖投影方式。mconv.gs—計(jì)算水氣輻合。panels_demo.gs,panel.gsf—演示在一個(gè)物理頁上輸出多幅圖的方法，_vpg.panel”openc:/pcgrads/sample/model.ctl””set”openc:/pcgrads/sample/model.ctl””setlon50140””setlat0 90””setlev1000100”*”definet225=”pinterp(t,lev,225)”sett 1””setz1”與上例一樣,這里總有這樣一驚異的維數(shù)設(shè)定?！眃”pinterp(z,lev,225) 225hPa高度?！眃efinet225=”pinterp(t,lev,225)”dt225” 225hPa溫度以下要把以下要把pinterp.gs文件的內(nèi)容包括內(nèi)。functionpinterp(field,pgrid,plev)．．md18.gs”openc:/pcgrads/sample/model.ctl””setlon140””setlat30”*definethevariablevirticalscale”setlev1000100””sett 1””querydims”rec=sublin(result,4)ztype=subwrd(rec,3)if(ztype=“varying“)z1=subwrd(rec,11)z2=subwrd(rec,13)elsesay“youmustsetlevvariable...quit“returnendif”definetc=t-273.16””definetd=”q2td(q,tc,lev,z1,z2)*changem/stoknot”definewsp=mag(u,v)*1.943844””definedir=atan2(v,u)*180/3.1415926”*extradefinethevirticalscaleforthet_skewmap”setlev1050200”*plotskew(-1,-1,-1,-1)plotskew(tc,td,wsp,dir)functionq2td(q,tc,plevs,z1,z2)*functiontoconvertspecifichumidityintodewpointtemperature*qkg/kg(afunctionofp)*tccelsius(afunctionofp)*plevshPa*tdcelsius(returndewpointtemperature)”definet=”tc”+273.16””definep=”plevs”*100”rd=287.05;rv=461.51;tmelt=273.16;c1es=610.78;c2es=rd/rv;c3les=17.269;c4les=35.86;c3ies=21.875;c4ies=7.66;rvd=rv-rd*sayrd””rv””tmelt””c1es””c2es””c3les””c4les””c3ies””c4ies””rvd”definealfa=log(”rv”*p*q/(”c1es”*”rd”+q*”rvd”))”*”dalfa”if(”t”<tmelt)c3es=c3iesc4es=c4ies”definetdk=(”c3es”*”tmelt”-”c4es”*alfa)/(”c3es”-alfa)”elsec3es=c3lesc3es=c3lesc4es=c4les”definetdk=(”c3es”*”tmelt”-”c4es”*alfa)/(”c3es”-alfa)”endif“definetd=tdk-273.16“*”dtd”return(td);以下要把plotskew.gs文件的內(nèi)容包括內(nèi)。functionplotskew(sndtemp,snddewp,sndspd,snddir)。。script_math_demo.g—演示數(shù)學(xué)函樹的使用。其中很多math_開頭的在很多手冊(cè)中并未提到，但值得一。string.gs—寫字串。ga->stringstring 位置用鼠標(biāo)點(diǎn)擊圖上的位置即可。wxsym.gs—演示氣象符號(hào)”openc:/pcgrads/sample/model.ctl””setlon50”openc:/pcgrads/sample/model.ctl””setlon50140””setlat0 90””setlev1000100””sett 1””setz1””d”d”zinterp(lev,z,6000)”definet225=”zinterp(t,z,6000)6000米等高面上的氣壓場(chǎng)。6000米等高面上的溫度場(chǎng)?！眃t225”以下要把zinterp.gs文件的內(nèi)容包括內(nèi)。functionzinterp(field,zgrid,zlev)．．xanim.gs—增加的動(dòng)畫顯示工具。use.gs—增加的open命令.traj.gs—計(jì)算水平向前和向后的軌跡。sweat_index.gs—計(jì)算出汗指數(shù)。cress_section.gs—依據(jù)defval_demo.gs設(shè)計(jì)的畫剖面圖工具。cress_section.gs”openmodel.ctl””setimprundefault.gs”*getthedimensionmaximumlimilations.”queryfile”rec =sublin(result,5)_xsize=subwrd(rec ,3)_ysize=subwrd(rec ,6)_zsize=subwrd(rec ,9)*say_xsize””_ysize””_zsize”setx1”_xsize”sety1”_ysize”setz1”_zsizegetinformation”querygxinfo”rec =sublin(result,2)_pxsiz=subwrd(rec,4)_pysiz=subwrd(rec,6)*say“width=“_pxsiz“andhiegh=“_pysiz******************************** changebelow ************************ifyouwanttoprintatestdata,setthisto1orhaveonleonegraphicarea;youcanoutputthecrosssectionmapbesideahorizontialfielsbytosetprint_grid=0print_grid=0giveavariableandaleveltoindicatwherethecross_sectionwillbeonthemap.zlev=500”definedisplay=z(lev=”zlev”)/9.8”givethevariablestoindictwhichvariablesyouwanttomakethecross_sectmap.nowonlyworkwithonevariable,thisshouldbeupdattoservalvariables”definevarb =t-273.16”*”definevarb2=*”definevarb3=*definethelimitationfortheinterestingarea.lon1=70;lon2=140;lat1=15;lat2=60;lev1=1000;lev2=100*ifyouusetheworldcoordinatesgiveworld=1;forgridcoordinatesworld=0world=1******************************** changeabove ***************************???Problemwhitvarpasstofunction???????*cross(zlev,dpv,var,70,140,15,60,1000,100,1);*functioncross(zlev,display,varb,lon1,lon2,lat1,lat2,lev1,lev2,world)assign11(odd)gridpointsalongthecross_sectionx-axis.youcanchangethisvaluebyassignnp=newnumbertohavemoregridpoints.np=11np1=np-1restrictthedimensionenvironmenttotheinterestingareaif(world)converttheworldcoordinatestogridoneif(world)”setlon”lon1-0.5””lon2+0.5”setlat”lat1-0.5””lat2+0.5”setlev”lev1””lev2”querydims”ln2=sublin(result,2);ln3=sublin(result,3);ln4=sublin(result,4)lon1=subwrd(ln2,11);lon2=subwrd(ln2,13);lat1=subwrd(ln3,11)lat2=subwrd(ln3,13);lev1=subwrd(ln4,11);lev2=subwrd(ln4,13)lon1=math_int(lon1); lon2=math_int(lon2) ;lat1=math_int(lat1)lat2=math_int(lat2); lev1=math_int(lev1); lev2=math_int(lev2)”setx”lon1-0.5””lon2+0.5”sety”lat1-0.5””lat2+0.5”setz”lev1””lev2else”setx”lon1-0.5””lon2+0.5”sety”lat1-0.5””lat2+0.5”setz”lev1””lev2endifthe“defval“commandswillonlyworkwith2-Ddefinedvariables”set lev”zlev”definevar=”varbDisplaythevariablewithcoloredgridcellsandtheirvaluesif(print_grid)”setgxoutgrid””setgridoff””setmprojoff””setdigsiz0.18”else”setmprojscaled””setgxoutcontour”endifsettwoareatoshowthemaps,thelifeshowahorizontialfieldandcrosssectionlineandtherightshowthecoress_sectionmapsforlandscapeortop/bottomforportrait.if(!print_grid)if(_pysiz<_pxsiz)landscape“”setv0”_pxsiz/2” 0 ”_pysizelseportrait“”setv0”_pxsiz” ”_pysiz/2””_pysizendif”display”displayelse”d”var”setgxoutcontour””d”var ”printimcr_0.pngwhite””display”displayendif*”qdims”;sayresult*Usethemousetoclickonagridpoint.say”Clickonanygridpointtogetacrosssectionstartpoint””qpos”xscreen_s=subwrd(result,3)yscreen_s=subwrd(result,4)Convertvirtualpositionstoscreenpositions”qpp2xy”xscreen_s””yscreen_sxscreen_s=subwrd(result,3)yscreen_s=subwrd(result,6)Convertscreenpositionstogridcoordinates”qxy2gr”xscreen_s””yscreen_sxg.0=subwrd(result,3)yg.0=subwrd(result,6)Convertgridcoordinatestoworldcoordinates”qgr2w”xg.0””yg.0lon.0=subwrd(result,3)lat.0=subwrd(result,6)say”Clickonanygridpointtogetaendpoint””qpos”xscreen_e=subwrd(result,3)yscreen_e=subwrd(result,4)Convertvirtualpositionstoscreenpositions”qpp2xy”xscreen_e””yscreen_exscreen_e=subwrd(result,3)yscreen_e=subwrd(result,6)Convertscreenpositionstogridcoordinates”qxy2gr”xscreen_e””yscreen_exg.np1=subwrd(result,3)yg.np1=subwrd(result,6)Convertgridcoordinatestoworldcoordinates”qgr2w”xg.np1””yg.np1lon.np1=subwrd(result,3)lat.np1=subwrd(result,6)if(0)say”(xscreen_s,yscreen_s)=(”xscreen_s””yscreen_s”)”say”(xscreen_e,yscreen_e)=(”xscreen_e””yscreen_e”)”say”(xg.0,yg.0) =(”xg.0””yg.0”)”say”(xg.np1,yg.np1)=(”xg.np1””yg.np1”)”say”(lon.0,lar.0) =(”lon.0””lat.0”)”say”(lon.”np1”,lat.”np1”) =(”lon.np1””lat.np1”)”say””endif*drawalineasthecrosssectionline”drawline”xscreen_s