search

metno / emep-ctm

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

Chemical group question. #37

Open mvieno opened 6 years ago

mvieno commented 6 years ago

Hi, in CM_ChemGroup PM10 group is defined as follow:

integer, public, target, save, dimension(26) :: &
  PM10_GROUP = (/ SO4,NO3_F,NO3_C,NH4_F,ASH_F,ASH_C,POM_C_FFUEL,EC_F_WOOD_NEW,EC_F_WOOD_AGE,EC_C_WOOD,EC_F_FFUEL_NEW,EC_F_FFUEL_AGE,EC_C_FFUEL,REMPPM25,REMPPM_C,FFIRE_BC,FFIRE_REMPPM25,OM25_P,SEASALT_F,SEASALT_C,DUST_ROAD_F,DUST_ROAD_C,DUST_WB_F,DUST_WB_C,DUST_SAH_F,DUST_SAH_C /)

whereas in the DDEP_PM10_GROUP and WDEP_PM10_GROUP the ON25_P is missing? probably this is wanted.

thanks Massimo

avaldebe commented 6 years ago

Hi @mvieno, @mifads

As far as I understand OM25_P is not a "real variable", but the sum of all OM particulate contributions from the SOA species. @mifads am I correct? Should OM25_P contribute to PM10 dry/wet depositions?

mifads commented 6 years ago

OM25_p is not a 'real' species, but rather a very handy sum of all inert and condensed OM species. This is used as for example a compound like ASOC_ug1 has both gas and particle phase contributions. We don't want the gas-phase part counted as PM2.5 or PM10.

When calculating e.g. dry deposition, the code does not use OM25_p, but operates on each of the 'real' species separately. It calculates how much of e.g. ASOC_ug1 is gas and deposits that as gas, and the rest as particle.

mvieno commented 6 years ago

Hi, thanks. That's clear now.

I am still not clear in which specie, let say, the dry dep for the particulate part of ASOC_ug1, is then accounted for in the DDEP_PM10_GROUP. As both PM10_GROUP and DDEP_PM10_GROUP do include the same variables but OM25_p.

integer, public, target, save, dimension(25) :: &
  DDEP_PM10_GROUP = (/ SO4,NO3_F,NO3_C,NH4_F,ASH_F,ASH_C,POM_C_FFUEL,EC_F_WOOD_NEW,EC_F_WOOD_AGE,EC_C_WOOD,EC_F_FFUEL_NEW,EC_F_FFUEL_AGE,EC_C_FFUEL,REMPPM25,REMPPM_C,FFIRE_BC,FFIRE_REMPPM25,SEASALT_F,SEASALT_C,DUST_ROAD_F,DUST_ROAD_C,DUST_WB_F,DUST_WB_C,DUST_SAH_F,DUST_SAH_C /)

Massimo

avaldebe commented 5 years ago

@mifads ?

avaldebe commented 5 years ago

@mifads can you answer this question?

mifads commented 5 years ago

Looks like a bug for PM outputs, but probably a small one. We do account for dry dep of organics, though it is complicated by the gas-particle partitioning. As we never output deposition of PM ourselves I have never come across this. Can be added to the (long) post-summer to-do list.

avaldebe commented 3 years ago

@mifads

was this addressed on the code we're about to release?

Cheers, Á.

mifads commented 3 years ago

Looks like the bug is still there. In short, the model does deposit e.g. ASOC_ug1 with both its gas and particle phases, but these are not included in those output DDEP and WDEP PM groups. Hopefully we can fix it next time around!