calc_conductance_modflow_perLen: Calculate classic MODFLOW-style stream conductance (per unit...

In scantle/pbjr: Modflow-USG Polyline Boundary Junction (PBJ) R Package

Description

Calculate classic MODFLOW-style stream conductance (per unit length) Per Unit Length Conductance = width x streambed hydraulic conductivity / streambed thickness

Usage

calc_conductance_modflow_perLen(k_str, str_width, str_thick)

Arguments

k_str	numeric, streambed hydraulic conductivity
str_width	numeric, stream width
str_thick	numeric, streambed thickness (head loss zone)

Value

numeric, conductance

Author(s)

Leland Scantlebury

Examples

#-- Read in shapefiles
str <- read_sf(system.file("extdata", "MehlandHill2010_stream.shp", package = "pbjr"))
tri <- read_sf(system.file("extdata", "720_triangles.shp", package = "pbjr"))
vor <- read_sf(system.file("extdata", "720_voronoi.shp", package = "pbjr"))
str <- line_explode(str)

#-- Calculate barycentric weight DF
swdf <- calc_stream_voronoi_weights(stream = str, voronoi = vor, triangles = tri)

#-- Calculate distances
swdf <- stream_elev_from_slope(swdf = swdf, slope = 0.0015, initial_elev = 50)

#-- Calculate conductances
swdf$Conductance <- calc_conductance_modflow_perLen(k_str = 1, str_width = 1,
                                                    thickness = 0.5)

scantle/pbjr documentation built on Dec. 22, 2021, 10:19 p.m.