既有點源檔案

Table of contents

背景
直接選取TEDS11點源資料庫
- add_tzpp.py程式說明
add TZPP on BASE2 and BASE3 background

背景

基於點源在公版模式排放檔案中已被網格均化，除非是大型污染源獨佔網格空間，否則難以由背景其他污染源當中切割出來。
此處範例以中部某電廠機組為例，該廠因排放高度超過第1格範圍，周圍沒有其他污染源，因此辨識較為容易。

直接選取TEDS11點源資料庫

這個方案較前述新增點源增量方案單純一些，程式碼控制在100行之內

add_tzpp.py程式說明

完整程式碼詳add_tzpp.py
TEDS11資料庫(有管制編號CP_NO)：/nas1/cmaqruns/2019base/data/ptse/twn/目錄下之fortBE.413_teds11.ptse01.nc(CAMx nc file)
- 按照管制編號及煙囪高度2個變數來搜尋資料庫
常數與時變量檔案模版：目錄下之teds11.1901.timvar.nc、及teds11.1901.const.nc
變量檔
- 公版的啟動時間(0:ibeg、前月25日至月底)較長，需要填入暫用值。
- 回存時又遇到多維度批次篩選議題

kuang@DEVP /nas2/cmaq2019/download/input/201901/grid03/smoke
$ cat add_tzpp.py
...
var=np.zeros(shape=(nv,nt,nlay,l_tzpp,ncol))
for v in V[3]:
  var4=nc0[v][:]
  iv=V[3].index(v)
  var[iv,:,:,:,:]=var4[:,:,tzpp.index,:]
...
for v in V[3]:
  iv=V[3].index(v)
  nc[v][:ibeg,:,:,:]=var[iv,:ibeg,:,:,:]
  nc[v][ibeg:nt0,:,:,:]=var[iv,:(nt0-ibeg),:,:,:]
...  

常數檔
- 因各項煙囪參數並沒有時間、層數的維度，只有ROW的維度較為單純，此時ncf的indexing與np.array行為一樣，可以直接指定數值
- (YLOCA, XLOCA)、(ROW, COL)需要以新的座標系統校正

...
fname='teds11.1901.const.nc'
os.system(ncks+' -O -d ROW,1,'+str(l_tzpp)+' '+root+fname+' '+fname)
nc0 = netCDF4.Dataset(root+fname,'r')
nc = netCDF4.Dataset(fname,'r+')
V=[list(filter(lambda x:nc[x].ndim==j, [i for i in nc.variables])) for j in [1,2,3,4]]
nt,nlay,nrow,ncol=(nc[V[3][0]].shape[i] for i in range(4))
for v in V[3]:
  nc[v][0,0,:,0]=nc0[v][0,0,tzpp.index,0]
for i in atts:
  if i not in dir(nc00):continue
  exec('nc.'+i+'=nc00.'+i)
nc.NROWS=l_tzpp
pnyc = Proj(proj='lcc', datum='NAD83', lat_1=nc.P_ALP, lat_2=nc.P_BET,lat_0=nc.YCENT, lon_0=nc.XCENT, x_0=0, y_0=0.0)
lat,lon=nc['LATITUDE'][0,0,:,0],nc['LONGITUDE'][0,0,:,0]
x0,y0=pnyc(lon,lat, inverse=False)
nc['XLOCA'][0,0,:,0]=x0[:]
nc['YLOCA'][0,0,:,0]=y0[:]
nc['COL'][0,0,:,0]=np.array((x0[:]-nc.XORIG)/nc.XCELL,dtype=int)
nc['ROW'][0,0,:,0]=np.array((y0[:]-nc.YORIG)/nc.YCELL,dtype=int)
nc['TFLAG'][0,:,0]=nc.SDATE
nc['TFLAG'][0,:,1]=nc.STIME
nc.close()

add TZPP on BASE2 and BASE3 background

將同樣高空點源排放量(add_tzpp.py、加在2組不同背景條件下敏感性的比較

背景條件說明

氣象：2019年1月
BASE2:2組面源，分別是生物源及背景基準排放量(剔除特定高空點源)
- bio3taiwan：${cmaqproject}/smoke/b3gts_l.20181225.38.d4.ea2019_d4.ncf
- basetaiwan：${cmaqproject}/smoke/cmaq_cb06r3_ae7_aq.01-20181225.38.TW3-d4.BaseEms.ncf_0-8NoTZPP
BASE3：除了前述2者，再加上egts第3層網格排放內插之排放量
- bio3taiwan：${cmaqproject}/smoke/b3gts_l.20181225.38.d4.ea2019_d4.ncf
- basetaiwan：${cmaqproject}/smoke/cmaq_cb06r3_ae7_aq.01-20181225.38.TW3-d4.BaseEms.ncf
- d3_to_d4：${cmaqproject}/smoke/egts_l.20181225.38.d4.ea2019_d4.ncf

Background O3 Max Hr Comparisons


BASE2之O₃全月最大小時值	BASE3之O₃全月最大小時值

TZPP PM2.5 Increments

2019年1月最大日均值
BASE2環境背景TZPP之最大增量：3.41 μg/m³
BASE3環境背景TZPP之最大增量：3.59 μg/m³


BASE2之PM_2.5增量	BASE3之PM_2.5增量