str(gwl)
poi <- as.POSIXct(c("2000-01-01", "2001-01-01")) # period of interest
site_no <- "432700112470801" # well USGS 1
is <- gwl$lev_dt >= poi[1] & gwl$lev_dt < poi[2] & gwl$site_no == site_no
d <- gwl[is, ]
plotrix::plotCI(
x = d$lev_dt,
y = d$sl_lev_va,
li = d$sl_lev_va - d$sl_lev_acy_va,
ui = d$sl_lev_va + d$sl_lev_acy_va
)
Run the code above in your browser using DataLab