列印m3.nc的時間標籤

Table of contents

• 背景
• 程式使用
• pr_tflag.py 內容

---

背景

• nc檔案的時間標籤是程式能否執行的關鍵。如果時間不正確，CMAQ系列程式很快就會跳出錯誤訊息。
• m3.nc檔案的時間標籤變數為TFLAG，其維度為：TSTEP, VAR, DATE-TIME
• 雖然可以使用ncdump印出檔案內容逐一確認，但是每一個濃度項目都有TFLAG，ncdump似乎就不好用了。
• 此處使用簡單的python程式，將TFLAG打印至螢幕以供檢查

程式使用

• python pr_tflag.py m3.nc_file_name

pr_tflag.py 內容

#!/opt/anaconda3/envs/py37/bin/python
import numpy as np
import netCDF4
import os,sys
fname=sys.argv[1]
nc = netCDF4.Dataset(fname,'r')
v='TFLAG'
print(np.array(nc.variables[v][:,0,:]))