[BUG/ISSUE] Online MEGAN diagnostics not working properly

[yantosca]

Dylan Millet wrote:

I am using the online version of the megan emissions and I see that the standard HEMCO emission diagnostics do not work for this configuration.

There is a separate place in the megan module where emissions for some (non-standard) emissions are being saved out via Diagn_Create.

It seemed to me that the most straightforward way to save out the online MEGAN emissions for isoprene and other advected VOCs would be to add in additional calls to Diagn_Create.

However this is complicated by the fact that Diagn_Create is looking for single-precision input and many of the flux variables (e.g., Inst%FLUXISOP) are cast as double precision.

For my purposes I have kludged this by changing the relevant variables in megan_mod to single precision. But in general it might be worth creating an interface for Diagn_Create so that it can accept either sp or dp.

[stale]

This issue has been automatically marked as stale because it has not had recent activity. If there are no updates within 7 days it will be closed. You can add the "never stale" tag to prevent the Stale bot from closing this issue.

[drnimbusrain]

@yantosca @msulprizio

I have also been having some recent issues using the latest develop branch of the standalone HEMCO MEGANv2.1 extension (https://github.com/geoschem/HEMCO/blob/dev/src/Extensions/hcox_megan_mod.F90) and seemingly the diagnostic outputs for monoterpenes, e.g., MTPA, MTPO, LIMO, and SESQ.

In a test run using MERRA2 driving meteorology, all other MEGAN diagnostic variables (e.g., ISOP, ACET, etc.) seem to have good output data for my test run, but the four lumped monoterpene variables are always set to zero for all hours, with the following warnings in the my attached HEMCO log file:

HEMCO WARNING: Diagnostics counter is zero - return empty array: SESQ_MEGAN
--> LOCATION: DiagnCont_PrepareOutput (hco_diagn_mod.F90)
--- Added diagnostics: SESQ_MEGAN
HEMCO WARNING: Diagnostics counter is zero - return empty array: LIMO_MEGAN
--> LOCATION: DiagnCont_PrepareOutput (hco_diagn_mod.F90)
--- Added diagnostics: LIMO_MEGAN
HEMCO WARNING: Diagnostics counter is zero - return empty array: MTPO_MEGAN
--> LOCATION: DiagnCont_PrepareOutput (hco_diagn_mod.F90)
--- Added diagnostics: MTPO_MEGAN
HEMCO WARNING: Diagnostics counter is zero - return empty array: MTPA_MEGAN
--> LOCATION: DiagnCont_PrepareOutput (hco_diagn_mod.F90)

I also attach my HEMCO_sa MEGANv2.1 configuration files for this run.

HEMCO.log
HEMCO_Config.txt
HEMCO_sa_Diagn.txt
HEMCO_sa_Grid.txt
HEMCO_sa_Spec.txt
HEMCO_sa_Time.txt