# HG changeset patch # User vlabarre <Vincent.Labarre@oca.eu> # Date 1667469187 -3600 # Thu Nov 03 10:53:07 2022 +0100 # Node ID 036ff2425c1af43a11854b26bab302a75937f0b2 # Parent 2177ed75cf537c804372f96abd2232c2a13cf703 remove script save_ratio_EA_EK_vs_Fh.py diff --git a/2022strat_turb_polo/py/save_ratio_EA_EK_vs_Fh.py b/2022strat_turb_polo/py/save_ratio_EA_EK_vs_Fh.py deleted file mode 100644 --- a/2022strat_turb_polo/py/save_ratio_EA_EK_vs_Fh.py +++ /dev/null @@ -1,37 +0,0 @@ -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()