existing heating: adjust 2020 stock to ARIADNE database totals

Also adjust space heating reduction and DH shares to better values.

Author: nworbmot

config/config.yaml

@@ -4,7 +4,7 @@
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#run
 run:
-  name: "240112-365H-forceboilersconstraint2"
+  name: "240112-365H-adaptexistingheat"
   scenarios:
     enable: true
   shared_resources: true #stops recalculating
@@ -122,7 +122,7 @@ h2_import_max:
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#sector
 sector:
   district_heating:
-    potential: 0.6
+    potential: 0.4
     progress:
       2020: 0.0
       2030: 0.3
@@ -146,7 +146,16 @@ industry:
     2040: 0.4
     2045: 0.35
     2050: 0.3
-  #TBD what to include in config
+
+  # this needs to be taken from ariadne database
+  reduce_space_heat_exogenously_factor:
+    2020: 0.0
+    2025: 0.05
+    2030: 0.10
+    2035: 0.16
+    2040: 0.22
+    2045: 0.25
+    2050: 0.35
 
 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#solving
 solving:

workflow/Snakefile

@@ -55,7 +55,7 @@ rule retrieve_ariadne_scenario_data:
     params:
         iiasa=config["iiasa_database"]
     output:
-        data="resources/ariadne_scenario_data.csv"
+        data=RESOURCES + "ariadne_scenario_data.csv"
     log:
         "logs/retrieve_ariadne_scenario_data.log"
     resources:
@@ -87,4 +87,23 @@ use rule solve_sector_network_myopic from pypsaeur with:
         custom_extra_functionality=workflow.source_path("scripts/additional_functionality.py"),
     input:
         **rules.solve_sector_network_myopic.input,
-        co2_totals_name=RESOURCES + "co2_totals.csv",
+        co2_totals_name=RESOURCES + "co2_totals.csv",
+
+
+rule modify_existing_heating:
+    input:
+        ariadne=RESOURCES + "ariadne_scenario_data.csv",
+        existing_heating="data/existing_infrastructure/existing_heating_raw.csv",
+    output:
+        existing_heating=RESOURCES + "existing_heating.csv",
+    resources:
+        mem_mb=1000
+    script:
+        "scripts/modify_existing_heating.py"
+
+
+
+use rule add_existing_baseyear from pypsaeur with:
+    input:
+        **{k: v for k, v in rules.add_existing_baseyear.input.items() if k != "existing_heating"},
+	existing_heating=RESOURCES + "existing_heating.csv",

workflow/scripts/modify_existing_heating.py

@@ -0,0 +1,31 @@
+
+
+import pandas as pd
+
+
+existing_heating = pd.read_csv(snakemake.input.existing_heating,
+                               index_col=0)
+
+ariadne = pd.read_csv(snakemake.input.ariadne,
+                      index_col=0)
+
+print("before",existing_heating.loc["Germany"])
+
+mapping = {"gas boiler" : "Gas Boiler",
+                  "oil boiler" : "Oil Boiler",
+                  "air heat pump" : "Heat Pump|Electrical|Air",
+                  "ground heat pump" : "Heat Pump|Electrical|Ground",
+                  "biomass boiler" : "Biomass Boiler"}
+
+year = "2020"
+for tech in mapping:
+    stock = ariadne.at[f"Stock|Space Heating|{mapping[tech]}",year]
+    peak = stock*existing_heating.loc["Germany"].sum()/ariadne.at[f"Stock|Space Heating",year]
+    existing_heating.at["Germany",tech] = peak
+
+
+
+print("before",existing_heating.loc["Germany"])
+
+
+existing_heating.to_csv(snakemake.output.existing_heating)