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

leap_days need to be filtered for daily data too in heat

Browse Source
This commit is contained in:
Fabian Neumann 2022-07-29 20:29:35 +02:00
parent 88ec493b3a
commit ee0c69727f
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/build_heat_demand.py
View File

@ -25,9 +25,12 @@ if __name__ == '__main__':
else: else:
snapshots = snakemake.config['snapshots'] snapshots = snakemake.config['snapshots']
drop_leap_day = snakemake.config["atlite"].get("drop_leap_day", False)
time = pd.date_range(freq='h', **snapshots) time = pd.date_range(freq='h', **snapshots)
if snakemake.config["atlite"].get("drop_leap_day", False): daily = pd.date_range(freq='D', **snapshots)
if drop_leap_day:
time = time[~((time.month == 2) & (time.day == 29))] time = time[~((time.month == 2) & (time.day == 29))]
daily = daily[~((daily.month == 2) & (daily.day == 29))]
cutout = atlite.Cutout(cutout_name).sel(time=time) cutout = atlite.Cutout(cutout_name).sel(time=time)
@ -44,6 +47,6 @@ if __name__ == '__main__':
M = I.T.dot(np.diag(I.dot(stacked_pop))) M = I.T.dot(np.diag(I.dot(stacked_pop)))
heat_demand = cutout.heat_demand( heat_demand = cutout.heat_demand(
matrix=M.T, index=clustered_regions.index) matrix=M.T, index=clustered_regions.index).sel(time=daily)
heat_demand.to_netcdf(snakemake.output[f"heat_demand_{area}"]) heat_demand.to_netcdf(snakemake.output[f"heat_demand_{area}"])