Dac tutorial message models

message_ix/util/tutorial.py

@@ -3,17 +3,19 @@
 from functools import partial
 
 from message_ix import Scenario
-from message_ix.report import Key, Reporter, operator
+from message_ix.report import Key, Reporter, computations
 
 log = logging.getLogger(__name__)
 
 PLOTS = [
-    ("activity", operator.stacked_bar, "out:nl-t-ya", "GWa"),
-    ("capacity", operator.stacked_bar, "CAP:nl-t-ya", "GW"),
-    ("demand", operator.stacked_bar, "demand:n-c-y", "GWa"),
-    ("extraction", operator.stacked_bar, "EXT:n-c-g-y", "GW"),
-    ("new capacity", operator.stacked_bar, "CAP_NEW:nl-t-yv", "GWa"),
-    ("prices", operator.stacked_bar, "PRICE_COMMODITY:n-c-y", "¢/kW·h"),
+    ("activity", computations.stacked_bar, "out:nl-t-ya", "GWa"),
+    ("capacity", computations.stacked_bar, "CAP:nl-t-ya", "GW"),
+    ("removal capacity", computations.stacked_bar, "CAP:nl-t-ya", "tCO2/yr"),
+    ("demand", computations.stacked_bar, "demand:n-c-y", "GWa"),
+    ("emission", computations.stacked_bar, "emi:nl-t-ya", "tCO2"),
+    ("extraction", computations.stacked_bar, "EXT:n-c-g-y", "GW"),
+    ("new capacity", computations.stacked_bar, "CAP_NEW:nl-t-yv", "GWa"),
+    ("prices", computations.stacked_bar, "PRICE_COMMODITY:n-c-y", "¢/kW·h"),
 ]
 
 
@@ -27,6 +29,8 @@ def prepare_plots(rep: Reporter, input_costs="$/GWa") -> None:
     - ``plot extraction``
     - ``plot fossil supply curve``
     - ``plot capacity``
+    - ``plot removal capacity``
+    - ``plot emission``
     - ``plot new capacity``
     - ``plot prices``
 
@@ -47,7 +51,7 @@ def prepare_plots(rep: Reporter, input_costs="$/GWa") -> None:
     # Add one node to the reporter for each plot
     for title, func, key_str, units in PLOTS:
         # Convert the string to a Key object so as to reference its .dims
-        key = Key(key_str)
+        key = Key.from_str_or_key(key_str)
 
         # Operation for the reporter
         comp = partial(
@@ -68,7 +72,7 @@ def prepare_plots(rep: Reporter, input_costs="$/GWa") -> None:
         "plot fossil supply curve",
         (
             partial(
-                operator.plot_cumulative,
+                computations.plot_cumulative,
                 labels=("Fossil supply", "Resource volume", "Cost"),
             ),
             "resource_volume:n-g",

tutorial/README.rst

@@ -158,6 +158,9 @@ framework, such as used in global research applications of |MESSAGEix|.
    #. Modeling of a multi-node energy system and representing trade between nodes
       (:tut:`westeros/westeros_multinode_energy_trade.ipynb`).
 
+   #. Including carbon removal technologies
+      (:tut:`westeros/westeros_carbon_removal.ipynb`).
+
 #. Use other features of :mod:`message_ix` and :mod:`ixmp`:
 
    #. ⭐ After the MESSAGE model has solved, use the :mod:`.message_ix.report`

tutorial/westeros/data/westeros_carbon_removal_data.yaml

@@ -0,0 +1,98 @@
+daccs:
+  year_init: 700
+  inv_cost_: 
+    par_name: inv_cost
+    value: 100
+    unit: USD/(tCO2/yr)
+    node_loc:
+      Westeros: 1
+    year_vtg:
+      rate: 0
+  fix_cost_: 
+    par_name: fix_cost
+    value: 5
+    unit: USD/(tCO2/yr)/yr
+    node_loc:
+      Westeros: 1
+    year_vtg:
+      rate: 0
+    year_act:
+      rate: 0
+  var_cost_: 
+    par_name: var_cost
+    value: 5
+    unit: USD/tCO2
+    node_loc:
+      Westeros: 1
+    year_vtg:
+      rate: 0
+    year_act:
+      rate: 0
+  input_:
+    par_name: input
+    value: 0.0028
+    unit: '-'
+    node_loc:
+      Westeros: 1
+    mode:
+      standard: 1
+    commodity:
+      electricity: 1
+    level:
+      final: 1
+  output_:
+    par_name: output
+    value: 1
+    unit: tCO2
+    node_loc:
+      Westeros: 1
+    mode:
+      standard: 1
+    commodity:
+      CO2: 1
+    level:
+      final: 1
+  capacity_factor_:
+    par_name: capacity_factor
+    value: 0.913
+    unit: '-'
+    node_loc:
+      Westeros: 1
+  emission_factor_:
+    par_name: emission_factor
+    value: -1
+    unit: tCO2/tCO2
+    node_loc:
+      Westeros: 1
+    mode:
+      standard: 1
+    emission:
+      CO2: 1
+  technical_lifetime_:
+    par_name: technical_lifetime
+    value: 25
+    unit: y
+    node_loc:
+      Westeros: 1
+  initial_new_capacity_up_: 
+    par_name: initial_new_capacity_up
+    value: 0.5
+    unit: Mt CO2/yr
+    node_loc:
+      Westeros: 1
+    year_vtg:
+      rate: 0
+    time:
+      year: 1
+  growth_new_capacity_up_: 
+    par_name: growth_new_capacity_up
+    value: 0.05
+    unit: '-'
+    node_loc:
+      Westeros: 1
+    year_vtg:
+      rate: 0
+    time:
+      year: 1
+
+