Example aas etal read
import pyaerocom as pya
pya.__version__
NOTE: needs version 0.8.1.dev3 or later!
pya.const.BASEDIR = '/home/notebook/shared-ns1000k/inputs/pyaerocom-testdata/'
pya.const.OBS_IDS
NOTE: these are all observation data IDs that are theoretically available, but some of them are not provided in the pyaerocom-testdata database. You can check using pya.browse_database
read_factory = pya.io.ReadUngridded('GAWTADsubsetAasEtAl')
read_factory.get_reader().PROVIDES_VARIABLES
obsdata = read_factory.read('GAWTADsubsetAasEtAl', 'concso4')
obsdata
obsdata.plot_station_coordinates();
len(obsdata.metadata)
print(obsdata.unique_station_names)
station = obsdata.to_station_data('Tana*')
print(station)
station.plot_timeseries('concso4');
help(station.compute_trend)
trend_info = station.compute_trend('concso4', start_year=2001, stop_year=2009)
trend_info
trend_info = station.compute_trend('concso4', start_year=2011, stop_year=2016)
trend_info
help(station.trends['concso4'].plot)
station.trends['concso4'].plot(period='2001-2009');
station.trends['concso4'].plot(period='2011-2016');
station.compute_trend('concso4', 2000, 2018)
station.trends['concso4'].plot(period='2000-2018');