calc_WP: Calculate Wetted Perimeter

In jschap1/WSEw: River Geometry Package for Studying Water Surface Elevation-Width Relationships

Description

Calculates wetted perimeter for a particular water surface elevation, ranging from empty channel to bankfull conditions

Usage

calc_WP(x, b, WSE, method = "wide")

Arguments

x	distance from horizontal datum
b	channel bed elevation
WSE	water surface elevation
method	"neal" for Jeff Neal's approximation in Neal et al. (2015), "wide" for wide-channel approximation

Details

Uses a numerical method to approximate the integral of the difference between WSE and bed elevation The numerical method only counts area where the WSE is higher than the bed elevation Needs more work to figure out how to implement arbitrary spatial discretization Performs poorly for a small number of data points or WSE close to bed elevation, see the example

Value

WP wetted perimeter

Examples

calc_WP(x = cross_sections$x[[1]], b = cross_sections$b[[1]], WSE = 138)
x <- c(1,1,1,2,3,4,4,4)
b <- c(3,2,1,1,1,1,2,3)
WP <- calc_WP(x, b, WSE = 3)
WP <- calc_WP(x, b, WSE = 2)

jschap1/WSEw documentation built on May 12, 2019, 10:50 a.m.