pypsa-eur/scripts/build_ammonia_production.py

"""Build ammonia production."""

import pandas as pd

country_to_alpha2 = {
    "Austriae": "AT",
    "Bulgaria": "BG",
    "Belgiume": "BE",
    "Croatia": "HR",
    "Czechia": "CZ",
    "Estonia": "EE",
    "Finland": "FI",
    "France": "FR",
    "Germany": "DE",
    "Greece": "GR",
    "Hungarye": "HU",
    "Italye": "IT",
    "Lithuania": "LT",
    "Netherlands": "NL",
    "Norwaye": "NO",
    "Poland": "PL",
    "Romania": "RO",
    "Serbia": "RS",
    "Slovakia": "SK",
    "Spain": "ES",
    "Switzerland": "CH",
    "United Kingdom": "GB",
}

if __name__ == '__main__':
    if 'snakemake' not in globals():
        from helper import mock_snakemake
        snakemake = mock_snakemake('build_ammonia_production')

    ammonia = pd.read_excel(snakemake.input.usgs,
                            sheet_name="T12",
                            skiprows=5,
                            header=0,
                            index_col=0,
                            skipfooter=19)

    ammonia.rename(country_to_alpha2, inplace=True)

    years = [str(i) for i in range(2013, 2018)]
    countries = country_to_alpha2.values()
    ammonia = ammonia.loc[countries, years].astype(float)

    # convert from ktonN to ktonNH3
    ammonia *= 17 / 14

    ammonia.index.name = "ktonNH3/a"

    ammonia.to_csv(snakemake.output.ammonia_production)