if (len(aperture_names) == len(aperture_values)): nb_apertures = len(aperture_names) else: raise ValueError("aperture_values and aperture_names do not have the same length") # We generate the values for the distances delta_distance = (d_max - d_min) / float(nb_distance-1) distances = np.array([d_min + i * delta_distance for i in range(nb_distance)]) #~ # We generate the values for the focal #~ delta_focal = (focal_max - focal_min) / float(nb_focals-1) #~ focals = np.array([focal_min + i * delta_focal for i in range(nb_focals)]) # We prepare the plot fig = pl.figure(1) nb_subplot = autiwa.get_subplot_shape(len(focals)) fig.subplots_adjust(left=0.12, bottom=0.1, right=0.96, top=0.95, wspace=0.26, hspace=0.26) #~ myxfmt = ScalarFormatter(useOffset=True) #~ myxfmt._set_offset(1e5) #~ myxfmt.set_scientific(True) #~ myxfmt.set_powerlimits((-3, 3)) myxfmt = FormatStrFormatter('%.0f') myminorxfmt = FormatStrFormatter('%.1f') # We generate a list of colors colors = [ '#'+li for li in autiwa.colorList(len(dofs))] for focal_length in focals:
phi[too_low] = phi[too_low] + 360. phi[too_high] = phi[too_high] - 360. delta_longitude = long_of_peri_outer - long_of_peri_inner delta_longitude = delta_longitude%(360.) too_low = delta_longitude < ANGLE_MIN too_high = delta_longitude > ANGLE_MAX delta_longitude[too_low] = delta_longitude[too_low] + 360. delta_longitude[too_high] = delta_longitude[too_high] - 360. # We chose the number of plot in x and y axis for the p.multi # environment in order to plot ALL the resonant angles and have x and # y numbers as close on from another as possible. There are q+1 # resonant angles, the period ratio and w1 - w2, i.e q+3 plots (nb_lines, nb_rows) = autiwa.get_subplot_shape(q+3) # on trace les plots myxfmt = ScalarFormatter(useOffset=True) myxfmt._set_offset(1e5) myxfmt.set_scientific(True) myxfmt.set_powerlimits((-3, 3)) sys.stdout.write("Displaying %5.1f %% \r" % (75.)) sys.stdout.flush() fig = pl.figure(1) pl.clf() fig.subplots_adjust(left=0.12, bottom=0.1, right=0.96, top=0.95, wspace=0.26, hspace=0.26) # We create a variable to store the index of the subplot