请教一下SMOKE的输入数据问题

HarryChang

目前在做CMAQ，需要用SMOKE处理排放数据，关于输入数据有一些疑问，想向用过的大神们请教一下。
我现在只考虑AREA部分的数据，暂不处理Mobile和Point，在进行smkinven的时候就出问题了，查看日志文件没看明白是怎么回事。
我的输入数据：
（1）ge_dat目录下自己制作的内容有：
国省县查找表costcy：根据实际情况手动制作
网格描述文件GRIDDESC：从WRF中输出
节假日holidays：根据实际情况手动制作
时间分配权重因子文件amptro_monthly/daily/hourly：根据历史数据制作
查找表amptref_tref：根据时间分配权重文件制作
（2）inventory目录下的部分：
area排放清单ARINV：根据MEIC排放数据制作
invtable：直接cp的测试数据nctox数据里面的invtable

我现在想先做一下smkinven和temporal，就暂时没有制作空间分配的静态文件，没想到清单读入阶段就报了错。想请教一下是哪一方面的问题。是我的输入数据不够还是什么问题呢？
下面是日志里面给的报错信息，可以看出来是某个中间文件没有生成。

HarryChang

文件读入，找到相应的路劲了，应该没有问题
     File "ARINV" opened for input on unit:  98
     /home/harry.chang/SMOKE/realdata/inventory/meic2016/area/ARINV_2016_zhangjiakou_1degree_orl.txt

     Value for USE_EXP_GEO_CODES not defined;returning default:   FALSE

     File "COSTCY" opened for input on unit:  97
     /home/harry.chang/SMOKE/realdata/ge_dat/costcy.txt


     File "INVTABLE" opened for input on unit:  96
     /home/harry.chang/SMOKE/realdata/inventory/meic2016/other/invtable_hapcap_cb05soa.txt


然后报了几个警告信息，不知道什么意思，是我的污染物关键字NOX没有被认出来的意思？
     WARNING: No kept pollutants found at line      12 . The source will be dropped.
     WARNING: No kept pollutants found at line      13 . The source will be dropped.
     WARNING: No kept pollutants found at line      14 . The source will be dropped.
     WARNING: No kept pollutants found at line      15 . The source will be dropped.
     WARNING: No kept pollutants found at line      16 . The source will be dropped.
     WARNING: No kept pollutants found at line      17 . The source will be dropped.
     WARNING: No kept pollutants found at line      18 . The source will be dropped.
     WARNING: No kept pollutants found at line      19 . The source will be dropped.
     WARNING: No kept pollutants found at line      20 . The source will be dropped.
          First 10 lines of current inventory:
          "03713","2199001000","02","NOX","9.316457"
          "03711","2199001000","02","NOX","94.05709"
          "03711","2199001000","02","NOX","49.197"
          "03711","2199001000","02","NOX","46.37614"
          "03712","2199001000","02","NOX","35.23962"
          "03713","2199001000","02","NOX","15.67034"
          "03714","2199001000","02","NOX","97.5354"
          "03714","2199001000","02","NOX","41.51996"
          "03711","2199001000","02","NOX","38.32697"
          "03708","2199001000","02","NOX","128.362"
（后面还有一串）
     WARNING: No kept pollutants found at line      21 . The source will be dropped.
     WARNING: No kept pollutants found at line      22 . The source will be dropped.
     WARNING: No kept pollutants found at line      23 . The source will be dropped.

然后又报了一堆告警，说CAS不在清单里面？
我看nctox数据里面并没有用POLID而是用的POLCODE，我在清单里面就没有弄CAS码，直接敲的污染物名

然后下面就是报错的部分了
File "ASCC" opened for output on unit:  94
     /home/harry.chang/SMOKE/realdata/inventory/zjk/ASCC.area.zjk.txt

     Writing out ACTUAL SCC file...
     Generating unique lists from inventory data...
     Value for SMKINVEN_FORMULA defined but empty; returning default:  ' '
     Value for SMKINVEN_FORMULA defined but empty; returning default:  ' '
     Value for SMK_BASEYR_OVERRIDE:  0
     Value for INVNAME1:  'area.zjk'
     Value for INVNAME1:  'area.zjk'
     Value for INVNAME2:  'areasrc.zjk'
     Value for INVNAME2:  'areasrc.zjk'

     File "AREA" opened for output on unit:  93
     /home/harry.chang/SMOKE/realdata/inventory/zjk/area.map.area.zjk.txt

     Value for AREA:  '/home/harry.chang/SMOKE/realdata/inventory/zjk/area.map.area.zjk.txt'
     Value for AREA:  '/home/harry.chang/SMOKE/realdata/inventory/zjk/area.map.area.zjk.txt'
     NOTE: The number of inventory years was      0  and the inventory base
          year was set to the most frequent year, ******
     Number of variables per file array is not allocated for file set IOAPI_INV;
     using default of  2048 variables per file
     
     Value for IOAPI_CHECK_HEADERS not defined;returning default:   FALSE
     Value for IOAPI_OFFSET_64 not defined;returning default:   TRUE
     Value for IOAPI_CFMETA not defined;returning default:   FALSE
     Value for IOAPI_CMAQMETA not defined; returning defaultval ':  'NONE'
     Value for IOAPI_CMAQMETA not defined; returning defaultval ':  'NONE'
     Value for IOAPI_SMOKEMETA not defined; returning defaultval ':  'NONE'
     Value for IOAPI_SMOKEMETA not defined; returning defaultval ':  'NONE'
     Value for IOAPI_TEXTMETA not defined; returning defaultval ':  'NONE'
     Value for IOAPI_TEXTMETA not defined; returning defaultval ':  'NONE'
     Error creating netCDF variable TZONES
     netCDF error number  -47  processing file "IOAPI_INV"
     NetCDF: NC_UNLIMITED in the wrong index
     ERROR: Could not open I/O API inventory file for file name:
          /home/harry.chang/SMOKE/realdata/inventory/zjk/area.zjk.ncf
没看懂什么意思，NC_UNLIMITED in the wrong index，是我的netCDF库有问题还是因为上面清单的问题导致排放信息没有被读入，从而导致了文件无法生成呢？
报错的程序：*** ERROR ABORT in subroutine OPENINVOUT

HarryChang

有没有懂的前辈赐教，指点一下应该怎么修改，或者是哪个地方出了问题。
感激不尽！

如果有人跟我一样是新手，也在试着接触SMOKE的话，欢迎讨论。

lwkraul2019

你有安装ioapi么，提示的错误都是这个

HarryChang

lwkraul2019 发表于 2019-6-28 15:10
你有安装ioapi么，提示的错误都是这个

安装了的，可能是软链接没指对地方，我再看一下

L0725

smoke楼主弄通了吗？关于smoke学习可否给些建议，按照英文manual吗？看了好几天了，内容太多了，一头雾水。

哈尔的移动城堡

L0725 发表于 2019-7-27 11:45
smoke楼主弄通了吗？关于smoke学习可否给些建议，按照英文manual吗？看了好几天了，内容太多了，一头雾水。

同样在学习smoke，测试数据一直有问题，说文件缺失，请问您遇到这种了吗

Kevin_Zhang

查看下清单的格式是否正确！

HarryChang

Kevin_Zhang 发表于 2019-8-11 13:38
查看下清单的格式是否正确！

谢谢回复，现在转战ISAT工具了

HarryChang

哈尔的移动城堡 发表于 2019-8-9 11:36
同样在学习smoke，测试数据一直有问题，说文件缺失，请问您遇到这种了吗

我现在改用ISAT了