2 Read Utah monitoring locations

sites=readWQP(type="sites", statecode="US:49")
## 
  |                                                                       
  |                                                                 |   0%
  |                                                                       
  |=================================================================| 100%
plot(LatitudeMeasure~LongitudeMeasure, sites)

Some wrong coordinates are present - get rid of those.

sites=subset(sites, LongitudeMeasure< -100)
plot(LatitudeMeasure~LongitudeMeasure, sites)