change the profile script

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