# NOT RUN {
## get data from CDEC
# define station of interest
s <- 'SPW'
# get metadata
s.info <- CDEC_StationInfo(s)
# format title for cumulative PPT
title.text <- sprintf("%s [%s]", s.info$site.meta$Name, s)
# get data
x <- CDECquery(id=s, sensor=45, interval='D', start='1900-01-01', end='2030-01-01')
## NOTE: requires sharpshootR >= 1.4.02
# plot
PCP_plot(x, ylab='Cumulative PPT (inches)', main=title.text)
# }
Run the code above in your browser using DataLab