diff --git a/testcases/CBC_ExBshear/fort_00.90 b/testcases/CBC_ExBshear/fort_00.90 new file mode 100644 index 0000000000000000000000000000000000000000..a65bd0441c9e8059f9f1ceb7544df753f45c87a0 --- /dev/null +++ b/testcases/CBC_ExBshear/fort_00.90 @@ -0,0 +1,97 @@ +&BASIC + nrun = 99999999 + dt = 0.01 + tmax = 50 + maxruntime = 72000 + job2load = -1 +/ +&GRID + pmax = 2 + jmax = 1 + Nx = 64 + Lx = 120 + Ny = 48 + Ly = 120 + Nz = 16 + SG = .f. + Nexc = 0 +/ +&GEOMETRY + geom = 's-alpha' + !geom = 'miller' + q0 = 1.4 + shear = 0.8 + eps = 0.18 + kappa = 1.0 + s_kappa= 0.0 + delta = 0.0 + s_delta= 0.0 + zeta = 0.0 + s_zeta = 0.0 + parallel_bc = 'dirichlet' + shift_y= 0.0 +/ +&DIAGNOSTICS + dtsave_0d = 1 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 1 + dtsave_5d = 20 + write_doubleprecision = .f. + write_gamma = .t. + write_hf = .t. + write_phi = .t. + write_Na00 = .t. + write_Napj = .t. + write_dens = .t. + write_fvel = .t. + write_temp = .t. +/ +&MODEL + LINEARITY = 'nonlinear' + Na = 1 ! number of species + mu_x = 0.0 + mu_y = 0.0 + N_HD = 4 + mu_z = 1.0 + HYP_V = 'hypcoll' + mu_p = 0.0 + mu_j = 0.0 + nu = 0.05 + beta = 0.0 + ADIAB_E = .t. + ExBrate = 0.01 +/ +&CLOSURE + hierarchy_closure='truncation' + dmax = -1 + nonlinear_closure='truncation' + nmax = 0 +/ +&SPECIES + ! ions + name_ = 'ions' + tau_ = 1.0 + sigma_= 1.0 + q_ = 1.0 + k_N_ = 2.22 + k_T_ = 6.96 +/ + +&COLLISION + collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau) + GK_CO = .f. + INTERSPECIES = .true. + mat_file = 'gk_sugama_P_20_J_10_N_150_kpm_8.0.h5' +/ +&INITIAL + INIT_OPT = 'blob' + ACT_ON_MODES = 'donothing' + init_background = 0.0 + init_noiselvl = 0.005 + iseed = 42 +/ +&TIME_INTEGRATION + !numerical_scheme = 'RK4' + numerical_scheme = 'SSP_RK3' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_00.90 b/testcases/ExB_shear_atomic_tests/fort_00.90 new file mode 100644 index 0000000000000000000000000000000000000000..29ec2972c9c8fbb690570215f5a4446f09f8ed53 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_00.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.00075 + tmax = 2.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 186 + Lx = 100 + Ny = 186 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.01 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.01 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 1 + mu_y = 1 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 1 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ikxZF = 2 + ZFamp = 400 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =0 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'ricci' + init_amp = 0.1 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_01.90 b/testcases/ExB_shear_atomic_tests/fort_01.90 new file mode 100644 index 0000000000000000000000000000000000000000..2aa6249f0ade57f2e74e98b317b2a05928776589 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_01.90 @@ -0,0 +1,98 @@ +&BASIC + nrun = 100000000 + dt = 0.00075 + tmax = 0.5 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 186 + Lx = 100 + Ny = 186 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.01 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.01 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'linear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 1 + mu_y = 1 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 1 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 5.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =0 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'ricci' +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_02.90 b/testcases/ExB_shear_atomic_tests/fort_02.90 new file mode 100644 index 0000000000000000000000000000000000000000..898315bbe614b7ce923ea8569153332c14447bd1 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_02.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'linear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 3.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 3 + iky_init = 3 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_03.90 b/testcases/ExB_shear_atomic_tests/fort_03.90 new file mode 100644 index 0000000000000000000000000000000000000000..d1cbbdfaec4c99e3bc455b25bd6669bcc6be77c7 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_03.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 2 + iky_init = 2 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_04.90 b/testcases/ExB_shear_atomic_tests/fort_04.90 new file mode 100644 index 0000000000000000000000000000000000000000..05e8113b1088666dfe5a2bb5b6ec523d64a89775 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_04.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 3 + iky_init = 3 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_05.90 b/testcases/ExB_shear_atomic_tests/fort_05.90 new file mode 100644 index 0000000000000000000000000000000000000000..9fcce56d077c0dce70c9d61caf1b940c8129ed5d --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_05.90 @@ -0,0 +1,100 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 120 + Ny = 64 + Ly = 90 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 3.0 + ExB_NL_CORRECTION = .true. + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/