update scripts

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

testcases/smallest_problem/fort_01.90

+&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'
+/

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;

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/';