# coding: utf-8
#!/usr/bin/env python
import matplotlib.pyplot as plt
import numpy as np

data = np.array([
    [1997, 42105],
    [1998, 41800],
    [1999, 41772],
    [2000, 41905],
    [2001, 42246],
    [2002, 42665],
    [2003, 42581],
    [2004, 42430],
    [2005, 42183],
    [2006, 41945],
    [2007, 42043],
    [2008, 41966],
    [2009, 41881],
    [2010, 41843],
    [2011, 41745],
    [2012, 41570],
    [2013, 41759],
    [2014, 42014],
    [2015, 42731],
    [2016, 43134],
    [2017, 43478],
    [2018, 43893],
    [2019, 44398],
    [2020, 44662],
    [2021, 45118]])

# Prepare axis
fig = plt.figure(figsize=(8, 3))
ax = fig.add_subplot(111)

# Plot data
years = data[:, 0]
people = data[:, 1]
baseline = 40000
ax.fill_between(years, people, y2=baseline, color='#000066', edgecolor='black',
                alpha=0.6, lw=2)

# Cosmetics
plt.xlim(min(years), max(years))
ax.minorticks_on()
fig.subplots_adjust(bottom=0.15)
plt.grid()
plt.title(u'Bevölkerungsentwicklung in Kaufbeuren')
plt.xlabel('Jahr')
plt.ylabel('Einwohner')

# Save figure
plt.savefig('KaufbeurenBevoelkerung.svg', transparent=True)