ESCOMP / ESCOMP/CAM

HEMCO interface for CAM-chem missing surface emission history fields (SF*)

Open
#1,413 0 comments 0 reactions 0 assignees View on GitHub
bug
Dominant language
No language data
Stars
91
Forks
183
Avg merge
6d 2m
Merged PRs (30d)
9

Description

### What happened?

Reported by Behrooz Roozitalab (ACOM)

When running CAM-chem/MUSICA with HEMCO enabled the surface emission history fields (starting with `SF`) are unavailable.

This is due to `chemistry/mozart/chemistry.F90` initializing the diagnostics based on a `has_emis` array of logical flags in `mo_srf_emissions`:
```fortran
if ( n>0 ) then
if (has_emis(m) .or. aero_has_emis(solsym(m)) .or. ocean_emis_species(solsym(m)) .or. srf_emis_diag(n)) then
srf_emis_diag(n) = .true.
```

However, the HEMCO compatibility layer for surface emissions in `hco_cc_emissions.F90` does not set this flag at the time of `chem_init`. **This does not affect the actual simulation results** but the diagnostics are unavailable. An alternative is to use the `HCO_` history output fields.

### What are the steps to reproduce the bug?

`use_hemco = .true.` in a CAM-chem compset and check field list

### What CAM tag were you using?

cam6_4_122

### What machine were you running CAM on?

CISL machine (e.g. cheyenne)

### What compiler were you using?

Intel

### Path to a case directory, if applicable

_No response_

### Will you be addressing this bug yourself?

Any CAM SE can do this

### Extra info

Copying @lizziel

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.