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

2022strat_turb_toro/py/save_spectra_1R.py: plot versus k/kb

parent adf9c5f6
No related branches found
No related tags found
No related merge requests found
......@@ -28,9 +28,9 @@
data = sim.output.spectra.load1d_mean(tmin)
k = data["kx"]
line, = ax0.plot(
k,
(line,) = ax0.plot(
k / kb,
data["spectra_E_kx"] * k ** (5 / 3),
label=rf"$F_h={mean_values['Fh']:.2f}$",
)
......@@ -33,7 +33,7 @@
data["spectra_E_kx"] * k ** (5 / 3),
label=rf"$F_h={mean_values['Fh']:.2f}$",
)
ax0.axvline(kb, linestyle=":", color=line.get_color())
# ax0.axvline(kb, linestyle=":", color=line.get_color())
k = data["kz"]
......@@ -38,5 +38,5 @@
k = data["kz"]
ax1.plot(k, data["spectra_E_kz"] * k ** (5 / 3))
ax1.axvline(kb, linestyle=":", color=line.get_color())
ax1.plot(k / kb, data["spectra_E_kz"] * k ** (5 / 3))
# ax1.axvline(kb, linestyle=":", color=line.get_color())
......@@ -42,6 +42,6 @@
ax0.set_xlabel("$k_h$")
ax1.set_xlabel("$k_z$")
ax0.set_xlabel("$k_h / k_b$")
ax1.set_xlabel("$k_z / k_b$")
ax0.set_ylabel("$E_K(k_h) {k_h}^{5/3}$")
ax1.set_ylabel("$E_K(k_z) {k_z}^{5/3}$")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment