Link Search Menu Expand Document
Focus on Air Quality

產生D1範圍之CMAQ初始濃度檔案序列

Table of contents

背景

  • 全球模式模擬結果是以逐6小時儲存，不但檔案很大，難以管理，也不能同步處理，耗費處理的時間。
  • 因此需要有一個程式按其時間進行拆解，拆解前需要準備好各結果檔案的模版。此為本程式的目的。
  • 拆解後全球模擬結果之水平的內插、空品項目的對照等程序，便可按個別檔案同時進行，最後再以ncrcat、按批次需要的日期予以整併即可。

程式說明

程式執行

引數

  • 年月(4碼)

輸入檔

  1. 全球模式模擬結果全月檔案(經ncrcat及垂直內插處理)
    • 檔名規則：'moz_41_20'+yrmn+'.nc'
    • 只讀取時間標籤
  2. D1範圍CMAQ濃度檔模版：ICON_tmp.d1
    • 除了可以用做初始檔，此檔亦將作為是邊界條件之數據來源。
    • 規格如下
kuang@node03 /nas1/cmaqruns/2016base/data/bcon
$ nc=ICON_tmp.d1
(py37)
kuang@node03 /nas1/cmaqruns/2016base/data/bcon
$ ncdump -h $nc|H
netcdf ICON_tmp {
dimensions:
        TSTEP = UNLIMITED ; // (1 currently)
        LAY = 40 ;
        ROW = 57 ;
        COL = 57 ;
        VAR = 248 ;
        DATE-TIME = 2 ;
variables:
        float AACD(TSTEP, LAY, ROW, COL) ;

輸出檔

  • 每個時間D1範圍CMAQ之濃度檔
  • 檔案命名規則：ICON_20YYJJJHH.d1，YY：年代、JJJ：Julian day、HH：小時

分段說明

  • 引用模組及讀取引數
import numpy as np
import netCDF4
import os,sys,subprocess

if (len(sys.argv) != 2):
  print ('usage: '+sys.argv[0]+' YYMM(1601)')
yrmn=sys.argv[1]
  • 程式會使用到系統的ncks、ncatted 2支程式
ncks=subprocess.check_output('which ncks',shell=True).decode('utf8').strip('\n')
ncatted=subprocess.check_output('which ncatted',shell=True).decode('utf8').strip('\n')
  • 開啟舊檔，讀取時間標籤tflag
fname='moz_41_20'+yrmn+'.nc'
nc = netCDF4.Dataset(fname,'r')
tflag=nc.variables['TFLAG'][:,0,:]
nt,dt=tflag.shape
  • 將每個時間分別寫成CMAQ初始濃度檔案備用
fname='ICON_tmp.d1'
for j in range(nt):
  fnamej=fname.replace('tmp',str(tflag[j,0])+'{:02d}'.format(int(tflag[j,1]/10000)))
#  if '201600418' not in fnamej:continue
  os.system('cp '+fname+' '+fnamej) 
  nc = netCDF4.Dataset(fnamej,'r+')
  nc.variables['TFLAG'][0,:,0]=[tflag[j,0] for v in range(nc.NVARS)]
  nc.variables['TFLAG'][0,:,1]=[tflag[j,1] for v in range(nc.NVARS)]
  nc.SDATE=tflag[j,0]
  nc.STIME=tflag[j,1]
  nc.close()
  os.system(ncatted+' -O -a TSTEP,global,o,i,60000 '+fnamej)
  a=fnamej+'.tmp'
  os.system(ncks+' -O --mk_rec_dmn TSTEP '+fnamej+' '+a+';mv -f '+a+' '+fnamej)

程式下載

Reference

  • sinotec2, NCKS 在空品模式中的應用, FAQ, Dec 10 2021
  • sinotec2, MOZARD/WACCM模式輸出轉成CMAQ初始條件_垂直對照, FAQ, Dec 15 2021