]
fig, axs = plt.subplots(3, 1, sharex=True)
fig.subplots_adjust(hspace=0)
for i in range(len(strtype)):
pathstring.append(
'../alsi/polymer/quenchedpoly/{:s}'.format(str(chargeratio[i])) +
'/{:s}/'.format(str(strtype[i])))
plotlines = []
for j in range(len(pathstring)):
ph = []
for p in sorted(
Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.list')):
print(p)
d = mu.monAnalysis(p, 'iotube', 'tt')
e = mu.monAnalysis(p, 'itube', 'tt')
f = mu.monAnalysis(p, 'otube', 'tt')
labelstring = 'pH = {:d}'.format(int(p.parts[-2][2]))
print(p.parts)
print(p.parts[-2])
print(p.parts[-2][2])
ph.append(p.parts[-2][2])
d[:, 1] = e[:, 1] + f[:, 1]
for k in range(len(d)):
d[k][2] = max(e[k][2], f[k][2])
]
linecolor2 = [
'powderblue', 'lightblue', 'lightskyblue', 'royalblue', 'blue',
'mediumblue', 'darkblue'
]
for i in range(len(strtype)):
pathstring.append(
'alsi/polymer/quenchedpoly/{:s}'.format(str(chargeratio[i])) +
'/{:s}/'.format(str(strtype[i])))
for j in range(len(pathstring)):
for p in sorted(
Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.list')):
print(p)
d = mu.monAnalysis(p, 'iotube', 'tt') #choose itube/otube/iotube
e = mu.monAnalysis(p, 'itube', 'tt')
f = mu.monAnalysis(p, 'otube', 'tt')
titlestring = '{:s}'.format(p.parts[-4])
labelstring = 'ph = {:d}'.format(int(p.parts[-2][2]))
print(p.parts)
print(p.parts[-2])
print(p.parts[-2][2])
ph.append(p.parts[-2][2])
#reetemp = d
plt.figure(1)
#plt.subplot(311)
plt.errorbar(d[:, 0],