omega = sys.argv[2]
scan = {}
scan['\\Lambda']=Lambda
scan['\\omega']=omega
# load in 2d grid of data in data2d for each observable at the
# specified gamma,k24 pair.
for i,k24 in enumerate(k24s):
scan['k_{24}']=str(k24)
obsfwd = ObservableData(["\\gamma_s"],scan_dir='scanforward',scan=scan,
loadsuf=loadsuf,savesuf=savesuf)
obsbkwd = ObservableData(["\\gamma_s"],scan_dir='scanbackward',scan=scan,
loadsuf=loadsuf,savesuf=savesuf)
bkwdexists = os.path.isfile(obsbkwd.observables_fname())
if not bkwdexists:
# save a copy of scanforward data under the scanbackward filename,
# as if the file doesn't exist it's assumed that the scanbackward
# and scanforward files are the same (below the critical point region
# in gamma,k24 space)
np.savetxt(obsbkwd.observables_fname(),obsfwd.data,
fmt='\t'.join(["%13.6e"] + ["%15.8e"]*5))
scan = {}
scan['\\Lambda']=sys.argv[1]
scan['\\omega']=sys.argv[2]
# load in 2d grid of data in data2d for each observable at the
# specified gamma,k24 pair.
for k24 in k24s:
scan['k_{24}']=str(k24)
obsfwd = ObservableData(["\\gamma_s"],scan_dir='scanforward',scan=scan,
loadsuf=loadsuf,savesuf=savesuf)
obsbkwd = ObservableData(["\\gamma_s"],scan_dir='scanbackward',scan=scan,
loadsuf=loadsuf,savesuf=savesuf)
obsbkwd.sort_observables()
end_bkwd = obsbkwd.data.shape[0]
datanew = np.concatenate((obsbkwd.data,obsfwd.data[end_bkwd:,:]),axis=0)
np.savetxt(obsbkwd.observables_fname(),datanew,
fmt='\t'.join(["%13.6e"]+["%15.8e"]*5))