|
| 1 | +################################ |
| 2 | +####### GENERAL PARAMETERS ###### |
| 3 | +################################# |
| 4 | +max_step = 250 |
| 5 | +amr.n_cell = 32 32 32 |
| 6 | +amr.max_grid_size = 32 |
| 7 | +amr.blocking_factor = 32 |
| 8 | +geometry.dims = 3 |
| 9 | +geometry.prob_lo = -4.e-6 -4.e-6 -4.e-6 |
| 10 | +geometry.prob_hi = 4.e-6 4.e-6 4.e-6 |
| 11 | +amr.max_level = 0 |
| 12 | +boundary.field_lo = periodic periodic periodic |
| 13 | +boundary.field_hi = periodic periodic periodic |
| 14 | + |
| 15 | +################################# |
| 16 | +############ NUMERICS ########### |
| 17 | +################################# |
| 18 | +warpx.verbose = 1 |
| 19 | +warpx.use_filter = 0 |
| 20 | +warpx.cfl = 0.9 |
| 21 | + |
| 22 | +algo.em_solver_medium = macroscopic # vacuum/macroscopic |
| 23 | + |
| 24 | +algo.macroscopic_sigma_method = laxwendroff # laxwendroff or backwardeuler |
| 25 | + |
| 26 | +macroscopic.sigma_function(x,y,z) = "sigma_0" |
| 27 | + |
| 28 | +macroscopic.epsilon_function(x,y,z) = "epsilon_0" |
| 29 | + |
| 30 | +macroscopic.mu_function(x,y,z) = "mu_0 " |
| 31 | + |
| 32 | +################################# |
| 33 | +############ FIELDS ############# |
| 34 | +################################# |
| 35 | + |
| 36 | +my_constants.pi = 3.14159265359 |
| 37 | +my_constants.L = 50.e-6 |
| 38 | +my_constants.c = 299792458. |
| 39 | +my_constants.wavelength = 4.e-6 |
| 40 | + |
| 41 | +my_constants.sigma_0 = 0.0 |
| 42 | +my_constants.epsilon_0 = 8.8541878128e-12 |
| 43 | +my_constants.mu_0 = 1.25663706212e-06 |
| 44 | + |
| 45 | +warpx.E_ext_grid_init_style = parse_E_ext_grid_function |
| 46 | +warpx.Ex_external_grid_function(x,y,z) = "1.e5*exp(-(y)**2/(L/32)**2)*cos(2*pi*(y)/wavelength)" |
| 47 | +warpx.Ey_external_grid_function(x,y,z) = "0." |
| 48 | +warpx.Ez_external_grid_function(x,y,z) = "0." |
| 49 | + |
| 50 | +warpx.B_ext_grid_init_style = parse_B_ext_grid_function |
| 51 | +warpx.Bx_external_grid_function(x,y,z) = "0." |
| 52 | +warpx.By_external_grid_function(x,y,z) = "0." |
| 53 | +warpx.Bz_external_grid_function(x,y,z) = "(-1.e5*exp(-(y)**2/(L/32)**2)*cos(2*pi*(y)/wavelength)/c)" |
| 54 | + |
| 55 | +warpx.E_excitation_on_grid_style = parse_E_excitation_grid_function |
| 56 | + |
| 57 | +warpx.Ex_excitation_flag_function(x,y,z) = "0." |
| 58 | +warpx.Ey_excitation_flag_function(x,y,z) = "0." |
| 59 | +warpx.Ez_excitation_flag_function(x,y,z) = "0." |
| 60 | + |
| 61 | +warpx.Ex_excitation_grid_function(x,y,z,t) = "0." |
| 62 | +warpx.Ey_excitation_grid_function(x,y,z,t) = "0." |
| 63 | +warpx.Ez_excitation_grid_function(x,y,z,t) = "0." |
| 64 | + |
| 65 | +warpx.B_excitation_on_grid_style = parse_B_excitation_grid_function |
| 66 | + |
| 67 | +warpx.Bx_excitation_flag_function(x,y,z) = "0." |
| 68 | +warpx.By_excitation_flag_function(x,y,z) = "0." |
| 69 | +warpx.Bz_excitation_flag_function(x,y,z) = "0." |
| 70 | + |
| 71 | +warpx.Bx_excitation_grid_function(x,y,z,t) = "0." |
| 72 | +warpx.By_excitation_grid_function(x,y,z,t) = "0." |
| 73 | +warpx.Bz_excitation_grid_function(x,y,z,t) = "0." |
| 74 | + |
| 75 | +# Diagnostics |
| 76 | +diagnostics.diags_names = plt |
| 77 | +plt.intervals = 1 |
| 78 | +plt.fields_to_plot = Ex Ey Ez Bx By Bz |
| 79 | +plt.diag_type = Full |
| 80 | +plt.plot_raw_fields = 1 |
0 commit comments