Skip to content
Snippets Groups Projects
Commit 367f790a authored by Pierre Augier's avatar Pierre Augier
Browse files

Fix cond 2023sed-parametric/tuto_sed_async.py

parent 1df28a12
No related branches found
No related tags found
1 merge request!100Topic/default/prepare 0.0.7
......@@ -59,6 +59,8 @@
return tau[:, 3], np.gradient(tau[:, 3], y), field.time
cond_height = y > 1.6 * bed_height
nb_saved_times = 1
tau, grad_tau, t_now = get_tau()
......@@ -71,7 +73,7 @@
ax1.set_xlabel("$t$")
ax1.set_ylabel("residuals p_rbgh")
ax0.plot(gradp * (y[-1] - y[y > 0.1]), y[y > 0.1], c="r")
ax0.plot(gradp * (y[-1] - y[cond_height]), y[cond_height], c="r")
(line,) = ax0.plot(tau, y)
ax0.set_title(f"t = {t_now}")
......@@ -101,7 +103,6 @@
ax0.set_title(f"t = {t_now}")
fig.canvas.draw()
cond_height = (y > 1.1 * bed_height) & (y < 1.5 * bed_height)
percentage = 100 * abs(grad_tau[cond_height] + gradp).mean() / gradp
execution_time = timedelta(seconds=time() - time_start)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment