Skip to content
Snippets Groups Projects
Commit 036ff2425c1a authored by vlabarre's avatar vlabarre
Browse files

remove script save_ratio_EA_EK_vs_Fh.py

parent 2177ed75cf53
No related branches found
No related tags found
1 merge request!62022strat_turb_polo!
import numpy as np
import matplotlib.pyplot as plt
from util import save_fig, plot, Fh_limit
from util_dataframe import df
tmp = df[df.Re > 2000].copy()
tmp["E_A/E_K"] = tmp.EA / (tmp.EKh)
ax = plot(
tmp, "Fh", "E_A/E_K", c=np.log10(tmp["R2"]), vmin=0.5, vmax=2, logy=True
)
ax.set_xlabel("$F_h$", fontsize=12)
ax.set_ylabel(r"$E_A / E_K$", fontsize=12)
ax.axvline(Fh_limit, linestyle=":")
ax.axvline(1.0, linestyle=":")
xs = np.linspace(1.5e-1, 5e-1, 4)
ax.plot(xs, 5e-2 * xs**-1)
ax.text(0.16, 0.12, "${F_h}^{-1}$")
xs = np.linspace(6e-1, 4, 4)
ax.plot(xs, 5e-2 * xs**-2)
ax.text(1.2, 0.05, "${F_h}^{-2}$")
fig = ax.figure
fig.text(0.84, 0.07, r"$\log_{10}(\mathcal{R}_2)$", fontsize=12)
fig.tight_layout()
save_fig(fig, "fig_ratio_EA_EK_vs_Fh.png")
if __name__ == "__main__":
plt.show()
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