r = np.linspace(-20.0, 20.0, 1000) cs = 1.0 U = -cs / abs(r) U[U < Umin] = np.NaN laser = -r / 30.0 Ubent = U + laser lr = r[-1] - r[0] lU = Umax - Umin r0 = abs(cs / Xe_Z0_Ip) fig = plot.figure() ax1, ax1b = plot.get_axis_two_scales(fig, scale_y=cst.Eh_to_eV, ax2_ylabel="Energy [eV]") # Ip ax1.plot([-r0, r0], [-Xe_Z0_Ip, -Xe_Z0_Ip], "-m") # Threshold ax1.axhline(0.0, color="k", ls="-", alpha=0.5) # U(r) ax1.plot(r, U, "--k", alpha=0.5, label="Unperturbed ion") # Laser ax1.plot(r, laser, "-r", label="Laser") # U(r) + laser ax1.plot(r, Ubent, "-k", label="Effective")
#ax1, ax1b = plot.get_axis_two_scales(fig1, #scale_y = cst.Eh_to_eV, #ax2_ylabel = 'Energy [eV]') #ax2, ax2b = plot.get_axis_two_scales(fig2, #scale_y = cst.Eh_to_eV, #ax2_ylabel = 'Energy [eV]') #ax3, ax3b = plot.get_axis_two_scales(fig3, #scale_y = cst.Eh_to_eV, #ax2_ylabel = 'Energy [eV]') #ax4, ax4b = plot.get_axis_two_scales(fig4, #scale_y = cst.Eh_to_eV, #ax2_ylabel = 'Energy [eV]') fig = plot.figure() ax1, ax1b = plot.get_axis_two_scales(fig, scale_y = cst.Eh_to_eV, ax2_ylabel = 'Energy [eV]', subplot = 221) ax2, ax2b = plot.get_axis_two_scales(fig, scale_y = cst.Eh_to_eV, ax2_ylabel = 'Energy [eV]', sharex = ax1, sharey = ax1, subplot = 222) ax3, ax3b = plot.get_axis_two_scales(fig, scale_y = cst.Eh_to_eV, ax2_ylabel = 'Energy [eV]', sharex = ax1, sharey = ax1, subplot = 223) ax4, ax4b = plot.get_axis_two_scales(fig, scale_y = cst.Eh_to_eV, ax2_ylabel = 'Energy [eV]', sharex = ax1, sharey = ax1,