Merge pull request #239 from iiasa/ssp-dev-mm

Add scripts to modify "Low" and "Very Low" ScenarioMIP/SSP scenarios

Author: measrainsey

message_ix_models/data/ue-shares/trp_gas0.05.xlsx

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:20766b4dda378f9bc44879ff70b4f5fa9c4d6e0c0e467e157b0e8b00812f44f3
+size 17460

message_ix_models/data/ue-shares/trp_gas0.05_elec0.9.xlsx

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:44fd66b543a316fb3a837d4a70d64262a416074fd230dc85e5ecf4b62614f516
+size 17546

message_ix_models/project/ssp/script/__init__.py

@@ -0,0 +1,96 @@
+import ixmp  # type: ignore
+import message_ix  # type: ignore
+
+from message_ix_models.project.ssp.script.util.functions import (
+    add_balance_equality,
+    add_steel_sector_nze,
+    modify_steel_growth,
+    modify_steel_initial,
+    remove_bof_steel_lower,
+)
+
+# selections
+sel_scen = "SSP1"
+scen_suffix = ""
+rem_bof_steel = True
+mod_growth_steel = True
+mod_initial_steel = True
+add_steel_target = False
+
+# parameters
+trp_year_start = 2035
+mult_price = 5.5
+rc_years = [2060, 2070, 2080, 2090, 2100, 2110]
+steel_years = [2030, 2035, 2040, 2045, 2050, 2055, 2060, 2070, 2080, 2090, 2100, 2110]
+steel_growth = 0.075
+steel_inital = 1.0
+nze_targets = [
+    4.2,
+    3.7,
+    3.2,
+    2.7,
+]
+
+# model and scenario names
+snames = {"SSP1": "SSP1 - Low Emissions", "SSP2": "SSP2 - Low Emissions"}
+svers = {"SSP1": 1, "SSP2": 1}
+model_orig = "SSP_" + sel_scen + "_v1.0"
+scenario_orig = snames[sel_scen]
+
+if rem_bof_steel:
+    scen_suffix += "_bof"
+if mod_growth_steel:
+    scen_suffix += "_growth"
+if mod_initial_steel:
+    scen_suffix += "_initial"
+if add_steel_target:
+    scen_suffix += "_nzsteel"
+
+# target scenario
+model_target = "MM_ScenarioMIP"
+scenario_target = "Low_" + sel_scen + scen_suffix  # + "_v" + str(scen_vers)
+
+# connect to database
+mp = ixmp.Platform("ixmp_dev")
+
+# load scenario
+s_orig = message_ix.Scenario(
+    mp, model=model_orig, scenario=scenario_orig, version=svers[sel_scen]
+)
+
+# clone scenario
+s_tar = s_orig.clone(model_target, scenario_target, keep_solution=False)
+s_tar.set_as_default()
+
+# modify steel sector
+if rem_bof_steel:
+    remove_bof_steel_lower(s_tar, steel_years)
+
+if mod_growth_steel:
+    modify_steel_growth(
+        s_tar,
+        ["dri_gas_steel", "dri_h2_steel", "eaf_steel"],
+        steel_years,
+        steel_growth,
+    )
+
+if mod_initial_steel:
+    modify_steel_initial(
+        s_tar,
+        ["dri_gas_steel", "dri_h2_steel"],
+        steel_years,
+        steel_inital,
+    )
+
+if add_steel_target:
+    add_steel_sector_nze(s_tar, nze_targets)
+
+# add balance equality
+add_balance_equality(s_tar)
+
+solve_typ = "MESSAGE-MACRO"
+solve_args = dict(model=solve_typ)
+s_tar.solve(**solve_args)
+s_tar.set_as_default()
+
+mp.close_db()

message_ix_models/project/ssp/script/scenarios/low.py

