In [3]:

    

%load_ext autoreload
%autoreload 2


    

    




The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [15]:

    

barriers = []
energies = []


    

In [16]:

    

loc = []
with open("outputs/output_1.dat","r") as f1:
    temp = f1.readline()
    for line in f1:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(1)
energies.append(loc)


    

In [17]:

    

loc = []
with open("outputs/output_2.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(2)
energies.append(loc)


    

In [18]:

    

loc = []
with open("outputs/output_3.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(3)
energies.append(loc)


    

In [19]:

    

loc = []
with open("outputs/output_4.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(4)
energies.append(loc)


    

In [20]:

    

loc = []
with open("outputs/output_5.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(5)
energies.append(loc)


    

In [21]:

    

loc = []
with open("outputs/output_6.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(6)
energies.append(loc)


    

In [22]:

    

loc = []
with open("outputs/output_7.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(7)
energies.append(loc)


    

In [23]:

    

loc = []
with open("outputs/output_8.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(8)
energies.append(loc)


    

In [24]:

    

loc = []
with open("outputs/output_9.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(9)
energies.append(loc)


    

In [25]:

    

loc = []
with open("outputs/output_10.dat","r") as f2:
    temp = f2.readline()
    for line in f2:
        x = float(line.strip().split()[0])
        if x < 100:
            loc.append(x)
barriers.append(10)
energies.append(loc)


    

In [26]:

    

import matplotlib.pyplot as plt


    

In [28]:

    

for i in range(10):
    for j in range(len(energies[i])):
        plt.scatter(barriers[i],energies[i][j])
plt.show()


    

In [ ]: