grads軟件系列學(xué)習(xí)第3章基本操作命令

第3章基本操作命令

GrADS系統(tǒng)的運(yùn)行方式，是采用命令行輸入的方式交互式地顯示圖形。3.1基本命令

啟動(dòng)GrADS系統(tǒng)后，在ga->提示符后輸入命令。在交互式環(huán)境內(nèi)常用基本命令如下：

open、

set、

display、

clear……

下面逐一介紹。

open<路徑>數(shù)據(jù)描述文件名

open命令用于打開(kāi)GrADS的數(shù)據(jù)文件，啟動(dòng)GrADS后首先需要打開(kāi)至少一個(gè)數(shù)據(jù)描述文件，命令如下：

openfilename主要用于打開(kāi)數(shù)據(jù)描述文件或者數(shù)據(jù)文件,命令形式為ga->opendata.ctlga->sdfopenu.ncga->xdfopenwind.ddf關(guān)閉數(shù)據(jù)文件(close)ga->close2ga->close1注意要先關(guān)閉后打開(kāi)的文件。set各類(lèi)選項(xiàng)

set命令用于設(shè)置各種環(huán)境參數(shù)，包括維數(shù)環(huán)境、圖形類(lèi)型、圖形要素、屏幕顯示等。如設(shè)置維數(shù)環(huán)境(set)ga->setlon60160ga->setlat2050ga->settimeJul1983ga->setlev200lon、lat、time、lev為地圖坐標(biāo),也可以換為格點(diǎn)坐標(biāo)(x\y\z\t)表示。設(shè)定的維數(shù)環(huán)境決定了隨后參與繪圖或者計(jì)算的數(shù)據(jù)范圍。

display(或者d)為顯示命令

display命令是對(duì)表達(dá)式處理后進(jìn)行屏幕圖形顯示。最簡(jiǎn)單的表達(dá)式是變量名的縮寫(xiě)。如變量名為v則顯示命令:dv

clear(或c)

clear是清屏命令，清除圖形窗口的內(nèi)容。

quit

quit命令用于退出GrADS系統(tǒng)。

query(或q)選項(xiàng)

query是系統(tǒng)環(huán)境設(shè)置的查詢(xún)命令。主要用于查詢(xún)數(shù)據(jù)資料的時(shí)空結(jié)構(gòu)、變量等信息。ga->qfile2

如querydefine：可知道定義了那些變量

dims：當(dāng)前的維數(shù)環(huán)境

filen：查詢(xún)第n號(hào)描述文件的內(nèi)容

files：打開(kāi)n個(gè)文件的次序

gxinfo：用在d之后，告訴用戶(hù)圖的一些信息

shades：用在d之后，告訴用戶(hù)一些陰影的信息

time：時(shí)間設(shè)置信息定義新變量的命令(define/modify)

define命令用于定義新的變量，所定義的新變量可以用于隨后的表達(dá)式中。ga->defineuave=aave(u,t=1,t=50)ga->modifyuaveseasonal(uave被定義為氣候變量)初始化環(huán)境(Reinit/reset)注意Reinit和reset的區(qū)別,前者是關(guān)閉所有文件,重新初始化,后者則是重置初始環(huán)境,但不關(guān)閉數(shù)據(jù)文件,reinit常作為GrADS命令批處理文件中的第一個(gè)命令?；A(chǔ)繪圖指令(draw)

draw命令均用在d之后,用于對(duì)圖形進(jìn)行標(biāo)注字符串或者做標(biāo)記、連線等run*.gs

run命令用于執(zhí)行文件*.gs中定義的操作。

GrADS中幾個(gè)重要gs文件(1)生成圖形(.gmf)文件,并在(85°E,30°N)處標(biāo)字符G′reinit′′opend:\llp\data.ctl′′enableprintd:\llp\h500.gmf′′dh′′qw2xy8530′x1=subwrd(result,3);y1=sbuwrd(result,6)′drawstring′x1′′y1′G′′print′′disableprint′draw命令都要用在display之后。第二種方法（printim命令）：形式：printim<路徑>filenameoptionprintim命令在GrADSv1.8以上版本有效，可以在批處理文件中使用。例：輸出1000x800的PNG圖象：

printimout.pngx1000y800輸出800x600白色背景的GIF圖象：printimgifimage.outgifx800y600white(2)netCDF數(shù)據(jù)的讀寫(xiě)并生成新的數(shù)據(jù)文件′reinit′′sdfopenh.nc′′opendata.ctl′′setgxoutfwrite′′setfwrited:\llp\h500.grd′*(設(shè)定數(shù)據(jù)時(shí)空范圍)′settimejan1982dec1983′′setlon70100′′setlat2040′′setlev500′′dh′′disablefwrite′此例中,若資料為netCDF格式,就用sdfopen所在命令行;否則用open命令行。(3)繪制陰影圖形,0線加粗′reinit′′opendata.ctl′′setgxoutshaded′′setcsmoothon′′dh′′setgxoutcontour′′dh′′setclevls0′′setcthick8′′dh′

當(dāng)陰影圖形和等值線圖形在一起畫(huà),要先畫(huà)陰影后畫(huà)等直線,這一例實(shí)現(xiàn)了3幅圖的疊加,陰影圖、等值線圖、單獨(dú)加粗的0線圖形。(4)熟悉Define命令和modify命令的使用技巧。特別注意,使用define命令前所設(shè)定的維數(shù)環(huán)境是新變量的維數(shù)環(huán)境,而不是原變量的維數(shù)環(huán)境,但是,新的GrADS版本中,define命令可以省略?！鋙pendata.ctl′*(設(shè)定了氣候變量hcli的維數(shù)環(huán)境)′sett112′′setlev500′*定義新變量hcli′′definehcli=ave(h,t+0,t=48,12)′′modifyhcliseasonal′*(設(shè)定了異常變量hano的維數(shù)環(huán)境)′sett148′′definehano=h-hcli′′settimejul1983′′dhano′(5)熟悉GrADS腳本語(yǔ)言中的循環(huán)語(yǔ)句“While.....endwhile”,下面一段程序可以連續(xù)顯示多個(gè)時(shí)次500hPa高度場(chǎng)圖形?！鋙pendata.ctl′t=1while(t<=48)′c′′sett′t′′′dh′pulldummyt=t+1endwhile1.數(shù)據(jù)說(shuō)明

nc.dat是二進(jìn)制格點(diǎn)數(shù)據(jù)文件，為一組模式輸出的全球5天的數(shù)值預(yù)報(bào)結(jié)果。有1個(gè)層次：1000，多個(gè)要素：h85,t85,u85…，時(shí)間從1948年1月開(kāi)始，間隔1月，已知X方向間隔為144個(gè)單位，Y方向間隔為73個(gè)單位。

關(guān)于數(shù)據(jù)文件nc.dat的具體說(shuō)明，