Skip to content

Build cop profiles

Snakemake rule script to calculate the heat pump coefficient of performance with atlite

build_cop_profiles(pop_map, cutout, output_path)

Build COP time profiles with atlite and write outputs to output_path as hf5

Parameters:

Name Type Description Default
pop_map DataFrame

the population map (node resolution)

 required
cutout cutout

the atlite cutout (weather data)

 required
output_path PathLike

the path to write the output to as hdf5

 required
Source code in workflow/scripts/build_cop_profiles.py 
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
def build_cop_profiles(pop_map: pd.DateOffset, cutout: atlite.cutout, output_path: os.PathLike):
    """Build COP time profiles with atlite and write outputs to output_path as hf5

    Args:
        pop_map (pd.DataFrame): the population map (node resolution)
        cutout (atlite.cutout): the atlite cutout (weather data)
        output_path (os.PathLike): the path to write the output to as hdf5
    """

    pop_matrix = sp.sparse.csr_matrix(pop_map.T)
    index = pop_map.columns
    index.name = "provinces"

    soil_temp = cutout.soil_temperature(matrix=pop_matrix, index=index)
    soil_temp["time"] = (
        pd.DatetimeIndex(soil_temp["time"].values, tz="UTC")
        .tz_convert(TIMEZONE)
        .tz_localize(None)
        .values
    )

    with pd.HDFStore(snakemake.input.temp, mode="r") as store:
        temp = store["temperature"]

    source_T = temp
    source_soil_T = soil_temp.to_pandas().divide(pop_map.sum())

    # quadratic regression based on Staffell et al. (2012)
    # https://doi.org/10.1039/C2EE22653G

    sink_T = 55.0  # Based on DTU / large area radiators

    delta_T = sink_T - source_T

    # TODO make this user set and document
    # For ASHP
    def ashp_cop(d):
        return 6.81 - 0.121 * d + 0.000630 * d**2

    cop = ashp_cop(delta_T)

    delta_soil_T = sink_T - source_soil_T

    # TODO make this user set and document
    # For GSHP
    def gshp_cop(d):
        return 8.77 - 0.150 * d + 0.000734 * d**2

    cop_soil = gshp_cop(delta_soil_T)

    with pd.HDFStore(output_path, mode="w", complevel=4) as store:
        store["ashp_cop_profiles"] = cop
        store["gshp_cop_profiles"] = cop_soil