namelist.input的emiss_inpt_opt、emiss_opt和chem_opt之间是什么关系？应该如何设...

jazzyl

问题1）这三个参数都各自管的是什么？

问题2）在一个有排放数据的WRF/Chem模拟中，这三个参数在WRF/Chem处理排放中被读取的顺序是什么，分别起到了什么作用？其背后的运行逻辑是什么？

请大神们和坛友指导。涉及代码的讨论也可以。下面是我目前对这两个问题的一些理解。

问题1）这三个参数都各自管的是什么？

chem_opt设置了WRF/Chem模式的化学机制。

emiss_opt设置the correct emission input option (比如emiss_opt=3 for NEI emissions, or emiss_opt=5 for RETRO/EDGAR) for the anthropogenic emissions，还有其他的数字，比如6表示GOCART simple emissions，7表示MOZART emissions.

官方tutorial说The choice of emiss_opt depends upon the emissions data arrays that are going to beread in. 使用合适的emiss_opt将能保证所有的排放化学组分被读取。

emiss_inpt_opt设置的是emission speciation的类型，0表示no emissions data read，1表示emissions are speciation for RADM2/SORGAM. Recommended when using the NEI-05 or EDGAR/RETRO emissions speciated for RADM2 chemical mechanism, 还有其他数字表示不同化学机制下的emissionspeciation的类型。

但是如果已经设置了chem_opt（化学机制）和emiss_opt（排放数据的类型），那其实就已经决定了emission speciation，还要emiss_inpt_opt这个参数干嘛呢？

问题2）在一个有排放数据的WRF/Chem模拟中，这三个参数在WRF/Chem处理排放中被读取的顺序是什么，分别起到了什么作用？其背后的运行逻辑是什么？

第一步，首先使用convert_emiss.exe生成排放数据文件时程序中根据emiss_opt的设置以不同方式读取binary emission dataset，然后输出wrfchemi数据文件。这一步骤可以在chem/convert_emiss.F代码中发现：

1. if(config_flags%emiss_opt == ecptec)
   
2. …
   
3. grid%emis_ant(ips:ipe,kps,jps:jpe,p_e_pm_25)=dumc1(ips:ipe,jps:jpe)
   
4. …
   
5. end if

复制代码

第二步，wrf.exe在运行过程中调用share/mediation_integrate.F中的med_read_wrf_chem_emissions函数读取wrfchemi数据文件。

1. CALL construct_filename1 (inpname, 'wrfchemi_00z' , grid%id , 2 )
   
2. WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Open file ',TRIM(inpname)
   
3. …
   
4. CALL open_r_dataset ( grid%auxinput5_oid, TRIM(inpname) , grid , config_flags, & "DATASET=AUXINPUT5", ierr )

复制代码

第三步，wrf.exe在运行过程中调用chem/module_emissions_anthropogenics.F中subroutine add_anthropogenics 函数开始根据chem_opt的不同对不同化学成分加上emission中各个变量进行计算。

1. if( config_flags%chem_opt == MOZCART_KPP ) then
   
2.              chem(i,k,j,p_p10) = chem(i,k,j,p_p10) + conv_rho_aer*emis_ant(i,k,j,p_e_pm_10)
   
3.              chem(i,k,j,p_p25) = chem(i,k,j,p_p25) + conv_rho_aer*emis_ant(i,k,j,p_e_pm_25)
   
4. end if

复制代码

我感觉这三个步骤还很不完善，中间有很多缺失的地方，而且我也没有在代码中找到emiss_inpt_opt这个参数。欢迎各位加入讨论！

天下第贰

讲的很好，我也一直不懂这个emiss_inpt_opt是干什么的

jazzyl

@清风 @andrewsoong  
请大神们不吝指导。

jazzyl

@astiny
大牛，能否一起来讨论和指导一下。

astiny

呜呜。我都两年没摸wrfchem了。而且我也不是wrfchem巨子。
我是这么理解的：首先chem_opt肯定就是wrfchem里面运行时候的化学机制，我相信这个应该是没有问题的，因为如果你其它选项都不变，只改这个，那么输出的变量都会变的（比如mosaic-8bin变成4bin之类）。
其次，emiss_opt好像是convert_emiss.exe生成排放数据文件的选项，这也没有问题。
然而更多的时候，排放源文件并不一定是convert_emiss.exe生成的。有时候会加上自己的一些东西（点源、日变化等等），这些东西很多都是自己写程序，生成nc文件的。这时候就不需要emiss_opt这个选项，有他没他都一样。
然后，emiss_inpt_opt这个选项是决定了你怎么读排放源文件的。这个就取决于你的排放源文件具体是怎么做的了。这个选项才是真实的，最终决定你怎么读取排放源文件的选项。如果你不知道怎么做排放源的时候，你就可以参考emiss_inpt_opt这段程序关联的数据读取方法，定制化自己的排放源文件，里面的变量名、单位等等信息。

mico2008

涨知识了。感谢

kermit

涨知识了。感谢

planck

一般来讲emiss_inpt_opt都会和emiss_opt对应

xlilium

astiny 发表于 2020-6-2 11:17
呜呜。我都两年没摸wrfchem了。而且我也不是wrfchem巨子。
我是这么理解的：首先chem_opt肯定就是wrfchem ...

感谢大神的讲解，有个问题是emiss_opt都有按照saprc99排放的选项，但是emiss_inpt_opt里面并没有诶，请问如果要用saprc99机制的话emiss_inpt_opt要选哪个呢？谢谢！

mohouwa

xlilium 发表于 2020-6-15 16:44
感谢大神的讲解，有个问题是emiss_opt都有按照saprc99排放的选项，但是emiss_inpt_opt里面并没有诶，请问 ...

你好！我最近也是在用SAPRC99的机制进行模拟，目前也是不大清楚这一步怎么选择，请问您现在解决了这个问题了吗？