# 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) + \