write.DRNpackage: Writes Modflow-USG Drain Package

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

Description

Writes Modflow-USG Drain Package

Usage

write.DRNpackage(drndf, filename, nSPs, IDRNCB, SPwarnings = T)

Arguments

drndf	DataFrame/Geometry of drain nodes as returned by calc_stream_drn with conductances and elevation columns required for stress periods as 'cond[sp]' and 'elev[sp]'
filename	character, name/location of output file
nSPs	integer, number of stress periods in simulation
IDRNCB	integer, CBB flow flag. See details (or PBJ package manual) for more info
SPwarnings	T/F (optional) turn on (True) or off (False) warnings about reused or missing SP data

Author(s)

Leland Scantlebury

See Also

calc_stream_drn write.PBJpackage

Examples

#-- Read in shapefiles
stream <- read_sf(system.file("extdata", "straight_river.shp", package = "pbjr"))
tri <- read_sf(system.file("extdata", "straight_triangles.shp", package = "pbjr"))
vor <- read_sf(system.file("extdata", "straight_voronoi.shp", package = "pbjr"))

#-- Explode polyline
stream <- line_explode(stream)

#-- Create DRNDF
drndf <- calc_stream_drn(stream = stream, voronoi = vor)

#-- Add conductance
drndf$cond1 <- calc_conductance_modflow(drndf, k_streambed = 1,
                                        str_width = 1, thickness = 0.5)
#-- Add elevation
drndf$elev1 <- 90.0

#-- Write file
write.DRNpackage(drndf, filename = paste0(tempdir(),'/straight_model.drn'), nSPs=2, IDRNCB=50)

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