capacity_past.Oil = smooth(capacity_past.Oil)
capacity_past["Solar"] = 0
capacity_past["Import"] = 0
capacity_past["PumpedStorage"] = 0
capacity_past["CoalCCS"] = 0
capacity_past["GasCCS"] = 0
capacity_past["BioCCS"] = 0
capacity_past["BigHydro"] = (capacity_past.LargeHydro +
capacity_past.InterHydro)
capacity_past["Hydro"] = (capacity_past.BigHydro + capacity_past.SmallHydro)
addcol_Renewable(capacity_past)
addcol_Renewable4(capacity_past)
show("""
Vietnam historical capacity additions by fuel type (MW)
Source: Capacities listed in EVN activity report 2016, dated by internet search
""")
show(capacity_past[fuels])
show()
show("""
Vietnam historical generation capacity by fuel type (MW)
Source: Capacities listed in EVN activity report 2016, dated by internet search
""")
show(capacity_past[fuels].cumsum())
show()
CAPACITY_PAST.Oil = smooth(CAPACITY_PAST.Oil)
CAPACITY_PAST["Solar"] = 0
CAPACITY_PAST["Import"] = 0
CAPACITY_PAST["PumpedStorage"] = 0
CAPACITY_PAST["CoalCCS"] = 0
CAPACITY_PAST["GasCCS"] = 0
CAPACITY_PAST["BioCCS"] = 0
CAPACITY_PAST["BigHydro"] = (CAPACITY_PAST.LargeHydro +
CAPACITY_PAST.InterHydro)
CAPACITY_PAST["Hydro"] = (CAPACITY_PAST.BigHydro + CAPACITY_PAST.SmallHydro)
addcol_Renewable(CAPACITY_PAST)
addcol_Renewable4(CAPACITY_PAST)
show("""
Vietnam historical capacity additions by fuel type (MW)
Source: Capacities listed in EVN activity report 2016, dated by internet search
""")
show(CAPACITY_PAST[PLANT_TYPE])
show()
show("""
Vietnam historical generation capacity by fuel type (MW)
Source: Capacities listed in EVN activity report 2016, dated by internet search
""")
show(CAPACITY_PAST[PLANT_TYPE].cumsum())
show()
alternative.plot_plan(sys.argv[2])
#%% Validation: compares to PDP7A
show("""
*****************************************************
*** Comparing our baseline with PDP7 ***
*****************************************************
""")
compared = ["Coal", "Gas", "BigHydro", "Renewable4", "Nuclear"]
tocompare = alternative.capacities.loc[[2020, 2025, 2030]]
tocompare["Gas"] = tocompare.Gas + tocompare.Oil
tocompare["Nuclear"] = 0
addcol_Renewable4(tocompare)
tocompare = tocompare[compared]
relerror = (tocompare - capacities_PDP7A) / capacities_PDP7A
relerror = relerror[compared]
show("PDP7A")
show(capacities_PDP7A[compared])
show()
show("Baseline scenario, Gas includes Oil")
show(tocompare)
show()
show("Relative error")
show(relerror)
show(
"Note: Gas 2030 is larger in baseline because we replace nuclear with gas")