Ejemplo n.º 1
0
            "1d-" + root + "-crunch.in", root + ".dbs", "aqueous" + ".dbs",
            input_file
        ], path_to_amanzi)
        x_amanziU_crunch, c_amanziU_crunch = GetXY_AmanziU_1D(
            path_to_amanzi, root, comp, 1)
        unstruct_crunch = len(x_amanziU_crunch)

    except:
        unstruct_crunch = False

    # parflow + pflotran
    try:
        path_to_parflow = os.path.join(PF_DIR, root + '_1d', 'pflotran')

        compPF = "tritium_pf.out.PrimaryMobile.00.Tritium.00005.txt"
        x_parflow_pflo, c_parflow_pflo = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        parflow_pflo = len(x_parflow_pflo)

    except:
        parflow_pflo = 0

    # parflow + crunch
    try:
        path_to_parflow = os.path.join(PF_DIR, root + '_1d', 'crunch')

        compPF = "tritium_pf.out.PrimaryMobile.00.Tritium.00005.txt"
        x_parflow_crunch, c_parflow_crunch = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        parflow_crunch = len(x_parflow_crunch)

    except:
Ejemplo n.º 2
0
        compS = "A_Sorbed_Concentration"
        x_amanziS_crunch, v_amanziS_crunch = GetXY_AmanziS_1D(
            path_to_amanzi, root_amanziS, compS, 1)
        struct_c = len(x_amanziS_crunch)
    except:
        struct_c = 0

    # parflow + pflotran
    try:
        path_to_parflow = os.path.join(PF_DIR, root + '_1d', 'pflotran')

        c_parflow_pflo = [[] for x in range(len(components))]
        v_parflow_pflo = [[] for x in range(len(components))]

        compPF = "isotherms_pf_pflotran.out.PrimaryMobile.00.A.00005.txt"
        x_parflow_pflo, c_parflow_pflo[0] = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        compPF = "isotherms_pf_pflotran.out.PrimaryMobile.01.B.00005.txt"
        x_parflow_pflo, c_parflow_pflo[1] = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        compPF = "isotherms_pf_pflotran.out.PrimaryMobile.02.C.00005.txt"
        x_parflow_pflo, c_parflow_pflo[2] = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        parflow_pflo = len(x_parflow_pflo)

        compPF = "isotherms_pf_pflotran.out.PrimarySorbed.00.A.00005.txt"
        x_parflow_pflo, v_parflow_pflo[0] = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        compPF = "isotherms_pf_pflotran.out.PrimarySorbed.01.B.00005.txt"
        x_parflow_pflo, v_parflow_pflo[1] = GetXY_ParFlow_1D_100(
            compPF, path=path_to_parflow)
        compPF = "isotherms_pf_pflotran.out.PrimarySorbed.02.C.00005.txt"