From aa820fa0464480ec5e11c1782d2ce1791ad9f370 Mon Sep 17 00:00:00 2001 From: lisazeyen Date: Mon, 29 Jul 2024 16:00:37 +0200 Subject: [PATCH] change industry reference year to 2019 --- config/config.default.yaml | 2 +- scripts/build_industrial_production_per_country.py | 8 ++++++-- 2 files changed, 7 insertions(+), 3 deletions(-) diff --git a/config/config.default.yaml b/config/config.default.yaml index 373f4a6c..521013b0 100644 --- a/config/config.default.yaml +++ b/config/config.default.yaml @@ -716,7 +716,7 @@ industry: MWh_CH4_per_tMeOH: 10.25 MWh_MeOH_per_tMeOH: 5.528 hotmaps_locate_missing: false - reference_year: 2021 + reference_year: 2019 # docs in https://pypsa-eur.readthedocs.io/en/latest/configuration.html#costs diff --git a/scripts/build_industrial_production_per_country.py b/scripts/build_industrial_production_per_country.py index 1bb75271..332627dd 100644 --- a/scripts/build_industrial_production_per_country.py +++ b/scripts/build_industrial_production_per_country.py @@ -293,6 +293,7 @@ def separate_basic_chemicals(demand, year): """ Separate basic chemicals into ammonia, chlorine, methanol and HVC. """ + # ammonia data from 2017-2021 ammonia = pd.read_csv(snakemake.input.ammonia_production, index_col=0) there = ammonia.index.intersection(demand.index) @@ -301,8 +302,11 @@ def separate_basic_chemicals(demand, year): logger.info(f"Following countries have no ammonia demand: {missing.tolist()}") demand["Ammonia"] = 0.0 - - demand.loc[there, "Ammonia"] = ammonia.loc[there, str(year)] + + year_to_use = min(max(year, 2017), 2021) + if year_to_use != year: + logger.info(f"Using data from {year_to_use} for ammonia production.") + demand.loc[there, "Ammonia"] = ammonia.loc[there, str(year_to_use)] demand["Basic chemicals"] -= demand["Ammonia"]