SMOKE 处理MEIC清单Temporal子程序运行出错

Alpha_su

我使用SMOKE（4.5版）进行分月份的MEIC（2017数据）多部门排放清单预处理操作，预备将结果输入CMAQ进行分析。但是运行出现报错，出错位置位于Temporal子程序，报错内容截取如下（来自于运行输出的temporal.meic.CN_36km.20190124.log文件)：

1. <p><p>     "IOAPI_DAT" opened as OLD:READ-ONLY   
   
2.      File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/PM2_5.ncf"
   
3.      File type GRDDED3
   
4.      Execution ID "????????????????"
   
5.      Grid name ""
   
6.      Dimensions: 61147 rows, 1 cols, 1 lays, 7 vbles
   
7.      NetCDF ID:        11  opened as READONLY            
   
8.      Time-independent data.
   
9.      Closing file IOAPI_DAT
   
10.      ERROR: Missing or negative emission value(s) in inventory for "PM2_5".
    
11.       
    
12.      "IOAPI_DAT" opened as OLD:READ-ONLY   
    
13.      File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/SO2.ncf"
    
14.      File type GRDDED3
    
15.      Execution ID "????????????????"
    
16.      Grid name ""
    
17.      Dimensions: 59334 rows, 1 cols, 1 lays, 7 vbles
    
18.      NetCDF ID:        11  opened as READONLY            
    
19.      Time-independent data.
    
20.      Closing file IOAPI_DAT
    
21.      ERROR: Missing or negative emission value(s) in inventory for "SO2".
    
22.       
    
23.      "IOAPI_DAT" opened as OLD:READ-ONLY   
    
24.      File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/VOC.ncf"
    
25.      File type GRDDED3
    
26.      Execution ID "????????????????"
    
27.      Grid name ""
    
28.      Dimensions: 61342 rows, 1 cols, 1 lays, 7 vbles
    
29.      NetCDF ID:        11  opened as READONLY            
    
30.      Time-independent data.
    
31.      Closing file IOAPI_DAT
    
32.      ERROR: Missing or negative emission value(s) in inventory for "VOC".
    
33.       
    
34.      "IOAPI_DAT" opened as OLD:READ-ONLY   
    
35.      File name "/home/Alpha_su/SMOKE/install_SMK2/data/inventory/CN_36km/meic.CN_36km_dat/PMC.ncf"
    
36.      File type GRDDED3
    
37.      Execution ID "????????????????"
    
38.      Grid name ""
    
39.      Dimensions: 237 rows, 1 cols, 1 lays, 7 vbles
    
40.      NetCDF ID:        11  opened as READONLY            
    
41.      Time-independent data.
    
42.      Closing file IOAPI_DAT</p><p>     *** **ERROR** ABORT in subroutine TEMPORAL
    
43.      Problem with input data.
    
44. </p>

复制代码

受限于篇幅，报错内容仅选取PM2.5、PMC、VOC以及SO2四个变量。当我使用ncdump -v命令检查各个变量所对应的ncf文件时，发现这些文件中的确都为负值（-9.999e+36）.
而有意思的是，当我将清单的时间尺度从“逐月”变换为“逐年”时（即，在arinv.meic.lst清单列表文件中的Month字段设置为0），程序能够输出结果而不报错。
请问各位有遇到过这种情况吗？
附上自己的assign、run以及清单列表文件。http://bbs.06climate.com/forum.php?mod=attachment&aid=MTA1MzAyfDc0MmFmMzc2ZmVjNDljNzQ4YzhjNzZjMjMxYzU0MjA3fDE3NDQ5ODg0ODc%3D&request=yes&_f=.txt

ydhsola

你好，在你的month文件中，12个月的分配系数都设置为0.1，这样加在一起为1.2，超过了1，我猜是不是这个情况？
此外呢，想了解一下你的GRIDMASK文件，即assign里的GRIDMASK_n4c.nc是如何制作的。因为我一直用的edgar默认的GRIDMASK_EDGAR.ncf。感谢。

Alpha_su

ydhsola，谢谢你的回复！
你说得对，month文件那里的确是设置错了，后来我已修改成12个月份相加为1 （每个月为0.083333），但是还是没有解决问题。后来我继续排查发现，问题应该出在Smkinv程序执行环节，因为该子程序生成/home/Alpha_su/SMOKE/install_SMK2/data/inventory/China_36km/meic.China_36km_dat/路径下的排放清单。但是为什么会导致ARINV清单列表中Month为1-12时所有清单值都是负值（-9.999e+36），而当月份设置为0的时候则不是负值（程序正常运行）？
关于GRIDMASK文件那里，由于之前我使用EDGAR对应的原始GRIDMASK文件运行MEIC清单数据会报错（主要原因可能是两个清单空间分辨率不同），然后参考CMAS论坛自己制作了一个适用于MEIC的GRIDMASK（参考网址：在CMAS论坛中搜索“ERROR ABORT in subroutine OPENMRGIN”主题就能找到）

ydhsola

Alpha_su 发表于 2022-9-8 09:23
ydhsola，谢谢你的回复！
你说得对，month文件那里的确是设置错了，后来我已修改成12个月份相加为1 （每个 ...

十分感谢