comp, 3)
unstruct = len(x_amanziU)
except:
unstruct = 0
# Amanzi: structured
try:
input_file = "amanzi_infiltration_clay_sand_1d-s.xml"
path_to_amanzi = "output_2c-s"
root_amanzi = "plot"
run_amanzi_standard.run_amanzi(input_file, 1, [input_file],
path_to_amanzi)
comp = "water_Pressure"
x_amanziS, c_amanziS = GetXY_AmanziS_1D(path_to_amanzi, root_amanzi,
comp, 2)
struct = len(x_amanziS)
except:
struct = 0
# Amanzi: analytic
try:
path_to_amanzi = "golden_output"
root_amanzi = 'case_2c_plot'
comp = 'pressure.cell.0'
x_amanziU_gold, c_amanziU_gold = GetXY_AmanziU_1D(
path_to_amanzi, root_amanzi, comp, 3)
unstruct_gold = len(x_amanziU_gold)
except:
unstruct_gold = 0
alq_crunch = len(x_amanzi_alquimia_crunch)
except:
alq_crunch = 0
# AmanziS + Alquimia + PFloTran chemistry
try:
input_file = os.path.join("amanzi-s-1d-" + root + "-alq-pflo.xml")
path_to_amanzi = "output-s-alq-pflo"
run_amanzi_standard.run_amanzi(
input_file, 1, ["1d-" + root + ".in", root + ".dat", input_file],
path_to_amanzi)
root_amanziS = "plt"
compS = "tracer_water_Concentration"
x_amanziS, c_amanziS = GetXY_AmanziS_1D(path_to_amanzi, root_amanziS,
compS, 1)
struct = len(x_amanziS)
except:
struct = 0
# AmanziS + Alquimia + CrunchFlow chemistry
try:
input_file = os.path.join("amanzi-s-1d-" + root + "-alq-crunch.xml")
path_to_amanziS = "output-s-alq-crunch"
run_amanzi_standard.run_amanzi(
input_file, 1,
["1d-" + root + "-crunch.in", root + ".dbs", input_file],
path_to_amanziS)
root_amanziS = "plt"
compS = "tracer_water_Concentration"
x_amanziS_crunch, c_amanziS_crunch = GetXY_AmanziS_1D(
