1
Fork 0

Switched dep to watl

This commit is contained in:
Edgar P. Burkhart 2022-03-04 10:52:28 +01:00
parent 071fde81b5
commit 1be7a5285b
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227

View file

@ -36,7 +36,7 @@ dt = config.getfloat("post", "dt")
f = 1 / dt
log.info(f"Computing reflection coefficient at x={x0}")
eta = data["dep"][t > t0, x0] - data["botl"][x0]
eta = data["watl"][t > t0, x0]
u = data["vel"][t > t0, 0, x0]
phi_eta = np.abs(sgl.csd(eta, eta, f))
@ -50,13 +50,13 @@ R = np.sqrt(
# Plotting
log.info("Plotting results")
fig, (ax_dep, ax_vel) = plt.subplots(2)
fig, (ax_watl, ax_vel) = plt.subplots(2)
ax_dep.plot(t, data["dep"][:, x0] - data["botl"][x0], label="dep")
ax_dep.set(xlabel="t (s)", ylabel="z (m)")
ax_dep.autoscale(axis="x", tight=True)
ax_dep.grid()
ax_dep.axvline(t0, c="k", alpha=0.2)
ax_watl.plot(t, data["watl"][:, x0], label="watl")
ax_watl.set(xlabel="t (s)", ylabel="z (m)")
ax_watl.autoscale(axis="x", tight=True)
ax_watl.grid()
ax_watl.axvline(t0, c="k", alpha=0.2)
ax_vel.plot(t, data["vel"][:, 0, x0], label="vel")
ax_vel.set(xlabel="t (s)", ylabel="U (m/s)")
@ -75,7 +75,7 @@ ax_r.grid()
fig_x, ax_x = plt.subplots()
ax_x.plot(-data["botl"], color="k")
ax_x.plot(data["dep"][np.argmin(np.abs(t - t0)), :] - data["botl"])
ax_x.plot(data["watl"][np.argmin(np.abs(t - t0)), :])
ax_x.axvline(x0, c="k", alpha=0.2)
ax_x.set(xlabel="x (m)", ylabel="z (m)")
ax_x.autoscale(axis="x", tight=True)