예제 #1
0
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")
예제 #2
0
파일: aci.py 프로젝트: nbigaouette/thesis
#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,