"Gas"]
additions.loc[RETROFIT_PERIOD, "Gas"] = 0
# Ramp up some BioCCS, quadratically
BIOCCS_TREND = 10 # The increase of annual capacity installed (MW / yr)
bioCCS_2050 = BIOCCS_TREND * len(RETROFIT_PERIOD)
bioCCS_ramp = pd.Series(range(0, bioCCS_2050, BIOCCS_TREND), RETROFIT_PERIOD)
additions.loc[RETROFIT_PERIOD, "BioCCS"] = bioCCS_ramp
# Save a bit on Gas CCS, keeping the END_YEAR generation unchanged
savedGasCCS = (bioCCS_ramp * baseline.capacity_factor.at[END_YEAR, "BioCCS"] /
baseline.capacity_factor.at[END_YEAR, "GasCCS"])
additions.loc[RETROFIT_PERIOD, "GasCCS"] -= savedGasCCS
withCCS = Plan(additions, retirement[technologies], baseline.capacity_factor,
baseline.net_import)
withCCS.__doc__ = "With CCS"
if __name__ == '__main__':
if (len(sys.argv) == 2) and (sys.argv[1] == "summarize"):
withCCS.summarize()
if (len(sys.argv) == 3) and (sys.argv[1] == "plot"):
withCCS.plot_plan(sys.argv[2])
#print(withCCS)
#
#withCCS.summarize()
#print(repr(withCCS))
#withCCS.plot_plan("withCCS.pdf")
The electricity produced per year is the less than in the baseline scenario,
because capacity factor is lower and plant lifetime is shorter for gas than for coal.
"""
import sys
from Plan import Plan
from plan_baseline import (additions as baseline_additions, retirement as
baseline_retirement, CAPACITY_FACTOR, net_import)
# %%
additions = baseline_additions.copy()
additions["Coal"] = (additions["Coal"] + additions["Gas"]) / 2
additions["Gas"] = additions["Coal"]
retirement = baseline_retirement.copy()
retirement["Coal"] = (retirement["Coal"] + retirement["Gas"]) / 2
retirement["Gas"] = retirement["Coal"]
#%% Main statement
moreGas = Plan(additions, retirement, CAPACITY_FACTOR, net_import)
moreGas.__doc__ = "Install as much gas power new capacity as coal power new capacity."
if __name__ == '__main__':
if (len(sys.argv) == 2) and (sys.argv[1] == "summarize"):
moreGas.summarize()
if (len(sys.argv) == 3) and (sys.argv[1] == "plot"):
moreGas.plot_plan(sys.argv[2])
CAPACITY_FACTOR["Wind"] = extend("Renewable", 0.3, "Wind")
CAPACITY_FACTOR["Solar"] = extend("Renewable", 0.23, "Solar")
CAPACITY_FACTOR["CoalCCS"] = CAPACITY_FACTOR["Coal"]
CAPACITY_FACTOR["GasCCS"] = CAPACITY_FACTOR["Gas"]
CAPACITY_FACTOR["BioCCS"] = CAPACITY_FACTOR["Biomass"]
CAPACITY_FACTOR = CAPACITY_FACTOR.where(CAPACITY_FACTOR < 1)
net_import = extend("Import", 7000, "Import", PRODUCTION_PAST,
production_PDP7A)
# %% Main statement
baseline = Plan(additions, retirement, CAPACITY_FACTOR, net_import)
baseline.__doc__ = "Baseline - PDP7A extended"
if __name__ == '__main__':
if (len(sys.argv) == 2) and (sys.argv[1] == "summarize"):
baseline.summarize()
if (len(sys.argv) == 3) and (sys.argv[1] == "plot"):
baseline.plot_plan(sys.argv[2])
# %% Validation: compares to PDP7A
show("""
*****************************************************
*** Comparing our baseline with PDP7 ***
*****************************************************
""")
