vertical lines

NCAR · Aug 4, 2023 · 7bc1822 · 7bc1822
1 parent e9df05f
commit 7bc1822
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/MPAS_true_failure_testing/post_ect_script.py b/MPAS_true_failure_testing/post_ect_script.py
@@ -113,6 +113,7 @@ def main(argv):
         plt.xscale("symlog")
         plt.plot(perturbations, avg_pca_fails/PCA_dims)
         plt.xticks(perturbations)
+        plt.axvline(x=0, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT PCA Fails vs\n {var_name} Perturbation")
         plt.ylabel("PCA Fail Percent")
         plt.xlabel("Perturbation Factor")
@@ -123,6 +124,7 @@ def main(argv):
         plt.xscale("symlog")
         plt.plot(test_vals, avg_pca_fails/PCA_dims)
         plt.xticks(test_vals)
+        plt.axvline(x=default_var_value, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT PCA Fails vs\n {var_name} Value")
         plt.ylabel("PCA Fail Percent")
         plt.xlabel(f"{var_name}")
@@ -133,6 +135,7 @@ def main(argv):
         plt.xscale("linear")
         plt.plot(test_vals, avg_pca_fails/PCA_dims)
         plt.xticks(test_vals)
+        plt.axvline(x=default_var_value, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT PCA Fails vs\n {var_name} Value")
         plt.ylabel("PCA Fail Percent")
         plt.xlabel(f"{var_name}")
@@ -143,6 +146,7 @@ def main(argv):
         plt.xscale("symlog")
         plt.plot(perturbations, avg_eet_fails)
         plt.xticks(perturbations)
+        plt.axvline(x=0, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT EET Rate vs\n {var_name} Perturbation")
         plt.ylabel("EET Fail Percent")
         plt.xlabel(f"Perturbation Factor")
@@ -153,6 +157,7 @@ def main(argv):
         plt.xscale("symlog")
         plt.plot(test_vals, avg_eet_fails)
         plt.xticks(test_vals)
+        plt.axvline(x=default_var_value, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT EET Rate vs\n {var_name} Value")
         plt.ylabel("EET Fail Percent")
         plt.xlabel(f"{var_name}")
@@ -163,6 +168,7 @@ def main(argv):
         plt.xscale("linear")
         plt.plot(test_vals, avg_eet_fails)
         plt.xticks(test_vals)
+        plt.axvline(x=default_var_value, color='r', label='axvline - full height')
         plt.title(f"MPAS UF-ECT EET Rate vs\n {var_name} Value")
         plt.ylabel("EET Fail Percent")
         plt.xlabel(f"{var_name}")