Hi there,
I am using CMAQv5.4 with enabling MEGAN, and how to output the MEGAN mass emissions diagnostic files similar to enabling BEIS modules?
Any recommendations for how to solve this or investigate further?
Thanks!
Hi,
You can use the DESID diagnostic output framework described in section B.3.7.2 here: CMAQ/CMAQ_UG_appendixB_emissions_control.md at main · USEPA/CMAQ · GitHub
For the Stream name, use ‘MIOG’ for MEGAN. For BEIS, use ‘BIOG’.
Best wishes,
Ben Murphy
3 Likes
Thanks, @Ben_Murphy for responding so quickly!
To add to this response, please note that the diagnostic emission output files generated through the DESID framework will be mole based rather than mass based (unless I am missing an option in DESID to control that), so if you want to compare MEGAN diagnostic mole-based emissions from the DESID1_MIOG file to existing non-DESID B3GTS mass-based BEIS diagnostic emissions generated using the B3GTS_DIAG flag, you’ll have to do some unit conversions using molecular weights, though for some lumped mechanism species there may be some issues with that. In that case, it’d be cleaner to also generate BEIS mole-based diagnostic emission files through DESID.
Here is an example of the modification to CMAQ_Control_DESID.nml needed to implement MEGAN diagnostic mole-based emission files through DESID, as documented in the link posted by @Ben_Murphy above:
> diff CMAQ_Control_DESID.nml CMAQ_Control_DESID_MEGAN.nml
> 188c188
> < Desid_N_Diag_Rules = 0 ! Exact Number of Diagnostic Rules Below
> ---
> > Desid_N_Diag_Rules = 1 ! Exact Number of Diagnostic Rules Below
> 196c196,200
> <
> ---
> >
> > Desid_Diag_Streams_Nml(1,:)= 'MIOG'
> > Desid_Diag_Fmt_Nml(1) = 'COLSUM' ! Options: 2D, 2DCOL, 3D
> > Desid_Diag_Spec_Nml(1,:) = 'ALL'
> >
3 Likes
Thx. Based on your suggestions, CCTM_DESID1_MIOG file can be created.
Hi, @hogrefe.christian.
I still cannot create the CCTM_DESID1_MIOG files. And the CTM files show that the name of first variable is always mistaken code whether I choose ‘ALL’ or ‘ISOP’, 'TERP'.
The figure shows that log when I set Desid_Diag_Spec_Nml(1, :) = 'ALL'.
Yikes. I’m sorry for this frustration. I’ve tried to reproduce this error using our Benchmark scenario and everything seems to be working ok. Can you post your version of the CMAQ_Control_DESID.nml file that is prescribing the diagnostic output?
Thx, I have solved this problem by using another computer.
Hm, I’m not sure why the behavior would go away on a new system. Perhaps the reinstall/recompile step ironed out potential problems. If you end up looking further into this a figuring out system-to-system issues or differences, please report them if possible.
Thank you and happy coding!
Ben
Hi, Hogrefe
I have run CMAQ5.4-MEGAN with saprc07tc_ae6_aq mechanism to produce CCTM_DESID1_MIOG files. However, I couldn’t output Emission Rate of ISOPRENE from MIOG emissions no matter setting the Desid_Diag_Spec_Nml(1, ：) = ‘ALL’
or Desid_Diag_Spec_Nml(1,：) =‘NO’,‘ISOPRENE’,‘TERP’,‘PROPENE’,‘SESQ’,‘APIN’. Here’s the variable list log file of CCTM_DESID1_MIOG files.
DESID_ALL.txt (16.2 KB)
DESID_specific.txt (6.8 KB)
Given that you’re getting output for some but not all species in the CCTM_DESID1_MIOG files even when specifying ‘ALL’, I’m wondering whether this points to some issue with how the MEGAN-internal emissions for the missing species (like isoprene) are mapped to saprc07tc_ae6_aq. I’ll reach out to some of my colleagues who know more about this.
Thanks for your response. I’m looking forward to solve the issue.
Hi Hogrefe,
The issue might be attributed to settings in my CMAQ_Control_DESID* file. I have commented out following line:
!‘EVERYWHERE’, ‘ALL’ ,‘ISOPRENE’ ,‘ISOPRENE’ ,‘GAS’ ,1. ,‘UNIT’,‘a’,
And I rewrote new rules to output ISOPRENE emission rate.
‘EVERYWHERE’, ‘MIOG’ ,‘ISOPRENE’ ,‘ISOPRENE’ ,‘GAS’ ,1. ,‘UNIT’,‘a’,
This is certainly unexpected behavior. Would you mind sharing one of the CTM_LOG files for each of the cases:
(1) the run with the default DESID rules and diagnostic output set to ‘ALL’
(2) the run with the default DESID rules and diagnostic output set to ‘NO’,‘ISOPRENE’,‘TERP’,‘PROPENE’,‘SESQ’,‘APIN’
(3) the run with the DESID rule changed to specify MIOG for the ISOPRENE scaling?
Best wishes,
Ben
Hi Ben,
Thanks for your kind reply. I think it might be caused by my carelessness. Here’s the emission checking information and the DESID file before and after modification:
- Before
|> Checking for unused Emissions Variables: ============================================ GRIDDED_EMIS | Gridded Area Emissions File 1: 19 unused variables. ARO1 BC CH4 CO2 ETHE IPRD ISOPRENE N20 NONR NOX NROG NVOL OC PHEN PM2_5 PMFINE PRD2 TRP1 VOC MIOG | Megan Biogenic Emissions: 19 unused variables. CH4 ISOPRENE TERP BENZ ARO1 ARO2MN PRD2 CCOOOH IPRD HCN CH3BR CH3CL CH3I TRS N2O HCHO_PRIMARY CCHO_PRIMARY OCS PYRYACID SEASPRAY | Sea Spray Aerosol Emissions: 18 unused variables. PMFINE_CR_VI PMFINE_NI PMFINE_AS PMFINE_BE PMFINE_PHG PMFINE_PB PMFINE_CD PMFINE_MN_HAPS PMFINE_BR PMCOARSE_CR_VI PMCOARSE_NI PMCOARSE_AS PMCOARSE_BE PMCOARSE_PHG PMCOARSE_PB PMCOARSE_CD PMCOARSE_MN_HAPS PMCOARSE_BR NOTE: Some Emissions Variables are unused by your current emission control configuration. You may want to include these emissions if they are relevant to your application.
I commented out ‘ISOPRENE’ for ‘ALL' emission streams because there is ISOPRENE emission calculated by MEGAN V2.1 in my gridded emission files (varlist as bellow), thus no ISOPRENE emission information was output. Lately, I modified the rules to
‘EVERYWHERE’, ‘MIOG’ ,‘ISOPRENE’ ,‘ISOPRENE’ ,‘GAS’ ,1. ,‘UNIT’,‘a’,
and I got the ISOPRENE emission rate.
CMAQ_Control_DESID_saprc07tc_ae6_aq_megan.txt (15.5 KB)
emis_vars.txt (58.7 KB)
- After
CMAQ_Control_DESID_saprc07tc_ae6_aq_megan_modified.txt (16.1 KB)
1 Like
Ah ok. I understand better now that you have Isoprene on your gridded file that you would like to ignore while also using Isoprene from MEGAN. So is everything lining up now with your expectations? Are there parts of the issue that are still unanswered or need to be addressed?
Best,
Ben
Hi Ben，
Yes, the problem has been completely resolved. Thank you again for your kind attention. ::
1 Like