AOD johan
import pyaerocom as pya
import matplotlib.pyplot as plt
import pandas as pd
from warnings import filterwarnings
filterwarnings('ignore')
pya.change_verbosity('critical', log=pya.const.print_log) # don't output warnings
import seaborn as sns
# Use seaborn style defaults and set the default figure size
#sns.set(rc={'figure.figsize':(11, 4)})
pya.__version__
pya.const.BASEDIR = './shared-ns1000k/inputs/pyaerocom-testdata'
pya.const.BASEDIR
obs_reader = pya.io.ReadUngridded('AeronetSDAV3Lev2.daily')
print(obs_reader)
obs_data = obs_reader.read(vars_to_retrieve='od550aer')
obs_data #displays data type
DATA_ID = 'TM5_AP3-CTRL2019'
_model_data = pya.io.ReadGridded(DATA_ID)
model_data = _model_data.read_var('od550aer', start=2010)
display(model_data)
# Example accessing data
axes = obs_data['Leipzig']['od550aer']['2015':'2015']
axes.plot(style='.')
col_data = pya.colocation.colocate_gridded_ungridded(model_data,
obs_data,
ts_type='monthly',
start=2010)