From 783e026ebfde8669d40675eab18872b243f523e8 Mon Sep 17 00:00:00 2001 From: Tom Brown Date: Mon, 7 Sep 2020 18:41:07 +0200 Subject: [PATCH] Separate ammonia from basic chemicals in today's industry energy --- Snakefile | 2 ++ ...ustrial_energy_demand_per_country_today.py | 21 ++++++++++++++++++- 2 files changed, 22 insertions(+), 1 deletion(-) diff --git a/Snakefile b/Snakefile index 0b153428..07868043 100644 --- a/Snakefile +++ b/Snakefile @@ -186,6 +186,8 @@ rule build_industrial_production_per_country_tomorrow: script: 'scripts/build_industrial_production_per_country_tomorrow.py' rule build_industrial_energy_demand_per_country_today: + input: + ammonia_production="resources/ammonia_production.csv" output: industrial_energy_demand_per_country_today="resources/industrial_energy_demand_per_country_today.csv" threads: 1 diff --git a/scripts/build_industrial_energy_demand_per_country_today.py b/scripts/build_industrial_energy_demand_per_country_today.py index 91b39576..7e91636f 100644 --- a/scripts/build_industrial_energy_demand_per_country_today.py +++ b/scripts/build_industrial_energy_demand_per_country_today.py @@ -71,6 +71,11 @@ summaries = {} ois_subs = ['Mining and quarrying','Construction','Non-specified'] +#MtNH3/a +ammonia = pd.read_csv(snakemake.input.ammonia_production, + index_col=0)/1e3 + + for ct in eu28: print(ct) @@ -98,7 +103,21 @@ for ct in eu28: summary['Other Industrial Sectors'] = summary[ois_subs].sum(axis=1) summary.drop(columns=ois_subs,inplace=True) - summaries[ct] = summary*ktoe_to_twh + summary.drop(index=['all'],inplace=True) + + summary *= ktoe_to_twh + + summary['Basic chemicals'] += summary['Basic chemicals feedstock'] + summary.drop(columns=['Basic chemicals feedstock'], inplace=True) + + summary['Ammonia'] = 0. + summary.at['gas','Ammonia'] = snakemake.config['industry']['MWh_CH4_per_tNH3_SMR']*ammonia[str(year)].get(ct,0.) + summary.at['electricity','Ammonia'] = snakemake.config['industry']['MWh_elec_per_tNH3_SMR']*ammonia[str(year)].get(ct,0.) + summary['Basic chemicals (without ammonia)'] = summary['Basic chemicals'] - summary['Ammonia'] + summary.loc[summary['Basic chemicals (without ammonia)'] < 0, 'Basic chemicals (without ammonia)'] = 0. + summary.drop(columns=['Basic chemicals'], inplace=True) + + summaries[ct] = summary final_summary = pd.concat(summaries,axis=1)