Home

ncdump應用範例

背景 netCDF是地球科學領域重要的檔案格式。然而因其二進位檔案格式的限制,不方便以文字檔案方式檢閱,官網提供了ncdump讓使用者可以輕鬆讀取其內容。 由於netCDF檔案使用了階層的資訊技術,除了ncdump,使用python或其他高階的軟體(或程式庫),可以直接讀取其內容,不必再以文字檔模式循序讀取。 以下以ioapi協定之內容說明CDL(Common Data format Language)的內容順序如下表所示。 CDL的內容順序 項次 名稱 內容 說明 1 dimensions 所有變數應用到的維度名稱、長度、是否受限 ...

Read more

hsinda3G.py程式說明

背景 這支程式讀取傳統的CAMx點源point_source格式檔案,從中抽取開發計畫的排放數據,填入CAMx6的nc檔案中。 CAMx6點源nc檔案的內容、以及與CAMx7的差異詳見比較表CAMx67。 程式之執行 檔案準備 CAMx6點源檔案模版 檔案位置:/nas1/camxruns/2016_v7/ptse/XindaG3/template.nc 檔案維度內容 kuang@master /nas1/camxruns/2016_v7/ptse/XindaG3 $ ncdump -h $nc|H netcdf template { dimensions: TSTEP = UNLIMITED ; // (1 currently) ...

Read more

pt_const程式說明

說明 整體工作流程詳見CMAQ點源常數檔案之準備([[2022-07-06-pt_const]]) 目的:產生CMAQ所需的點源排放量中常數、不隨時間變化的部分 相關程式見[[2022-07-07-pt_timvar.py]]及[[2022-07-06-pt_timvar]]說明 執行 引數:CAMx分月高空點源,MM=01~12 for MM in {01..12}; do python pt_constLL.py fortBE.413_teds11.ptse$MM.nc done 輸入檔(ptse檔案):ptsEnHRBE.py港區計算結果([[2022-07-15-ptsEnHRBE.py程式說明]]) 結果:teds11.YYMM.c...

Read more

EMISSCTRL_NML之設定

背景 $EMISSCTRL_NML是CMAQ主程式執行時的一個環境變數,指向排放源檔案的物種對照、與其加權乘數。 可以避免反應機制之標準污染物名稱與排放檔案有出入時,仍然可以順利進行模擬。 如果只有單純的增減量方案或測試,可以藉由乘數的加權來達成,不必另外再進行繁瑣的排放檔前處理過程。 基本上,$EMISSCTRL_NML是個對照的文字檔。 與模式選用的反應機制密切相關 雖然對照(dict)不是Fortran程式的強項,CMAQ還是提供了類似功能的作法。 因與排放修正的增減量有關,在執行時需注意妥善管理檔案及版本。 以下以EmissCtrl_cb6r3_ae7_...

Read more

pt_timvar程式說明

說明 有關pt_timvar的程式背景、策略構思、後處理及運作流程等等,詳見CMAQ點源變數檔案之準備。 版本說明 版本: TWN_3X3範圍、資料來源為TEDS資料庫,點源座標系統為LL EAsia_81K範圍,除了TEDS之外,還讀取REAS處理結果 粗細網格pt_timvar程式版本之比較 項目 EAsia_81K TWN_3X3 說明 資料庫 REAS+TEDS TEDS 前者範圍較大需要REAS資料庫 引數及主要檔 逐月CAMx...

Read more

CMAQ點源變數檔案之準備

前言 CMAQ所需的點源排放量檔案有2個,一者為煙囪參數、座標、指標等常數項,一者為排放量會隨時間改變。此處將討論後者之準備。 由於CMAQ執行與檔案系統都是以日期為主體,因此點源隨時間變化部分最後也要整理切割成以日期為單位的檔案。(<=v5.3) 準備CMAQ點源之前,需先進行REAS電廠以及TEDS點源之處理,整理成CAMx點源檔案。 CMAQ點源相關緣起、與CAMx點源檔案之比較、策略面之考量等等,參見CMAQ點源常數檔案之準備之討論。 就時間項而言,CAMx保持較大的時間範圍,雖然檔案較大,處理起來有些不方便,但檔案管理輕鬆很多。適合系統性的改變。CMAQ逐日檔案如果要進行系統性的修改,還是必須從源頭、系統性來改,在逐日檔案內一一修改將事倍功半。 CM...

Read more

CMAQ點源常數檔案之準備

前言 顧名思義,這支程式就是為了產生CMAQ所需的點源排放量中常數、不隨時間變化的部分。 因為是常數,沒有必要切割成(複製)成逐日檔。但至少在執行批次範圍內所有日期的變數檔可以共用此常數檔。簡言之,就是執行批次範圍內的煙囪個數、順序、必須保持一致,這使得很多未操作的污染源及煙道,必須仍然保留在排放檔中。 執行變數檔案轉換之前,必須先執行常數檔案之轉換,主要是REAS點源排放量的處理是在這一階段執行。 緣起 在CAMx(與傳統的UAM)中,點源檔案與其他檔案最大的不同就是有常數部分的檔頭,其他檔案可能只有4筆header,點源因有煙囪位置、排放條件等,這些不會隨著時間有所改變,因此放在檔頭。 即使CAMx 6/7等最新版nc格式,仍然採取此一策略。煙...

Read more

REASv3.1地面排放檔案之處理

背景 目標 全球或東亞地區之空氣污染排放量來源有很多,此處以日本國立環境研究所2019/12/4公開的REAS3.1 (Regional Emission inventory in ASia)為對象,由於早先此一數據庫乃提供區域酸雨研究所用,近年來雖然酸雨議題不再那麼嚴重,該研究所仍然持續發展更新,以應用在地區空氣污染的研究。  本項作業的目標就是將其0.25度解析度之數據庫,轉成網格模式(CAMx、CMAQ等)各層範圍的排放源,重要轉換包括: 經緯度系統之網格 → 正交、直角座標系統 網格範圍解析度 → d1、d2、等指定範圍解析度 污染物質名稱(NMV、PM) → 按反應機制之物質名稱 架構 REAS資料的架構是以污染物_年代來區分壓縮檔,共有12種污染項...

Read more