Some online MEGAN species cannot be saved to diagnostics in standalone HEMCO

[lwcugb]

Name and Institution (Required)

Name: Wei Li
Institution: NOAA ARL

Description of your issue or question

I am running HEMCOv3.8 standalone to save out natural emissions from MEGAN and soil NOx. The issue is I cannot save out MOH, CH2O, and C2H4 to the diagnostics. I have turned them on in the HEMCO_Diagn.rc file like below. It works for me in an old version v3.6.2.

InvMEGAN_MOH         MOH    108  -1  -1  2   kg/m2/s  Methanol_emission_flux_from_MEGAN_extension
InvMEGAN_CH2O       CH2O    108  -1  -1  2   kg/m2/s  Formaldehyde_emission_flux_from_MEGAN_extension
InvMEGAN_C2H4       C2H4    108  -1  -1  2   kg/m2/s  C2H4_biogenic_emission_flux_from_MEGAN_extension

I am using offline soil NOx for now because the online version needs 'DRYCOEFF' parameter to be read from the Olson_Land_Map_201203/Olson_2001_Drydep_Inputs.nc file. I am not sure how I should set it up to read the file in the standalone simulation. Any advice would be greatly appreciated!

Thanks,
Wei

[yantosca]

Thanks for writing @lwcugb. Would you be able to attach your entire HEMCO_sa_Config.rc, HEMCO_Config.rc, and HEMCO_Diagn.rc config files to this issue? That would help us diagnose the problem.

[lwcugb]

Sorry about that. I attached them below. I do not use HEMCO_sa_Config.rc.
HEMCO_Config.rc.megan.txt
HEMCO_Diagn.rc.txt

[yantosca]

Thanks for the config files @lwcugb. Try the following.

1. CH2O and C2H2 are not listed as MEGAN species. You'll need to add these in HEMCO_Config.rc:

108     MEGAN                  : on    ISOP/ACET/PRPE/C2H4/ALD2/MOH/EOH/MTPA/MTPO/LIMO/SESQ/SOAP/SOAS/C2H2/CH2O

2. Add these lines to HEMCO_Diagn.rc

EmisCH2O_Biogenic  CH2O   108    -1  -1   2   kg/m2/s  CH2O_emission_flux_from_biogenic_sources
EmisC2H2_Biogenic  CH2O   108    -1  -1   2   kg/m2/s  C2H2_emission_flux_from_biogenic_sources
EmisMOH_Biogenic   MOH    108    -1  -1   2   kg/m2/s  MOH_emission_flux_from_biogenic_sources

This should save out the diagnostics properly. If not, then let me know!

[lwcugb]

It does not seem to work. I added CH2O to the species list and the model crashed with an error message:

HEMCO ERROR: Invalid species names: CH2O

I do not know why this happened because I do see CH2O in the hcox_megan_mod.F90 file.

Thanks,
Wei

[yantosca]

@lwcugb: You might need to add CH2O to the HEMCO_sa_Spec.rc file. That tells HEMCO that CH2O is a species. Also you'd need to copy the name and molecular weight from the GEOS-Chem species_database.yml file.

[lwcugb]

I added C2H4 and MOH to the HEMCO_sa_Spec.rc file and now they can be saved. CH2O is already in the file. I guess CH2O may not be a MEGAN emission species in HEMCO. I will not save it.

[yantosca]

Thanks for the feedback! CH2O is definitely an emission species in HEMCO. You should be able to archive it.

[lwcugb]

I am only running biogenic emissions, and I mean CH2O is not emitted from MEGAN. Although there is emission factor of CH2O in the MEGAN extension codes, I did not find its flux array.

HEMCO/src/Extensions/hcox_megan_mod.F90

Lines 198 to 229 in 4e81881

	! Emission arrays (kg/m2/s)
	REAL(hp), POINTER :: FLUXISOP(:,:)
	REAL(hp), POINTER :: FLUXMONO(:,:)
	REAL(hp), POINTER :: FLUXACET(:,:)
	REAL(hp), POINTER :: FLUXPRPE(:,:)
	REAL(hp), POINTER :: FLUXC2H4(:,:)
	REAL(hp), POINTER :: FLUXLIMO(:,:)
	REAL(hp), POINTER :: FLUXMTPA(:,:)
	REAL(hp), POINTER :: FLUXMTPO(:,:)
	REAL(hp), POINTER :: FLUXSESQ(:,:)
	REAL(hp), POINTER :: FLUXSOAP(:,:)
	REAL(hp), POINTER :: FLUXSOAS(:,:)
	REAL(hp), POINTER :: FLUXALD2(:,:)
	REAL(hp), POINTER :: FLUXMOH(:,:)
	REAL(hp), POINTER :: FLUXEOH(:,:)
	! Emission arrays for use in diagnostics (kg/m2/s)
	REAL(sp), POINTER :: FLUXACETmb(:,:)
	REAL(sp), POINTER :: FLUXACETbg(:,:)
	REAL(sp), POINTER :: FLUXAPIN(:,:)
	REAL(sp), POINTER :: FLUXBPIN(:,:)
	REAL(sp), POINTER :: FLUXSABI(:,:)
	REAL(sp), POINTER :: FLUXMYRC(:,:)
	REAL(sp), POINTER :: FLUXCARE(:,:)
	REAL(sp), POINTER :: FLUXOCIM(:,:)
	REAL(sp), POINTER :: FLUXOMON(:,:)
	REAL(sp), POINTER :: FLUXFARN(:,:)
	REAL(sp), POINTER :: FLUXBCAR(:,:)
	REAL(sp), POINTER :: FLUXOSQT(:,:)
	REAL(sp), POINTER :: FLUXMBOX(:,:)
	REAL(sp), POINTER :: FLUXFAXX(:,:)
	REAL(sp), POINTER :: FLUXAAXX(:,:)

.