def study_height_synopsis(in_file):
#in_file = '/space/tohekorh/Spiral/bender_input/calc/consts/in_8-16_H=73.0.txt'
params = parse_input(in_file)
ir = params[("ir")]
rmin,rmax = params[("rmin")],params[("rmax")]
data_array, syn_dir = read_synopsis(in_file) # heights, E_bs, E_ss, consts
heights_hb, energies_hb = read_energies(ir)[:2]
plot_energies(ir, rmin, rmax, data_array[("heights")], data_array[("E_bs")], data_array[("E_ss")], \
heights_hb, energies_hb, syn_dir)
plot_amplitude(ir, rmin, rmax, data_array[("heights")], data_array[("consts")], syn_dir)
energies = []
for opt in opts:
data = read_h_consts(ir, opt)
all_data.append(data)
heights.append(data[('height')])
energies.append(data[('energies')])
heights_hb, energies_hb = read_energies(ir)[:2]
H = heights[0][ih]
hangle = H/2/pi
#energies = deform_energies(ue)
plot_energies(ir, heights, energies, heights_hb, energies_hb)
for all_data_sys in all_data:
for ida, data in enumerate(all_data_sys):
if ida == n or n == None:
H = data[0]
consts = data[1]
[nr, rmin, rmax] = data[3]
[nphi, phimin, phimax] = data[4]
hangle = H/2/pi
rs = np.linspace(rmin, rmax, nr)
phis = np.linspace(phimin, phimax, nphi)