列印m3.nc的時間標籤

Table of contents

背景

  • nc檔案的時間標籤是程式能否執行的關鍵。如果時間不正確,CMAQ系列程式很快就會跳出錯誤訊息。
  • m3.nc檔案的時間標籤變數為TFLAG,其維度為:TSTEP, VAR, DATE-TIME
  • 雖然可以使用ncdump印出檔案內容逐一確認,但是每一個濃度項目都有TFLAGncdump似乎就不好用了。
  • 此處使用簡單的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,:]))