diff --git a/PlotProfiles_for_paper.py b/PlotProfiles_for_paper.py index 706a0d9600b7d6ef8fc8a1c47dfb79af985df394..d44daacf4a00a5c59f75c3a3609dad03e1526623 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