[RegCNET] time disaggregation of yearly/monthly emissions

[Zakey Ashraf_Saber]

Hello Peter,

    In order to preform a time aggregation to yearly or monthly emission
data, you should processes the following part of the code:



MG_YEAR_TO_G_HOUR = 1.0D6/365.0/24.0 !1.0D6 - ktones

  DO IT = 1, TIME_LEN  ! TIME_LEN = 24
   DO L = 1, NLAT
    DO K = 1, NLON
     DO J = 1, NRCON
      DO I = 1, NREMIS
         TIME_DEP_FAC = MONTH_VAR(J,C_MM(IT)) *                      &
                      & WEEK_VAR(J,C_WEEK(IT)) * DAY_VAR(J,C_HH(IT))
         EMIS(IT,I,J,K,L) = FLOAT(SNAP(J))*                          &
                          & TIME_DEP_FAC* EMIS(I,J,K,L)*          &
                          & MG_YEAR_TO_G_HOUR
      ENDDO       !DO I = 1, NREMIS
     ENDDO       !DO J = 1, NRCON
    ENDDO       !DO K = 1, NLON
   ENDDO       !DO L = 1, NLAT
  ENDDO       !DO IT = 0, TIME_LEN+1


The SNAP code should be something like:


Snap codes:
  SNAP     =  (/   &
               1   & ! 1.  Power generation
              ,1   & ! 2.  Residential, commercial and other combustion
              ,1   & ! 3.  Industrial combustion
              ,1   & ! 4.  Industrial processes
              ,1   & ! 5.  Extraction distribution of fossil fuels
              ,1   & ! 6.  Solvent use
              ,1   & ! 7.  Road transport
              ,1   & ! 8.  Other mobile sources
              ,1   & ! 9.  Waste treatment and disposal
              ,1   & ! 10. Agriculture
              ,1  /) ! 11. Other sources and sinks

I hope this will work with you

If you need some further discussion contact me.


Best wishes
Ashraf


> Dear RegCM users,
>
> does RegCM perform time disaggregation (to create e.g. hourly
> emissions) of yearly/monthly emission data supplied in AERO.dat and if
> so, what monthly/daily and hourly disaggregation factors he uses?
>
> Thanks for the answer
>
> Peter
> Charles University, Prague
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>