khiA1_mixed_ion_kinetic_electron.py File Reference

Namespaces
namespace	khiA1_mixed_ion_kinetic_electron

Functions
def	khiA1_mixed_ion_kinetic_electron.str2float (string_array)

Variables
str	khiA1_mixed_ion_kinetic_electron.torder = 'SSPRK3'
str	khiA1_mixed_ion_kinetic_electron.sorder = 'thirdOrder'
float	khiA1_mixed_ion_kinetic_electron.cfl = 0.209
float	khiA1_mixed_ion_kinetic_electron.gamma = 2.0
str	khiA1_mixed_ion_kinetic_electron.numerical_flux_type = 'Rusanov'
int	khiA1_mixed_ion_kinetic_electron.nout = 100
float	khiA1_mixed_ion_kinetic_electron.target_tws = 60.0
float	khiA1_mixed_ion_kinetic_electron.dt_sim = 0.01
float	khiA1_mixed_ion_kinetic_electron.omega_s_over_omega_p = 0.2625
int	khiA1_mixed_ion_kinetic_electron.diag_writer_interval = 1e-2
bool	khiA1_mixed_ion_kinetic_electron.start_from_previous = False
	if starting from previous sim
warpy	khiA1_mixed_ion_kinetic_electron.configs = warpy.load_configs()
str	khiA1_mixed_ion_kinetic_electron.orig_loc = '/Users/imandatta/Documents/warpxm_runs/hyak_mox_outputs/hybrid/2d2v/khi/2021_02_19/'
str	khiA1_mixed_ion_kinetic_electron.orig_sim = 'khiA1_hybrid_pert_3rd_tws60_32x8x32x32_mixed'
warpy	khiA1_mixed_ion_kinetic_electron.orig_files = configs['test_rundir'](orig_loc, orig_sim)
int	khiA1_mixed_ion_kinetic_electron.frame_start_from_previous = 10
os	khiA1_mixed_ion_kinetic_electron.source_file = os.path.join(orig_files,'data', '%s_%d.h5'%(orig_sim,frame_start_from_previous))
list	khiA1_mixed_ion_kinetic_electron.tf_i_sds = ['left', 'right']
list	khiA1_mixed_ion_kinetic_electron.tf_i_vds = ['middle']
None	khiA1_mixed_ion_kinetic_electron.tf_e_sds = None
None	khiA1_mixed_ion_kinetic_electron.tf_e_vds = None
list	khiA1_mixed_ion_kinetic_electron.kinetic_i_sds = ['middle']
list	khiA1_mixed_ion_kinetic_electron.kinetic_i_vds = ['left', 'right']
list	khiA1_mixed_ion_kinetic_electron.kinetic_e_sds = ['left', 'middle', 'right']
None	khiA1_mixed_ion_kinetic_electron.kinetic_e_vds = None
float	khiA1_mixed_ion_kinetic_electron.qe = 1.602177e-19
float	khiA1_mixed_ion_kinetic_electron.mp = 1.672622e-27
float	khiA1_mixed_ion_kinetic_electron.mu0 = 1.256637e-6
float	khiA1_mixed_ion_kinetic_electron.eps0 = 8.854187e-12
int	khiA1_mixed_ion_kinetic_electron.c_light = 1/(mu0*eps0)**.5
float	khiA1_mixed_ion_kinetic_electron.omega_p_tau = 1.00
float	khiA1_mixed_ion_kinetic_electron.omega_c_tau = 1.00
float	khiA1_mixed_ion_kinetic_electron.skin_depth_norm = 1.0 / (omega_c_tau)
float	khiA1_mixed_ion_kinetic_electron.n0 = 2.00e20
math	khiA1_mixed_ion_kinetic_electron.omega_p = math.sqrt((qe**2 * n0)/(eps0 * mp))
float	khiA1_mixed_ion_kinetic_electron.tau = omega_p_tau / omega_p
float	khiA1_mixed_ion_kinetic_electron.omega_c = omega_c_tau / tau
float	khiA1_mixed_ion_kinetic_electron.B0 = omega_c * mp/qe
float	khiA1_mixed_ion_kinetic_electron.v0 = B0/math.sqrt(mp*n0*mu0)
float	khiA1_mixed_ion_kinetic_electron.L = v0*tau
float	khiA1_mixed_ion_kinetic_electron.p0 = B0**2/mu0
float	khiA1_mixed_ion_kinetic_electron.T0 = p0/n0
float	khiA1_mixed_ion_kinetic_electron.T0_ev = T0 / qe
math	khiA1_mixed_ion_kinetic_electron.v_th0 = math.sqrt(T0/mp)
math	khiA1_mixed_ion_kinetic_electron.rl0 = v_th0/omega_c
int	khiA1_mixed_ion_kinetic_electron.dp0 = c_light / omega_p
tuple	khiA1_mixed_ion_kinetic_electron.L_debye0 = (eps0*T0/(n0*qe**2))**0.5
float	khiA1_mixed_ion_kinetic_electron.mass_ratio = 25.0
float	khiA1_mixed_ion_kinetic_electron.Ai = 1.0
float	khiA1_mixed_ion_kinetic_electron.Ae = Ai / mass_ratio
float	khiA1_mixed_ion_kinetic_electron.Zi = 1.0
float	khiA1_mixed_ion_kinetic_electron.Ze = -1.0
float	khiA1_mixed_ion_kinetic_electron.d_width_sim = 1.0/20.0
float	khiA1_mixed_ion_kinetic_electron.b_sim = -10.0
float	khiA1_mixed_ion_kinetic_electron.T_sim = 6.25e-4
float	khiA1_mixed_ion_kinetic_electron.E0_sim = 2.00E-2
float	khiA1_mixed_ion_kinetic_electron.Ex0_sim = E0_sim / omega_c_tau
float	khiA1_mixed_ion_kinetic_electron.Bz_x0_sim = 1.0
float	khiA1_mixed_ion_kinetic_electron.Tx_sim = T_sim / 2.0
float	khiA1_mixed_ion_kinetic_electron.ky_sim_d_sim = 0.4
float	khiA1_mixed_ion_kinetic_electron.Lx_sim = 1.0
float	khiA1_mixed_ion_kinetic_electron.Ly_sim = 2.0*math.pi*d_width_sim/ky_sim_d_sim
list	khiA1_mixed_ion_kinetic_electron.center_sim = [0,0]
float	khiA1_mixed_ion_kinetic_electron.perturbation_sim = 2.0E-4
float	khiA1_mixed_ion_kinetic_electron.fluid_ion_perturbation_sim = 0.0
float	khiA1_mixed_ion_kinetic_electron.fluid_electron_perturbation_sim = perturbation_sim
float	khiA1_mixed_ion_kinetic_electron.omega_ci_sim = Ai*Zi/Bz_x0_sim
float	khiA1_mixed_ion_kinetic_electron.omega_ce_sim = Ae*Ze/Bz_x0_sim
math	khiA1_mixed_ion_kinetic_electron.rli_sim_over_d_sim = math.sqrt(T_sim/Ai)/omega_ci_sim / d_width_sim
float	khiA1_mixed_ion_kinetic_electron.n_sim_ref = 1.0
math	khiA1_mixed_ion_kinetic_electron.Ld_sim_over_d_sim = math.sqrt(T_sim / n_sim_ref /Ze**2) / d_width_sim
int	khiA1_mixed_ion_kinetic_electron.beta_sim_left = 2 * (2 * n_sim_ref)*T_sim/Bz_x0_sim**2
int	khiA1_mixed_ion_kinetic_electron.number_density_min = 1e-12
int	khiA1_mixed_ion_kinetic_electron.pressure_min = number_density_min*T_sim
	khiA1_mixed_ion_kinetic_electron.ni_aux_fit_xvec = str2float(ni_aux_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.ni_aux_fit_a = str2float(ni_aux_fit_a)
	khiA1_mixed_ion_kinetic_electron.ni_aux_fit_b = ni_aux_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.ni_aux_reader = csv.reader(ni_aux_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.phi_fit_xvec = str2float(phi_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.phi_fit_a = str2float(phi_fit_a)
	khiA1_mixed_ion_kinetic_electron.phi_fit_b = phi_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.phi_reader = csv.reader(phi_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.Bz_fit_xvec = str2float(Bz_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.Bz_fit_a = str2float(Bz_fit_a)
	khiA1_mixed_ion_kinetic_electron.Bz_fit_b = Bz_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.Bz_reader = csv.reader(Bz_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.Ex_fit_xvec = str2float(Ex_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.Ex_fit_a = str2float(Ex_fit_a)
	khiA1_mixed_ion_kinetic_electron.Ex_fit_b = Ex_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.Ex_reader = csv.reader(Ex_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.Ay_fit_xvec = str2float(Ay_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.Ay_fit_a = str2float(Ay_fit_a)
	khiA1_mixed_ion_kinetic_electron.Ay_fit_b = Ay_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.Ay_reader = csv.reader(Ay_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.ni_fit_xvec = str2float(ni_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.ni_fit_a = str2float(ni_fit_a)
	khiA1_mixed_ion_kinetic_electron.ni_fit_b = ni_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.ni_reader = csv.reader(ni_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.ne_fit_xvec = str2float(ne_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.ne_fit_a = str2float(ne_fit_a)
	khiA1_mixed_ion_kinetic_electron.ne_fit_b = ne_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.ne_reader = csv.reader(ne_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.vyi_fit_xvec = str2float(vyi_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.vyi_fit_a = str2float(vyi_fit_a)
	khiA1_mixed_ion_kinetic_electron.vyi_fit_b = vyi_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.vyi_reader = csv.reader(vyi_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.vye_fit_xvec = str2float(vye_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.vye_fit_a = str2float(vye_fit_a)
	khiA1_mixed_ion_kinetic_electron.vye_fit_b = vye_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.vye_reader = csv.reader(vye_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.Tyi_fit_xvec = str2float(Tyi_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.Tyi_fit_a = str2float(Tyi_fit_a)
	khiA1_mixed_ion_kinetic_electron.Tyi_fit_b = Tyi_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.Tyi_reader = csv.reader(Tyi_file, delimiter=',')
	khiA1_mixed_ion_kinetic_electron.Tye_fit_xvec = str2float(Tye_fit_xvec)
	khiA1_mixed_ion_kinetic_electron.Tye_fit_a = str2float(Tye_fit_a)
	khiA1_mixed_ion_kinetic_electron.Tye_fit_b = Tye_fit_b[0:2]
csv	khiA1_mixed_ion_kinetic_electron.Tye_reader = csv.reader(Tye_file, delimiter=',')
float	khiA1_mixed_ion_kinetic_electron.min_vx_ion = -0.2
float	khiA1_mixed_ion_kinetic_electron.min_vy_ion = -0.2
float	khiA1_mixed_ion_kinetic_electron.max_vx_ion = +0.2
float	khiA1_mixed_ion_kinetic_electron.max_vy_ion = +0.2
float	khiA1_mixed_ion_kinetic_electron.min_vx_electron = -1.0
float	khiA1_mixed_ion_kinetic_electron.min_vy_electron = -1.0
float	khiA1_mixed_ion_kinetic_electron.max_vx_electron = +1.0
float	khiA1_mixed_ion_kinetic_electron.max_vy_electron = +1.0
int	khiA1_mixed_ion_kinetic_electron.v_space_num_vx = 32
int	khiA1_mixed_ion_kinetic_electron.v_space_num_vy = 32
list	khiA1_mixed_ion_kinetic_electron.v_space_num_elements = [v_space_num_vx, v_space_num_vy]
list	khiA1_mixed_ion_kinetic_electron.min_v_ion = [min_vx_ion, min_vy_ion]
list	khiA1_mixed_ion_kinetic_electron.dv_ion = [(max_vx_ion-min_vx_ion)/v_space_num_vx, (max_vy_ion-min_vy_ion)/v_space_num_vy]
list	khiA1_mixed_ion_kinetic_electron.min_v_electron = [min_vx_electron, min_vy_electron]
list	khiA1_mixed_ion_kinetic_electron.dv_electron = [(max_vx_electron-min_vx_electron)/v_space_num_vx, (max_vy_electron-min_vy_electron)/v_space_num_vy]
str	khiA1_mixed_ion_kinetic_electron.meshfile = '../../../../../../user_runs/andrew_ho/meshes/khi_xl8m16r8_y8.inp'
warpy	khiA1_mixed_ion_kinetic_electron.mesh
list	khiA1_mixed_ion_kinetic_electron.fluid_comps = ['rho', 'px', 'py', 'pz', 'e']
warpy	khiA1_mixed_ion_kinetic_electron.ion_fluid
list	khiA1_mixed_ion_kinetic_electron.field_comps = ['Ex', 'Ey', 'Ez', 'Bx', 'By', 'Bz']
warpy	khiA1_mixed_ion_kinetic_electron.field
warpy	khiA1_mixed_ion_kinetic_electron.ion_distribution
warpy	khiA1_mixed_ion_kinetic_electron.electron_distribution
warpy	khiA1_mixed_ion_kinetic_electron.kinetic_ion_fluid
warpy	khiA1_mixed_ion_kinetic_electron.kinetic_electron_fluid
warpy	khiA1_mixed_ion_kinetic_electron.chi_metric_i
warpy	khiA1_mixed_ion_kinetic_electron.chi_metric_e
warpy	khiA1_mixed_ion_kinetic_electron.ion_fluid_ic
warpy	khiA1_mixed_ion_kinetic_electron.ion_distribution_ic
warpy	khiA1_mixed_ion_kinetic_electron.electron_distribution_ic
warpy	khiA1_mixed_ion_kinetic_electron.field_ic
warpy	khiA1_mixed_ion_kinetic_electron.fluid_ic_va
warpy	khiA1_mixed_ion_kinetic_electron.ion_distribution_ic_va
warpy	khiA1_mixed_ion_kinetic_electron.electron_distribution_ic_va
warpy	khiA1_mixed_ion_kinetic_electron.field_ic_va
warpy	khiA1_mixed_ion_kinetic_electron.tf_ic_ha
	load from file
warpy	khiA1_mixed_ion_kinetic_electron.ion_distribution_ic_ha
warpy	khiA1_mixed_ion_kinetic_electron.electron_distribution_ic_ha
warpy	khiA1_mixed_ion_kinetic_electron.field_ic_ha
list	khiA1_mixed_ion_kinetic_electron.ic_ha = [tf_ic_ha, ion_distribution_ic_ha, electron_distribution_ic_ha, field_ic_ha]
list	khiA1_mixed_ion_kinetic_electron.variable_adjusters_tf = []
list	khiA1_mixed_ion_kinetic_electron.variable_adjusters_field = []
list	khiA1_mixed_ion_kinetic_electron.variable_adjusters_kinetic = []
list	khiA1_mixed_ion_kinetic_electron.applications_moments_ions = []
list	khiA1_mixed_ion_kinetic_electron.applications_2nd_moments_ions = []
list	khiA1_mixed_ion_kinetic_electron.applications_moments_electrons = []
list	khiA1_mixed_ion_kinetic_electron.applications_2nd_moments_electrons = []
list	khiA1_mixed_ion_kinetic_electron.applications_maxwellian_metric_ions = []
list	khiA1_mixed_ion_kinetic_electron.applications_maxwellian_metric_electrons = []
list	khiA1_mixed_ion_kinetic_electron.kinetic_to_tf_va = []
list	khiA1_mixed_ion_kinetic_electron.tf_to_kinetic_va = []
list	khiA1_mixed_ion_kinetic_electron.virtual_bc_apps_tf = []
list	khiA1_mixed_ion_kinetic_electron.virtual_bc_apps_kinetic_ions = []
warpy	khiA1_mixed_ion_kinetic_electron.ion_timestep_start_moment_taker
warpy	khiA1_mixed_ion_kinetic_electron.ion_timestep_start_second_moment_taker
warpy	khiA1_mixed_ion_kinetic_electron.ion_timestep_start_chi_metric_moment_taker
warpy	khiA1_mixed_ion_kinetic_electron.electron_timestep_start_moment_taker
warpy	khiA1_mixed_ion_kinetic_electron.electron_timestep_start_second_moment_taker
warpy	khiA1_mixed_ion_kinetic_electron.electron_timestep_start_chi_metric_moment_taker
list	khiA1_mixed_ion_kinetic_electron.extra_timestep_start_writer
list	khiA1_mixed_ion_kinetic_electron.solver_apps_tf = []
list	khiA1_mixed_ion_kinetic_electron.spatial_solver_tf = []
list	khiA1_mixed_ion_kinetic_electron.spatial_solver_kinetic = []
list	khiA1_mixed_ion_kinetic_electron.solver_apps_kinetic_ions
list	khiA1_mixed_ion_kinetic_electron.solver_apps_kinetic_electrons
list	khiA1_mixed_ion_kinetic_electron.solver_apps_maxwell_flux
list	khiA1_mixed_ion_kinetic_electron.solver_apps_maxwell_tf_fluid_source
list	khiA1_mixed_ion_kinetic_electron.solver_apps_maxwell_kinetic_ion_fluid_source
list	khiA1_mixed_ion_kinetic_electron.solver_apps_maxwell_kinetic_electron_fluid_source
list	khiA1_mixed_ion_kinetic_electron.spatial_solver_field = []
list	khiA1_mixed_ion_kinetic_electron.ss_s = spatial_solver_tf + spatial_solver_kinetic + spatial_solver_field
list	khiA1_mixed_ion_kinetic_electron.va_s = variable_adjusters_tf + variable_adjusters_kinetic + variable_adjusters_field + kinetic_to_tf_va + tf_to_kinetic_va
warpy	khiA1_mixed_ion_kinetic_electron.time_integrator
list	khiA1_mixed_ion_kinetic_electron.time_sim = [0, target_tws / omega_p_tau / omega_s_over_omega_p ]
warpy	khiA1_mixed_ion_kinetic_electron.dt_controller = warpy.dt_calc.stability_dt(init_dt = dt_sim)
warpy	khiA1_mixed_ion_kinetic_electron.ion_probe_fluid = kinetic_ion_fluid
list	khiA1_mixed_ion_kinetic_electron.ion_probe_fluid_components = ['ni_0', 'nvxi_0', 'nvyi_0', 'nvzi_0', 'nTi_0']
warpy	khiA1_mixed_ion_kinetic_electron.ions_probe0
warpy	khiA1_mixed_ion_kinetic_electron.integ_Ex2
warpy	khiA1_mixed_ion_kinetic_electron.integ_Ey2
warpy	khiA1_mixed_ion_kinetic_electron.integ_Ez2
warpy	khiA1_mixed_ion_kinetic_electron.integ_Bx2
warpy	khiA1_mixed_ion_kinetic_electron.integ_By2
warpy	khiA1_mixed_ion_kinetic_electron.integ_Bz2
warpy	khiA1_mixed_ion_kinetic_electron.diag_writer
list	khiA1_mixed_ion_kinetic_electron.write_vars
warpy	khiA1_mixed_ion_kinetic_electron.writer
datetime	khiA1_mixed_ion_kinetic_electron.now = datetime.datetime.now().strftime("%Y_%m_%d")
f	khiA1_mixed_ion_kinetic_electron.sim_name = f'khiA1_mixed_i_kin_e_pert_3rd_tws60_32x8x32x32'
warpy	khiA1_mixed_ion_kinetic_electron.config_args = warpy.load_configs()
warpy	khiA1_mixed_ion_kinetic_electron.case_loc = config_args['test_rundir'](f'hybrid/2d2v/khi/{now}', '')
f	khiA1_mixed_ion_kinetic_electron.sim_param
warpy	khiA1_mixed_ion_kinetic_electron.sim
	khiA1_mixed_ion_kinetic_electron.gen_xdmf
	Run Simulation ... now = datetime.datetime.now() sim.run('hybrid/2d2v/khi/'+now.strftime("%Y_%m_%d"),.
	khiA1_mixed_ion_kinetic_electron.detect_nonscalar