예제 #1
0
        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],
예제 #2
0
    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
예제 #3
0
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,
예제 #4
0
                         '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