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)
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
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,