sharpshootR (version 1.6)

PCP_plot: Percentiles of Cumulative Precipitation

Description

Generate a plot representing percentiles of cumulative precipitation, given a historic record, and criteria for selecting a year of data for comparison.

Usage

PCP_plot(x, this.year, method = "exemplar", 
q.color = "RoyalBlue", c.color = "firebrick", ...)

Arguments

x

result from CDECquery() for now, will need to generalize to other sources

this.year

a single water year, e.g. 2020

method

'exemplar' or 'daily', currently 'exemplar' is the only method available

q.color

color of percentiles cumulative precipitation

c.color

color of selected year

…

addtional arguments to plot()

Value

Currently nothing is returned.

Details

This is very much a work in progress. Further examples at http://ncss-tech.github.io/AQP/sharpshootR/CDEC.html

See Also

waterDayYear

Examples

# NOT RUN {
if(requireNamespace("curl") &
   curl::has_internet()
   ) {
  
  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='2000-01-01', end='2030-01-01')

## NOTE: requires sharpshootR >= 1.4.02
# plot
par(mar=c(4.5, 4.5, 2.5, 1.5))
PCP_plot(x[1:(nrow(x)-60), ], ylab='Cumulative PPT (inches)', main=title.text, this.year = 2020)
  
}

# }