1
Fork 0

Post-process output images

This commit is contained in:
Edgar P. Burkhart 2022-03-03 11:16:40 +01:00
parent 420b9cfd4f
commit 7061abe284
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227
2 changed files with 27 additions and 11 deletions

View file

@ -13,4 +13,7 @@ out=out_data
[swash]
input=sws/INPUT.sws
swashrun=/opt/swash/swash/swashrun
out=out_swash
out=out
[post]
out=out_post

View file

@ -25,32 +25,45 @@ vel = read_nohead_vect(root.joinpath("vel.dat"), n)
n_t = botl.shape[0]
# Cospectral calculations
pos_x = n//10
pos_x = n // 10
f = 1 / 0.25
eta = (dep - botl)[n_t//2:, pos_x]
u = vel[n_t//2:, 0, pos_x]
eta = (dep - botl)[n_t // 2 :, pos_x]
u = vel[n_t // 2 :, 0, pos_x]
phi_eta = np.abs(sgl.csd(eta, eta)[1])
phi_u = np.abs(sgl.csd(u, u)[1])
phi_eta_u = np.abs(sgl.csd(eta, u)[1])
phi_eta = np.abs(sgl.csd(eta, eta, f))
phi_u = np.abs(sgl.csd(u, u, f))
phi_eta_u = np.abs(sgl.csd(eta, u, f))
R = np.sqrt((phi_eta + phi_u - 2*phi_eta_u)/(phi_eta + phi_u + 2*phi_eta_u))
R = np.sqrt(
(phi_eta[1] + phi_u[1] - 2 * phi_eta_u[1])
/ (phi_eta[1] + phi_u[1] + 2 * phi_eta_u[1])
)
# Plotting
fig, (ax_dep, ax_vel) = plt.subplots(2)
ax_dep.plot((dep - botl)[:, pos_x], label="dep", color="#0066ff")
ax_dep.set(xlabel="t (s)", ylabel="z (m)")
ax_dep.autoscale(axis="x", tight=True)
ax_dep.grid()
ax_vel.plot(vel[:, 0, pos_x], label="vel")
ax_vel.set(xlabel="t (s)", ylabel="U (m/s)")
ax_vel.autoscale(axis="x", tight=True)
ax_vel.grid()
fig.tight_layout()
fig_r, ax_r = plt.subplots()
ax_r.plot(R)
ax_r.set(ylim=(0,1))
ax_r.plot(phi_eta[0], R)
ax_r.autoscale(axis="x", tight=True)
ax_r.set(ylim=(0, 1), xlabel="f (Hz)", ylabel="R")
ax_r.grid()
plt.show(block=True)
out = pathlib.Path(config.get("post", "out"))
out.mkdir(exist_ok=True)
fig.savefig(out.joinpath(f"{pos_x}.png"))
fig_r.savefig(out.joinpath(f"R{pos_x}.png"))