Example EBAS artic SO4 concentrations
import pyaerocom as pya
pya.const.BASEDIR = '/home/notebook/shared-ns1000k/inputs/pyaerocom-testdata/'
obsdata = pya.io.ReadUngridded().read('EBASMC', 'concso4', lat_range=[66, 90])
obsdata.plot_station_coordinates(var_name='concso4', start=2010, markersize=50);
obsdata.unique_station_names
station = obsdata.to_station_data('Zeppelin mountain (Ny-Ålesund)')
ax = station.plot_timeseries('concso4')
ax = station.plot_timeseries('concso4', freq='monthly', lw=4, ax=ax)