calc_conductance_modflow: Calculate classic MODFLOW-style stream conductance...

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

Description

Calculate classic MODFLOW-style stream conductance Conductance = length x width x streambed hydraulic conductivity / streambed thickness

Usage

calc_conductance_modflow(swdf, k_str, str_width, str_thick)

Arguments

swdf	DataFrame/Geometry of node barycentric coordinates as returned by calc_stream_voronoi_weights
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(swdf, k_str = 1,
                                             str_width = 1, thickness = 0.5)

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