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)

path <- system.file("extdata/example1-4pools.xlsx", package = "mocaredd")

cs <- read_xlsx(path = path, sheet = "c_stocks", na = "NA")
ad <- read_xlsx(path = path, sheet = "AD_lu_transitions", na = "NA")
usr <- read_xlsx(path = path, sheet = "user_inputs", na = "NA")
time <- read_xlsx(path = path, 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)
#> Error in pmap(combi, function(period, lu) {    c_sub <- filter(.cs, .data$c_period == period, .data$c_lu_id ==         lu)    c_check <- fct_check_pool(.c_sub = c_sub, .c_unit = .usr$c_unit,         .c_fraction = .usr$c_fraction)}): could not find function "pmap"

sim_FREL <- fct_combine_mcs_P(
  .data = sim_trans,
  .time = time_clean,
  .period_type = "REF",
  .ad_annual = usr$ad_annual
)
#> Error in purrr::map(unique(time_ref$period_type), function(x) {    dplyr::summarise(dplyr::group_by(dplyr::filter(dplyr::left_join(.data,         time_ref, by = c(time_period = "period_no")), .data$period_type ==         x), .data$sim_no, .data$period_type), E_sim = round(sum(.data$E_sim *         .data$nb_years)/length_ref, 0), .groups = "drop")}): ℹ In index: 1.
#> Caused by error:
#> ! object 'sim_trans' not found

hist(sim_FREL$E_sim)
#> Error: object 'sim_FREL' not found
round(median(sim_FREL$E_sim))
#> Error: object 'sim_FREL' not found