# pairs = [(0.22, 0.0), (0.0, 0.22), (-0.22, 0.0), (0.0, -0.22), \
# (0.44, 0.0), (0.0, 0.44), (-0.44, 0.0), (0.0, -0.44)]
pairs = [(0.33, 0.0), (0.0, 0.33), (-0.33, 0.0), (0.0, -0.33)]
pl.figure()
for pair in pairs:
dataset = grps[pair].sort_index(0, "Z (m)")
pl.plot(dataset["Temperature (K)"], dataset["Z (m)"], label = str(pair))
index = pl.argmin(abs(means[k]["Time"] - (filling_levels[j] / massflows[k])))
pl.plot(means[k].loc[index][1:][0:].values, heights, label = "Mean")
h_kolbenfuss = hoehe_kolbenfuss.hoehe_kolbenfuss(filling_levels[j] / 1000.0)
h_speicher = 2.335
pl.plot([60.0, 60.0, 60.0 + dt[k], 60.0 + dt[k]], \
[0.0, h_speicher - h_kolbenfuss, h_speicher - h_kolbenfuss, h_speicher], \
label = "Plug-flow")
pl.xlim([58.0, 86.0])
pl.ylim([0.0, 2.4])
pl.grid(True)
pl.ylabel("Height (m)")
pl.xlabel("Temperature ($^{\circ}$C)")
pl.legend(loc="lower right")
pl.title("Temperature distribution in " + sim + \
os.chdir("..")
# -- fig_comp_plug_flow --
for level in filling_levels:
pl.figure()
for k, sim in enumerate(sims):
index = pl.argmin(abs(means[k]["Time"] - (level / massflows[k])))
pl.plot(means[k].loc[index][1:][0:].values, heights, label = sim)
h_kolbenfuss = hoehe_kolbenfuss.hoehe_kolbenfuss(level / 1000.0)
h_speicher = 2.335
pl.plot([60.0, 60.0, 60.0 + dt[0], 60.0 + dt[0]], \
[0.0, h_speicher - h_kolbenfuss, h_speicher - h_kolbenfuss, h_speicher], \
label = "Plug-flow ")
pl.xlim([58.0, 86.0])
pl.ylim([0.2, 2.4])
pl.grid(True)
pl.ylabel("Height (m)")
pl.xlabel("Temperature ($^{\circ}$C)")
pl.legend(loc="lower right")
pl.title("Mean temperature per stratum (filling level = " + str(level) + \