else:
p1 = [(northing[0]), 0]
p2 = [(northing[0]), 20e3]
x_axis = slicedir[1] # direction of p1p2 profil
SI = 2
pEXP.plot_line(xp, yp, gravity, p1, p2, interp=True, Xaxis=x_axis)
pEXP.plot_field(xp, yp, gravity, shape_grid)
#%%
# Upward continuation of the field data with discretisation in altitude controlled by the number of layers (nlay) and the maximum elevation desired (max_elevation)
mesh, label_prop = dEXP.upwc(xp,
yp,
zp,
gravity,
shape_grid,
zmin=0,
zmax=max_elevation,
nlayers=nlay,
qorder=qorder)
plt, cmap = pEXP.plot_xy(mesh,
label=label_prop,
Xaxis=x_axis,
p1p2=np.array([p1, p2]))
plt.colorbar(cmap)
#%%
xderiv = transform.derivx(xp, yp, gravity, shape_grid, order=qorder)
yderiv = transform.derivy(xp, yp, gravity, shape_grid, order=qorder)
xp, yp, U, shape = dEXP.pad_edges(xp, yp, U, shape, pad_type=0) # reflexion=5
p1 = [min(yp), 0]
p2 = [max(yp), 0]
x_axis = 'y'
pEXP.plot_line(xp, yp, U, p1, p2, interp=interp, Xaxis=x_axis)
#%%
# Upward continuation of the field data
mesh, label_prop = dEXP.upwc(xp,
yp,
zp,
U,
shape,
zmin=0,
zmax=max_elevation,
nlayers=nlay,
qorder=qorder)
plt, cmap = pEXP.plot_xy(mesh, label=label_prop, Xaxis=x_axis)
plt.colorbar(cmap)
#%%
# ridges identification
# dEXP.ridges_minmax_plot(xp, yp, mesh, p1, p2,
# label=label_prop,
# fix_peak_nb=2,
# method_peak='find_peaks')
x_resolution=interp_size,
savefig=False,
interp=interp,
smooth=smooth,
Xaxis=x_axis)
plt.savefig('zderiv' + str(file) + '.png', dpi=450)
#%% ------- upward continuation of the field data
p = [p1_s, p2_s]
mesh, label_prop = dEXP.upwc(XFs,
YFs,
zp,
UF,
shape,
zmin=0,
zmax=max_elevation,
nlayers=nlay,
qorder=qorder)
plt, cmap = pEXP.plot_xy(mesh,
label=label_prop,
Xaxis=x_axis,
Vminmax=[0, 0.0125],
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)
