Get intersection probability

Get intersection probability

get_intersection_probability(
  wells_array,
  theta_range = c(-pi, pi),
  alpha_range = c(0, 100),
  self_treat = FALSE,
  return_option = 1,
  show_progress = TRUE
)

Arguments

wells_array	Wells object prepared with get_septic_well_array
theta_range	Vector describing the min and max of the uniform distribution for the mean lateral direction of flow
alpha_range	Vector describing the min and max of distance to the groundwater divide
self_treat	Logical. If TRUE, any well at (x = 0, y = 0) will be removed
return_option	Either 1 or 2. See Return, below
show_progress	Logical that determines if progress bar is shown

Value

If return_option = 1: Probability of contamination If return_option = 2: Returns a list containing the well array after clipping and shifting and probabilities of contamination

Details

Note that theta_range must not cross -pi or pi.

Examples

library(units)
wells_array <- get_septic_well_array(hh_grid_example, "septic")
prob <- get_intersection_probability(wells_array, show_progress = FALSE)