sharpshootR (version 1.8)

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, this.day = NULL, 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

this.day

optional integer representing days since start of selected water year

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 https://ncss-tech.github.io/AQP/sharpshootR/CDEC.html, and https://ncss-tech.github.io/AQP/sharpshootR/cumulative-PPT.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.6.1
# 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)
  
}

# }