add coal and coke demand for integrated steelworks (closes #717)

scripts/prepare_sector_network.py

@@ -2883,6 +2883,25 @@ def add_industry(n, costs):
             p_set=p_set,
         )
 
+    primary_steel = snakemake.config["industry"]["St_primary_fraction"]
+    dri_steel = snakemake.config["industry"]["DRI_fraction"]
+    bof_steel = primary_steel - dri_steel
+
+    if bof_steel > 0:
+
+        add_carrier_buses(n, "coal")
+
+        mwh_coal_per_mwh_coke = 1.366 # from eurostat energy balance
+        p_set = (industrial_demand["coal"].sum() + mwh_coal_per_mwh_coke * industrial_demand["coke"].sum()) / nhours
+
+        n.madd(
+            "Load",
+            "coal for industry",
+            bus=spatial.coal.nodes,
+            carrier="coal for industry",
+            p_set=p_set,
+        )
+
 
 def add_waste_heat(n):
     # TODO options?