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

idees eurostat rescaling also for period 1990-1999

Browse Source
This commit is contained in:
Fabian Neumann 2024-03-13 13:40:23 +01:00
parent a429b48874
commit 2eb159af7c
1 changed files with 42 additions and 34 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
scripts/build_energy_totals.py
View File

@ -856,20 +856,25 @@ def rescale_idees_from_eurostat(
]
for country in idees_countries:
slicer_target = idx[country, 2016:2021, :, :]
slicer_source = idx[country, 2015, :, :]
filling_years = [(2015, slice(2016, 2021)), (2000, slice(1990, 1999))]
for source_year, target_years in filling_years:
slicer_source = idx[country, source_year, :, :]
slicer_target = idx[country, target_years, :, :]
for sector, mapping in mappings.items():
sector_ratio = ratio.loc[
(country, slice(None), slice(None), sector)
].droplevel("lvl2")
energy.loc[slicer_target, mapping["total"]] = cartesian(
sector_ratio.loc[2016:2021, "total"],
energy.loc[slicer_source, mapping["total"]].squeeze(),
sector_ratio.loc[target_years, "total"],
energy.loc[slicer_source, mapping["total"]].squeeze(axis=0),
).values
energy.loc[slicer_target, mapping["elec"]] = cartesian(
sector_ratio.loc[2016:2021, "ele"],
energy.loc[slicer_source, mapping["elec"]].squeeze(),
sector_ratio.loc[target_years, "ele"],
energy.loc[slicer_source, mapping["elec"]].squeeze(axis=0),
).values
level_drops = ["country", "lvl2", "lvl3"]
@ -884,15 +889,18 @@ def rescale_idees_from_eurostat(
nav = ratio.loc[slicer, "total"].droplevel(level_drops)
energy.loc[slicer_target, avia_inter] = cartesian(
avi_i.loc[2016:2021], energy.loc[slicer_source, avia_inter].squeeze()
avi_i.loc[target_years],
energy.loc[slicer_source, avia_inter].squeeze(axis=0),
).values
energy.loc[slicer_target, avia_domestic] = cartesian(
avi_d.loc[2016:2021], energy.loc[slicer_source, avia_domestic].squeeze()
avi_d.loc[target_years],
energy.loc[slicer_source, avia_domestic].squeeze(axis=0),
).values
energy.loc[slicer_target, navigation] = cartesian(
nav.loc[2016:2021], energy.loc[slicer_source, navigation]
nav.loc[target_years],
energy.loc[slicer_source, navigation].squeeze(axis=0),
).values
return energy