Description Usage Arguments Value Examples
View source: R/anem_preparation.R
Define an aquifer as confined or unconfined and with saturated hydraulic conductivity. Other optional parameters can be defined including boundaries and recharge.
1 | define_aquifer(aquifer_type, Ksat, ...)
|
aquifer_type |
"confined" or "unconfined" |
Ksat |
Saturated hydraulic conductivity |
... |
Optional parameters including boundaries, recharge_params, and aquifer thickness |
This function returns an S3 "aquifer" object that behaves as a list and contains the named items:
aquifer_type: "confined" or "unconfined"
Ksat: Saturated hydraulic conductivity
h0: Hydraulic head of the undisturbed aquifer (i.e., no pumping)
z0: Thickness of the aquifer (applies only to confined aquifers)
bounds: Boundaries of the aquifer, if any (generated using define_boundaries()
)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | (aquifer <- define_aquifer("confined",1e-4))
(aquifer <- define_aquifer("confined",1e-4,h0=100,z0=10))
bounds_df1 <- data.frame(bound_type=c("CH","NF","NF","NF"),
x1=c(0,10,13,1),y1=c(0,10,9,-1),x2=c(10,13,1,0),y2=c(10,9,-1,0))
aquifer_confined <- define_aquifer("confined",1e-3,bounds=bounds_df1,h0=100,z0=10)
print(aquifer_confined)
bounds_df2 <- data.frame(bound_type=c("CH","CH","NF","NF"),
x1=c(0,0,10,10),y1=c(0,10,10,0),x2=c(0,10,10,0),y2=c(10,10,0,0))
aquifer_confined <- define_aquifer("unconfined",1e-3,bounds=bounds_df2,h0=100)
aquifer_confined
recharge_params <- list(recharge_type="F",recharge_vector=c(0,0,3,3),flow=1,x0=3,y0=3)
aquifer <- define_aquifer("confined",1e-3,h0=50,z0=10,recharge=recharge_params)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.