WRF-Chem中CLM陆面模式耦合的MEGAN2.1疑问

17732591679

各位大佬大家好，近期刚刚接触WRF-Chem，MEGAN新手，想使用里面的MEGAN2.1进行模拟现遇到两个疑问：

1. WRF-Chem中CLM耦合的MEGAN2.1需要一个排放因子输入文件，如何获取该文件
2. 如何让WRF-Chem读取该排放因子文件

问题1补充：
git的更新日志上写了需要读入一个排放因子文件，但我没能找到这个文件的相关信息，也没有找到示例文件。
github.com/wrf-model/WRF/pull/1425

Among the key changes to module_sf_clm.F are:

• a MEGANv2.1 emissions factor file containing standard emissions as a function of PFT and 20 classes of compounds is read in
• new modules for handling the MEGAN categories are created (and the use of these modules by other modules necessitated the re-ordering of the subroutines in the source file)
• the environmental emissions factors are tracked through running averaged variables
• PFTs can now be determined from a geographical mapping in the wrfinput file rather than from an internal vegetation type-based table - the code is looking for a variable of name 'PCT_PFT' in the wrfinput file; it is up to the user to use external software to insert the variable containing the percentages for each CLM PFT
• mappings of MEGAN species to WRF-Chem species packages are included.
  

问题2补充：
在WRF phys下的源代码module_sf_clm.F中有关于MEGAN2.1的namelist描述，我将他添加到了WRF运行路径下的namelist.input中，WRF似乎仍未尝试读取对应路径下的文件
  !-------------------------------------------------------------------------
  !
  ! This reads the megan_emis_nl namelist group in drv_flds_in and parses the
  ! namelist information for the driver, CLM, and CAM.
  !
  ! Namelist variables:
  !   megan_specifier, megan_mapped_emisfctrs, megan_factors_file
  !
  ! megan_specifier is a series of strings where each string contains one
  !  CAM chemistry constituent name (left of = sign) and one or more MEGAN
  !  compounds (seperated by + sign if more than one).  The specification of
  !  the MEGAN compounds to the right of the = signs tells the MEGAN VOC
  !  model within CLM how to construct the VOC fluxes using the factors in
  !  megan_factors_file and land surface state.
  !
  ! megan_factors_file read by CLM contains valid MEGAN compound names,
  !  MEGAN class groupings and scalar emission factors
  !
  ! megan_mapped_emisfctrs switch is used to tell the MEGAN model to use
  !  mapped emission factors read in from the CLM surface data input file
  !  rather than the scalar factors from megan_factors_file
  !
  ! Example:
  ! &megan_emis_nl
  !  megan_specifier = 'ISOP = isoprene',
  !     'C10H16 = myrcene + sabinene + limonene + carene_3 + ocimene_t_b + pinene_b + ...',
  !     'CH3OH = methanol',
  !     'C2H5OH = ethanol',
  !     'CH2O = formaldehyde',
  !     'CH3CHO = acetaldehyde',
  ! ...
  !  megan_factors_file = '$datapath/megan_emis_factors.nc'
  ! /
  !-------------------------------------------------------------------------

WRF-Chem运行报错信息：
  module_start: phy_init: Before call to bl_init
  clminit: beginning for domain =  1
megan_specifier diagnostic in clminit
megan_specifier(1)
isopr=isoprene+2.4*MBO_2m3e2ol+1.25*MBO_3m2e1ol+1.25*MBO_3m3e1ol
megan_specifier(2)
c10h16=myrcene+sabinene+limonene+carene_3+ocimene_t_b+pinene_b+pinene_a+2met_styrene+phellandrene_a+thujene_a+terpinene_a+terpinene_g+terpinolene+phellandrene_b+camphene+bornene+fenchene_a+ocimene_al+ocimene_c_b+tricyclene+estragole+piperitone+linalool+terpineol_4+terpineol_a+1.25*linalool_OXD_c+1.25*linalool_OXD_t+1.4*ionone_b+1.5*farnescene_a+1.5*caryophyllene_b+1.5*acoradiene+1.5*aromadendrene+1.5*bergamotene_a+1.5*bergamotene_b+1.5*bisabolene_a+1.5*bisabolene_b+1.5*bourbonene_b+1.5*cadinene_d+1.5*cadinene_g+1.5*cedrene_a+1.5*copaene_a+1.5*cubebene_a+1.5*cubebene_b+1.5*elemene_b+1.5*farnescene_b+1.5*germacrene_B+1.5*germacrene_D+1.5*gurjunene_b+1.5*humulene_a+1.5*humulene_g+1.5*isolongifolene+1.5*longifolene+1.5*longipinene+1.5*muurolene_a+1.5*muurolene_g+1.5*selinene_b+1.5*selinene_d+1.5*nerolidol_c+1.5*nerolidol_t+1.4*geranyl_acetone+1.4*terpinyl_ACT_a+1.1*2met_nonatriene+1.2*jasmone+1.6*met_jasmonate
megan_specifier(10)
c2h5oh=ethanol
megan_specifier(20)
so2=1.53*diallyl_2s+2met_2s+0.5*hydrogen_s+0.75*met_mercaptan+1.8*met_propenyl_2s+2.3*PPPP_2s
megan_factors_file

clminit: B4 call to megan_factors_init
megan_factors_init: failed to open file
Invalid argument
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

Mrmake1

请问这个问题您解决了吗，可不可以加个好友探讨一下这个问题