GrADS繪圖軟件使用手冊

GrADS繪圖軟件實(shí)用手冊2002年1月

目錄GrADS繪圖軟件概述GrADS繪圖軟件簡介Internet上的GrADS資源GrADS在Internet上的主頁windows環(huán)境下GrADS資源GrADS繪圖軟件的安裝（windows環(huán)境）在windows環(huán)境下安裝GrADS軟件包Xserver的安裝GrADS繪圖模板GrADS示例演示1.1啟動(dòng)GrADS1.2退出GrADS1.3示例演示GrADS命令的使用GrADS繪圖模板GrADS模板的高級(jí)應(yīng)用GrADS描述語言GrADS高級(jí)模板的應(yīng)用GrADS數(shù)據(jù)格式格點(diǎn)數(shù)據(jù)描述文件1.1數(shù)據(jù)描述文件各項(xiàng)解釋1.2生成和文件的程序代碼片段站點(diǎn)數(shù)據(jù)的格式附錄如何精確控制圖形輸出的尺寸—Landscape紙型臺(tái)站資料的顯示Linux環(huán)境下的安裝

GrADS繪圖軟件概述GrADS繪圖軟件簡介TheGridAnalysisandDisplaySystem(GrADS)是一套應(yīng)用廣泛、使用方便的科學(xué)數(shù)據(jù)繪圖軟件包。其主要特點(diǎn)：GrADS屬于自由軟件，可以從Internet上免費(fèi)獲得。可運(yùn)行于各種Windows和Unix工作平臺(tái)。GrADS可用于4D數(shù)據(jù)的分析。既經(jīng)度、緯度、層（氣壓層、高度層等）和時(shí)間/xyzt4維。數(shù)據(jù)可以是格點(diǎn)化的數(shù)據(jù)或離散點(diǎn)數(shù)據(jù)。GrADS特別適用于氣象類數(shù)據(jù)的分析。但也完全可以用于更廣泛類型的數(shù)據(jù)分析。GrADS有多種顯示方式：等值線、流線、矢量圖、風(fēng)矢量圖、站點(diǎn)填圖、折線圖、直方圖等多種兩維圖形?？商幚矶喾N數(shù)據(jù)格式的數(shù)據(jù)。GRIB、NetCDF、HDF-SDS等通用數(shù)據(jù)格式和系統(tǒng)自定義的一種二進(jìn)制數(shù)據(jù)格式。采用命令行輸入的方式交互式地顯示圖形。并有多種命令對數(shù)據(jù)進(jìn)行再加工。如求平均；計(jì)算渦度、散度、垂直積分、計(jì)算差分等。圖形可以按多種格式存儲(chǔ)：ps、png、jpg、tiff、gif、windowsmetafile等。自身以metafile格式存儲(chǔ)。提供多種輔助工具軟件。如看圖、打印、圖形格式轉(zhuǎn)換（gv）等工具。Internet上的GrADS資源2.1GrADS在Internet上的主頁GrADS主頁地址：rg/grads從GrADS主頁上可以找到預(yù)編譯好了的適合于windows環(huán)境下的GrADS軟件包，和適合于各種UNIX環(huán)境下的GrADS軟件包。關(guān)于GrADS在UNIX環(huán)境下的資源請參考附錄。 2.2windows環(huán)境下GrADS資源G——GrADS軟件包安裝執(zhí)行程序。目前在windows環(huán)境下最新的版本為GrADS第1.8版。下載GrADS演示數(shù)據(jù)：從rg/grads/sprite/tutorial下載：model.le.dat和文件；或者GRIB碼格式的數(shù)據(jù)：model.grb、model.gmp和model.ctlGrADS繪圖軟件的安裝（windows環(huán)境）在windows下運(yùn)行GrADS的條件：安裝GrADS軟件包安裝XSERVER軟件包在windows環(huán)境下安裝GrADS軟件包運(yùn)行上述可執(zhí)行文件(GrADS1.8sl8.win32.exe)進(jìn)入第一個(gè)畫面：用鼠標(biāo)選擇“Install”鍵，進(jìn)入第二步，顯示GrADS的一般信息，選擇“OK”進(jìn)入第三步，顯示license信息如下：選擇“IAcceptTems”鍵：進(jìn)入第四步，要確定最終安裝目錄。系統(tǒng)顯示缺省安裝在:“c:\ProgramFile\pcGrADS”目錄下。但為了在Dos環(huán)境下運(yùn)行方便，最好不要安裝在目錄名中有空隔的目錄下—如“ProgramFile”。因此這里可改為“c:\pcGrADS”或其它。修改完后用鼠標(biāo)點(diǎn)擊“OK”鍵安裝程序開始安裝。安裝程序?qū)⑿薷哪愕腶utoexec.bat文件中的setpath=（…原path加上);c:\pcgrads\win32(若安裝在ProgramFile目錄下，c:\progra~1\pcgrads\win32)。安裝完后可從“開始”看到如下畫面：表示GrADS軟件包已安裝成功，現(xiàn)在可以開始使用GrADS軟件了。Getting_Started和tutorial是入門和初學(xué)指南；users是在線幫助手冊。 在開始使用前，先來考察一下安裝程序在GrADS的目標(biāo)目錄（如上例中的“c:\pcGrADS”）下安裝的內(nèi)容：win32下安裝有GrADS的所有可執(zhí)行程序；dat目錄下有運(yùn)行GrADS所需的字庫和地圖數(shù)據(jù)文件；doc目錄下存放幫助和說明文件；lib目錄下存放了一些有趣的模板。GrADS常用命令列表Grads經(jīng)典GrADS，支持Grib,NetCDF,IEEE格式的數(shù)據(jù)。GradsC壓縮模式的GrADS，占內(nèi)存最小。GradsGUI（假）圖形界面的。GradsHDF與經(jīng)典GrADS相當(dāng)，多增加了HDF-SDS格式。GradsNC與經(jīng)典GrADS相同。gv32Windows下看圖（metafile格式）及格式轉(zhuǎn)換工具。rxvt在Windows下打開一個(gè)命名行窗口。 3.2 Xserver的安裝 在windows環(huán)境下運(yùn)行GrADS軟件包，需要有一個(gè)能在PC機(jī)上運(yùn)行的XServer系統(tǒng)支持。這部分軟件并不屬于GrADS支持范圍，因此在GrADS主頁上并不能下載。但I(xiàn)nternet也有許多這方面的免費(fèi)軟件，如：x-win32、omni、Xmanager等軟件。這些軟件安裝并不復(fù)雜，并配有說明文檔。本手冊以Xmanager作為GrADS的XServer支持系統(tǒng)。Xmgr139e.exe——既為Xmanager的安裝執(zhí)行程序，序列號(hào)：018112-136905-999104。安裝后即在桌面上生成Xmanager目錄圖標(biāo)，用鼠標(biāo)雙擊打開該目錄：在開始啟動(dòng)GrADS之前，用鼠標(biāo)雙擊圖標(biāo)，運(yùn)行XServer。也可以從開始程序XXmanager來啟動(dòng)XServer。GrADS繪圖模板初學(xué)者可以通過“初學(xué)指南”既參考GrADS軟件的“Tutorial”聯(lián)機(jī)幫助文件開始學(xué)習(xí)。該文檔可以以鼠標(biāo)點(diǎn)擊“開始”—>“程序”—>“win32GrADS”—>“tutorial”打開Tutorial”聯(lián)機(jī)幫助文件。運(yùn)行tutorial的條件需要有一組演示數(shù)據(jù)：Model.le.dat(數(shù)據(jù)文件—二進(jìn)制)，Model.le.ctl(描述文件—ASCII碼)或Model.grb，Model.ctl，Model.gmp(GRIB碼數(shù)據(jù))第一組（）為模式輸出的五天的預(yù)報(bào)結(jié)果；第二組（Model.grb）為同一結(jié)果按GRIB格式存儲(chǔ)的數(shù)據(jù)，而*.ctl文件（ASCII碼）是對應(yīng)數(shù)據(jù)文件的描述文件。GrADS并不直接使用“數(shù)據(jù)文件”，而是通過“描述文件”間接使用“數(shù)據(jù)文件”。GrADS中“打開一個(gè)數(shù)據(jù)文件”既是指打開一個(gè)數(shù)據(jù)描述文件。關(guān)于GrADS數(shù)據(jù)格式和用戶如何生成該格式的數(shù)據(jù)，將有專門章節(jié)講述；關(guān)于GRIB、NETCDF等數(shù)據(jù)格式的使用也將有專門論述。以下以第一組（或第二組）數(shù)據(jù)說明GrADS的使用。在使用數(shù)據(jù)前，先將數(shù)據(jù)，和文件存于硬盤上的某一目錄下，如c:\pcgrads\sample。首先來了解一下數(shù)據(jù)的內(nèi)容。是一組模式輸出的全球5天數(shù)值預(yù)報(bào)結(jié)果，包括了多個(gè)要素、多層，按經(jīng)緯度網(wǎng)格存放的數(shù)據(jù)。以下是清單：GrADS示例演示目的：學(xué)習(xí)GrADS的基本使用方法。介紹GrADS的常用命令。1.1啟動(dòng)GrADS啟動(dòng)xmanager打開一個(gè)MS-DOS窗口。在MS-DOS窗口輸入命令：cdc:\pcgrads\sample再輸入命令:grads，顯示如下提示：首先顯示GrADS的一般信息，最底行提示是用“Landscape”模式（11X8.5英寸）或portrait（8.5X11英寸—GrADS中長度和大小單位都用英寸）模式顯示圖形輸出窗口。輸入L回車或回車是“Landscape”模式；輸入P回車是用portrait模式。之后進(jìn)入GrADS的命令交互模式，等待用戶輸入命令。此時(shí)你的桌面應(yīng)如下圖所示的特征。若想啟動(dòng)GrADS直接進(jìn)入命令交互模式，輸入命令：grads–l，直接進(jìn)入“Landscape”模式或grads–p，直接進(jìn)入portrait模式。要獲得GrADS命令行幫助，請運(yùn)行g(shù)radsgrads-b以批處理方式運(yùn)行，不顯示圖形輸出grads-l以“Landscape”模式運(yùn)行。grads-p以“portrait”模式運(yùn)行。grads-c進(jìn)入GrADS后，執(zhí)行批處理命令文件，顯示圖形輸出。grads-cl以上選項(xiàng)可聯(lián)合使用。左上角一個(gè)窗口是圖形顯示區(qū)，GrADS的所有圖形輸出結(jié)果在此窗口下顯示；下部為原打開的MS-DOS窗口，此時(shí)正在運(yùn)行GrADS，提示顯示“ga->”，說明正處于GrADS命令等待狀態(tài)。GrADS的所有命令都只能通過該窗口輸出，所有文字信息也都由此窗口輸出。GrADS命令采用以下格式：ga->命令<參數(shù)<參數(shù)，，>，，>輸入的參數(shù)不包括“<>”符號(hào)，“<>”內(nèi)的部份可以省略。所為的命令交互模式既是在GrADS命令提示符下，一步步輸入各種GrADS命令產(chǎn)生各種圖形。1.2退出GrADSga->quit1.3示例演示GrADS命令的使用。 例1ga->openmodel.le.ctl進(jìn)入的第一步，打開一個(gè)描述文件（擴(kuò)展名ctl可省略），GrADS并不直接使用“數(shù)據(jù)文件”，而是通過“描述文件”間接使用“數(shù)據(jù)文件”。ga->qfile顯示GrADS打開數(shù)據(jù)的內(nèi)容，顯示結(jié)果應(yīng)與描述文件說明一致。(可以直接運(yùn)行g(shù)a->q顯示q命令所代的參數(shù)及多項(xiàng)功能解釋。)ga->dps顯示地面氣壓（hr）如上圖?！癲”是“顯示”命令，“ps”是文件中定義的變量名稱，代表地面氣壓。例2緊接上例輸入以下命令：ga->c清除圖面。如不清除圖面，GrADS后續(xù)顯示的圖形將與已存在的圖形產(chǎn)生疊加。ga->setlat40ga->setlon-90ga->setlev500ga->sett1GrADS視所有物理量為四維數(shù)據(jù)，以lon|lat|lev|time—稱為“世界坐標(biāo)“或x|y|z|t—稱為“網(wǎng)格坐標(biāo)”。兩種方式表示。缺省采用第一種方式，水平（lon|lat）自動(dòng)認(rèn)為是經(jīng)度/緯度坐標(biāo)。以x|y|z|t方式表示維數(shù)時(shí)，每一維都是一組從1開始的序列數(shù)，對應(yīng)網(wǎng)格的序號(hào)。本例中：x從1到73（代表從00到3600共72個(gè)格點(diǎn)，第73個(gè)格點(diǎn)即第1個(gè)格點(diǎn)）；y從1到46（代表從-900到900共46個(gè)格點(diǎn)）;z從1到7（從1000hPa開始，共7層）;t從1到5（5天）運(yùn)行“openmodel.le.ctl”命令后，lon|lat或x|y是可變的（lon:00~3600;lat:-900~900）。而垂直和時(shí)間維取固定值，lev=1000hPa或z=1;hr或t=1?！皊et維數(shù)參數(shù)數(shù)值1<數(shù)值2>”命令改變當(dāng)前維數(shù)設(shè)定值。當(dāng)取“數(shù)值1數(shù)值2”時(shí)，表示該維是變化的；而只取一個(gè)值時(shí)，表示該維取固定值，對于網(wǎng)格坐標(biāo)可以是分?jǐn)?shù)值，系統(tǒng)會(huì)自動(dòng)作內(nèi)插。維數(shù)參數(shù)一直保持到再次設(shè)定時(shí)都有效ga->qdimsga->dz顯示位于(400N,900W)，500hPa層，1hr一點(diǎn)的位勢高度。例3緊接上例輸入以下命令：ga->c清除圖面。ga->setlon-1800ga->dz顯示沿400N,，1800W至0度，500hPa層，1hr位勢高度剖面圖。例4ga->c清除圖面。ga->setlat090ga->dz顯示西北半球500hPa層，1hr位勢高度。Lon、lev、t設(shè)定值采用前例使用值不變。例5ga->c清除圖面。ga->sett15gaga->dz動(dòng)畫顯示西北半球500hPa層，1hr至1hr位勢高度。例6ga->c清除圖面。ga->setlat-9090ga->setlon-90ga->setlev1000100ga->sett1ga->dt顯示沿900W，1000-100hPa溫度剖面圖。ga->du在上圖基礎(chǔ)上再疊加上東西風(fēng)分量。例7ga->c清除圖面。ga->setlat40ga->setlon-1800ga->setlev500ga->sett15ga->dz顯示500hPa沿400N，高度的時(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)顯示500hPa高度5天平均。ga->dz-ave(z,t=1,t=5)1hr高度與平均的偏差。例10ga->c清除圖面。ga->dz-ave(z,x=1,x=72)高度與緯向平均值的偏差。例11ga->c清除圖面。ga->dz(t=2)-z(t=1)兩個(gè)時(shí)刻高度的差?；驁?zhí)ga->dz(t+1)-z例12ga->c清除圖面。ga->dz(lev=500)-z(lev=700)兩個(gè)高度間的厚度?；驁?zhí)ga->setz3ga->dz(z+1)-z例13ga->openmodel.ctlga->openmodel.le.ctl同時(shí)打開兩個(gè)文件。ga->dz.2(lev=500)-z.1(lev=700)兩個(gè)高度間的厚度。或執(zhí)ga->setz3ga->dz.1(z+1)-z.2結(jié)果顯示如例12所示在GrADS中一般只簡單使用變量名就可以了，但變量名完全定義格式如下：變量名·文件序號(hào)（維數(shù)參數(shù)+/-/=某一數(shù)值，…）“變量名”是指用戶在“·ctl”文件中定義的變量名稱。GrADS可以同時(shí)打開多個(gè)文件，并為每個(gè)打開的文件編一個(gè)序列號(hào)（從1開始）。上例中，序列號(hào)1指model.ctl文件，序列號(hào)2指。缺省情況下序列號(hào)為1，并可省略。維數(shù)參數(shù)指lat/lon/lev/time或x/y/z/t。下畫線部分可省或部分省略。對于有兩個(gè)不同文件中的數(shù)據(jù)作運(yùn)算時(shí)，兩種數(shù)據(jù)網(wǎng)格要一致。例14ga->reinit刪除所有設(shè)置重新回到剛進(jìn)入GrADS狀態(tài)。ga->openmodel.ctlga->setlat090ga->setlon-1800ga->dvint(ps,q,275)作q的垂直質(zhì)量積分，計(jì)算可降水量。vint(ps,q,top)=,ps和top：hPa。例15ga->reset刪除open命令后的所有設(shè)置。ga->setlat090ga->setlon-1800ga->setlev500ga->dhcurl(u,v)由風(fēng)場導(dǎo)出渦度場。例16ga->c清除圖面。ga->setgxoutshaded以分色圖形方式輸出。ga->dhcurl(u,v)由風(fēng)場導(dǎo)出渦度場。ga->’cbarn10’畫圖例。1:相對長短，>1放大；0:水平；1:垂直。ga->setgxoutcontour以等值線方式輸出（缺省方式）。ga->dz疊加500hPa高度場。注意疊加次序。ga->drawtitle500hPaHeightsandVorticity寫圖標(biāo)題。例17ga->c清除圖面。ga->gxoutvector以箭頭方式表示矢量場。ga->du;v;q風(fēng)矢量場。顯示矢量時(shí)，dx分量；y分量<；標(biāo)量>?！?lt;>”內(nèi)的部分只起標(biāo)顏色的作用。例18ga->c清除圖面。ga->setgxoutstream以流線方式表示矢量場。ga->du;v;q風(fēng)流線場。例19ga->c清除圖面。ga->setgxoutbarb以WMO風(fēng)標(biāo)方式表示矢量場。ga->du;v;q風(fēng)標(biāo)。例20ga->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ù)值。例21ga->c清除圖面。ga->setgxoutcontourga->setlat1580ga->setlon-140-40ga->setlev500ga->setmpdrawon如為off，不畫地圖背景（非經(jīng)緯度數(shù)據(jù)需此項(xiàng)）ga->setpolion如為off不畫國界省界等。ga->setmapauto由系統(tǒng)自定地圖背景的顏色、線型和粗細(xì)。而setmapcolor<line_style<thickness>>由用戶自定義。ga->setmpvals-140-401580極射投影經(jīng)緯度范圍。ga->setmprojnps北半球極射投影(其它取值：latlon:等比經(jīng)緯度投影（缺?。?；scaled：經(jīng)緯度投影并沖滿整個(gè)繪圖區(qū)；sps:南半球極射投影；Lambert:lambert方式投影;robinson:robinson式投影；off:關(guān)閉地圖投影，缺省是打開的。)ga->setmpdsetmres取值：lowres/mres/hires分別代表可使用低分辨率（缺省）、中分辨率和高分辨率的地圖數(shù)據(jù)庫。ga->dz例22ga->c清除圖面。ga->gxoutcontourga->setlat-9090ga->setlon-180180ga->setlev500ga->mprojscaledga->mpdsetmresga->dz例23ga->c清除圖面。ga->gxoutcontourga->setlat-9090ga->setlon-180180ga->setlev500ga->mprojrobinsonga->mpdsetmresga->dz例24ga->c清除圖面。ga->setgxoutcontourga->setlat090ga->setlon-180180ga->setlev500ga->mpvals-1801800-90ga->mprojnpsga->mpdsetmresga->dzga->setgxoutvectorga->du;v;q結(jié)果如封面所示。GrADS繪圖模板上述使用方法是GrADS的基本使用方法，其特點(diǎn)是在GrADS系統(tǒng)提示符下，用戶需要輸入一系列的GrADS命令來完成繪畫。但很容易造成輸入錯(cuò)誤，特別是有些命令設(shè)定后，如果不再重新設(shè)置，是永久有效，而有些命令只是一次有效。如果用戶對初次繪畫的效果不滿意，要增加一些命令反復(fù)繪制時(shí)，按GrADS的基本方法使用會(huì)有許多重復(fù)命令，因此效率會(huì)很低。增加效率的一種方法是編制“繪圖模板”，所為“繪圖模板”就是把繪制命令預(yù)先編到一個(gè)文件中，在GrADS下以批處理方式執(zhí)行。這種文件被稱為“模板”或叫“描述語言”文件—“scripts”。對照上例編制一些簡單的模板設(shè)計(jì)。模板1md01.gs文件清單（對應(yīng)“例1”）’openmodel.le.ctl’‘dps’；﹡本例中，維數(shù)參數(shù)采用打開數(shù)據(jù)文件后的缺省值。﹡模板文件中以“﹡”開始的行是注解行。﹡命令必須用單引號(hào)括起來，行尾一定要有回車。﹡作為一種良好習(xí)慣，在結(jié)尾單獨(dú)有一行“；”執(zhí)行模板1：在DOS命令提示符下輸入命令：c:\pcgrads\sample\grads–clsample01.gsc:\pcgrads\sample\grads–cl“sample01參數(shù)1參數(shù)2…”(.gs擴(kuò)張名可以省略)c:\pcgrads\sample\grads–bclsample01.gs或先進(jìn)入GrADS，在GrADS命令提示符下輸入命令ga->runsample01.gs參數(shù)1參數(shù)2…或ga->sample01參數(shù)1參數(shù)2…（在GrADS命令提示符下，可反復(fù)執(zhí)行run命令。但考慮到命令間相互有影響，在GrADS命令提示符下運(yùn)行多個(gè)模板時(shí)，之間可以先運(yùn)行“ga->reinit/reset”命令再運(yùn)行run命令。）最后退出GrADS在退出前，如果滿意，可把圖形存于文件。有多種方法：ga->printimfile.pngwhite存于名為“file.png”的文件(png格式)。printim命令格式：ga->printimfile選項(xiàng)1，選項(xiàng)2…選項(xiàng)：gif—為gif存于“file”文件（缺省為。png格式）white—白色背景。black—黑色背景。xnnnynnn—輸出圖形的水平（xnnn）和垂直（ynnn）大?。c(diǎn)陣）。如：ga->printimfile.gifgifwhitex800y600另一種方式是：ga->enableprintfile-name.gmfga->print執(zhí)行輸出，結(jié)果存于file-name.gmf—用戶指定的文件，metafile格式,最好代gmf后綴，ga->disableprint只有執(zhí)行了disable命令后，print命令的結(jié)果才真正存于文件。metafile文件可用gv32查看或轉(zhuǎn)換成wmf格式，或用gxps/gxeps轉(zhuǎn)成ps/eps格式。metafile和ps/eps格式圖象的印刷效果是最佳的。最后退出GrADSga->quit模板2md02.gs文件清單（對應(yīng)“例3”‘reinit’’openc:/pcgrads/sample/model.le.ctl’*注意這里路徑的使用與DOS不同。也可使用DOS式的路徑。‘setlat40’‘setlon-1800’‘setlev500’*‘sett1’*‘setccolor0’*設(shè)定顏色[注1]。缺省取前景色。*‘setcmark3’*設(shè)定折線圖節(jié)點(diǎn)標(biāo)記[注2]。缺省取2。*‘setcstyle1’*設(shè)定線型[注3]。缺省取1。*‘setcthick1’*設(shè)定線粗細(xì)[注4]。缺省取1。*‘setgridon33’*設(shè)定是否畫網(wǎng)格線[注5]。缺省gridon?！甦z’‘printimsp02.pngwhite’*把圖形存于文件[注1]：顏色取值：0：黑；1：白；2：紅；3：綠；4：藍(lán)；5：；6：；7：黃；8：桔黃；9：紫；10：黃綠；11：中藍(lán)；12：深藍(lán)；13：；14：深紫；15：灰色。[注2]：標(biāo)記取值：0：不作標(biāo)記；1：＋；2：○；3：●；4：□；5：■；6：X；7：

