a = np.loadtxt("#distance.xvg.1#")
plt.figure(figsize=(10,10))
plt.subplot(211)
x=a[:,1]
n, bins, patches = plt.hist(x, 50, facecolor='m', alpha=0.75)
plt.xlabel('Distance')
plt.subplot(212)
x=a[:,0]
y=a[:,1]
plt.scatter(x, y)
plt.xlabel('time, ps')
plt.ylabel('dist')
plt.savefig('plot_vr.png')
plt.show()
a = np.loadtxt("distance.xvg")
plt.figure(figsize=(10,10))
plt.subplot(211)
x=a[:,1]
n, bins, patches = plt.hist(x, 50, facecolor='m', alpha=0.75)
plt.xlabel('Distance')
plt.subplot(212)
x=a[:,0]
y=a[:,1]
plt.scatter(x, y)
plt.xlabel('time, ps')
plt.ylabel('dist')
plt.savefig('plot_be.png')
plt.show()
a = np.loadtxt("#distance.xvg.2#")
plt.figure(figsize=(10,10))
plt.subplot(211)
x=a[:,1]
n, bins, patches = plt.hist(x, 50, facecolor='m', alpha=0.75)
plt.xlabel('Distance')
plt.subplot(212)
x=a[:,0]
y=a[:,1]
plt.scatter(x, y)
plt.xlabel('time, ps')
plt.ylabel('dist')
plt.savefig('plot_nh.png')
plt.show()
a = np.loadtxt("#distance.xvg.3#")
plt.figure(figsize=(10,10))
plt.subplot(211)
x=a[:,1]
n, bins, patches = plt.hist(x, 50, facecolor='m', alpha=0.75)
plt.xlabel('Distance')
plt.subplot(212)
x=a[:,0]
y=a[:,1]
plt.scatter(x, y)
plt.xlabel('time, ps')
plt.ylabel('dist')
plt.savefig('plot_an.png')
plt.show()
a = np.loadtxt("#distance.xvg.4#")
plt.figure(figsize=(10,10))
plt.subplot(211)
x=a[:,1]
n, bins, patches = plt.hist(x, 50, facecolor='m', alpha=0.75)
plt.xlabel('Distance')
plt.subplot(212)
x=a[:,0]
y=a[:,1]
plt.scatter(x, y)
plt.xlabel('time, ps')
plt.ylabel('dist')
plt.savefig('plot_sd.png')
plt.show()