create a demonstration folder

demo/ExB_shear/2D_Z_pinch/GS2_method/fort_01.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.05
+  tmax       = 200
+  maxruntime = 72000
+  job2load   = 0
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = 0.5
+  dtsave_3d = 0.5
+  dtsave_5d = 100
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 1.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.1
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+  ExBrate = 0.5
+/
+&CLOSURE
+  hierarchy_closure='truncation'
+  !hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/ExB_shear/2D_Z_pinch/background_Er/fort_01.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.025
+  tmax       = 200
+  maxruntime = 72000
+  job2load   = 0
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = 0.5
+  dtsave_3d = 0.5
+  dtsave_5d = 100
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 1.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.1
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+  ZFrate  = 1.0
+  ikxZF   = 1
+/
+&CLOSURE
+  hierarchy_closure='truncation'
+  !hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/ExB_shear/2D_Z_pinch/control/fort_01.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.05
+  tmax       = 200
+  maxruntime = 72000
+  job2load   = 0
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = 0.5
+  dtsave_3d = 0.5
+  dtsave_5d = 100
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 1.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.1
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+  ExBrate = 0
+/
+&CLOSURE
+  hierarchy_closure='truncation'
+  !hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/ExB_shear/2D_Z_pinch/fort_00.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.05
+  tmax       = 100
+  maxruntime = 72000
+  job2load   = -1
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = 0.5
+  dtsave_3d = 0.5
+  dtsave_5d = 100
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 1.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.1
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+/
+&CLOSURE
+  hierarchy_closure='truncation'
+  !hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/ExB_shear/GS2_method/fast_analysis.m

+gyacomodir = '../../../';
+addpath(genpath([gyacomodir,'matlab'])) % ... add
+addpath(genpath([gyacomodir,'matlab/plot'])) % ... add
+addpath(genpath([gyacomodir,'matlab/compute'])) % ... add
+addpath(genpath([gyacomodir,'matlab/load'])) % ... add
+default_plots_options
+
+J0 = 00; J1 = 00;
+
+% Load basic info (grids and time traces)
+DATADIR = [pwd,'/'];
+data    = {};
+data    = compile_results_low_mem(data,DATADIR,J0,J1);
+[data.Na00, data.Ts3D] = compile_results_3Da(DATADIR,J0,J1,'Na00');
+data.Ni00 = reshape(data.Na00(1,:,:,:,:),data.grids.Nky,data.grids.Nkx,data.grids.Nz,numel(data.Ts3D));
+
+% field snapshots
+options.INTERP    = 0;
+options.POLARPLOT = 0;
+options.AXISEQUAL = 1;
+options.NORMALIZE = 0;
+options.LOGSCALE  = 0;
+options.CLIMAUTO  = 1;
+options.NAME      = ['N_i^{00}'];
+options.PLAN      = 'xy'; options.COMP =floor(data.grids.Nz/2)+1;
+options.TIME      = [0 0.5 1.0];
+options.RESOLUTION = 256;
+options.BWR       = 0; % bluewhitered plot or gray
+fig = photomaton(data,options);
+colormap(gray)
+clim('auto')
+data.FIGDIR = DATADIR;
+save_figure(data,fig,'.png');

demo/ExB_shear/GS2_method/fort_00.90

+&BASIC
+  nrun       = 100000000
+  dt         = 0.001
+  tmax       = 1.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   = 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'
+/

demo/ExB_shear/background_Er/fast_analysis.m

+gyacomodir = '../../../';
+addpath(genpath([gyacomodir,'matlab'])) % ... add
+addpath(genpath([gyacomodir,'matlab/plot'])) % ... add
+addpath(genpath([gyacomodir,'matlab/compute'])) % ... add
+addpath(genpath([gyacomodir,'matlab/load'])) % ... add
+default_plots_options
+
+J0 = 00; J1 = 00;
+
+% Load basic info (grids and time traces)
+DATADIR = [pwd,'/'];
+data    = {};
+data    = compile_results_low_mem(data,DATADIR,J0,J1);
+[data.Na00, data.Ts3D] = compile_results_3Da(DATADIR,J0,J1,'Na00');
+data.Ni00 = reshape(data.Na00(1,:,:,:,:),data.grids.Nky,data.grids.Nkx,data.grids.Nz,numel(data.Ts3D));
+
+% field snapshots
+options.INTERP    = 0;
+options.POLARPLOT = 0;
+options.AXISEQUAL = 1;
+options.NORMALIZE = 0;
+options.LOGSCALE  = 0;
+options.CLIMAUTO  = 1;
+options.NAME      = ['N_i^{00}'];
+options.PLAN      = 'xy'; options.COMP =floor(data.grids.Nz/2)+1;
+options.TIME      = [0 2.0 4.0];
+options.RESOLUTION = 256;
+options.BWR       = 0; % bluewhitered plot or gray
+fig = photomaton(data,options);
+colormap(gray)
+clim('auto')
+data.FIGDIR = DATADIR;
+% save_figure(data,fig,'.png');
+
+
+if 0
+%% MOVIES %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+options.INTERP    = 1;
+options.POLARPLOT = 0;
+options.BWR       = 0; % bluewhitered plot or gray
+options.CLIMAUTO  = 1; % adjust the colormap auto
+options.NAME     = ['N_i^{00}'];
+% options.NAME     = ['N_i^{00}'];
+options.PLAN      = ['xy'];
+options.COMP      = 1;
+options.TIME      =  data.Ts3D(1:1:end);
+data.EPS          = 0.1;
+data.a = data.EPS * 2000;
+options.RESOLUTION = 64;
+options.FPS       = 12;
+options.RMAXIS    = 1;
+create_film(data,options,'.gif')
+end

