get_mirror_point: Mirror a well across a boundary

In gopalpenny/anem: Analytical element models for simple aquifers

Description

Mirror a well across a boundary, specified as constant head (CH) or no flow (NF) and defined by the line m x + b)

Usage

get_mirror_point(well, boundary, new_wID = NA)

Arguments

well	single well as a list (or sf object single well feature) containing coordinates (well$x, well$y), pumping rate (well$Q), and well_image ("Actual" or "Image (+Q)" or "Image (-Q)"). which is positive for injection, negative for pumping. Should also have: diam, path, orig_wID.
boundary	single line as a list (or sf object single line feature) containing slope (as boundary$m) and intercept (as boundary$b), as well as the boundary type (bound_type) as "CH" for constant head or "NF" for no flow
new_wID	integer ID for the new well to be created

Value

a row (same class as well input) containing the new imaged well and the following columns: x, y, Q, diam, path, orig_wID, transform (boundary type), source_bound

Examples

## Not run: 
well <- define_wells(wID=1,x=0,y=0,path=1,orig_wID=1)
boundary <- data.frame(bID=1,m=-1,b=1,bound_type="NF")
get_mirror_point(well, boundary, new_wID=2)

well <- define_wells(wID=1,x=0,y=0,Q=0.5,diam=0.75,path=1,orig_wID=1,well_group="a")
boundary <- data.frame(bID=1,m=-1,b=1,bound_type="NF")
get_mirror_point(well, boundary, new_wID=2)

well <- define_wells(wID=1,x=2,y=2,Q=0.5,diam=0.75,path=1,orig_wID=1,well_group="a")
boundary <- data.frame(bID=1,m=Inf,b=1,bound_type="NF")
get_mirror_point(well, boundary, new_wID=2) %>% select(x,y)


well <- define_wells(wID=1,x=2,y=2,Q=0.5,diam=0.75,path=1,orig_wID=1,well_group="a")
boundary <- data.frame(bID=1,m=0,b=1,bound_type="NF")
get_mirror_point(well, boundary, new_wID=2) %>% select(x,y)

## End(Not run)

gopalpenny/anem documentation built on Dec. 20, 2020, 5:27 a.m.