ncdump應用範例
背景
netCDF是地球科學領域重要的檔案格式。然而因其二進位檔案格式的限制,不方便以文字檔案方式檢閱,官網提供了ncdump讓使用者可以輕鬆讀取其內容。
由於netCDF檔案使用了階層的資訊技術,除了ncdump,使用python或其他高階的軟體(或程式庫),可以直接讀取其內容,不必再以文字檔模式循序讀取。
以下以ioapi協定之內容說明CDL(Common Data format Language)的內容順序如下表所示。
CDL的內容順序
項次
名稱
內容
說明
1
dimensions
所有變數應用到的維度名稱、長度、是否受限
...
>
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)
...
>
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...
>
EMISSCTRL_NML之設定
背景
$EMISSCTRL_NML是CMAQ主程式執行時的一個環境變數,指向排放源檔案的物種對照、與其加權乘數。
可以避免反應機制之標準污染物名稱與排放檔案有出入時,仍然可以順利進行模擬。
如果只有單純的增減量方案或測試,可以藉由乘數的加權來達成,不必另外再進行繁瑣的排放檔前處理過程。
基本上,$EMISSCTRL_NML是個對照的文字檔。
與模式選用的反應機制密切相關
雖然對照(dict)不是Fortran程式的強項,CMAQ還是提供了類似功能的作法。
因與排放修正的增減量有關,在執行時需注意妥善管理檔案及版本。
以下以EmissCtrl_cb6r3_ae7_...
>
pt_timvar程式說明
說明
有關pt_timvar的程式背景、策略構思、後處理及運作流程等等,詳見CMAQ點源變數檔案之準備。
版本說明
版本:
TWN_3X3範圍、資料來源為TEDS資料庫,點源座標系統為LL
EAsia_81K範圍,除了TEDS之外,還讀取REAS處理結果
粗細網格pt_timvar程式版本之比較
項目
EAsia_81K
TWN_3X3
說明
資料庫
REAS+TEDS
TEDS
前者範圍較大需要REAS資料庫
引數及主要檔
逐月CAMx...
>
CMAQ點源變數檔案之準備
前言
CMAQ所需的點源排放量檔案有2個,一者為煙囪參數、座標、指標等常數項,一者為排放量會隨時間改變。此處將討論後者之準備。
由於CMAQ執行與檔案系統都是以日期為主體,因此點源隨時間變化部分最後也要整理切割成以日期為單位的檔案。(<=v5.3)
準備CMAQ點源之前,需先進行REAS電廠以及TEDS點源之處理,整理成CAMx點源檔案。
CMAQ點源相關緣起、與CAMx點源檔案之比較、策略面之考量等等,參見CMAQ點源常數檔案之準備之討論。
就時間項而言,CAMx保持較大的時間範圍,雖然檔案較大,處理起來有些不方便,但檔案管理輕鬆很多。適合系統性的改變。CMAQ逐日檔案如果要進行系統性的修改,還是必須從源頭、系統性來改,在逐日檔案內一一修改將事倍功半。
CM...
>
CMAQ點源常數檔案之準備
前言
顧名思義,這支程式就是為了產生CMAQ所需的點源排放量中常數、不隨時間變化的部分。
因為是常數,沒有必要切割成(複製)成逐日檔。但至少在執行批次範圍內所有日期的變數檔可以共用此常數檔。簡言之,就是執行批次範圍內的煙囪個數、順序、必須保持一致,這使得很多未操作的污染源及煙道,必須仍然保留在排放檔中。
執行變數檔案轉換之前,必須先執行常數檔案之轉換,主要是REAS點源排放量的處理是在這一階段執行。
緣起
在CAMx(與傳統的UAM)中,點源檔案與其他檔案最大的不同就是有常數部分的檔頭,其他檔案可能只有4筆header,點源因有煙囪位置、排放條件等,這些不會隨著時間有所改變,因此放在檔頭。
即使CAMx 6/7等最新版nc格式,仍然採取此一策略。煙...
>
REASv3.1地面排放檔案之處理
背景
目標
全球或東亞地區之空氣污染排放量來源有很多,此處以日本國立環境研究所2019/12/4公開的REAS3.1 (Regional Emission inventory in ASia)為對象,由於早先此一數據庫乃提供區域酸雨研究所用,近年來雖然酸雨議題不再那麼嚴重,該研究所仍然持續發展更新,以應用在地區空氣污染的研究。
本項作業的目標就是將其0.25度解析度之數據庫,轉成網格模式(CAMx、CMAQ等)各層範圍的排放源,重要轉換包括:
經緯度系統之網格 → 正交、直角座標系統
網格範圍解析度 → d1、d2、等指定範圍解析度
污染物質名稱(NMV、PM) → 按反應機制之物質名稱
架構
REAS資料的架構是以污染物_年代來區分壓縮檔,共有12種污染項...
>
156 post articles, 20 pages.