# hardwired for 1d-tritium: time and comp
time = 'Time: 5.00000E+01 y'
comp = 'Total_' + root.title() + ' [M]'
x_pflotran, c_pflotran = GetXY_PFloTran_1D(path_to_pflotran, root_pflotran,
time, comp)
# crunchflow GIMRT
path_to_crunchflow = "crunchflow"
# hardwired for 1d-tritium-crunch.in: time and comp
times_CF = 'totcon5.out'
comp = 0
ignore = 4
x_crunchflow, c_crunchflow = GetXY_CrunchFlow_1D(path_to_crunchflow, root,
times_CF, comp, ignore)
CWD = os.getcwd()
local_path = ""
# subplots
fig, ax = plt.subplots(figsize=(10, 8))
# AmanziS + Alqumia + 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 + "-trim.in", root + ".dat", input_file],
path_to_amanzi)
pH_pflotran = [[] for x in range(len(timespflo))]
comp = 'pH'
for i, time in enumerate(timespflo):
x_pflotran, c_pflotran = GetXY_PFloTran_1D(path_to_pflotran,root_pflotran,time,comp)
pH_pflotran[i] = c_pflotran
# read crunchflow results --->
path_to_crunch = "crunchflow"
try:
u_crunchflow = [ [ [] for x in range(len(compcrunch)) ] for x in range(len(times_CF)) ]
ignore = 4
for i, time in enumerate(times_CF):
for j, comp in enumerate(compcrunch):
x_crunchflow, c_crunchflow = GetXY_CrunchFlow_1D(path_to_crunch,root,time,comp,ignore)
u_crunchflow[i][j] = c_crunchflow
v_crunchflow = [ [ [] for x in range(len(compcrunch)) ] for x in range(len(times_CF_surf)) ]
ignore = 4
for i, time in enumerate(times_CF_surf):
for j, comp in enumerate(compcrunch):
x_crunchflow, c_crunchflow = GetXY_CrunchFlow_1D(path_to_crunch,root,time,comp,ignore)
v_crunchflow[i][j] = c_crunchflow
pH_crunchflow = [ [] for x in range(len(times_CF_pH)) ]
ignore = 4
comp = 0
for i, time in enumerate(times_CF_pH):
y_crunchflow, c_crunchflow = GetXY_CrunchFlow_1D(path_to_crunch,root,time,comp,ignore)
pH_crunchflow[i] = c_crunchflow