diff --git a/swash/processing/post.py b/swash/processing/post.py index 5194ca5..f3ed9e9 100644 --- a/swash/processing/post.py +++ b/swash/processing/post.py @@ -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)