From 27019b9e1667059feb895b139dca82e21aadc9a9 Mon Sep 17 00:00:00 2001
From: Antoine <antoine.hoffmann@epfl.ch>
Date: Fri, 22 Sep 2023 10:08:46 +0200
Subject: [PATCH] improve the test ExB

---
 testcases/ExB_shear_atomic_tests/fort_02.90 | 11 ++++++-----
 testcases/ExB_shear_atomic_tests/fort_03.90 | 13 +++++++------
 testcases/ExB_shear_atomic_tests/fort_04.90 |  9 +++++----
 wk/test_ExB_shear.m                         | 21 +++++++++++----------
 4 files changed, 29 insertions(+), 25 deletions(-)

diff --git a/testcases/ExB_shear_atomic_tests/fort_02.90 b/testcases/ExB_shear_atomic_tests/fort_02.90
index 898315bb..f4070b56 100644
--- a/testcases/ExB_shear_atomic_tests/fort_02.90
+++ b/testcases/ExB_shear_atomic_tests/fort_02.90
@@ -1,7 +1,7 @@
 &BASIC
   nrun       = 100000000
   dt         = 0.001
-  tmax       = 5.0
+  tmax       = 2.0
   maxruntime = 356400
   job2load   = -1
 /
@@ -36,7 +36,7 @@
   dtsave_0d = 0.05
   dtsave_1d = -1
   dtsave_2d = -1
-  dtsave_3d = 0.05
+  dtsave_3d = 0.02
   dtsave_5d = 100
   write_doubleprecision = .true.
   write_gamma = .true.
@@ -48,7 +48,7 @@
   write_temp  = .true.
 /
 &MODEL
-LINEARITY = 'linear'
+LINEARITY = 'nonlinear'
 RM_LD_T_EQ= .false.
   Na      = 1
   mu_x    = 0
@@ -63,7 +63,8 @@ RM_LD_T_EQ= .false.
   k_cB    = 0
   lambdaD = 0
   beta    = 0.0
-  ExBrate = 3.0
+  ExBrate = 0.0
+  ExB_NL_CORRECTION = .false.
   ADIAB_E = .true.
   ADIAB_I = .false.
   tau_i   = 1
@@ -92,7 +93,7 @@ RM_LD_T_EQ= .false.
 /
 &INITIAL
   INIT_OPT      = 'mom00_single_mode'
-  init_amp      = 100
+  init_amp      = 1000
   ikx_init      = 3
   iky_init      = 3
 /
diff --git a/testcases/ExB_shear_atomic_tests/fort_03.90 b/testcases/ExB_shear_atomic_tests/fort_03.90
index d1cbbdfa..6a943378 100644
--- a/testcases/ExB_shear_atomic_tests/fort_03.90
+++ b/testcases/ExB_shear_atomic_tests/fort_03.90
@@ -1,7 +1,7 @@
 &BASIC
   nrun       = 100000000
   dt         = 0.001
-  tmax       = 5.0
+  tmax       = 2.0
   maxruntime = 356400
   job2load   = -1
 /
@@ -36,7 +36,7 @@
   dtsave_0d = 0.05
   dtsave_1d = -1
   dtsave_2d = -1
-  dtsave_3d = 0.05
+  dtsave_3d = 0.02
   dtsave_5d = 100
   write_doubleprecision = .true.
   write_gamma = .true.
@@ -63,7 +63,8 @@ RM_LD_T_EQ= .false.
   k_cB    = 0
   lambdaD = 0
   beta    = 0.0
-  ExBrate = 0.0
+  ExBrate = 2.0
+  ExB_NL_CORRECTION = .false.
   ADIAB_E = .true.
   ADIAB_I = .false.
   tau_i   = 1
@@ -92,9 +93,9 @@ RM_LD_T_EQ= .false.
 /
 &INITIAL
   INIT_OPT      = 'mom00_single_mode'
-  init_amp      = 100
-  ikx_init      = 2
-  iky_init      = 2
+  init_amp      = 1000
+  ikx_init      = 3
+  iky_init      = 3
 /
 &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
index 05e8113b..3506310b 100644
--- a/testcases/ExB_shear_atomic_tests/fort_04.90
+++ b/testcases/ExB_shear_atomic_tests/fort_04.90
@@ -1,7 +1,7 @@
 &BASIC
   nrun       = 100000000
   dt         = 0.001
-  tmax       = 5.0
+  tmax       = 2.0
   maxruntime = 356400
   job2load   = -1
 /
@@ -36,7 +36,7 @@
   dtsave_0d = 0.05
   dtsave_1d = -1
   dtsave_2d = -1
-  dtsave_3d = 0.05
+  dtsave_3d = 0.02
   dtsave_5d = 100
   write_doubleprecision = .true.
   write_gamma = .true.
@@ -63,7 +63,8 @@ RM_LD_T_EQ= .false.
   k_cB    = 0
   lambdaD = 0
   beta    = 0.0
-  ExBrate = 0.0
+  ExBrate = 2.0
+  ExB_NL_CORRECTION = .true.
   ADIAB_E = .true.
   ADIAB_I = .false.
   tau_i   = 1
@@ -92,7 +93,7 @@ RM_LD_T_EQ= .false.
 /
 &INITIAL
   INIT_OPT      = 'mom00_single_mode'
-  init_amp      = 100
+  init_amp      = 1000
   ikx_init      = 3
   iky_init      = 3
 /
diff --git a/wk/test_ExB_shear.m b/wk/test_ExB_shear.m
index 9958861c..26484c7b 100644
--- a/wk/test_ExB_shear.m
+++ b/wk/test_ExB_shear.m
@@ -3,12 +3,13 @@ MVIN  = ['cd ',resdir];
 curdir = pwd;
 MVOUT = ['cd ',curdir];
 %------
-% MPIRUN = '/opt/homebrew/bin/mpirun';
-MPIRUN   = 'mpirun';
+MPIRUN = '/opt/homebrew/bin/mpirun';
+% MPIRUN   = 'mpirun';
 %------
-EXECDIR  = '/home/ahoffman/gyacomo/bin/';
-EXECNAME = 'gyacomo23_dp_O1';
-% EXECNAME = 'gyacomo23_dp';
+% EXECDIR  = '/home/ahoffman/gyacomo/bin/';
+EXECDIR  = '/Users/ahoffmann/gyacomo/bin/';
+% EXECNAME = 'gyacomo23_dp_O1';
+EXECNAME = 'gyacomo23_dp';
 % EXECNAME = 'gyacomo23_sp';
 % EXECNAME = 'gyacomo23_debug';
 %------
@@ -20,12 +21,12 @@ NP = '6'; PARA = '1 6 1';
 % INNAME = '0';
 % Picture, bckg ExB shear
 % INNAME = '1';
-% Linear with parallel modes initialization + ExB shear
+% NL with padrallel modes initialization no ExB shear
 % INNAME = '2';
-% NL with parallel modes initialization no ExB shear
-INNAME = '3';
-% Ultra reduced NL with parallel modes and 8x8 grid no shear
-% INNAME = '4';
+% NL with parallel modes initialization with ExB shear, no with factor
+% INNAME = '3';
+% NL with parallel modes initialization with ExB shear, with factor
+INNAME = '4';
 % NL with parallel modes initialization + ExB shear
 % INNAME = '5';
 
-- 
GitLab