From 2c84811291d7b41a8d2df346bf0209740239a744 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Tue, 28 Feb 2023 16:22:44 +0100
Subject: [PATCH] update scripts
---
testcases/smallest_problem/fort_00.90 | 22 +++----
testcases/smallest_problem/fort_01.90 | 83 +++++++++++++++++++++++++++
wk/analysis_gyacomo.m | 6 +-
wk/header_3D_results.m | 5 +-
4 files changed, 101 insertions(+), 15 deletions(-)
create mode 100644 testcases/smallest_problem/fort_01.90
diff --git a/testcases/smallest_problem/fort_00.90 b/testcases/smallest_problem/fort_00.90
index e2108230..7a0a83bf 100644
--- a/testcases/smallest_problem/fort_00.90
+++ b/testcases/smallest_problem/fort_00.90
@@ -44,35 +44,35 @@
! Collisionality
CLOS = 0
NL_CLOS = -1
- LINEARITY = 'nonlinear'
- KIN_E = .t.
- mu_x = 1.0
- mu_y = 1.0
+ LINEARITY = 'linear'
+ KIN_E = .f.
+ mu_x = 0.0
+ mu_y = 0.0
N_HD = 4
mu_z = 0.0
mu_p = 0
mu_j = 0
- nu = 0.1
+ nu = 1
tau_e = 1
tau_i = 1
sigma_e = 0.023338
sigma_i = 1
q_e = -1
q_i = 1
- K_Ne = 2.0
- K_Te = 0.4
- K_Ni = 2.0
- K_Ti = 0.4
+ K_Ne = 6.96
+ K_Te = 2.22
+ K_Ni = 6.96
+ K_Ti = 2.22
beta = 0
/
&COLLISION_PAR
collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
- gyrokin_CO = .true.
+ gyrokin_CO = .f.
interspecies = .true.
!mat_file = 'gk_sugama_P_20_J_10_N_150_kpm_8.0.h5'
/
&INITIAL_CON
- INIT_OPT = 'phi'
+ INIT_OPT = 'blob'
ACT_ON_MODES = 'donothing'
init_background = 0
init_noiselvl = 0.005
diff --git a/testcases/smallest_problem/fort_01.90 b/testcases/smallest_problem/fort_01.90
new file mode 100644
index 00000000..f5a7d560
--- /dev/null
+++ b/testcases/smallest_problem/fort_01.90
@@ -0,0 +1,83 @@
+&BASIC
+ nrun = 1
+ dt = 0.01
+ tmax = 1
+ maxruntime = 356400
+/
+&GRID
+ pmaxe = 2
+ jmaxe = 1
+ pmaxi = 2
+ jmaxi = 1
+ Nx = 2
+ Lx = 200
+ Ny = 2
+ Ly = 60
+ Nz = 4
+ SG = .f.
+/
+&GEOMETRY
+ geom = 's-alpha'
+ q0 = 1.4
+ shear = 0.8
+ eps = 0.18
+ parallel_bc = 'dirichlet'
+/
+&OUTPUT_PAR
+ nsave_0d = 1
+ nsave_1d = -1
+ nsave_2d = -1
+ nsave_3d = 1
+ nsave_5d = 1
+ write_doubleprecision = .f.
+ write_gamma = .t.
+ write_hf = .t.
+ write_phi = .t.
+ write_Na00 = .f.
+ write_Napj = .t.
+ write_Sapj = .f.
+ write_dens = .t.
+ write_temp = .t.
+ job2load = -1
+/
+&MODEL_PAR
+ ! Collisionality
+ CLOS = 0
+ NL_CLOS = -1
+ LINEARITY = 'linear'
+ KIN_E = .f.
+ mu_x = 0.0
+ mu_y = 0.0
+ N_HD = 4
+ mu_z = 0.0
+ mu_p = 0
+ mu_j = 0
+ nu = 1
+ tau_e = 1
+ tau_i = 1
+ sigma_e = 0.023338
+ sigma_i = 1
+ q_e = -1
+ q_i = 1
+ K_Ne = 6.96
+ K_Te = 2.22
+ K_Ni = 6.96
+ K_Ti = 2.22
+ beta = 0
+/
+&COLLISION_PAR
+ collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+ gyrokin_CO = .t.
+ 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
+ init_noiselvl = 0.005
+ iseed = 42
+/
+&TIME_INTEGRATION_PAR
+ numerical_scheme = 'RK4'
+/
diff --git a/wk/analysis_gyacomo.m b/wk/analysis_gyacomo.m
index a41a72ef..855f6ea9 100644
--- a/wk/analysis_gyacomo.m
+++ b/wk/analysis_gyacomo.m
@@ -110,9 +110,9 @@ options.NAME = 'N_i^{00}';
% options.NAME = 's_{Ex}';
% options.NAME = 'Q_x';
% options.NAME = 'k^2n_e';
-options.PLAN = 'kxky';
-options.COMP = 24;
-options.TIME = [600 700 800 1000];
+options.PLAN = 'xz';
+options.COMP = 'avg';
+options.TIME = [500 1000 1500];
options.RESOLUTION = 256;
data.a = data.EPS * 2e3;
diff --git a/wk/header_3D_results.m b/wk/header_3D_results.m
index 32541c52..58375962 100644
--- a/wk/header_3D_results.m
+++ b/wk/header_3D_results.m
@@ -9,7 +9,7 @@ PARTITION = '/misc/gyacomo_outputs/';
% resdir = 'paper_2_nonlinear/kT_6.96/3x2x128x64x24';
% resdir = 'paper_2_nonlinear/kT_6.96/5x3x128x64x24'; %+ diff study
% resdir = 'paper_2_nonlinear/kT_6.96/5x3x128x64x24_dv4_diff';
-resdir = 'paper_2_nonlinear/kT_6.96/optimal_muz_nu_5x3x128x64x24';
+% resdir = 'paper_2_nonlinear/kT_6.96/optimal_muz_nu_5x3x128x64x24';
% resdir = 'paper_2_nonlinear/kT_6.96/7x4x128x64x24';
% resdir = 'paper_2_nonlinear/kT_6.96/9x5x128x64x24';
% low resolution (Cyclic)
@@ -51,9 +51,12 @@ resdir = 'paper_2_nonlinear/kT_6.96/optimal_muz_nu_5x3x128x64x24';
% resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x24_MUxy_0';
% resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x24_NL_-1';
% resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x24_nuDG_0.01';
+% resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x64';
% resdir = 'paper_2_nonlinear/kT_5.3/7x4x192x96x64';
% resdir = 'paper_2_nonlinear/kT_5.3/9x5x128x64x24';
+% resdir = 'paper_2_nonlinear/kT_5.3/9x5x128x64x64';
% resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x24';
+resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x64';
%% Old stuff
% resdir = 'CBC/kT_4.5_128x64x16x13x7/';
--
GitLab