comp = "Tritium_water_Concentration"
x_amanziS_crunch, c_amanziS_crunch = GetXY_AmanziS_1D(
path_to_amanzi, root_amanzi, comp, 1)
struct_crunch = len(x_amanziS_crunch)
except:
struct_crunch = 0
# AmanziU + Native chemistry
try:
comp = 'total_component_concentration.cell.Tritium conc'
input_file = os.path.join("amanzi-u-1d-" + root + ".xml")
path_to_amanzi = "output-u"
run_amanzi_standard.run_amanzi(input_file, 1, [input_file],
path_to_amanzi)
x_amanziU_native, c_amanziU_native = GetXY_AmanziU_1D(
path_to_amanzi, root, comp, 1)
native = len(x_amanziU_native)
except:
native = 0
# HARDWIRE to exclude native
native = 0
# AmanziU + Alqumia + PFloTran chemistry
try:
comp = 'total_component_concentration.cell.Tritium conc'
input_file = os.path.join("amanzi-u-1d-" + root + "-alq-pflo.xml")
path_to_amanzi = "output-u-alq-pflo"
run_amanzi_standard.run_amanzi(
input_file, 1,
["1d-" + root + "-trim.in", root + ".dat", input_file],
local_path = ""
# Amanzi + Native chemistry
try:
input_file = os.path.join("amanzi-u-1d-" + root + ".xml")
path_to_amanzi = "output-u"
run_amanzi_standard.run_amanzi(input_file, 1,
[root + ".bgd", input_file],
path_to_amanzi)
time = timesama[0]
# tot conc
u_amanzi_native = [[] for x in range(len(totcama))]
for j, comp in enumerate(totcama):
x_amanzi_native, c_amanzi_native = GetXY_AmanziU_1D(
path_to_amanzi, 'plot', comp, 1)
u_amanzi_native[j] = c_amanzi_native
# sorb conc
v_amanzi_native = [[] for x in range(len(sorbama))]
for j, sorb in enumerate(sorbama):
x_amanzi_native, c_amanzi_native = GetXY_AmanziU_1D(
path_to_amanzi, 'plot', sorb, 1)
v_amanzi_native[j] = c_amanzi_native
# mineral volume fraction
w_amanzi_native = [[] for x in range(len(vfama))]
for j, vf in enumerate(vfama):
x_amanzi_native, c_amanzi_native = GetXY_AmanziU_1D(
path_to_amanzi, 'plot', vf, 1)
w_amanzi_native[j] = c_amanzi_native
from compare_field_results import GetXY_AmanziU_1D
from compare_field_results import GetXY_AmanziS_1D
if __name__ == "__main__":
# Amanzi: unstructured
try:
input_file = "amanzi_infiltration_sand_loam_1d-u.xml"
path_to_amanzi = "output_2b-u"
root_amanzi = 'case_2b_plot'
run_amanzi_standard.run_amanzi(input_file, 1, [input_file],
path_to_amanzi)
comp = 'pressure.cell.0'
x_amanziU, c_amanziU = GetXY_AmanziU_1D(path_to_amanzi, root_amanzi,
comp, 3)
unstruct = len(x_amanziU)
except:
unstruct = 0
# Amanzi: structured
try:
input_file = "amanzi_infiltration_sand_loam_1d-s.xml"
path_to_amanzi = "output_2b-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,
'Mg++_Sorbed_Concentration', \
'Cl-_Sorbed_Concentration']
# AmanziU + Native chemistry
try:
input_file = os.path.join("amanzi-u-1d-" + root + ".xml")
path_to_amanzi = "output-u"
run_amanzi_standard.run_amanzi(input_file, 1,
[root + ".bgd", input_file],
path_to_amanzi)
u_amanzi_native = [[[] for x in range(len(amanzi_components))]
for x in range(len(times))]
for i, time in enumerate(times):
for j, comp in enumerate(amanzi_components):
x_amanzi_native, c_amanzi_native = GetXY_AmanziU_1D(
path_to_amanzi, root, comp, 1)
u_amanzi_native[i][j] = c_amanzi_native
v_amanzi_native = [[[] for x in range(len(amanzi_sorbed))]
for x in range(len(times))]
for i, time in enumerate(times):
for j, comp in enumerate(amanzi_sorbed):
x_amanzi_native, c_amanzi_native = GetXY_AmanziU_1D(
path_to_amanzi, root, comp, 1)
v_amanzi_native[i][j] = c_amanzi_native
native = True
except:
native = False