fbff32dcfc
Previously the DH share was being multiplied by the population weighting, reducing the DH share with multiple nodes.
35 lines
1.1 KiB
Python
35 lines
1.1 KiB
Python
# -*- coding: utf-8 -*-
|
|
# SPDX-FileCopyrightText: : 2020-2023 The PyPSA-Eur Authors
|
|
#
|
|
# SPDX-License-Identifier: MIT
|
|
"""
|
|
Distribute country-level energy demands by population.
|
|
"""
|
|
|
|
import pandas as pd
|
|
|
|
if __name__ == "__main__":
|
|
if "snakemake" not in globals():
|
|
from _helpers import mock_snakemake
|
|
|
|
snakemake = mock_snakemake(
|
|
"build_population_weighted_energy_totals",
|
|
simpl="",
|
|
clusters=48,
|
|
)
|
|
|
|
pop_layout = pd.read_csv(snakemake.input.clustered_pop_layout, index_col=0)
|
|
|
|
energy_totals = pd.read_csv(snakemake.input.energy_totals, index_col=0)
|
|
|
|
nodal_energy_totals = energy_totals.loc[pop_layout.ct].fillna(0.0)
|
|
nodal_energy_totals.index = pop_layout.index
|
|
nodal_energy_totals = nodal_energy_totals.multiply(pop_layout.fraction, axis=0)
|
|
|
|
# district heating share should not be divided by population fraction
|
|
dh_share = energy_totals["district heat share"].loc[pop_layout.ct].fillna(0.0)
|
|
dh_share.index = pop_layout.index
|
|
nodal_energy_totals["district heat share"] = dh_share
|
|
|
|
nodal_energy_totals.to_csv(snakemake.output[0])
|