DailyEmisOut = T about emep-ctm HOT 4 OPEN

In https://emep-ctm.readthedocs.io/en/latest/Output.html it is described that to output emissions the option DailyEmisOut = T need to be set.

The documentation refers to a new feature on the next release (hopefully before the end of March 2022).

HourlyEmisOut = T, which is available on the current open source version of the model, will enable hourly emission output.
If you also want daily emission output, and can not wait for the next release, you need to modify Derived_mod.f90 line 695 as follows:

  do  ind = 1, size(EMIS_FILE)
    dname = "Emis_mgm2_" // trim(EMIS_FILE(ind))
    if(HourlyEmisOut)then
       call AddNewDeriv( dname, "TotEmis", "-", "-", "mg/m2", &
-           ind , -99, T,  1.0e6,  F,  'YMH' )
+           ind , -99, T,  1.0e6,  F,  'YMDH' )
    else
       call AddNewDeriv( dname, "TotEmis", "-", "-", "mg/m2", &
            ind , -99, T,  1.0e6,  F,  'YM' )
    endif
  end do ! ind

Is it possible to output other biogenic emissions but "Emis_mgm2_BioNatC5H8" and "Emis_mgm2_BioNatTERP"?

As far as I can tell all EMIS_BioNat are written out.
If you want daily and hourly EMIS_BioNat output you need to modify Derived_mod.f90 line 683 as follows:

  do  ind = 1, NEMIS_BioNat
    if(EMIS_BioNat(ind)(1:5)=="ASH_L")cycle   ! skip ASH_LxxByy for AshInversion
    dname = "Emis_mgm2_BioNat" // trim(EMIS_BioNat(ind) )
    if(MasterProc) write(*,'(a,i4,a)') dtxt//'NatEmis ', ind, trim(dname)
    call AddNewDeriv( dname, "NatEmis", "-", "-", "mg/m2", &
-                ind , -99, T ,    1.0e6,     F, 'YM' )
+                ind , -99, T ,    1.0e6,     F, 'YMDH' )
  end do

from emep-ctm.

Hi, thank you for your reply. I tried setting "HourlyEmisOut = T", but no hourly emissions were output in the " *hourInst.nc file ".

from emep-ctm.

Hi, thank you for your reply. I tried setting "HourlyEmisOut = T", but no hourly emissions were output in the " *hourInst.nc file ".

Did you also made the changes I wrote on my previous message?

from emep-ctm.

Hi, sorry - I understood your message that way, that I can use "HourlyEmisOut = T" already in the current version without adapting the source code and recompiling. I thought the recompiliation is only needed for Daily output. So I tried without adapting the source code and recompiling.

from emep-ctm.