Welcome to sulfure dioxide!

Cumbrevieja is a volcano located in Palma island. It is an active volcano producing sulfur dioxide gases in the atmosphere.

The weekend 25 to 26 of Semptembre 2021 polution from the volcano arrived to Catalonia as you can see in the following graph.

The graph was created using mathplotlib library in python and the code is below

Data are obtained from the air pollution website of Generalitat de Catalunya.


  Import matplotlib.pyplot as plt  # it makes accesible all the instructions and codes contained in matplotlib n order to create grapghs

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]  
# x data is 48 hours data obtained from Generalitat date base corresponding to SO2 in the weekend of interest
y = [ [4,4,4,4,4,4,4,4,4,4,4,3,3,3,3,3,2,3,2,2,3,2,2,2,2,2,2,2,2,2,2,2,2,3,3,6,6,6,6,5,5,6,5,6,5,6,5,5],
      [6,6,5,6,6,5,5,5,5,6,5,6,7,6,7,7,8,7,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,4,5,6,7,5,5,4,2,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,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] ]
      # y data are SO2 levels hour by hour in the 3 cities indicates in labels
      
labels=['Rubí', 'Sant Seloni', 'Sant Vicenç']
colors=['r','g','b'] 
# colours of the  3 cities are  r:red, g: green and b: blue

# 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()