Dr. Kuang's Utilities

D4範圍地面排放檔案之轉換

 
  • Jul 04, 2022
  • views
>

前言

此程式將適用新版CAMx D4範圍地面排放源的netCDF格式檔案，轉成cmaq可以使用的IOAPI格式。2個模式其他檔案格式的比較與互相轉換，可以詳見 CMAQ/CAMx排放量檔案之轉換的說明。

近年來CAMx模式在IO格式的選擇逐漸向CMAQ系統靠攏，除了傳統的uamiv格式之外，也容許nc格式的IO，在camx7版本所使用的nc檔案已經與CMAQ非常接近了。然而畢竟是2個不同的模式，檔案內容有些微的差異。可能因為差異不大，2個模式的官網似乎也沒有強調，但使用者還是必須注意、處理之。

官網提供的CAMx2IOAPI是轉換uamiv格式，並不能適用新版CAMx的nc格式。

除了D4範圍，同樣程式(不同模版)也可以應用在東亞地區粗網格排放量檔案的轉換，詳見CMAQ模式筆記。

CAMx地面排放源檔案的準備

D4範圍CAMx所需檔案{PTH}包括：

  • 面源(area)、
  • 生物源(biog)、
  • 線源(line)、
  • 地面點源(ptseG)、
  • 地面船舶(ship)等5大類，

處理成逐月檔案、檔名原則為傳統之fortBE.413_teds.{PTH}{MM}.nc，詳見TEDS Emission Manipulation。

主要轉換項目

  • 時間項LST改成UTC：CAMx模擬是用LST，CMAQ是用UTC。
  • 單位轉換：CAMx模擬是用gmole/Hr，CMAQ是用gmole/s。
  • 污染物名稱：假設2個模式選用一樣的反應機制，在此不做轉換。

接續執行

CMAQ需要逐日檔案，因此還需要將每月檔案分成逐日儲存在月份目錄下備用（brk_day2.cs(*), see  CMAQ/CAMx排放量檔案之轉換)。

YY=19
for p in area biog line ptse ship; do for m in {01..12};do brk_day2.cs ${p}_TWN_3X3.${YY}${m}.nc;done;done

camx2cmaq.py執行方式

需先準備

  • CMAQ_d4範圍排放檔案之模版template.nc
  • fortBE檔案位置（home）
  • 同步執行各個分類(pathway)
for i in {0..4};do sub python camx2cmaq.py $i;done

程式碼

  • 時間項LST改成UTC
    • jul2=[dt2jul(i+datetime.timedelta(hours=-8)) for i in sdatetime]
    • 時間日期的轉換詳見dtconvertor@FAQ
  • 單位轉換：nc[v][:]=nc0[v][:]/3600.

參考資料與連結

  • netCDF格式
  • cmaq IOAPI格式
  • Here:D4範圍地面排放檔案之轉換
  • Parent:
    • CMAQ/CAMx排放量檔案之轉換
  • Relatives:
    • TEDS Emission Manipulation
    • Dr. Kuang's Evernotes
Search