search

metno / emep-ctm

Open Source EMEP/MSC-W model
GNU General Public License v3.0
29 stars 19 forks source link

How to output short lived gas concentrations #24

Closed PaulHamer closed 6 years ago

PaulHamer commented 7 years ago

I am trying to diagnose the behaviour of the chemistry in EMEP. Specifically I want to create some diagnostics offline that given information on the reaction budgets for some of the important chemical species in the model. In order to do this, I would like to be able to output the concentrations of the different species beyond the default list in the config_emep.nml input file. I have tried a few tests trying to output OH concentrations, and I have gotten as far as testing this in the OutputConcs namelist:

'OH' ,'molec/cm3' ,'2d','AIR_CONCS','SHL' ,'YMDI',

However, when I try this I get the following error message after EMEP crashes:

NetCDF: Name contains illegal characters ERRMSG: def2d:SURF_molec/cm3_OH STOP-ALL ERROR: Error in netcdf routine

Please can anyone explain how to output chemical species like OH, and the ones in the list below? Thanks in advance.

In the end, I would like to output the following list of chemical species from the chemical mechanism used by EMEP: NO3, HO2, CH3O2, C2H5O2, CH3COO2, SECC4H9O2, MEKO2, ETHRO2, PRRO2, OXYO2, MALO2, ISRO2, MVKO2, MACRO2, MACO3, ISNIR, ISONO3, OP, and (as previously mentioned) OH.

I should clarify, I am using EMEP version rv4.10.

avaldebe commented 7 years ago

Unit conversion is not supported for short lived species (SHL). Therefore, the unit you give will be used for meta data and for the variable name. The runtime error message you get is from the netCDF library telling you that you can not name a variable SURF_molec/cm3_OH. Try with the following definition:

'OH' ,'mcm3' ,'2d','AIR_CONCS','SHL' ,'YMDI',
PaulHamer commented 7 years ago

Thanks, that seems to work.