ridge_nb=None)
    if x_axis == 'x':
        plt.scatter(easting, upward, marker=(5, 1), c='red')
        plt.annotate(str(masses), [easting, upward])
    else:
        plt.scatter(northing, upward, marker=(5, 1), c='red')
        plt.annotate(str(masses), [northing, upward])

    #%%

    qratio = [1, 0]
    mesh_dexp, label_dexp = dEXP.dEXP_ratio(xp,
                                            yp,
                                            zp,
                                            gravity,
                                            shape_grid,
                                            zmin=0,
                                            zmax=max_elevation,
                                            nlayers=nlay,
                                            qorders=qratio)
    fig = plt.figure()
    ax = plt.gca()

    plt, cmap = pEXP.plot_xy(mesh_dexp,
                             label=label_dexp,
                             markerMax=True,
                             qratio=str(qratio),
                             p1p2=np.array([p1, p2]),
                             ax=ax,
                             Xaxis=x_axis)  #, ldg=)
    plt.colorbar(cmap)
Esempio n. 2
0
        max(x_raw) / 2 - ano_prop['HWD'][1] / 2,
        max(x_raw) / 2 + ano_prop['HWD'][1] / 2,
        ano_prop['HWD'][2] + ano_prop['HWD'][0] / 2,
        ano_prop['HWD'][2] - ano_prop['HWD'][0] / 2
    ]
    xxzz = [x1, x2, z1, z2]
    CT = ano_prop['SoilR'] / ano_prop['AnoR']

    #%%

    qratio = [1, 0]
    mesh_ratio, label_ratio = dEXP.dEXP_ratio(xp,
                                              yp,
                                              zp,
                                              U,
                                              shape,
                                              zmin=0,
                                              zmax=max_elevation,
                                              nlayers=nlay,
                                              qorders=qratio)

    #%%
    # save data loop

    MESHratio.append(mesh_ratio)
    LABELratio.append(label_ratio)
    XXZZ.append(xxzz)
    CTm.append(CT)

#%%
# Plot the results
Esempio n. 3
0
                             Xaxis=x_axis,
                             Vminmax=[0, 0.35],
                             p1p2=p)
    cbar = plt.colorbar(cmap, shrink=0.25, pad=0.04)
    cbar.set_label('upwc voltage (V)')
    plt.tight_layout()
    plt.savefig('upwc voltage' + str(file) + '.png', dpi=450)

    #%% DEXP ratio
    # x_axis = 'x'
    qratio = [1, 0]
    mesh_dexp, label_dexp = dEXP.dEXP_ratio(Xs,
                                            Ys,
                                            zp,
                                            U,
                                            shape,
                                            zmin=0,
                                            zmax=max_elevation,
                                            nlayers=nlay,
                                            qorders=qratio)
    fig = plt.figure()
    ax = plt.gca()
    plt, cmap = pEXP.plot_xy(mesh_dexp,
                             label=label_dexp,
                             markerMax=True,
                             qratio=str(qratio),
                             Vminmax=[0, 0.075],
                             p1p2=np.array([p1_s, p2_s]),
                             ax=ax,
                             Xaxis=x_axis)  #, ldg=)
    # plt, cmap = pEXP.plot_xy(mesh_dexp, label=label_dexp,