From daa121201bcd04266932848ba996d908d5f4eae5 Mon Sep 17 00:00:00 2001
From: weisse02 <andreas.weissenbrunner@ptb.de>
Date: Mon, 20 Mar 2023 01:34:37 +0100
Subject: [PATCH] change the profile script
---
PlotProfiles_for_paper.py | 32 ++++++++++++++++++++++++--------
1 file changed, 24 insertions(+), 8 deletions(-)
diff --git a/PlotProfiles_for_paper.py b/PlotProfiles_for_paper.py
index 706a0d9..d44daac 100644
--- a/PlotProfiles_for_paper.py
+++ b/PlotProfiles_for_paper.py
@@ -75,19 +75,35 @@ def get_fully_SA(r):
# --------------------------------------------------------------
+plt.rcParams['figure.figsize'] = [5,5]
-
-flow = Elbow_profile('Double_S_Elbow')
+flow = Elbow_profile('SingleElbow')
+if flow.case == 'DoubleElbow':
+ case = "DE"
+elif flow.case == "Double_S_Elbow":
+ case = "DSE"
+elif flow.case == "SingleElbow":
+ case = "SE"
+#
ufully = get_fully_SA(flow.r.ravel())
-Rks = [0.51,1.67375,3.38375]
-rk = 1.17604257
-dls = [0,5,10]
+# Rks = [0.51,1.5981882,5.51] #for the DE and DSE cases
+Rks = [0.51,1.10375, 1.67375, 5.165, 10.01]
+rk = 1.5981882
+dls = [0.01,5,10]
+dl = 0.01
dist = 6
-for dl in dls:
+for rk in Rks:
u = flow.get_profile(Rk=rk,dl=dl,dist=dist,addfully=False)
u = u.reshape(3,-1)
u[2] += ufully
- nicecontour(flow.x,flow.y,u,minval=0.48,maxval=1.13,scalefac=0.3)
+ nicecontour(flow.x,flow.y,u,minval=0.48,maxval=1.13,scalefac=0.3,withcolorbar=False)
+
+ plt.xlabel("x")
+ plt.ylabel("y")
+ #plt.title(case + " $R_c$= "+ "{:.2f}".format(rk) + " $d_l$= " + "{:.2f}".format(dl) )
+ plt.title(case + " $R_c$= "+ "{:.2f}".format(rk))
plt.tight_layout()
- plt.savefig("C:/Users/weisse02/PTBbox/VirtMet/Auswertung_paper/Profile_plots/" + flow.case + "_rk_" + str(rk) + "_dist_" + str(dist) + "_dl_" + str(dl) + ".pdf")
\ No newline at end of file
+ plt.savefig("C:/Users/weisse02/PTBbox/VirtMet/Auswertung_paper/Profile_plots/" + flow.case + "_rk_" + str(rk) + "_dist_" + str(dist) + "_dl_" + str(dl) + ".pdf")
+#
+#
\ No newline at end of file
--
GitLab