氣象繪圖DIY

今天示範用 GRADS 繪圖軟件畫圖. 我敢告訴大家, 這個步驟是大多數做氣象研究的學生開工後第一個星期內必學的軟件.

先到 GRADS 的網站下載適用於 windows 環境下的版本 ftp://grads.iges.org/grads/1.8/grads-1.8sl11-win32e.exe 並執行安裝.

到美方 CDC 網站 http://www.cdc.noaa.gov/data/gridded/, 可以找到很多格點資料. 請點選 GPCP V2 PRECIPITATION. 進入下一頁後, 會發現有四個檔案可供下載. 請下載第一個 precip.mon.mean.nc, 並放到 C:\ 目錄下.

操作 GRADS 主要有兩個方式, 一個是 interactive, 另一個是行 script. 先從程式集中選取 Grads, 順利打開軟件並選取顯示模式後總共會有兩個視窗, 一個是輸入名令, 另一個顯示圖.

對於 *.nc 類的文件, 開檔的名令是: sdfopen c:\precip.mon.mean.nc

如能開檔, 己成功了一大半. 開檔後常做的是看看裏面有那些變數: q file 1 (q 是 query 的意思)
剛才下載的檔案只有一個變數 precip

輸入: d precip (d 是 display/顯示的意思)
可以顯示雨量, 由於並沒有其他設定值, 顯示的是 1979 年 1 月的全球雨量.

要選擇其他時間和繪畫範圍, 可輸入以下命令:
1) c (clear 的意思, 將剛才的畫面清除)
2) set lat 0 40 (將緯度設為 0 - 40)
3) set lon 100 140 (將經度設為 100 - 140)
4) set t 8 (由於是 79 年 1 月開始的月資料, t=8 是 1979 年 8 月)
5) d precip
6) quit (離開 grads)


要用 script 來執行繪圖(例如執行一些重覆或自動的工作), 先將以下數行連所有標點抄到新檔 c:\test.gs
'sdfopen c:\precip.mon.mean.nc'
'set lat 0 40'
'set lon 100 140'
'set t 8'
'set grads off'
'set gxout shaded'
'd precip'
'run cbarn'
'printim c:\test.png white'
'quit'


再次進入 grads 後, 只需輸入: run c:\test 便可繪圖並將結果輸出至 c:\test.png

在 unix/linux 環境下只需在 shell 執行(windows 好像不行): grads -bc "run c:\test.gs" 便可.