Skip to content
Snippets Groups Projects
Commit ce1c4dfc authored by Anna Lanteri's avatar Anna Lanteri
Browse files

adding inventories to emissions

parent 514ad49d
No related branches found
No related tags found
No related merge requests found
%% Cell type:markdown id:14c2e348-a778-4c62-914f-740d8248493b tags:
# Historical emissions
%% Cell type:code id:b5954e9e-dfcf-45d6-bd71-275eb1bbb845 tags:
``` python
import xarray as xr
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import glob
from IPython.display import Image
```
%% Cell type:code id:26750783-2379-4532-8520-736187905daf tags:
``` python
path_hist = '/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/*DECK_road_*.nc'
files_hist = glob.glob(path_hist)
files_hist
```
%% Output
['/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_192001-192912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_192001-192912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_176001-176912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_201001-201012.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_197001-197912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_195001-195912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_183001-183912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_200001-200912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_198001-198912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_193001-193912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_180001-180912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_185001-185912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_201401-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_183001-183912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_201301-201312.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_193001-193912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_176001-176912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_175001-175912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_199001-199912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_194001-194912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_175001-175912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_186001-186912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_201101-201112.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_195001-195912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_185001-185912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_201401-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_186001-186912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_198001-198912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_187001-187912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_177001-177912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_197001-197912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_184001-184912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_186001-186912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_189001-189912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_200001-200912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_195001-195912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_201001-201012.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_181001-181912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_180001-180912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_194001-194912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_187001-187912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_183001-183912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_179001-179912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_178001-178912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_200001-200912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_180001-180912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_191001-191912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_190001-190912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_175001-175912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_178001-178912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_187001-187912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_175001-175912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_184001-184912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_201101-201112.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_188001-188912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_190001-190912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_175001-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_175001-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_201301-201312.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_201201-201212.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_177001-177912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_191001-191912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_201201-201212.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_191001-191912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_185001-185912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_188001-188912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_200001-200912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_183001-183912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_198001-198912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_196001-196912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_184001-184912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_179001-179912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_181001-181912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_182001-182912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_190001-190912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_178001-178912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_188001-188912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_184001-184912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_199001-199912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_182001-182912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_197001-197912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_201401-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_201201-201212.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_177001-177912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_188001-188912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_201001-201012.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_181001-181912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_191001-191912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_199001-199912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_175001-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_182001-182912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_175001-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_193001-193912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_196001-196912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_189001-189912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_198001-198912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_201301-201312.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_189001-189912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_201001-201012.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_196001-196912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_190001-190912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_189001-189912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_196001-196912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_179001-179912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_178001-178912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_201101-201112.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_192001-192912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_180001-180912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_193001-193912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_201301-201312.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_186001-186912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_177001-177912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_195001-195912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_201101-201112.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_199001-199912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_197001-197912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_192001-192912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_176001-176912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_182001-182912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_187001-187912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_179001-179912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_194001-194912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMHC_176001-176912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_201201-201212.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_185001-185912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_181001-181912.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_aerosol_201401-201412.nc',
'/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_NMVOC_194001-194912.nc']
%% Cell type:code id:4eeb225e-79d7-4670-a0ba-74dd587fc28b tags:
``` python
path_road_misc = '/pool/data/MESSY/DATA/MESSy2/raw/offemis/CMIP6/CMIP6v6.1_DLR1.0_DECK_road_MISC_175001-201412.nc'
```
%% Cell type:code id:1a680b77-aeae-4651-b969-8f3ad062c545 tags:
``` python
seconds_in_year = 31536000.
area_earth = 5.1006447295*pow(10,14)
k_to_tera = pow(10, -9)
```
%% Cell type:code id:1f204678-097f-4519-96f3-22a871e553fe tags:
``` python
data_road_misc = xr.open_dataset(path_road_misc)
```
%% Cell type:code id:0d8a2c7f-c613-46f6-ac2c-a27f9b2bb5bf tags:
``` python
data_road_misc
```
%% Output
<xarray.Dataset>
Dimensions: (lon: 720, bnds: 2, lat: 360, lev: 1, time: 3180)
Coordinates:
* lon (lon) float64 -179.8 -179.2 -178.8 -178.2 ... 178.8 179.2 179.8
* lat (lat) float64 -89.75 -89.25 -88.75 -88.25 ... 88.75 89.25 89.75
* lev (lev) float64 45.0
* time (time) object 1750-01-16 00:00:00 ... 2014-12-16 00:00:00
Dimensions without coordinates: bnds
Data variables:
lon_bnds (lon, bnds) float64 ...
lat_bnds (lat, bnds) float64 ...
time_bnds (time, bnds) object ...
CO_flux (time, lev, lat, lon) float32 ...
NH3_flux (time, lev, lat, lon) float32 ...
NOx_flux (time, lev, lat, lon) float32 ...
SO2_flux (time, lev, lat, lon) float32 ...
Attributes: (12/39)
CDI: Climate Data Interface version 1.7.0 (http:/...
Conventions: CF-1.6
history: Tue Aug 28 23:02:02 2018: cdo -O -Q copy /sc...
source: CEDS-2017-08-30: Community Emissions Data Sy...
institution: Pacific Northwest National Laboratory - Join...
reference1: Lamarque et al.(2010), doi:10.5194/acp-10-70...
... ...
data_usage_tips: Note that these are monthly average fluxes.
reporting_unit: Mass flux of BC, reported as carbon mass
nominal_resolution: 50 km
source_id: CEDS-2017-08-30
tracking_id: hdl:21.14100/649247a5-afc6-4e25-b777-9d7a77a...
CDO: Climate Data Operators version 1.7.0 (http:/...
%% Cell type:code id:4cd472c4-a52b-4d43-94d0-24b6d3e42e3d tags:
``` python
co_flux_road_misc = data_road_misc['CO_flux'] # kg m-2 s-1
nh3_flux_road_misc = data_road_misc['NH3_flux'] # kg m-2 s-1
nox_flux_road_misc = data_road_misc['NOx_flux'] # kg(NO2) m-2 s-1
so2_flux_road_misc = data_road_misc['SO2_flux'] # kg m-2 s-1
```
%% Cell type:code id:83759f6f-2c64-4c0a-a392-aa915093f3d5 tags:
``` python
def zonal_average(emission):
# Average over lev
emission = emission.mean('lev')
# Create a weight array for the zonal average
weight_emission = np.cos(np.deg2rad(emission.lat))
weight_emission = weight_emission / weight_emission.sum()
# Average over lon, zonal average over lan
emission = emission.mean('lon')
emission_weighted = []
for year in range(emission.shape[0]):
emission_year_weighted = np.sum(emission[year, :] * weight_emission).compute().data.item()
emission_weighted.append(emission_year_weighted)
return emission_weighted
```
%% Cell type:code id:cb911dad-93ab-4ad5-bbcf-372c1d17f7c9 tags:
``` python
def total_emissions(flux):
yearly_flux = zonal_average(flux.resample(time="1YS").mean(dim="time"))
return np.array(yearly_flux) * area_earth * seconds_in_year * k_to_tera
```
%% Cell type:code id:ade27428-d000-479d-b279-f48270fd50c1 tags:
``` python
co_yearly_flux_road_misc = zonal_average(co_flux_road_misc.resample(time="1YS").mean(dim="time"))
```
%% Cell type:code id:8a608640-13a1-46cc-9844-af3862a50f23 tags:
``` python
co_yearly_emissions_road_misc = np.array(co_yearly_flux_road_misc) * area_earth * seconds_in_year * k_to_tera
```
%% Cell type:code id:d9ab5edf-e5d7-4938-9426-45db65b20112 tags:
``` python
nox_yearly_emissions_road_misc = total_emissions(nox_flux_road_misc)
```
%% Cell type:code id:9a68eccf-408b-442e-87c3-d3ed948a2182 tags:
``` python
fig, axs = plt.subplots(1,2)
fig.set_figwidth(12)
N = 265
dates = pd.date_range("1/1/1750", periods=265, freq='Y')
axs[0].plot(dates, co_yearly_emissions_road_misc)
axs[0].set_title('CO emissions')
axs[0].set_xlabel('Year')
axs[0].set_ylabel('Tg/a')
axs[1].plot(dates, nox_yearly_emissions_road_misc)
axs[1].set_title('nox emissions')
axs[1].set_xlabel('Year')
axs[1].set_ylabel('Tg/a')
fig.tight_layout()
```
%% Output
%% Cell type:code id:07a908a1-5a7d-4dd8-8ae5-260bfba993aa tags:
``` python
N = 57
dates = pd.date_range("1/1/1958", periods=57, freq='Y')
plt.plot(dates, test)
plt.title('CO emissions')
plt.xlabel('Year')
plt.ylabel('Tg/a')
fig.tight_layout()
```
%% Output
%% Cell type:code id:d695b8f6-f060-469e-aa91-22d6cb705524 tags:
``` python
```
......
Source diff could not be displayed: it is too large. Options to address this: view the blob.
Source diff could not be displayed: it is too large. Options to address this: view the blob.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment