Núvol de SO2 del volcà de la Palma
Per analitzar si el dia 25 i 26 de setembre de 2021 va arribar el núvol de SO2 del volcà de la Palma analitzaràs les 48 dades horàries de les més de 40 poblacions catalanes amb mesuradors automàtics de la Generalitat de Catalunya a l'enllaç disponible a la pàgina mare d'aquesta pàgina (Medi ambient) però en comptes d'emprar el llenguatge R i la biblioteca openair, empraràs Python i la bilbioteca matplotlib
import matplotlib.pyplot as plt
x = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]
y = [ [1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,1,1,1,1,1,2,5,5,5,5,4,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1],
[1,1,1,1,1,1,1,1,1,1,1,1,1,5,1,1,1,1,1,1,1,1,1,1,2,5,5,1,3,3,3,2,3,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1],
[3,3,3,3,3,3,3,4,4,4,3,3,3,4,4,3,3,3,3,3,3,3,4,4,3,3,3,3,2,2,3,3,3,2,3,3,3,3,4,3,3,3,3,3,3,3,3,3] ]
labels=['Igualada', 'Montseny', 'Barcelona (Eixample)']
colors=['r','g','b']
# loop over data, labels and colors
for i in range(len(y)):
plt.plot(x,y[i],'o-',color=colors[i],label=labels[i])
plt.title("AIR SULFUR DIOXIDE IN CATALONIA", loc='left')
plt.xlabel("Hours of weekend 25-26/09/2021")
plt.ylabel("SO2 microgram/m3")
plt.legend()
plt.show()
