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(