endogenizing steel and cement production and add existing industry capacities

Author: toniseibold

config/config.default.yaml

@@ -29,19 +29,24 @@ run:
   use_shadow_directory: false
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#foresight
-foresight: overnight
+foresight: myopic
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#scenario
 # Wildcard docs in https://pypsa-eur.readthedocs.io/en/latest/wildcards.html
 scenario:
   clusters:
-  - 50
+  - 39
+  # - 128
+  # - 256
   opts:
   - ''
   sector_opts:
   - ''
   planning_horizons:
-  - 2050
+  # - 2020
+  - 2030
+  - 2040
+  # - 2050
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#countries
 countries:
@@ -728,7 +733,7 @@ sector:
   - nearshore      # within 50 km of sea
     # - offshore
   methanol:
-    regional_methanol_demand: false
+    regional_methanol_demand: true
     methanol_reforming: false
     methanol_reforming_cc: false
     methanol_to_kerosene: false
@@ -739,7 +744,7 @@ sector:
       allam: false
     biomass_to_methanol: true
     biomass_to_methanol_cc: false
-  ammonia: true
+  ammonia: regional
   min_part_load_electrolysis: 0
   min_part_load_fischer_tropsch: 0.5
   min_part_load_methanolisation: 0.3
@@ -827,6 +832,12 @@ sector:
       methanol: 121
       gas: 122
       oil: 125
+  endogenous_sectors:
+    - steel
+    - cement
+    - ammonia
+    - methanol
+  industry_relocation: false
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#industry
 industry:
@@ -969,7 +980,7 @@ clustering:
       ramp_limit_down: max
   temporal:
     resolution_elec: false
-    resolution_sector: 8760H
+    resolution_sector: 365H
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#adjustments
 adjustments:

config/plotting.default.yaml

@@ -586,3 +586,21 @@ plotting:
     import NH3: '#e2ed74'
     import oil: '#93eda2'
     import methanol: '#87d0e6'
+    EAF: '#586357'
+    H2 DRI: '#618ab0'
+    gas DRI: '#baaf68'
+    BOF: '#1c1f1b'
+    steel: '#705238'
+    steel emission: '#756252'
+    steel emission CC: '#756252'
+    cement: '#c98349'
+    cement emission CC: '#c98349'
+    hbi: '#87d0e6'
+    DRI: '#87d0e6'
+    grey methanol: '#6b7887'
+    blue methanol: '#496687'
+    cement emission: '#c98349'
+    cement kiln: '#c98349'
+    cement production: '#c98349'
+    cement finishing: '#c98349'
+    clinker: '#c98349'

rules/build_sector.smk

@@ -1558,6 +1558,9 @@ rule prepare_sector_network:
             if config_provider("sector", "district_heating", "ates", "enable")(w)
             else []
         ),
+        industry_sector_ratios=resources(
+            "industry_sector_ratios_{planning_horizons}.csv"
+        ),
     output:
         resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc"

rules/solve_myopic.smk

@@ -14,6 +14,7 @@ rule add_existing_baseyear:
         energy_totals_year=config_provider("energy", "energy_totals_year"),
         countries=config_provider("countries"),
         MWh_NH3_per_tNH3=config_provider("industry", "MWh_NH3_per_tNH3"),
+        MWh_MeOH_per_tMeOH=config_provider("industry", "MWh_MeOH_per_tMeOH"),
     input:
         network=resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}.nc"
@@ -33,6 +34,7 @@ rule add_existing_baseyear:
         regions_onshore=resources("regions_onshore_base_s_{clusters}.geojson"),
         ammonia="data/ammonia_plants.csv",
         isi_database="data/1-s2.0-S0196890424010586-mmc2.xlsx",
+        gem_gspt="data/Global-Cement-and-Concrete-Tracker_July-2025.xlsx",
     output:
         resources(
             "networks/base_s_{clusters}_{opts}_{sector_opts}_{planning_horizons}_brownfield.nc"