Home

earth套件讀取CWB_WRF數據

背景 動機與目的 中央氣象局雖然每天進行WRF數值預報,網路上卻沒有開發者使用來做氣流線等加值應用。雖然過去發展了高解析度軌跡近5日預報(see also [[2022-11-04-daily_traj]])、calpuff每日的預報等等,卻還沒有嘗試以earth.nullschool套件來加以延伸、展示。 除了每日、即時的展示之外,還有儲存、檢視的內部需求。相較windy,earth.nullschool可以按照使用者的需要,在url地址列簡單的指令,即可顯示過去特定時間的氣象與(或)濃度場。這個方案目前似乎是網路存取、檢視龐雜的氣象-空品模擬系統成果的最佳平台。 earth的發展與應用 earth套件是Cameron Beccario(cambecc)早期為東京...

Read more

等濃度圖之繪製

背景 等濃度圖為模式模擬結果展示的一種重要型態,此類圖形展現說明可以按照模擬結果本身與背景的型態分成4大類,如表所示: 以模式本身的表現為主,對於背景、影響之標的特性為輔。此類表現型式會模擬結果會以色塊(shaded plot/raster plot)為主,背景(basemap)則以向量線條、點狀方式表現。如,適用在模擬尺度較大的情況,如NCL、GrADS(安裝及資源詳參[[2022-07-21-grads]]1)、VERDI(raster)、MeteoInfo等。 第2類需要較完整的底圖,如google或者是openstreet系列,臺灣內政部底圖雖然也有提供網路圖磚在GIS上應用,局部版本之圖磚目前尚無經驗。NCL雖然也可以連線下載地形圖如...

Read more

GrADS繪製CAMx等濃度圖

背景 雖然這個作業方式停留在2016年版本,迄今沒有顯著的更新,然而作為單一個pdf圖檔的整合工具,目前其他工具似乎也還沒有可以完全取代的功能。 2016年的github wiki,有較簡明的整體描述。此處將進一步詳細說明。 有關等值線的選擇 GrADS之2維等值圖 包括等值線(Contour)與等值色塊圖(Shaded-plotting), 前者的濃度值標示在等值線上,後者則需要另外標示色標(cbar)。 模式模擬規範並未指定等值線的形式,一般報告會傾向使用後者、具有一目瞭然的效果,學術期刊則以前者較為清楚。 matplotlib為python提供的繪圖模組,其效果會較 GrADS具有更多的...

Read more

濃度時序圖之GrADs腳本spec.gs

背景 此一腳本將在一頁A4版面內繪製8個測站的濃度時序圖,並集合所有頁面在同一pdf檔案內。 各污染項目的製圖腳本是以模版形式儲存。 模版內可更動的變數為第三行的spe:污染物名稱 可能的污染物名稱:包括在epa.ctl檔案中(so2 cmo ozn pmt nox p25 no2 voc) 控制修改模版與執行grads的(上層)程式:sss.cs NMHC只有部分測站有觀測值。作圖時仍然保留沒有測值之測站模擬值。 CMAQ類似的功能詳見[[2022-07-21-wsite]]1 GrADS安裝及資源詳參[[2022-07-21-grads]]2 分段說明 讀取epa.ctl檔案 由所在目錄下的...

Read more

從COMBINE結果中讀取測站位置之濃度值(wsite)

背景 wsite是CMAQ系統中的後處理程式之一,作為測站比對的重要程序之一,整體流程參考[CMAQ後處理程序圖][cmaq_pp]。 類似的程式功能在CAMx模式為aok(average vs observatoin, see[[2022-07-20-aok]]) CMAQ後處理程序圖 graph LR A(CCMS模擬結果 ACONC, APMDIAG) -- combine --> B(COMBINE_ACONC) B(COMBINE_ACONC) -- wsite --> C(MDL.csv) D(TZFILE) -- wsite --> C(MDL.csv) E(SITE_FILE) -- wsite --&...

Read more

測站模擬比較之時序圖

背景 測站模擬比較之時序圖是個簡單的XY-plot,用excel就可以輕鬆產出。何以會變得複雜,主要是因為污染物、空品站、月份等維度所造成。 環保署公版模式提供有python版本的後製工具,此處介紹GrADS版本的作業方式。 公版的後製工具可以產生個別測站項目之圖檔, 然而對撰寫模擬報告的作者而言,產生個別比較檔案還算容易,困難的是如將一個個的圖檔整併在頁面上,排列整齊、充分使用版面、還能保持各圖仍有一定的清晰度。這也是ovm2grads的工作重點。 ovm2grads的整體工作如下圖所示。 CAMx的後處理詳流程圖 GrADS安裝及資源詳參grads(also[[2022-07-21-grads]]) 其他CA...

Read more

GrADS筆記

背景 在python matplotlib還沒有流行以前,繪製報告水準等級的數據圖可以使用GrADS、matlab、GrapherTM等軟體來繪製。 而這三者中自動化、版權、支援等等考量條件下,以GrADS為最合理的方案。 GrADS是喬治梅森大學海洋、土地及大氣研究中心(The Center for Ocean-Land-Atmosphere Studies, COLA)持續發展維護的公開軟體,目前仍然持續進版,且擁有許多愛用者,詳見官網之介紹。 GrADS 為遵循GNU原則的開放與免費軟體,適用在個人電腦、工作站等跨平台作業系統,可以直接讀取netCDF、grib等檔格式,作為顯示與繪圖工具。 GrADS亦可直接讀取MM5, ...

Read more

CAMx模擬結果之比對(aok)

背景 AvrgvsHourlyObs.f(aok)這支程式讀取CAMx模擬結果(壓縮-[[2022-07-19-shk過後之地面濃度)、測站當月或指定時間範圍的測值,進行逐時之比較、同時也進行模式之性能評估abi_camx?.f系列程式,以及後續的製圖作業。 CAMx模式後處理整體流程如下圖所示。 aok之前一步驟為shk 下一步驟為ovm2gr.cs及sss.cs以grads繪製時間序列比較圖。 CAMx模式後處理整體流程 graph TD A(YYMMbase.avrg.grd02) -- shk --> B(YYMMbase.S.grd02) B --> H((aok & abi...

Read more