SMOKE TEMPORAL ERROR: Missing or negative emission value(s) in inventory for "VAR"

SMOKE Emissions Inventory
1

Hi everyone. I was trying to process monthly MEIC inventory dataset (200*320 pixels, 0.25 by 0.25 deg, 2017) with SMOKE. Everything is okay, except for following error when running the subroutine TEMPORAL:

     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/CO.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 61302 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "CO".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/NH3.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 48663 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "NH3".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/NOX.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 61302 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "NOX".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/PM10.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 59544 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "PM10".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/PM2_5.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 61147 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "PM2_5".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/SO2.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 59334 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "SO2".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/VOC.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 61342 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT
     ERROR: Missing or negative emission value(s) in inventory for "VOC".
      
     "IOAPI_DAT" opened as OLD:READ-ONLY   
     File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/PMC.ncf"
     File type GRDDED3 
     Execution ID "????????????????"
     Grid name ""
     Dimensions: 237 rows, 1 cols, 1 lays, 7 vbles
     NetCDF ID:        11  opened as READONLY            
     Time-independent data.
     Closing file IOAPI_DAT

     *** **ERROR** ABORT in subroutine TEMPORAL
     Problem with input data.

What I found interesting is that when I changed the “Month” field in ARINV file to 0 (which means using annual inventory), no error occurred. So can anyone help me out of this issue if i want to use monthly inventory ? Is this related to incorrect temporal allocation file?

All my log files for this run are listed as follows:
grdmat.meic.CN_36km.CH36km.txt (13.6 KB)
smkinven.meic.CN_36km.txt (29.1 KB)
smkmerge.meic.CN_36km.20190124.CH36km.txt (16.1 KB)
spcmat.meic.CN_36km.cmaq_cb05_soa.txt (19.3 KB)
temporal.meic.CN_36km.20190124.txt (19.3 KB)

And my ARINV file:
arinv.meic.lst.txt (3.8 KB)

Run Script file:
meic_smk_HEMI108k.csh (6.9 KB)

Assign file:
ASSIGNS.MEIC.cmaq.cb05_soa.China_108k.txt (18.8 KB)

temporal allocation file:
atpro_edgar_monthly.txt (1.0 KB)
atpro_edgar_weekly.txt (518 Bytes)
atref_edgar.txt (628 Bytes)
atpro_edgar_hourly.txt (1.4 KB)

2

Hi, Are you still having this issue?

3

Hi, thank you for your reply. Yes, I didn’t solve this issue yet. I had to use the annual inventory instead, which could be processed correctly.

Yi Su

4

The posted hemispheric scripts are designed to work with monthly datasets such as HTAP, whereas the dataset they are using is most likely an annual dataset.

The solution you came up with for their dataset is the correct solution.

5

Thank you for your reply! This helps me clarify my concern.

Yi Su