data-raw/costs/capacityfactors.R

In xavier-gutierrez/energy.markets:

calc.capacityfactors <- function(generators, plantgeneration)
{
  plantcapacity <- generators %>%
    select(yr, nameplate, utilcode, plntcode, overnightcategory, fuel.general) %>%
    rename(capacity = nameplate) %>%
    group_by(yr, utilcode, plntcode, overnightcategory, fuel.general) %>%
    summarise(capacity=sum(capacity)) %>%
    ungroup()

  plantgeneration <- plantgeneration %>%
    select(yr, generation, utilcode, plntcode, overnightcategory, fuel.general) %>%
    filter(generation > 0) # drop generators that used more energy than they produced

  cf <- plantcapacity %>%
    inner_join(plantgeneration, by=c("yr", "utilcode", "plntcode", "overnightcategory", "fuel.general") ) %>%
    mutate(potentialgeneration = capacity * 8760) %>%
    mutate(capacityfactor = generation/potentialgeneration) %>%
    select(yr, utilcode, plntcode, overnightcategory, fuel.general, capacityfactor)

  cf.clamp <- cf %>%
    mutate(capacityfactor = ifelse(capacityfactor > 1, 1, capacityfactor))
  # CF > 1 is data error. See figs/filterbyCF for analysis how filter(CF < 1) affected data

  data <- list(cf, cf.clamp)
  names(data) <- c("cf", "cf.clamp")
  # cf used in later calculations

  data
}