Add pollutantmean.R function