；8：△；9：不作標(biāo)記；10：①。[注3]：線型：1：實(shí)線；2：長斷線；3：短斷線；4：一長一短斷線；5：點(diǎn)線；6：點(diǎn)劃線；7：點(diǎn)點(diǎn)劃線。[注4]：線粗細(xì)取值：1-10（個(gè)點(diǎn)陣）。取值1-5（含）在屏幕上不顯示變化，5以上才顯示。因此，對于Printim輸出，5以上才顯示變化。[注5]：gridon<線型值><顏色值>;畫網(wǎng)格，并可指定線型和顏色（之一或全部或按缺?。ridhorizontal<線型值><顏色值>;只畫水平網(wǎng)格。gridvertical<線型值><顏色值>;只畫垂直網(wǎng)格。gridoff不畫網(wǎng)格。[注6]以“*”開使的命名行表示用戶可選擇是否使用。許多設(shè)置可首先考由GrADS自行決定，如不滿意，再由用戶自定義，以免畫蛇添足。模板3md03.gs文件清單（對應(yīng)“例4”）’openmodel.le.ctl’‘setlat09‘setlon-1800’‘setlev500’*‘setcint8’指定等值線間隔。*‘setcterpon’樣條平滑“ON”或“OFF”。*‘setcsmoothon’*‘setrgb1615622233’用戶自定義顏色[注1]。*‘setrgb171562341*‘setccolor16’*‘setclab%.1f’等值線標(biāo)記方式[注2]。*‘setclskip2每隔一條等值線標(biāo)記數(shù)值。*‘setclopts30.10.2’‘dz/9.8’‘printimmd03.gufgifwhite’；[注1]：‘setrgb用戶自定義顏色序號(hào)（16到99—1到15系統(tǒng)定義）紅綠藍(lán)‘?！凹t綠藍(lán)”取值范圍0-255[注2]：setclabon（等值線標(biāo)數(shù)值）/off（不標(biāo)記）/forded（強(qiáng)制標(biāo)記）/%gK（在標(biāo)記后加上字符K）/%g%%（加上字符%）/%·2f(保留2位小數(shù))/%03·1f（整數(shù)部分保留3為，若不足3位，前部用0補(bǔ)齊）模板4md04.gs文件清單（對應(yīng)“例4”）’openmodel.le.ctl’‘setlat09‘setlon-1800’‘setlev500’*畫指定值的等值線?！畇etclevs495523534556560564572584588’‘setccols123’可以分別指定以上每一條線的顏色。‘dz/9.8’；模板5md05.gs文件清單：’openmodel.le.ctl’‘setlat09‘setlon-1800’‘setlev500’‘setblack-0.10.1’跳過-0.1到0.1（指0線）不畫。*‘setcmin0’不畫低于0（含）以下的等值線。*’setcmax0’不畫超過0（不含）以上的等值線?！甦hdivg(u,v)*1.e5’散度*‘setblackoff’black設(shè)置在運(yùn)行“d”命令后自動(dòng)設(shè)為“off”?！畇etclevs0’特別處理“0”‘setccols1‘setcthick8‘dhdivg(u,v)*1.e5’；模板6md06.gs文件清單（對應(yīng)“例3”’openc:/pcgrads/sample/model.le.ctl’‘setlat40’‘setlon-1800’‘setlev500’*‘sett1’‘setccolor3’*‘setgridon33’‘setgxoutbar’以直方圖方式輸出?！畇etbargap20’直方圖間隔（20%）?！畇etbaroptsfilled’在直方圖中的矩形中填色；“outline”只畫矩形框而不填色?！畇etbarbasebottom’直方圖中的矩形從底部向上畫；“top”從頂部向下畫；“給一數(shù)值”，從這一數(shù)值開始，大于該值的向上畫，小于該值的向下畫。*GrADS坐標(biāo)軸設(shè)置。GrADS一般自動(dòng)設(shè)置這些值，但也可由用戶自定義。*“坐標(biāo)軸”是指輸出圖面上的2D坐標(biāo)，X軸指水平軸，Y軸指垂直軸。*與數(shù)據(jù)中定義的4D坐標(biāo)是兩個(gè)概念?！畇etxaxis-180010’設(shè)置X軸標(biāo)記范圍。[注1]‘setyaxis5150585050’設(shè)置Y軸標(biāo)記范圍。[注1]‘setylint200’設(shè)置Y軸標(biāo)記間隔。[注2]‘setxlint20’設(shè)置X軸標(biāo)記間隔。[注2]‘setxflipon’X軸翻轉(zhuǎn)。[注3]‘setyflipoff’X軸翻轉(zhuǎn)。[注3]‘setxlpos0b’X軸位置。[注4]‘setannot88’坐標(biāo)軸顏色和線粗細(xì)。[注5]‘setxlopts480.2’標(biāo)記特性。[注6]‘dz’；[注1]:setxaxis開始結(jié)束<間隔>/setyaxis開始結(jié)束<間隔>；設(shè)置X軸/Y軸標(biāo)記范圍，小心！??！“開始<結(jié)束”；開始—結(jié)束范圍要與維數(shù)參數(shù)的范圍一制（setlon-1800）。如不一制，圖形數(shù)據(jù)將沒有代表性，此時(shí)你給什么，它標(biāo)什么，而圖形不變。如不能預(yù)先定出范圍，可由先由GrADS自己決定或調(diào)整間隔（xlint/ylint）。使用xaxis/yaxis的情況是當(dāng)你處理非經(jīng)緯度數(shù)據(jù)時(shí)，你不想用經(jīng)緯度來標(biāo)記你的坐標(biāo)，而GrADS總是把水平坐標(biāo)處理成經(jīng)緯度來標(biāo)記。[注2]：setxlint/ylint間隔，如在此設(shè)置間隔，將取代xaxis/yaxis設(shè)置的間隔。[注3]：翻轉(zhuǎn)X/Y軸，同時(shí)圖像也作相應(yīng)的翻轉(zhuǎn)。[注4]：setxlposoffset（+/-英寸）b（底部）/t（頂部）setylposoffset（+/-英寸）l（左）/r（右）[注5]：setannot顏色〈粗細(xì)〉；邊框的顏色和線粗細(xì)。[注6]：setxlopts/ylopts顏色〈粗細(xì)〈大小—英寸〉〉；標(biāo)記數(shù)字的顏色、線粗細(xì)和字大小。模板7md07.gs文件清單：’openc:/pcgrads/sample/model.le.ctl’‘setlat160’‘setlon42‘setlev1000100’‘sett3’*分左右放置兩幅圖‘setvpage05.506’設(shè)置虛頁面大?。ㄓ⒋纾！畇etxyrevon’坐標(biāo)軸對調(diào)（‘d’命令后自動(dòng)失效）?！畇etannot88’坐標(biāo)軸特性將一直保持?！畇etxlopts410.15’‘setylopts1320.2’‘setzlogon’Z坐標(biāo)取對數(shù)坐標(biāo)（將一直保持）?！甦t’‘setvpage5.51106‘setxlint10‘setxlopts1310.15’‘setylopts420.2’‘dt’；[注1]：一個(gè)物理頁是指：水平11英寸X垂直8.5英寸(Landscape)或18.5X11(portrait)，其上可以畫分出多塊區(qū)域，實(shí)現(xiàn)多圖同時(shí)輸出。Setvpage水平開始結(jié)束垂直開始結(jié)束模板8md08.gs文件清單：’openc:/pcgrads/sample/model.le.ctl’‘setlat0160’‘setlon42‘setlev7‘sett1’*分左右放置兩幅圖。左圖中Y軸尺度由GrADS自定?！畇etvpage05.506’*‘setxyrevon’‘setannot88’‘setxlopts410.2‘setylopts120.2’‘dt’‘setylopts3’‘dz/9.8‘setvpage5.51106*‘setxyrevon’‘setmissonnon’連接斷點(diǎn)?！畇etylopts1’‘setvrange240330’‘dt’‘dz/9.8‘cbar_l-x2.5-y1.5–t“temp”“heigh”’標(biāo)圖例。-x2.5-y1.5位置；–t“temp”“heigh”標(biāo)名稱。模板9md09.gs文件清單：‘setvpage3803.4*‘setarrlabon’缺省為on‘setarrowhead-0.2*‘setarrscl0.1‘runsp17’‘setvpage383.46.8‘setstrmden3’調(diào)節(jié)流線密度（1稀—10疏）‘runsp18.gs’例18?！畇etvpage386.811*‘setdigsize0.04’‘runsp19.gs’例19?！畃rintimmd09.pngwhite’sp17.gs文件清單'openc:/pcgrads/sample/model.ctl''setlon-1800''setlat090''setlev500''setgxoutvector''du;v;q'*’dskip(u,2,2);skip(v,2,2);q’skip(變量名，X_skip,Y_skip),通過控制“X_skip,Y_skip”的數(shù)值來決定X和Y方向的取樣密度（取值1可省略不給）。上例表示u,v在X和Y方向每隔一點(diǎn)取值。對于vector或barb等圖形當(dāng)網(wǎng)格太密時(shí)，將無法畫，只有通過skip函數(shù)調(diào)整密度。sp18.gs文件清單'openc:/pcgrads/sample/model.ctl''setlon-1800''setlat090''setlev500''setgxoutstream''du;v;q'sp19.gs文件清單'openc:/pcgrads/sample/model.ctl''setlon-1800''setlat090''setlev500''setgxoutbarb''du;v;q'模板10md10.gs文件清單：‘’‘setlon0180‘setlat30‘setlev500'definezv=ave(z/9.8,lat=20,lat=40)'定義zv代表20~40N緯度帶間的平均位勢高度。'setgxoutlinefill'*setlfcols23'dz/9.8;zv'在z/9.8和zv兩條曲線間填色。當(dāng)z/9.8<zv時(shí)用紅色（2）；當(dāng)z/9.8>zv時(shí)用綠色（3）。模板11md11.gs文件清單：‘open’‘setlon0180‘setlat090‘setlev500‘setrbcols914411513310712826’設(shè)置彩虹色‘dz/9.8*以下給出GrADS一組基本繪圖指令的使用方法?！?第一組：‘setfont2’’setstring1c245’設(shè)置字串的顏色、位置、粗細(xì)、角度‘setstrsiz0.30.6’字水平大小、垂直大小‘drawstring83firtgroup’在（8,3）位置寫字符串”firtgroup”。位置可先運(yùn)行‘qpos’命令決定?！甦rawtitlefirstline\secondline’“\”起分行作用。*第二組：‘setline346’設(shè)置線的‘drawline3457’‘drawrec3456’畫矩形(3,4),(5,6)對角點(diǎn)‘drawrecf2133’‘drawmark10670.5’‘drawpolyf516264’‘drawwxsym2770.5-18’在（7,7）位置畫大小為0.5顏色為-1粗細(xì)為8的氣象符號(hào)‘setrbcols914411513310712826’設(shè)置彩虹色。GrADS中的多彩圖形都是以此彩虹色—rainbow輸出的。但用戶也可通過以上方式自定義。另外與此相關(guān)的命令還有：setccolorrainbow使用彩虹色;setccolorrerain使用反序的彩虹色;setrbrangelowhigh/setrbrangehighlowsetstrsiz字水平大小<垂直大小>setline線的顏色<線型<粗細(xì)>>setstring字串的顏色<位置<粗細(xì)<角度>>>位置取以下值：tl/tc/tr/l/c/r/bl/bc/br其中，c代表以起點(diǎn)位置居中；r從起點(diǎn)向右，等等。角度取正值逆時(shí)針，取負(fù)順時(shí)針轉(zhuǎn)。drawwxsym氣象符（1-41）xy大小<顏色<粗細(xì)>>顏色取負(fù)值時(shí)表示將由系統(tǒng)自定。預(yù)知每個(gè)符號(hào)的圖形，運(yùn)行runwxsym.gs，結(jié)果如下圖：GrADS模板的高級(jí)應(yīng)用模板12md12.gs文件清單—函數(shù)和調(diào)用在同一文件中：‘’由例21-24，設(shè)計(jì)通用地圖投影圖模板。lon1=-140;lon2=-40設(shè)置變量，用“；”號(hào)分隔或每行定義一個(gè)。lat1=15;lat2=80！??！不需用單引號(hào)括起來mp=nps*mp=sps*mp=lambert*mp=robinson‘setlev500‘sett1’‘setlon’lon1‘’lon2只有GrADS命令部分需要用單引號(hào)括起來?！畇etlat’lat1‘’lat2并注意留空格的方式?！甦efinezz=z/9.8’必須在setlon/lat后定義var=zz注意以上兩種“變量定義”方式的不同，不能直接定義var=z/9.8。zz代表了一種新數(shù)據(jù)，而var=z/9.8只是字符相等。*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ù)?！畇etmpdrawon’‘setpolioff’‘setmapauto’*’setmap556’*‘setmpvals’ln1‘’ln2‘’lt1‘’lt2‘setmproj’mpj‘setmpdsetmres’‘display’vars;模板13—函數(shù)和調(diào)用不在同一文件中md13.gs文件清單—調(diào)用部分：‘’由例21-24，設(shè)計(jì)通用地圖投影圖模板。lon1=-140;lon2=-40設(shè)置變量，用“；”號(hào)分隔或每行定義一個(gè)。lat1=15;lat2=80?。?！不需用單引號(hào)括起來mp=nps*mp=lambert‘setlev500‘sett1’‘setlon’lon1‘’lon2只有GrADS命令部分需要用單引號(hào)括起來?！畇etlat’lat1‘’lat2并注意留空格的方式?！甦efinezz=z/9.8’必須在setlon/lat后定義var=zz*‘md132’var‘’mp‘’lon1‘’lon2‘’lat1‘md132’var?。?！調(diào)用函數(shù)不再是以函數(shù)名，而是以函數(shù)所在的文件名。參數(shù)可以從后向前省略;md132.gs文件清單—函數(shù)部分：functionmaps(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)‘setmpdrawon’‘setpolioff’‘setmapauto’*’setmap556’*‘setmpvals’ln1‘’ln2‘’lt1‘’lt2‘setmproj’mpj‘setmpdsetmres’‘display’vardisplay可簡寫為“d”。;模板14—設(shè)計(jì)在一頁紙上輸出多幅圖功能的模板md14.gs文件清單：Functionall_in_one_page(args)*setgradsoff'setimprundefault.gs'在每次”d”命令前運(yùn)行default.gs模板。usage1="putseveraltempletsinonerealpage.arrangthemintwocollons/rows."usage2="md14templet1templet2...<-printim<file.png>>/<-print<file.gmf>>"*geteachsinglemapfromcommandlineargsi=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)breakendif*gettempletnamesmd.i=wrd*saymd.i'i='ii=i+1wrd=subwrd(args,i+1)endwhile*putthelasttempletaszeromd.i=0;num=math_nint(i/2)*2say"numberoftemplets"i"in"num"virtualpagesintwocollons/rows."say"outputoption-pt="pt"out_file="out_file*getpageinformation'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"pagewidth="xsiz"andpagehiegh="ysizif(ysiz<xsiz)page="landscape"elsepage="portrait"endif*setvpageandploteachsinglemapineachvirtualpage.if(page="landscape")*saypageif(num<4)i=0x1=0x2=xsiz/numwhile(i<num)y1=0y2=ysiz'setvpage'x1''x2''y1''y2'run'md.isay'setvpage'x1''x2''y1''y2x1=x2x2=x2+xsiz/numi=i+1endwhileif(0)*drawaboxoneachvirtulpage'setvpageoff'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)'setvpage'x1''x2''y1''y2*say'setvpage'x1''x2''y1''y2'run'md.k'setvpage'x1''x2''0''y1*say'setvpage'x1''x2''0''y1j=k+iif(md.j);'run'md.j;endif*say""x1=x2x2=x2+xsiz/ik=k+1endwhileif(0)*drawaboxoneachvirtulpage'setvpageoff''drawline0'ysiz/2''xsiz''ysiz/2k=1x1=xsiz/iwhile(k<i)'drawline'x1'0'x1''ysizk=k+1x1=x1+xsiz/iendwhileendifendifendifif(page="portrait")*saypageif(num<4)i=0y1=0y2=ysiz/numwhile(i<num)x1=0x2=xsiz'setvpage'x1''x2''y1''y2*say'setvpage'x1''x2''y1''y2j=num-1-i'run'md.jy1=y2y2=y2+ysiz/numi=i+1endwhileif(0)*drawaboxoneachvirtulpage'setvpageoff'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)'setvpage'x1''x2''y1''y2*say'setvpage'x1''x2''y1''y2l=j-1'run'md.l'setvpage'x2''xsiz''y1''y2*say'setvpage'x2''xsiz''y1''y2if(md.j);'run'md.j;endif*say""y1=y2y2=y2+ysiz/ij=j-2k=k+1endwhileif(0)*drawaboxoneachvirtulpage'setvpageoff''drawline'xsiz/2'0'xsiz/2''ysizk=1y1=ysiz/iwhile(k<i)'drawline0'y1''xsiz''y1y1=y1+ysiz/ik=k+1endwhileendifendifendif'setvpageoff''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;default.gs文件清單：‘setgradsoff’;執(zhí)行模板：dos>grads–cl“md14sp21sp22sp23sp24sp16”或在進(jìn)入grads后執(zhí)行：ga->md14sp12sp22sp23sp24sp16sp21sp22sp23sp24sp16即根據(jù)例21-24，16編寫的模板。例中有一處顯示出由于多個(gè)模板同時(shí)運(yùn)行，某些設(shè)置相互干擾（單獨(dú)運(yùn)行都是正常的）。改進(jìn)的方法是在每個(gè)模板中對相應(yīng)項(xiàng)顯式設(shè)置其參數(shù)，而不是用其缺省值。上圖中圖片大大小小，規(guī)格不能統(tǒng)一。如何精確控制圖片大小是個(gè)難題。setimprundefault.gs的作用是自動(dòng)在每次display命令前執(zhí)行一個(gè)模板（default.gs或其它）；因?yàn)樵诿看蝑isplay命令后，setgradson自動(dòng)打開。GrADS描述語言GrADS描述語言不用括起來。只有GrADS的命令要括起來。一行寫一條語句。多條語句寫在一行時(shí)，中間要用“；”分開。變量間需要用空格分開時(shí)（如GrADS命令的參數(shù)間）文本中輸入的空格是無效的，必需用‘‘—單/雙引號(hào)1個(gè)或多個(gè)空格單/雙引號(hào)的方式來強(qiáng)制加入空格。定義變量變量名（1~8字母數(shù)字）=數(shù)值/“字串”/‘字串’變量名2=變量名1使用變量。while(表達(dá)式)變量名·i=值定義數(shù)組，“i”取一數(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ù)過程體結(jié)束標(biāo)志）函數(shù)調(diào)用過程語句… name(variable1,variable2…)name()若沒有定義變量語句…目前GrADS要求函數(shù)過程的調(diào)用與函數(shù)過程體在同一文件中。方式2——函數(shù)過程體與函數(shù)過程的調(diào)用不在同一文件中定義函數(shù)過程函數(shù)過程文件清單functionname(variable) 緊接函數(shù)過程體 ；（函數(shù)過程體結(jié)束標(biāo)志）函數(shù)調(diào)用過程函數(shù)調(diào)用過程文件清單：語句… ‘函數(shù)過程文件名’參數(shù)1‘’參數(shù)2‘’…語句注意單引號(hào)和強(qiáng)制加入空格運(yùn)算符和表達(dá)式數(shù)學(xué)運(yùn)算：+，-，*，/邏輯運(yùn)算：=,!=,>,>=,<,<=,|,&,!,%(連接符)表達(dá)式即由運(yùn)算符連接起來的變量，其結(jié)果為0—假或非0—即真。控制語句if(express)語句elseendifif(express);語句;語句;endifwhile(express)語句endwhile內(nèi)部函數(shù)subwrd(string,nwd)從string字串中取第nwd（數(shù)字）個(gè)字。subline(string,nl)從string字串中取第nl（數(shù)字）行一整行。Substr(string,start,nl)從字串中取從start字開始nl長的字串。雜項(xiàng)輸出line=”PeterPan,theflyingone”;say“Shesaiditis“l(fā)ineprompt“Shesaiditis“l(fā)ine輸入pullvariable遇到pull語句后，程序等待用戶輸入，variable將代表用戶的輸入的值。query/q<參數(shù)>命令是一個(gè)非常重要的命令，用它可以獲得許多圖原信息。以下片段說明其一種用法及與描述語言其它部分信息交換的確方式。'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)（xsiz,ysiz）是圖形左上角的坐標(biāo)（英寸）。GrADS高級(jí)模板的應(yīng)用在pcgrads\lib目錄下存放有系統(tǒng)及愛好者開發(fā)的一些非常實(shí)用的模板，可以作為有用的工具和學(xué)習(xí)參考。由于該目錄是一個(gè)特殊目錄，因此這些工具可以在任何地方使用到。以下提到的各種模板，都為ascii碼文件，詳細(xì)使用方法可直接參考文件說明。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.asc此模板需要修改basemap.gs文件，在數(shù)據(jù)文件前加上路徑。 md15.gs'openc:/pcgrads/sample/model.ctl''setlon-180180''setlat-9090''setlev500''sett1''setmpdrawoff'不畫海陸線'setcmin9999'用basemap前要求畫任意圖，'dz/9.8'setcmin9999使等值線并不顯示出來。D命令后自動(dòng)矢效。'basemapl3'給陸地填綠色。'basemapo4'海洋填藍(lán)色。'setccolorrainbow'設(shè)為彩虹色，d命令后彩虹色自動(dòng)矢效。'dz/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>>>見例16，sf:相對長短，>1放大；vert=0:水平/1:垂直；xmid，ymid位置。cmap.gs—調(diào)整顏色gs->cmapnumberconnect_the_dots.gs—通過點(diǎn)擊鼠標(biāo)畫折線。define_colors.gs—定義新顏色值：21~79defval_demo.gs—演示“qdefval”和“setdefval…”的使用?？蓪狞c(diǎn)數(shù)據(jù)進(jìn)行修改。font.gs—演示字庫樣本。gs->font<0/1/2/3/4/5>isen.gs—等熵面圖。md16.gs'openc:/pcgrads/sample/model.ctl''setlon50140''setlat090''sett1''setlev1000100''defineue='isen(u,t,lev,320)320K等熵面層風(fēng)場ue,ve;u,v,t三維要素場.'defineve='isen(v,t,lev,320)lev氣壓層(hPa)'setz1''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說明了一個(gè)全局變量使用的例子。pinterp.gs—通過插值產(chǎn)生任意等壓面上的值。md17.gs'openc:/pcgrads/sample/model.ctl''setlon50140''setlat090''setlev1000100'*'definet225='pinterp(t,lev,225)'sett1''setz1'與上例相同,這里總有這樣一奇怪的維數(shù)設(shè)定。'd'pinterp(z,lev,225)225hPa高度。'definet225='pinterp(t,lev,225)'dt225'225hPa溫度以下要把pinterp.gs文件的內(nèi)容包括內(nèi)。functionpinterp(field,pgrid,plev)．．．plotskew.gs—T-logP圖md18.gs'openc:/pcgrads/sample/model.ctl''setlon140''setlat30'*definethevari