plot_bounds_behavior: Plot behavior on boundaries

In gopalpenny/anem: Analytical element models for simple aquifers

Description

Plot behavior on boundaries with and without well images

Usage

plot_bounds_behavior(wells, aquifer, length.out = 100)

Arguments

wells	wells object with each row containing rate Q [m^3/s], diam [m], radius of influence R [m], & coordinates x [m], y [m]
aquifer	Afuifer object containing aquifer_type, h0, Ksat, bounds, z0 (for confined case only)
length.out	The number of points to evaluate on each boundary

Value

Two ggplot objects that show behavior on the boundaries, one for head and one for flow

Examples

aquifer <- aquifer_confined_example
wells <- generate_image_wells(define_wells(wells_example), aquifer)
gg_list <- plot_bounds_behavior(wells,aquifer,length.out=20)

library(ggplot2)
p_domain <- ggplot() +
  geom_point(data=wells,aes(x,y,fill=Q),color="black",size=2,shape=21) +
  scale_fill_gradient2(low="blue",high="red",mid="gray") +
  geom_segment(data=aquifer$bounds,aes(x1,y1,xend=x2,yend=y2,linetype=bound_type)) +
  coord_equal()
gridExtra::grid.arrange(p_domain,gg_list$p_h,gg_list$p_f,nrow=1)
gg_list$table
gg_list$bounds_behavior

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