@@ -0,0 +1,127 @@
+import os
+
+import ixmp  # type: ignore
+import message_ix  # type: ignore
+
+from message_ix_models.project.ssp.script.util.functions import (
+    add_balance_equality,
+    add_steel_sector_nze,
+    modify_rc_bounds,
+    modify_steel_growth,
+    modify_steel_initial,
+    modify_tax_emission,
+    remove_bof_steel_lower,
+)
+from message_ix_models.project.ssp.script.util.shares import (
+    main as add_UE_share_constraints,
+)
+from message_ix_models.util import package_data_path
+
+# selections
+sel_scen = "LED"
+scen_suffix = "g05_e9"
+rem_bof_steel = True
+mod_growth_steel = True
+mod_initial_steel = True
+add_steel_target = False
+file_ue = "trp_gas0.05_elec0.9.xlsx"
+
+# parameters
+trp_year_start = 2035
+mult_price = 5.5
+rc_years = [2060, 2070, 2080, 2090, 2100, 2110]
+steel_years = [2030, 2035, 2040, 2045, 2050, 2055, 2060, 2070, 2080, 2090, 2100, 2110]
+steel_growth = 0.075
+steel_inital = 1.0
+nze_targets = [
+    0,
+    0,
+    0,
+    0,
+]
+
+# scenario names
+snames = {"SSP1": "SSP1 - Very Low Emissions", "LED": "SSP2 - Very Low Emissions"}
+svers = {"SSP1": 1, "LED": 2}
+model_orig = "SSP_" + sel_scen + "_v1.0"
+scenario_orig = snames[sel_scen]
+
+if rem_bof_steel:
+    scen_suffix += "_bof"
+if mod_growth_steel:
+    scen_suffix += "_growth"
+if mod_initial_steel:
+    scen_suffix += "_initial"
+if add_steel_target:
+    scen_suffix += "_nzsteel"
+
+# read UE share file
+path_ue = package_data_path("ue-shares")
+path_ue_file = os.path.join(path_ue, file_ue)
+
+# target scenario
+model_target = "MM_ScenarioMIP"
+scenario_target = "VL_" + sel_scen + "_" + scen_suffix
+
+# connect to database
+mp = ixmp.Platform("ixmp_dev")
+
+# load scenario
+s_orig = message_ix.Scenario(
+    mp, model=model_orig, scenario=scenario_orig, version=svers[sel_scen]
+)
+
+# clone scenario
+s_tar = s_orig.clone(model_target, scenario_target, keep_solution=False)
+s_tar.set_as_default()
+
+# modify bounds for some fuels in residential and commercial sector
+modify_rc_bounds(s_orig, s_tar, rc_years)
+
+# add UE share constraints to transport
+add_UE_share_constraints(
+    s_tar,  # scenario object
+    path_UE_share_input=path_ue_file,  # path
+    ssp=sel_scen,  # SSP-name i.e "LED"
+    start_year=trp_year_start,  # the year as of which a constraint should be added
+    calibration_year=2020,  # 2020
+    clean_relations=False,  # set to False
+    verbose=True,  # set to True
+)
+
+# modify tax/price emissions
+modify_tax_emission(s_orig, s_tar, mult_price)
+
+# modify steel sector
+if rem_bof_steel:
+    remove_bof_steel_lower(s_tar, steel_years)
+
+if mod_growth_steel:
+    modify_steel_growth(
+        s_tar,
+        ["dri_gas_steel", "dri_h2_steel", "eaf_steel"],
+        steel_years,
+        steel_growth,
+    )
+
+if mod_initial_steel:
+    modify_steel_initial(
+        s_tar,
+        ["dri_gas_steel", "dri_h2_steel"],
+        steel_years,
+        steel_inital,
+    )
+
+if add_steel_target:
+    add_steel_sector_nze(s_tar, nze_targets)
+
+# add balance equality
+add_balance_equality(s_tar)
+
+# solve
+solve_typ = "MESSAGE-MACRO"
+solve_args = dict(model=solve_typ)
+s_tar.solve(**solve_args)
+s_tar.set_as_default()
+
+mp.close_db()