Skip to content
Snippets Groups Projects
Commit 0b3c9a44 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

new test case for ExB shearing

parent 52e4fd8e
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 699 additions and 0 deletions
testcases/CBC_ExBshear/fort_00.90 0 → 100644
+ 97
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_00.90 0 → 100644
+ 101
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_01.90 0 → 100644
+ 98
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_02.90 0 → 100644
+ 101
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_03.90 0 → 100644
+ 101
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_04.90 0 → 100644
+ 101
− 0
View file @ 0b3c9a44
&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'
/
testcases/ExB_shear_atomic_tests/fort_05.90 0 → 100644
+ 100
− 0
View file @ 0b3c9a44
&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'
/
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