The Emissions Database for Global Atmospheric Research (EDGAR) is a project from the Joint Research Centre. This function reads the NetCDF and merge/aggregate into diferent chemical mechanisms
edgar_chem(path, chem)
path | Character; path to the NetCDF files from EDGAR. The directory must have one file for each of the following pollutants: "voc" from 1 to 25, "co", "nox", "nmvoc","so2", "nh3", "pm10", "pm2.5", "bc" and "oc" |
---|---|
chem | Character; chemical mechanism: "edgar", "radm", "radmsorg", "cbmz_mosaic", "cptec" and "ecb05_opt1".
|
RasterStack
Molecular weights were obtained from Development of Improved Chemical Speciation Database for Processing Emissions of Volatile Organic Compounds for Air Quality Models https://intra.engr.ucr.edu/~carter/emitdb/
Some mappings were obtained from:
Lopez-Norena, Ana and Fernandez, Rafael & Puliafito, SALVADOR. (2019). ESPECIACION DE INVENTARIOS DE EMISIONES DE AEROSOLES Y COMPUESTOS ORGANICOS VOLATILES PARA EL MODELO WRF-CHEM, APLICADO A LOS ESQUEMAS RADM-2, CBM-Z Y MOZART-4.
if (FALSE) { # do not run # Get EDGAR #### get_edgar(dataset = "v432_VOC_spec", destpath = "V50_432_AP/TOT/", sector = c( "TOTALS"), type = "nc", year = 2012) get_edgar(dataset = "v50_AP", destpath = "V50_432_AP/TOT", sector = c( "TOTALS"), type = "nc", year = 2014) get_edgar(dataset = "v432_VOC_spec", destpath = "V50_432_AP/TRO/", sector = c( "TRO"), type = "nc", year = 2012, ask = F) get_edgar(dataset = "v50_AP", destpath = "V50_432_AP/TRO", sector = c( "TRO_RES", "TRO_noRES"), type = "nc", year = 2014) totals <- list.files(path = "V50_432_AP/TOT/", full.names = TRUE, pattern = ".zip") lapply(totals, unzip, exdir = "V50_432_AP/TOT//") tros <- list.files(path = "V50_432_AP/TRO", full.names = TRUE, pattern = ".zip") lapply(tros, unzip, exdir = "V50_432_AP/TRO/") edgar_chem("V50_432_AP/TOT", "radm") }