bjoern/pypsa-eur
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

include review fneum build_industrial_production_per_country

Browse Source
This commit is contained in:
lisazeyen 2024-08-07 10:33:47 +02:00
parent 8794f53050
commit 6b3fc18386
1 changed files with 3 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/build_industrial_production_per_country.py
View File

@ -202,13 +202,12 @@ def get_energy_ratio(country, eurostat_dir, jrc_dir, year):
else: else:
ct_eurostat = country.replace("GB", "UK") ct_eurostat = country.replace("GB", "UK")
if ct_eurostat == "UK": if ct_eurostat == "UK":
year = 2019
logger.info("Assume Eurostat data for GB from 2019.") logger.info("Assume Eurostat data for GB from 2019.")
# estimate physical output, energy consumption in the sector and country # estimate physical output, energy consumption in the sector and country
fn = f"{eurostat_dir}/{ct_eurostat}-Energy-balance-sheets-April-2023-edition.xlsb" fn = f"{eurostat_dir}/{ct_eurostat}-Energy-balance-sheets-April-2023-edition.xlsb"
df = pd.read_excel( df = pd.read_excel(
fn, fn,
sheet_name=str(min(2021, year)), sheet_name=str(min(2019, year)),
index_col=2, index_col=2,
header=0, header=0,
skiprows=4, skiprows=4,
@ -292,7 +291,7 @@ def separate_basic_chemicals(demand, year):
""" """
Separate basic chemicals into ammonia, chlorine, methanol and HVC. Separate basic chemicals into ammonia, chlorine, methanol and HVC.
""" """
# ammonia data from 2017-2021 # ammonia data from 2018-2022
ammonia = pd.read_csv(snakemake.input.ammonia_production, index_col=0) ammonia = pd.read_csv(snakemake.input.ammonia_production, index_col=0)
there = ammonia.index.intersection(demand.index) there = ammonia.index.intersection(demand.index)
@ -304,7 +303,7 @@ def separate_basic_chemicals(demand, year):
year_to_use = min(max(year, 2018), 2022) year_to_use = min(max(year, 2018), 2022)
if year_to_use != year: if year_to_use != year:
logger.info(f"Using data from {year_to_use} for ammonia production.") logger.info(f"Year {year} outside data range. Using data from {year_to_use} for ammonia production.")
demand.loc[there, "Ammonia"] = ammonia.loc[there, str(year_to_use)] demand.loc[there, "Ammonia"] = ammonia.loc[there, str(year_to_use)]
demand["Basic chemicals"] -= demand["Ammonia"] demand["Basic chemicals"] -= demand["Ammonia"]