Added rough code to poluttantmean.R