#!/usr/bin/env python
# coding: utf-8

# In[1]:


import pyaerocom as pya
pya.__version__


# **NOTE**: needs version 0.8.1.dev3 or later!

# In[2]:


pya.const.BASEDIR = '/home/notebook/shared-ns1000k/inputs/pyaerocom-testdata/'


# In[3]:


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`

# In[4]:


read_factory = pya.io.ReadUngridded('GAWTADsubsetAasEtAl')
read_factory.get_reader().PROVIDES_VARIABLES


# In[5]:


obsdata = read_factory.read('GAWTADsubsetAasEtAl', 'concso4')
obsdata


# In[6]:


obsdata.plot_station_coordinates();


# In[7]:


len(obsdata.metadata)


# In[8]:


print(obsdata.unique_station_names)


# In[9]:


station = obsdata.to_station_data('Tana*')
print(station)


# In[10]:


station.plot_timeseries('concso4');


# In[11]:


help(station.compute_trend)


# In[12]:


trend_info = station.compute_trend('concso4', start_year=2001, stop_year=2009)
trend_info


# In[16]:


trend_info = station.compute_trend('concso4', start_year=2011, stop_year=2016)
trend_info


# In[18]:


help(station.trends['concso4'].plot)


# In[19]:


station.trends['concso4'].plot(period='2001-2009');


# In[20]:


station.trends['concso4'].plot(period='2011-2016');


# In[23]:


station.compute_trend('concso4', 2000, 2018)
station.trends['concso4'].plot(period='2000-2018');