Commit 160a21a5 authored by Mathieu RASSON's avatar Mathieu RASSON

Figures multi

parent 7731b49a
......@@ -81,7 +81,7 @@
"#tp = 1\n",
"Pcr = 1.8e9\n",
"#Pin = 5 * Pcr\n",
"Pin = 2*Pcr\n",
"Pin = 6*Pcr\n",
"\n",
"# Dispersion\n",
"k = 2 * np.pi / lambd\n",
......@@ -143,13 +143,13 @@
"metadata": {},
"outputs": [],
"source": [
"delta_z = 3e-4\n",
"z_max = 4\n",
"delta_z = 2e-4\n",
"z_max = 3\n",
"n_z = int(z_max/delta_z)\n",
"\n",
"crank = CrankNicolsonFilamentation()\n",
"\n",
"crank.set_grid(r_max=4e-3, n_r=300, z_min=0, z_max=z_max, n_z=n_z, t_min=-4e-13, t_max=4e-13, n_t=600)\n",
"crank.set_grid(r_max=3e-3, n_r=225, z_min=0, z_max=z_max, n_z=n_z, t_min=-3e-13, t_max=3e-13, n_t=450)\n",
"crank.set_parameters(Dr=diff_coeff, Dt=disp_coeff, V=potential, ofi=ofi, ava=ava ,K=K, sigma=sigma, \\\n",
" omega0=omega0, tau=tau, f=kerr_MPA)\n",
"\n",
......@@ -366,13 +366,13 @@
"# Field squared at given z,t along r\n",
"\n",
"plt.figure()\n",
"p1, = plt.plot(1000*np.real(crank.r_pts), 1/(c*mu)*np.abs(crank.E_matrix[7800,:,295])**2)\n",
"p2, = plt.plot(1000*np.real(crank.r_pts), 1/(c*mu)*np.abs(crank.E_matrix[7800,:,276])**2)\n",
"p1, = plt.plot(1000*np.real(crank.r_pts), 1/(c*mu)*np.abs(crank.E_matrix[9000,:,232])**2)\n",
"p2, = plt.plot(1000*np.real(crank.r_pts), 1/(c*mu)*np.abs(crank.E_matrix[9000,:,271])**2)\n",
"plt.xlim(-0.01, 1)\n",
"plt.xlabel('Radius (mm)')\n",
"plt.ylabel('Intensity ($W \\cdot m^{-2}$)')\n",
"plt.legend((p1, p2), ('Between pulses', 'On pulse'))\n",
"#plt.savefig('field_radius_2P.png', dpi=300)\n",
"#plt.savefig('field_radius3rd_6P.png', dpi=300)\n",
"plt.show()"
]
},
......@@ -386,14 +386,14 @@
"\n",
"plt.figure()\n",
"#plt.plot(np.log(np.abs(E[8300,0,:])**2))\n",
"p1, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[7500,0,:])**2)\n",
"p2, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[7900,0,:])**2)\n",
"p3, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[8300,0,:])**2)\n",
"plt.xlim(-100, 100)\n",
"p1, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[5000,0,:])**2)\n",
"p2, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[7500,0,:])**2)\n",
"p3, = plt.plot(1e15*np.real(crank.t_pts), 1/(c*mu)*np.abs(crank.E_matrix[9000,0,:])**2)\n",
"#plt.xlim(50, 75)\n",
"plt.xlabel('Time (fs)')\n",
"plt.ylabel('Intensity ($W \\cdot m^{-2}$)')\n",
"plt.legend((p1, p2, p3), ('First maximum', 'Minimum', 'Second maximum'))\n",
"#plt.savefig('field_time_2P.png', dpi=300)\n",
"plt.legend((p1, p2, p3), ('First filament', 'Second filament', 'Third filament'))\n",
"#plt.savefig('field_time_6P.png', dpi=300)\n",
"plt.show()"
]
},
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment