Combine MCS of emissions to a defined time period

Depending on how the period is defined and if the data are annualized or not, calculate the Emission Level for a reference or monitoring period for each simulation.

Usage

fct_combine_mcs_P(.data, .time, .period_type, .ad_annual)

Arguments

.data

a data frame containing the simulations

.time

the 'time' table from the tool input file (see template)

.period_type

"reference" or "monitoring"

.ad_annual

TRUE or FALSE, is the activity data annualized or not.

Value

A tibble with simulations at the final estimate per type of period.

Examples

library(readxl)
library(dplyr)
library(mocaredd)

cs <- read_xlsx(
  path = system.file("extdata/example1.xlsx", package = "mocaredd"),
  sheet = "c_stocks",
  na = "NA"
  )
ad <- read_xlsx(
  path = system.file("extdata/example1.xlsx", package = "mocaredd"),
  sheet = "AD_lu_transitions",
  na = "NA"
  )
usr <- read_xlsx(
  path = system.file("extdata/example1.xlsx", package = "mocaredd"),
  sheet = "user_inputs",
  na = "NA"
  )
time <- read_xlsx(
  path = system.file("extdata/example1.xlsx", package = "mocaredd"),
  sheet = "time_periods",
  na = "NA"
  )

ad_clean <- ad |> dplyr::filter(!is.na(trans_area) | !is.na(trans_pdf_a))
cs_clean <- cs |> dplyr::filter(!is.na(c_value) | !is.na(c_pdf_a))
time_clean <- time |> dplyr::mutate(nb_years = year_end - year_start + 1)

sim_trans <- fct_combine_mcs_E(.ad = ad_clean, .cs = cs_clean, .usr = usr)
#> Random simulations with seed: 93

sim_FREL <- fct_combine_mcs_P(
  .data = sim_trans,
  .time = time_clean,
  .period_type = "REF",
  .ad_annual = usr$ad_annual
)

hist(sim_FREL$E_sim)

round(median(sim_FREL$E_sim))
#> [1] 4591614