Industry workflow

Snakefile

@@ -315,6 +315,7 @@ rule build_population_layouts:
         pop_layout_total="resources/population_shares/pop_layout_total.nc",
         pop_layout_urban="resources/population_shares/pop_layout_urban.nc",
         pop_layout_rural="resources/population_shares/pop_layout_rural.nc",
+        gdp_layout="resources/gdp_shares/gdp_layout.nc",
     resources:
         mem_mb=20000,
     benchmark:
@@ -340,12 +341,14 @@ rule build_clustered_population_layouts:
         pop_layout_total="resources/population_shares/pop_layout_total.nc",
         pop_layout_urban="resources/population_shares/pop_layout_urban.nc",
         pop_layout_rural="resources/population_shares/pop_layout_rural.nc",
+        gdp_layout="resources/gdp_shares/gdp_layout.nc",
         regions_onshore=pypsaearth(
             "resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson"
         ),
         cutout=pypsaearth(CUTOUTS_PATH),
     output:
         clustered_pop_layout="resources/population_shares/pop_layout_elec_s{simpl}_{clusters}.csv",
+        clustered_gdp_layout="resources/gdp_shares/gdp_layout_elec_s{simpl}_{clusters}.csv",
     resources:
         mem_mb=10000,
     benchmark:
@@ -586,9 +589,12 @@ if config["custom_data"].get("industry_demand", False) == True:
                 "resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson"
             ),
             clustered_pop_layout="resources/population_shares/pop_layout_elec_s{simpl}_{clusters}.csv",
+            clustered_gdp_layout="resources/gdp_shares/gdp_layout_elec_s{simpl}_{clusters}.csv",
             industrial_database="resources/custom_data/industrial_database.csv",
-            #shapes_path=pypsaearth("resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson")
-            shapes_path=PYPSAEARTH_FOLDER + "/resources/shapes/MAR2.geojson",
+            shapes_path=pypsaearth(
+                "resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson"
+            ),
+            #shapes_path=PYPSAEARTH_FOLDER + "/resources/shapes/MAR2.geojson",
         output:
             industrial_distribution_key="resources/demand/industrial_distribution_key_elec_s{simpl}_{clusters}.csv",
         threads: 1
@@ -625,6 +631,7 @@ if config["custom_data"].get("industry_demand", False) == False:
                 "resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson"
             ),
             clustered_pop_layout="resources/population_shares/pop_layout_elec_s{simpl}_{clusters}.csv",
+            clustered_gdp_layout="resources/gdp_shares/gdp_layout_elec_s{simpl}_{clusters}.csv",
             industrial_database="data/industrial_database.csv",
             shapes_path=pypsaearth(
                 "resources/bus_regions/regions_onshore_elec_s{simpl}_{clusters}.geojson"
@@ -639,25 +646,29 @@ if config["custom_data"].get("industry_demand", False) == False:
         script:
             "scripts/build_industrial_distribution_key.py"
 
-    rule build_industrial_production_per_country_tomorrow:  #default data
+    rule build_base_industry_totals:  #default data
         input:
             industrial_production_per_country="data/industrial_production_per_country.csv",
+            #unsd_path="data/demand/unsd/data/",
+            energy_totals_base="data/energy_totals_base.csv",
         output:
-            industrial_production_per_country_tomorrow="resources/demand/industrial_production_per_country_tomorrow_{planning_horizons}_{demand}.csv",
+            base_industry_totals="resources/demand/base_industry_totals_{planning_horizons}_{demand}.csv",
         threads: 1
         resources:
             mem_mb=1000,
         benchmark:
-            "benchmarks/build_industrial_production_per_country_tomorrow_{planning_horizons}_{demand}"
+            "benchmarks/build_base_industry_totals_{planning_horizons}_{demand}"
         script:
-            "scripts/build_industrial_production_tomorrow.py"
+            "scripts/build_base_industry_totals.py"
 
     rule build_industry_demand:  #default data
         input:
             industry_sector_ratios="data/industry_sector_ratios.csv",
             industrial_distribution_key="resources/demand/industrial_distribution_key_elec_s{simpl}_{clusters}.csv",
-            industrial_production_per_country_tomorrow="resources/demand/industrial_production_per_country_tomorrow_{planning_horizons}_{demand}.csv",
-            industrial_production_per_country="data/industrial_production_per_country.csv",
+            #industrial_production_per_country_tomorrow="resources/demand/industrial_production_per_country_tomorrow_{planning_horizons}_{demand}.csv",
+            #industrial_production_per_country="data/industrial_production_per_country.csv",
+            base_industry_totals="resources/demand/base_industry_totals_{planning_horizons}_{demand}.csv",
+            industrial_database="data/industrial_database.csv",
             costs=CDIR
             + "costs_{}.csv".format(config["scenario"]["planning_horizons"][0]),
         output:

config.default.yaml

@@ -43,6 +43,8 @@ demand_data:
   update_data: true # if true, the workflow downloads the energy balances data saved in data/demand/unsd/data again. Turn on for the first run.
   base_year: 2019
 
+  other_industries: false # Whether or not to include industries that are not specified. some countries have has exageratted numbers, check carefully.
+
 H2_network: false
 H2_network_limit: 2000 #GWkm
 H2_repurposed_network: false

data/AL_production.csv

@@ -0,0 +1,34 @@
+country,production[ktons/a],Year,source
+CN,36000,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+IN,3700,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+RU,3600,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+CA,2900,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+AE,2700,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+AU,1580,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+BH,1400,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+NO,1300,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+US,1100,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+SA,932,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+IS,850,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+MY,760,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+ZA,714,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+BR,659,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+QA,616,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+MZ,571,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+DE,529,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+AR,438,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+FR,380,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+OM,380,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+IR,370,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+ES,350,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+NZ,340,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+EG,317,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+RO,283,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+KZ,258,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+ID,242,2018,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+SK,212,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+GR,184,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+SE,126,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+VE,108,2019,https://en.wikipedia.org/wiki/List_of_countries_by_aluminium_production
+NG,0,2019
+BJ,0,2019

data/demand/industry_growth_cagr.csv

@@ -0,0 +1,10 @@
+country,chemical and petrochemical,construction,food and tobacco,iron and steel,machinery,mining and quarrying,non-ferrous metals,non-metallic minerals,paper pulp and print,textile and leather,transport equipment,wood and wood products,other
+MA,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+NG,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+BJ,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+NA,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+BR,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+AE,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+EG,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+CL,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03
+BH,0.03,0.015,0.02,0.017,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03

data/export_ports.csv

@@ -9,3 +9,6 @@ dummy port 1,BJ,1,6.47,2.63
 dummy port 1,BR,1,-22.9068,-43.1729
 dummy port 1,NA,1,-21.05431,13.50664
 dummy port 1,AE,1,25.2048,55.2708
+AS Suways,EG, 1, 29.966667,32.55
+Damietta,EG,1,31.483333,31.75
+Port Said,EG,3,31.266667,32.3