Import Mads in Julia
import Mads
Change the working directory
cd(joinpath(Pkg.dir("Mads"), "examples", "contamination"))
Load Mads input file
md = Mads.loadmadsfile("w01.mads") # load Mads input file into Mads Dictionary
Mads.plotmadsproblem(md, keyword="all_wells")
There are 20 monitoring wells.
Each well has 2 measurement ports: shallow (3 m below the water table labeled a
) and deep (33 m below the water table labeled b
).
Contaminant concentrations are observed for 50 years at each well.
The contaminant transport is solved using the Anasol
package in Mads.
Analysis of the data from only 2 monitoring locations: w13a
and w20a
.
Mads.allwellsoff!(md) # turn off all wells
Mads.wellon!(md, "w13a") # use well w13a
Mads.wellon!(md, "w20a") # use well w20a
Mads.plotmadsproblem(md, keyword="w13a_w20a")
Mads.plotmatches(md, r"w13")
Mads.plotmatches(md, r"w20")
calib_param, calib_results = Mads.calibrate(md)
calib_predictions = Mads.forward(md, calib_param)
Mads.plotmatches(md, calib_predictions, r"w13")
Mads.plotmatches(md, calib_predictions, r"w20")