#! /usr/bin/env python
# -*- coding:utf8 -*-

import matplotlib.pyplot as plt
from math import *

plt.style.use('classic')
years = range(1918, 2019+1)
incidents = [12, 64, # 191x
 96,  98,  81,  51,  37,  56,  84, 111, 188, 240, # 192x
183, 163, 145, 153, 138, 181, 218, 163, 204, 239, # 193x
456, 331, 690, 809, 933, 733, 423, 370, 391, 349, # 194x
425, 491, 406, 335, 292, 212, 172, 213, 217, 196, # 195x
186, 211, 199, 158, 219, 257, 249, 311, 300, 320, # 196x
298, 271, 344, 333, 270, 316, 277, 340, 356, 328, # 197x
325, 272, 250, 238, 234, 261, 238, 277, 254, 265, # 198x
261, 240, 266, 275, 231, 266, 251, 232, 225, 221, # 199x
198, 211, 198, 201, 178, 194, 193, 171, 190, 163, # 200x
162, 154, 156, 139, 122, 122, 102, 101, 113, 125] # 201x

fsize = 20
fig = plt.figure(figsize=(600 / 90.0, 450 / 90.0), dpi=72)
plt.plot(years, incidents, 'bo-', linewidth=1, markersize=5)
plt.grid(True)
plt.xlim(floor(min(years)/10.)*10., ceil(max(years)/10.)*10.)
plt.ylim(0, 600)
plt.gca().ticklabel_format(useOffset=False)
plt.ylabel('Number of accidents')
plt.title('ACRO: air incidents %s-%s' % (min(years), max(years)))
plt.tight_layout()
plt.savefig('ACRO_incidents.svg')