Retrieve pollution data by station, in the original units, from the air quality server at Consulta de Concentraciones, or for earlier years use the archive files available from Contaminante, or Meteorología for meteorological data. There's a mistake in the 2016 wind speed data, so for this year, and only this year, the alternative Excel file was used.

get_station_data(criterion, pollutant, year, progress = interactive())

## Arguments

criterion Type of data to download. "HORARIOS" - Hourly data "MAXIMOS" - Daily maximums "MINIMOS" - Daily minimums The type of pollutant to download. "SO2" - Sulfur Dioxide (parts per billion) "CO" - Carbon Monoxide (parts per million) "NOX" - Nitrogen Oxides (parts per billion) "NO2" - Nitrogen Dioxide (parts per billion) "NO" - Nitric Oxide (parts per billion) "O3" - Ozone (parts per billion) "PM10" - Particulate matter 10 micrometers or less (micrograms per cubic meter) "PM25" - Particulate matter 2.5 micrometers or less (micrograms per cubic meter) "WSP" - Wind velocity (meters per second) "WDR" - Wind direction (degrees) "TMP" - Temperature (degrees Celsius) "RH" - Relative humidity (percentage) a numeric vector containing the years for which to download data (the earliest possible value is 1986) whether to display a progress bar (TRUE or FALSE). By default it will only display in an interactive session.

## Value

A data.frame with pollution data. When downloading "HORARIOS" the hours correspond to the Etc/GMT+6 timezone, with no daylight saving time

## Details

Temperature (TMP) archive values are correct to one decimal place, but the most recent data is only available rounded to the nearest integer.

## Warning

The data for the current month is in the process of being validated

stations for a data.frame with the location and names of all pollution measuring stations,

Other raw data functions: get_station_month_data

## Examples

# NOT RUN {
df2$mxc_time <- format(as.POSIXct(paste0(df2$date, " ", df2\$hour, ":00"),