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

cleaning

parent 7cfc2570
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 10 additions and 288 deletions
testcases/cyclone_example/cyclone_example.txt deleted 100644 → 0
+ 0
− 4
View file @ 7cfc2570
This is a testcase reproducing the cyclone base case of Dimits 2000
Adiabatic electrons, s-alpha geometry, gradlnN = 2.22, gradlnT = 6.96
With a small P,J=4,2 polynomial basis, one should observe the secondary instability (KHI) at t~50R/cs.
The saturated heat flux should be located around Qx ~ 30, i.e. Qx/QGB ~ 2 which is close to Dimits results.
testcases/cyclone_example/fort_00.90
+ 7
− 7
View file @ 32b5d860
&BASIC &BASIC
nrun = 99999999 nrun = 99999999
dt = 0.01 dt = 0.01
tmax = 50 tmax = 500
maxruntime = 356400 maxruntime = 72000
job2load = -1 job2load = -1
/ /
&GRID &GRID
...@@ -10,9 +10,9 @@ ...@@ -10,9 +10,9 @@
jmax = 2 jmax = 2
Nx = 128 Nx = 128
Lx = 120 Lx = 120
Ny = 48 Ny = 64
Ly = 120 Ly = 120
Nz = 16 Nz = 24
SG = .f. SG = .f.
Nexc = 0 Nexc = 0
/ /
...@@ -31,11 +31,11 @@ ...@@ -31,11 +31,11 @@
shift_y= 0.0 shift_y= 0.0
/ /
&OUTPUT_PAR &OUTPUT_PAR
dtsave_0d = 0.1 dtsave_0d = 1
dtsave_1d = -1 dtsave_1d = -1
dtsave_2d = -1 dtsave_2d = -1
dtsave_3d = 1 dtsave_3d = 5
dtsave_5d = 10 dtsave_5d = 20
write_doubleprecision = .f. write_doubleprecision = .f.
write_gamma = .t. write_gamma = .t.
write_hf = .t. write_hf = .t.
......
testcases/cyclone_example/fort_01.90 deleted 100644 → 0
+ 0
− 92
View file @ 7cfc2570
&BASIC
nrun = 99999999
dt = 0.01
tmax = 75
maxruntime = 356400
job2load = 0
/
&GRID
pmax = 4
jmax = 2
Nx = 128
Lx = 120
Ny = 48
Ly = 120
Nz = 6!16
SG = .f.
Nexc = 0
/
&GEOMETRY
geom = 's-alpha'
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
/
&OUTPUT_PAR
dtsave_0d = 0.1
dtsave_1d = -1
dtsave_2d = -1
dtsave_3d = 1
dtsave_5d = 10
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_PAR
! Collisionality
CLOS = 0
NL_CLOS = 0
LINEARITY = 'nonlinear'
Na = 1 ! number of species
mu_x = 1.0
mu_y = 1.0
N_HD = 4
mu_z = 2.0
HYP_V = 'hypcoll'
mu_p = 0.0
mu_j = 0.0
nu = 0.001
beta = 0.0
ADIAB_E = .t.
tau_e = 1.0
/
&SPECIES
! ions
name_ = 'ions'
tau_ = 1.0
sigma_= 1.0
q_ = 1.0
k_N_ = 2.22
k_T_ = 6.96
/
&COLLISION_PAR
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_CON
INIT_OPT = 'blob'
ACT_ON_MODES = 'donothing'
init_background = 0.0
init_noiselvl = 0.005
iseed = 42
/
&TIME_INTEGRATION_PAR
numerical_scheme = 'RK4'
/
testcases/cyclone_example/fort_02.90 deleted 100644 → 0
+ 0
− 92
View file @ 7cfc2570
&BASIC
nrun = 10000000
dt = 0.01
tmax = 150
maxruntime = 356400
job2load = 1
/
&GRID
pmax = 4
jmax = 2
Nx = 128
Lx = 120
Ny = 48
Ly = 120
Nz = 16
SG = .f.
Nexc = 0
/
&GEOMETRY
geom = 's-alpha'
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
/
&OUTPUT_PAR
dtsave_0d = 1
dtsave_1d = -1
dtsave_2d = -1
dtsave_3d = 1
dtsave_5d = 10
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_PAR
! Collisionality
CLOS = 0
NL_CLOS = 0
LINEARITY = 'nonlinear'
Na = 1 ! number of species
mu_x = 1.0
mu_y = 1.0
N_HD = 4
mu_z = 2.0
HYP_V = 'hypcoll'
mu_p = 0.0
mu_j = 0.0
nu = 0.001
beta = 0.0
ADIAB_E = .t.
tau_e = 1.0
/
&SPECIES
! ions
name_ = 'ions'
tau_ = 1.0
sigma_= 1.0
q_ = 1.0
k_N_ = 2.22
k_T_ = 6.96
/
&COLLISION_PAR
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_CON
INIT_OPT = 'blob'
ACT_ON_MODES = 'donothing'
init_background = 0.1
init_noiselvl = 0.00
iseed = 42
/
&TIME_INTEGRATION_PAR
numerical_scheme = 'RK4'
/
testcases/cyclone_example/fort_03.90 deleted 100644 → 0
+ 0
− 92
View file @ 7cfc2570
&BASIC
nrun = 10000000
dt = 0.01
tmax = 155
maxruntime = 356400
job2load = 2
/
&GRID
pmax = 4
jmax = 2
Nx = 128
Lx = 120
Ny = 48
Ly = 120
Nz = 16
SG = .f.
Nexc = 0
/
&GEOMETRY
geom = 's-alpha'
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
/
&OUTPUT_PAR
dtsave_0d = 0.01
dtsave_1d = -1
dtsave_2d = -1
dtsave_3d = 1
dtsave_5d = 10
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_PAR
! Collisionality
CLOS = 0
NL_CLOS = 0
LINEARITY = 'nonlinear'
Na = 1 ! number of species
mu_x = 1.0
mu_y = 1.0
N_HD = 4
mu_z = 2.0
HYP_V = 'hypcoll'
mu_p = 0.0
mu_j = 0.0
nu = 0.001
beta = 0.0
ADIAB_E = .t.
tau_e = 1.0
/
&SPECIES
! ions
name_ = 'ions'
tau_ = 1.0
sigma_= 1.0
q_ = 1.0
k_N_ = 2.22
k_T_ = 6.96
/
&COLLISION_PAR
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_CON
INIT_OPT = 'blob'
ACT_ON_MODES = 'donothing'
init_background = 0.1
init_noiselvl = 0.00
iseed = 42
/
&TIME_INTEGRATION_PAR
numerical_scheme = 'RK4'
/
testcases/cyclone_example/gyacomo23_debug 0 → 120000
+ 1
− 0
View file @ 32b5d860
../../bin/gyacomo23_debug
\ No newline at end of file
testcases/cyclone_example/gyacomo23_dp 0 → 120000
+ 1
− 0
View file @ 32b5d860
../../bin/gyacomo23_dp
\ No newline at end of file
testcases/cyclone_example/gyacomo23_sp 0 → 120000
+ 1
− 0
View file @ 32b5d860
../../bin/gyacomo23_sp
\ No newline at end of file
testcases/cyclone_example/gyacomo_debug deleted 120000 → 0
+ 0
− 1
View file @ 7cfc2570
../../bin/gyacomo_debug
\ No newline at end of file
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