except:
print('error: parflow directory not found in ', PF_DIR)
raise
# root name for problem
root = "tritium"
# pflotran
path_to_pflotran = "pflotran"
root_pflotran = "1d-" + root
# 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 = ""
path_to_pflotran = "pflotran"
root_pflotran = "1d-" + root
# hardwired for 1d-calcite: time and comp
times = [
'Time: 5.00000E+01 y',
]
components = ['Na+', 'Ca++', 'Mg++', 'Cl-']
try:
u_pflotran = [[[] for x in range(len(components))]
for x in range(len(times))]
for i, time in enumerate(times):
for j, comp in enumerate(components):
x_pflotran, c_pflotran = GetXY_PFloTran_1D(
path_to_pflotran, root_pflotran, time,
'Total_' + comp + ' [M]')
u_pflotran[i][j] = c_pflotran
v_pflotran = [[[] for x in range(len(components))]
for x in range(len(times))]
for i, time in enumerate(times):
for j, comp in enumerate(components):
x_pflotran, c_pflotran = GetXY_PFloTran_1D(
path_to_pflotran, root_pflotran, time,
'Total_Sorbed_' + comp + ' [mol_m^3]')
v_pflotran[i][j] = c_pflotran
pflotran = True
except: