GrADS繪圖軟件安裝及入門

GrADS繪圖軟件安裝于入門目 錄第一章 GrADS繪圖軟件概述GrADS繪圖軟件簡介GrADS繪圖軟件的安裝〔windows環(huán)境〕3.1windowsGrADS其次章 GrADS繪圖模板GrADS例如演示GrADSGrADS例如演示GrADS命令的使用其次章 GrADS繪圖軟件概述GrADS繪圖軟件簡介TheGridAnalysisandDisplaySystem〔GrADS〕的科學(xué)數(shù)據(jù)繪圖軟件包。其主要特點(diǎn):GrADSInternet上免費(fèi)獲得.WindowsUnix工作平臺(tái)。層等)和時(shí)間/xyzt 4維。數(shù)據(jù)可以是格點(diǎn)化的數(shù)據(jù)或離散點(diǎn)數(shù)據(jù).GrADS特別適用于氣象類數(shù)據(jù)的分析。但也完全可以用于更廣泛類型的數(shù)據(jù)分析。GrADS有多種顯示方式：等值線、流線、矢量圖、風(fēng)矢量圖、站點(diǎn)填圖、折線圖、直方圖等多種兩維圖形。式和系統(tǒng)自定義的一種二進(jìn)制數(shù)據(jù)格式。承受命令行輸入的方式交互式地顯示圖形.并有多種命令對數(shù)據(jù)進(jìn)展再加工。如求平均；計(jì)算渦度、散度、垂直積分、計(jì)算差分等。7ps、png、jpg、tiff、gif、windowsmetafile等。自身以metafile格式存儲(chǔ).(gv〕等工具。InternetGrADS資源GrADSInternetGrADS主頁地址：：///grads從GrADS主頁上可以找到預(yù)編譯好了的適合于windows環(huán)境下的GrADS軟件包，2。2 環(huán)境下GrADS資源GrADS第2。0版。下載GrADS演示數(shù)據(jù)：從ftp：///grads/sprite/tutorial下載：model。le。dat和model。le.ctl文件;GRIB碼格式的數(shù)據(jù)：model.grb、model.gmpmodel。ctlGrADS繪圖軟件的安裝〔windows環(huán)境〕windowsGrADS的條件：GrADS3.1在windowsGrADS進(jìn)入第一個(gè)畫面：OK”GrADSNext”進(jìn)入第三步，顯示信息如下：8選擇“IAcceptTems”鍵：進(jìn)入第四步，要確定最終安裝名目。系統(tǒng)顯示C:\OpenGrAD名目下。但為了在Dos最好不要安裝在名目名中有空隔的名目下-如“ProgramFile”。因此這里可改為“c：\OpenGrADS”或其它。OK“鍵安裝程序開頭安裝。安裝程序?qū)⑿薷哪愕腜rogramFile名目下，：～安裝完后可從“開頭“看到如下畫面：9表示GrADS軟件包已安裝成功，現(xiàn)在可以開頭使用GrADS軟件了.Getting_Startetutorial是入門和初學(xué)指南；users冊。GrADS的目標(biāo)名目〔如上例中的“c:\OpenGrADS”)下安裝的內(nèi)容：10需的字庫和地圖數(shù)據(jù)文件；doc名目下存放幫助和說明文件；lib名目下存放了一些好玩的模板.GrADS常用命令列表Grads 經(jīng)典GrADS，支持Grib,NetCDF,IEEE格式的數(shù)據(jù)。GradsC 式的GrADS,占內(nèi)存最小。GradsGUI 〔假〕圖形界面的.GradsHDF GradsNC GrADS一樣。gv32 Windows下看圖〔metafile格式〕及格式轉(zhuǎn)換工具。rxvt 在Windows下翻開一個(gè)命名行窗口。第三章GrADS第三章GrADS繪圖模板初學(xué)者可以通過“初學(xué)指南”既參考GrADS軟件的“Tutorial”聯(lián)機(jī)幫助文件開頭學(xué)習(xí)。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ù)格式的使用也將有特地論述.以下以第一組數(shù)據(jù)說明GrADS的使用在使用數(shù)據(jù)前先將數(shù)據(jù)和Model 。dat 文件存于硬盤上的某一目錄下，如C:\OpenGrADS\Contents\Resources\SampleDatasets.首先來了解一下數(shù)內(nèi)容。Model。le。dat是一組模式輸出的全球5天數(shù)值預(yù)報(bào)結(jié)果，包括了要素、多層,按經(jīng)緯度網(wǎng)格存放的數(shù)據(jù)。以下是Model。ctl清單：11ps01Surfacepressure［hPa］ua733Eastwardwind[m/s］va734Northwardwind［m/s］zg77Geopotentialheight[m］ta711AirTemperature［K］hus751Specifichumidity[kg/kg］ts011Surface(2m)airtemperature［K］pr059Totalprecipitationrate［kg/(m^2*s〕]GrADS例如演示目的：GrADS的根本使用方法。GrADS的常用命令。GrADSMS—DOS窗口。在 MS —DOS 窗口輸入命令:cd C ：\OpenGrADS\Contents\Resources\SampleDatasets再輸入命令:grads:12首先顯示GrADS的一般信息，最底行提示是用“Landscape”模式〔11X8.5英寸〕或portrait(8。5X11英寸—GrADS中長度和大小單位都用英寸〕模式顯示圖形輸出窗口。輸入 L回車或回車是模式；輸入P回車是用portrait模式。之后進(jìn)入GrADS的命令交互模式 ，等待用戶輸入命令.此時(shí)你的桌面應(yīng)如以以下圖所示的特征假設(shè)想啟動(dòng)GrADS直接進(jìn)入命令交互模式輸入命令“Landscape“模式或grads–p，直接進(jìn)入portrait模式。要獲命令行幫助，請運(yùn)行g(shù)rads–help。grads -b以批處理方式運(yùn)行，不顯示圖形輸出grads -l 以“Landscape“模式運(yùn)行。grads —p以“portrait”模式運(yùn)行。grads —cGrADS后，執(zhí)行批處理命令文件,顯示圖形輸出。grads -cl以上選項(xiàng)可聯(lián)合使用。13左上角一個(gè)窗口是圖形顯示區(qū)，GrADS的全部圖形輸出結(jié)果在此窗口下顯示;MS—DOS窗口，此時(shí)正在運(yùn)行GrADS，提示顯示“ga-〉”，說明正處于GrADSGrADS的全部命令都只能通過該窗口輸出，全部文字信息也都由此窗口輸出。GrADS命令承受以下格式:ga-〉命令〈參數(shù)〈參數(shù)，,>，〉>>”內(nèi)的部份可以省略。所為的命令交互模式既是在GrADS命令提示符下，一步步輸入各種GrADS命令產(chǎn)生各種圖形。GrADSga—〉quit1411。31ga-〉openmodel.le.ctl進(jìn)入的第一步,翻開一個(gè)描述文件〔擴(kuò)展名ctl可省略，GrADSga—>qfilega—〉dps顯示GrADS翻開數(shù)據(jù)的內(nèi)容，顯示結(jié)果應(yīng)與描q命令所代的參數(shù)及多項(xiàng)功能解釋。)的變量名稱，代表地面氣壓。2緊接上例輸入以下命令:ga—〉c 如不去除圖面，GrADS后續(xù)顯示的圖形將與已存在的圖形產(chǎn)生疊加。ga-〉setlat 40ga—〉setlon —90ga—>setlev 500ga-〉sett1GrADS視全部物理量為四維數(shù)據(jù)，以lon|lat|lev｜time-稱為“x｜y|z|t—稱為“網(wǎng)格坐標(biāo)“。兩種是經(jīng)度/緯度坐標(biāo).1開頭表從00到3600共7273個(gè)格點(diǎn)即第1個(gè)格點(diǎn)90046個(gè)格點(diǎn));z17〔1000hPa7層);t1〔5天命令后,lon|latx|y〔lon:00～3600;lat-90～900lev=1000hPaz=1;time=1987。1.2。0hrt=1。維數(shù)參數(shù)數(shù)值<;而只取一個(gè)值時(shí)，表示該維取固定值，對于可以是分?jǐn)?shù)值，系統(tǒng)會(huì)自動(dòng)作內(nèi)插.維數(shù)參數(shù)始終保持到再次設(shè)定時(shí)都有效.要觀看當(dāng)前維數(shù)狀況,運(yùn)行:ga—>q dimsga-dz 顯示位于(400N,900W),500hPa層,19hr一點(diǎn)的位勢15高度.3緊接上例輸入以下命令：ga->c 去除圖面。ga-〉setlon —1800gadz 顯示沿400，1800W至0度500hPa層1987.1.2。0hr位勢高度剖面圖.4ga->c 去除圖面。ga->setlat 0 90ga->dz 顯示西北半球500hPa層，1987.1.2。0hr位勢高度.Lon、lev、t設(shè)定值承受前例使用值不變。165ga->c 去除圖面。ga-〉sett 1 5gaga—〉dz 動(dòng)畫顯示西北半球500hPa層，1987。1。2。0hr至1987。1.2。0hr位勢高度.6ga->c 去除圖面。ga—〉setlat—9090ga-〉setlon-90ga-〉setlevga->sett11000100ga-〉dt 顯示沿900W,1000-100hPa溫度剖面圖。ga—>du 在上圖根底上再疊加上東西風(fēng)重量。177ga—>c 去除圖面.ga-〉setlat 40ga->setlon -180 0ga->setlev 500ga—〉sett1 5ga-〉d z 188ga->c 去除圖面。ga—〉d sqrt〔u*u+v*v)顯示500hPa全風(fēng)速值?；騡a—>dmag〔u,v〕GrADS內(nèi)部定義了多種函數(shù)。9ga-〉c 去除圖面。ga—>d ave(zt=1,t=5〕500hPa5天平均。ga->dz-ave〔z,t=1，t=5)1987.1。2.0hr高度與平均的偏差。1910ga-〉c 去除圖面.ga-〉dzave〔zx=1，x=72)高度與緯向平均值的偏差。11ga->c 去除圖面。ga->dz(t=2)-z〔t=1〕兩個(gè)時(shí)刻高度的差?；驁?zhí)ga->dz(t+1)—z2012ga->c 去除圖面.ga—〉dz(lev=500)—z(lev=700)兩個(gè)高度間的厚度。或執(zhí)ga—〉set z 3ga-〉dz〔z+1〕-z13ga->open model.ctlga—>open model。le.ctl 同時(shí)翻開兩個(gè)文件.21。2(lev=500〕—z.1〔lev=700〕兩個(gè)高度間的厚度?；驁?zhí)ga—〉setz 3ga—〉dz.1〔z+1〕-z.2 結(jié)果顯示如例12所示下：變量名·文件序號(維數(shù)參數(shù)+/-/=某一數(shù)值，…〕“變量名“是指用戶在“·ctl”文件中定義的變量名稱。GrADS可以同時(shí)翻開多個(gè)文件,并為每個(gè)翻開的文件編一個(gè)序列號〔從1開頭。上例中，序列號1指。1，并可省略。lat/lon/lev/timex/y/z/t。下畫線局部可省或局部省略。對于有兩個(gè)不同文件中的數(shù)據(jù)作運(yùn)算時(shí)，兩種數(shù)據(jù)網(wǎng)格要全都 。14ga-〉reinit 刪除全部設(shè)置重回到剛進(jìn)入GrADS狀態(tài)。ga->open model。ctlga->set lat 0 90ga-〉set lon -180 0ga->dvint〔ps,q,275)q的垂直質(zhì)量積分,計(jì)算可降水量。qdpvint〔ps,q,top)=1topqdpgps

， pstop：hPa。15ga->reset open命令后的全部設(shè)置。ga->set lat 0 90ga-〉set lon -180 0ga—>set lev 500ga-〉dhcurl〔u,v〕由風(fēng)場導(dǎo)出渦度場.2216ga-〉c 去除圖面。ga->setgxoutshaded 以分色圖形方式輸出。ga—〉dhcurl(u,v〕由風(fēng)場導(dǎo)出渦度場。ga->’cbarn 1 0’畫圖例。1:相對長短，>1放大；0:水平；1：垂直。ga—>setgxout contour 以等值線方式輸出(缺省方式。ga->d z疊加500hPa高度場。留意疊加次序 。ga—〉draw title 500hPaHeightsandVorticity 寫圖標(biāo)題。2317ga—>c 去除圖面.ga-〉gxoutvector 以箭頭方式表示矢量場.ga—>du;v;q風(fēng)矢量場.顯示矢量時(shí)，d x重量;y重量〈標(biāo)量.〈”內(nèi)的局部只起標(biāo)顏色的作用。18ga—>c 去除圖面。ga—〉setgxoutstream 以流線方式表示矢量場。ga—〉du;v；q 風(fēng)流線場.2419ga-〉c 去除圖面.ga-〉setgxoutbarb 以WMO風(fēng)標(biāo)方式表示矢量場。ga—>duvq 風(fēng)標(biāo)。20ga->c 去除圖面.ga—>setgxoutgrid 直接輸出網(wǎng)格點(diǎn)數(shù)值。ga—〉setdignum0 設(shè)置保存小數(shù)位數(shù)。ga—>setdigsize0.1 置數(shù)字大小(英寸。ga->du東西風(fēng)網(wǎng)格點(diǎn)數(shù)值。2521ga—〉c 去除圖面。ga->setgxoutcontourga-〉setlat 15 80ga—〉setlon —140 -40ga-〉setlev 500ga->set mpdraw on 〕ga-〉set poli on 如為off不畫國界省界等。ga-〉set map auto 由系統(tǒng)自定地圖背景的顏色、線型和粗細(xì)。而set ma