demo/ExB_shear/background_Er/fort.90

+&BASIC
+  nrun       = 100000000
+  dt         = 0.05
+  tmax       = 5.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.05
+  dtsave_1d = -1
+  dtsave_2d = -1
+  dtsave_3d = 0.025
+  dtsave_5d = -1
+  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      = 0
+  k_gB    = 0
+  k_cB    = 0
+  lambdaD = 0
+  beta    = 0.0
+  ExBrate = 0.0
+  ikxZF   = 1
+  ZFrate  =-0.5
+  ZF_ONLY = .t.
+  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 = 'mom00_mode'
+  Nmodes = 1
+/
+&MODE
+  I_   = 0
+  J_   = 1
+  amp_ = 1000
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/ExB_shear/background_Er/fort_00.90

+&BASIC
+  nrun       = 100000000
+  dt         = 0.05
+  tmax       = 10.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.05
+  dtsave_1d = -1
+  dtsave_2d = -1
+  dtsave_3d = 0.025
+  dtsave_5d = -1
+  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      = 0
+  k_gB    = 0
+  k_cB    = 0
+  lambdaD = 0
+  beta    = 0.0
+  ExBrate = 0
+  ikxZF   = 1
+  ZFrate  = 0.5
+  ZF_ONLY = .t.
+  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'
+/

demo/Laguerre_aliasing/2D_Z_pinch/Nmax_0/fort_00.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.02
+  tmax       = 500
+  maxruntime = 72000
+  job2load   = -1
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = -1
+  dtsave_3d = 2
+  dtsave_5d = 50
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 0.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.05
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+/
+&CLOSURE
+  !hierarchy_closure='truncation'
+  hierarchy_closure='max_degree'
+  dmax = 2
+  !nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nonlinear_closure='truncation'
+  nmax = 0
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/Laguerre_aliasing/2D_Z_pinch/anti_Laguerre_aliasing/fort_00.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.02
+  tmax       = 500
+  maxruntime = 72000
+  job2load   = -1
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = -1
+  dtsave_3d = 2
+  dtsave_5d = 50
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 0.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.05
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+/
+&CLOSURE
+  !hierarchy_closure='truncation'
+  hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/Laguerre_aliasing/2D_Z_pinch/fort_00.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.05
+  tmax       = 500
+  maxruntime = 72000
+  job2load   = -1
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = 0.5
+  dtsave_3d = 0.5
+  dtsave_5d = 100
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 0.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.05
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+/
+&CLOSURE
+  hierarchy_closure='truncation'
+  !hierarchy_closure='max_degree'
+  dmax = 2
+  nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/

demo/Laguerre_aliasing/2D_Z_pinch/maximized_sum/fort_00.90

+&BASIC
+  nrun       = 1e6
+  dt         = 0.02
+  tmax       = 500
+  maxruntime = 72000
+  job2load   = -1
+/
+&GRID
+  pmax   = 2
+  jmax   = 1
+  Nx     = 128
+  Lx     = 200
+  Ny     = 48
+  Ly     = 60
+  Nz     = 1
+  SG     = .f.
+  Nexc   = 1
+/
+&GEOMETRY
+  geom   = 'z-pinch'
+/
+&DIAGNOSTICS
+  dtsave_0d = 1
+  dtsave_1d = -1
+  dtsave_2d = -1
+  dtsave_3d = 2
+  dtsave_5d = 50
+  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      = 2 ! number of species
+  mu_x    = 1.0
+  mu_y    = 0.0
+  N_HD    = 4
+  mu_z    = 0.0
+  HYP_V   = 'hypcoll'
+  mu_p    = 0.0
+  mu_j    = 0.0
+  nu      = 0.05
+  beta    = 0.0
+  ADIAB_E = .f.
+  tau_i   = 1.0
+/
+&CLOSURE
+  !hierarchy_closure='truncation'
+  hierarchy_closure='max_degree'
+  dmax = 2
+  !nonlinear_closure='anti_laguerre_aliasing' !(truncation,full_sum,anti_laguerre_aliasing)
+  nonlinear_closure='full_sum'
+  nmax = -1
+/
+&SPECIES
+ ! ions
+ name_ = 'ions'
+ tau_  = 1.0
+ sigma_= 1.0
+ q_    = 1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&SPECIES
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau)
+  GK_CO           = .t.
+/
+&INITIAL
+  INIT_OPT         = 'phi' !(phi,blob)
+  init_background  = 0.0
+  init_noiselvl    = 0.005
+  iseed            = 42
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4'
+/