Beispiel #1
0
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)
Beispiel #2
0
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)