sharpshootR (version 1.8)

plotWB: Visualize Monthly Water Balance

Description

This function offers one possible visualization for the results of monthlyWB.

Usage

plotWB(
  WB,
  AWC,
  showAWC = "below",
  sw.col = "#377EB8",
  surplus.col = "#4DAF4A",
  et.col = "#E41A1C",
  deficit.col = "#FF7F00",
  pch = c("P", "E"),
  pt.cex = 0.85,
  lwd = 2,
  n.ticks = 8,
  grid.col = grey(0.65),
  month.cex = 1
)

Arguments

WB

output from monthlyWB

AWC

available water-holding capacity (mm)

showAWC

styling to show box for AWC below (default) or above zero point of axis

sw.col

color for soil water

surplus.col

color for surplus water

et.col

color for ET

deficit.col

color for deficit

pch

plotting character for PPT and PET points (c('P', 'E'))

pt.cex

character expansion factor for PPT and PET points

lwd

line width for PPT and PET curves

n.ticks

approximate number of tick marks on positive and negative y-axis

grid.col

horizontal grid line color

month.cex

character expansion factor for month labels (x-axis)

Examples

# NOT RUN {
if(requireNamespace('hydromad')) {

## 100mm (4") AWC
# monthly PPT and PET extracted from monthly estimates at Sonora, CA
AWC <- 100
PPT <- c(171, 151, 138, 71, 36, 7, 1, 2, 11, 48, 102, 145)
PET <- c(15.17, 18.26, 30.57, 42.95, 75.37, 108.05, 139.74, 128.9, 93.99, 59.84, 26.95, 14.2)

# water-year
# three years
x.wb <- monthlyWB(AWC, PPT, PET, S_init = 0, starting_month = 9, rep = 3)
x.wb[x.wb$mo == 'Sep', ]

# new-style representation, with AWC "above 0"
plotWB(x.wb, AWC, showAWC = 'below')

# old-style representation, with AWC "above 0"
plotWB(x.wb, AWC, showAWC = 'above')


# water year
# last iteration
x.wb <- monthlyWB(AWC, PPT, PET, S_init = 0, 
                  starting_month = 9, rep = 3, 
                  keep_last = TRUE
)

plotWB(x.wb, AWC)

}

# }