gammas = ['0.02','0.06','0.08','0.1','0.12']
for js,gamma in enumerate(gammas):
scan = {}
scan['\\gamma_s']=gamma
scan['k_{24}']=k24
scan['\\omega']=str(omega)
obsfwd = ObservableData(["\\Lambda"],scan_dir='scanforward',scan=scan,loadsuf=loadsuf,
savesuf=savesuf)
if not obsfwd.file_exists:
continue
obsfwd.sort_observables()
obsfwd.remove_duplicates()
Lambdas = obsfwd.data[:,0]
#obsbkwd = ObservableData(["\\Lambda"],scan_dir='scanbackward',scan=scan,loadsuf=loadsuf,
# savesuf=savesuf)
#obsbkwd.sort_observables()
ysfwd = [obsfwd.E(),obsfwd.R(),obsfwd.eta(),obsfwd.delta(),obsfwd.surfacetwist()]
#ysbkwd = [obsbkwd.E(),obsbkwd.R(),obsbkwd.eta(),obsbkwd.delta(),obsbkwd.surfacetwist()]
for i,observable in enumerate(observable_list):
scan_dir='scanforward',
scan=scan,
loadsuf=loadsuf,
savesuf=savesuf)
Lambdas = obsfwd.data[:, 0]
ysfwd = [obsfwd.R(), obsfwd.surfacetwist()]
if js == 1:
obsbkwd = ObservableData(["\\Lambda"],
scan_dir='scanbackward',
scan=scan,
loadsuf=loadsuf,
savesuf=savesuf)
obsbkwd.sort_observables()
ysbkwd = [obsbkwd.R(), obsbkwd.surfacetwist()]
for i, observable in enumerate(observable_list):
if observable == 'surfacetwist':
ylabel = r'$\psi(R)$'
else:
ylabel = fr'${observable}$'
if js == 0:
ax[observable].plot(Lambdas[1:],
ysfwd[i][1:],
'-',
color=colors[js],
label=rf"$(\gamma,k_{{24}})=({gamma},{k24})$",