Point source emissions not impacting output CMAQ 5.3.3

nderon · February 27, 2023, 2:11pm

Hello,

I am working on a CMAQ run using 2017 EQUATES emission data. I’ve successfully run the model, however the output did not match expectations. I reran the model with most diagnostic options turned on to find that the emission diagnostic files were all zero for my seven point source sectors.

However in the run logs (attached), the point source emissions appear to be read in correctly and scaled correctly. I’ve also confirmed that the point source files do contain correct emissions data, i.e. they are not zeroed out.

Feel like I’m missing something obvious here, but I cannot find anything amiss in the log files. Any help is greatly appreciated.

CTM_LOG_000.TXT (5.3 MB)
run_cctm.csh (36.2 KB)

GNU Fortran (Ubuntu 9.4.0-1ubuntu1~20.04.1) 9.4.0

lizadams · February 27, 2023, 3:27pm

Here is a link to a benchmark case that uses the EQUATES data and that contains build and run scripts. You could compare these scripts to what you are using to see if there are any differences that would explain your findings.

Here is a link to the DESID documentation for CMAQv5.3.3.

github.com

USEPA/CMAQ/blob/5.3.3/DOCS/Release_Notes/CMAQv5.3.1_DESID_families.md

# DESID Families
[Ben Murphy](mailto:murphy.ben@epa.gov), U.S. Environmental Protection Agency

## Brief Description
The Detailed Emission Scaling Isolation and Diagnostic (DESID) module is meant to make it trivial to modify emissions online in CMAQ through invoking complex scaling tasks with intuitive syntax for defining rules. For example, users have control over scaling individual compounds for each emission stream independently and over custom geographical areas of the model domain. One unattractive feature was that a lot of repetitive rules would be needed for complex, but easy-to-articulate tasks. For example, if a user wanted to scale NOx (NO + NO2) for 4 emission streams to 50%, then 8 rules would be required in the emission control file:  
```
! Region      | Stream Label  |Emission | CMAQ-        |Phase/|Scale |Basis |Op  
!  Label      |               |Surrogate| Species      |Mode  |Factor|      |
'EVERYWHERE'  , 'PT_EGU'      ,'NO2'    ,'NO2'         ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'GRIDDED'     ,'NO2'    ,'NO2'         ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'MOBILE'      ,'NO2'    ,'NO2'         ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'PT_NONEGU'   ,'NO2'    ,'NO2'         ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'PT_EGU'      ,'NO'     ,'NO'          ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'GRIDDED'     ,'NO'     ,'NO'          ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'MOBILE'      ,'NO'     ,'NO'          ,'GAS' ,0.50 ,'UNIT','m',
'EVERYWHERE'  , 'PT_NONEGU'   ,'NO'     ,'NO'          ,'GAS' ,0.50 ,'UNIT','m',
```  
This pull request implements family definitions for emission streams, chemical families, and geographical regions. So, if for example, the user defines the stream family "CONTROLLED_SOURCES" as the group made up by "PT_EGU","GRIDDED","MOBILE" and "PT_NONEGU" and the chemical family "NOX" made up by "NO" and "NO2", then the task above can be accomplished with one rule:
```
! Region      | Stream Label         |Emission | CMAQ-        |Phase/|Scale |Basis |Op

This file has been truncated. show original

Link to the tutorial for DESID for CMAQv5.3.3

github.com

USEPA/CMAQ/blob/5.3.3/DOCS/Users_Guide/Tutorials/CMAQ_UG_tutorial_emissions.md

## CMAQ Tutorial ##
### Prescribing Emissions Using the DESID, the Detailed Emissions Scaling, Isolation and Diagnostics Module ###
Purpose: This tutorial will guide users to utilizing the Emission Control namelist to perform some basic manipulation
of their emission streams. For additional questions, contact Ben Murphy (murphy.ben@epa.gov).

------------

### Definitions of Terms
- Stream: An Online emission source or a group of sources processed offline and read into CMAQ from one file. Common
examples of streams include biogenic VOCs, mobile sources, wind-blown dust, prescribed fires, electric-generating units,
residential heating, etc.  
- Species: A variable representing a chemical compound or group of compounds in CMAQ.  
- Surrogate: A variable representing a chemical compound or group of compounds on an emission Stream.  

### Important Notes to Remember
- Rules are applied in the order they are provided.

### Example Use Cases
- [1. Zero out emissions](#zero_out)  
- [2. Add emissions for a new tracer species](#add_emissions)

This file has been truncated. show original

Ben_Murphy · February 27, 2023, 3:53pm

Hi nderon,

Would you consider sending us also your EmissCtrl file so we can see the commands you’ve used to map the chemical species?

Best wishes,
Ben Murphy

nderon · February 27, 2023, 7:37pm

Find the EmisCtrl file attached. I couldn’t find any differences between the EQUATES benchmark scripts and my own that might explain this, but I’ll go back and take another pass at it.

EmissCtrl_cb6r3_ae7_aq.txt (31.0 KB)

Ben_Murphy · February 27, 2023, 8:35pm

You’re right, that EmissCtrl file looks fine to me.

Try setting the value of both EMISDIAG and EMISDIAG_SUM in your runscript to 2DSUM. Does that yield results in the point files more in line with what you expect?

-Ben

nderon · February 28, 2023, 1:42pm

Unfortunately I get the same output–zeroed out point sources.

Ben_Murphy · February 28, 2023, 3:30pm

Thanks for reporting back. Let’s inspect the contents of your point source files. Can you send a text file dump of the header information from one point source file (i.e. ncdump output)?

-Ben

nderon · February 28, 2023, 5:57pm

Sure thing. I’ve attached a file of an emissions header and the associated stack group header.

nonipm_stack_headers.txt (9.6 KB)
nonipm_headers.txt (20.2 KB)

nderon · March 8, 2023, 8:27pm

Any chance you had an opportunity to take a look at the headers provided? Any help is greatly appreciated.

Topic		Replies	Views
Modify model-ready emissions from EQUATES by conducting perturbation or adding noise CMAQ equates	1	432	December 14, 2023
Turning off emissions sources in CMAQ 5.3 v2 Run Time Errors and Issues	8	1035	November 9, 2020
Issues with CCTM output CMAQ	7	818	July 30, 2020
Pre-merged EQUATES emission data for 2019 CMAQ	3	27	October 7, 2024
Analysis of CMAQv5.3 Output Files CMAQ	1	314	July 20, 2020

Point source emissions not impacting output CMAQ 5.3.3

Related topics