better parameters handling

basic_scripts/parameters_Zpinch.in

+&BASIC
+  nrun       = 99999999 ! number of time step to perform
+  dt         = 0.01     ! time step (not adaptive)
+  tmax       = 100       ! maximal time [c_s/R]
+  maxruntime = 72000    ! maximal wallclock runtime [sec]
+  job2load   = -1       ! index of the previous run to restart (-1:start a new run)
+/
+&GRID
+  pmax   = 2            ! maximal degree of the Hermite basis (parallel velocity)
+  jmax   = 1            ! maximal degree of the Laguerre basis (magnetic moment)
+  Nx     = 64           ! number of points in the radial direction
+  Lx     = 200          ! size of the box in the radial direction [rho_s]
+  Ny     = 48           ! bumber of points in the binormal direction
+  Ly     = 60           ! size of the box in the binormal direction [rho_s]
+  Nz     = 1            ! number of points in the magnetic field direction
+  SG     = .f.          ! use a staggered grid in z (experimental)
+  Nexc   = -1           ! to fullfill the sheared boundary condition (set -1 for automatic)
+/
+&GEOMETRY
+  geom   = 'Z-pinch'    ! geometry model (Z-pinch,s-alpha,miller,circular)
+  q0     = 0.0          ! safety factor
+  shear  = 0.0          ! magnetic shear
+  eps    = 0.0          ! inverse aspect ratio
+  kappa  = 1.0          ! elongation
+  s_kappa= 0.0          ! elongation derivative
+  delta  = 0.0          ! triangularity
+  s_delta= 0.0          ! triangularity derivative
+  zeta   = 0.0          ! squareness
+  s_zeta = 0.0          ! squareness derivative
+  parallel_bc = 'dirichlet' ! to change the type of parallel boundary condition (experimental)
+  shift_y= 0.0          ! to add a shift in the parallel BC (experimental)
+  Npol   = 1.0          ! set the length of the z domain (-pi N_pol < z < pi N_pol)
+/
+&DIAGNOSTICS
+  dtsave_0d = 0.5       ! period of 0D diagnostics (time traces)
+  dtsave_1d = -1        ! period of 1D diagnostics (nothing)
+  dtsave_2d = -1        ! period of 2D diagnostics (nothing)
+  dtsave_3d = 1         ! period of 3D diagnostics (phi, Aparallel, fluid moments etc.)
+  dtsave_5d = 10        ! period of 5D diagnostics (full set of GMs)
+/
+&MODEL
+  LINEARITY = 'nonlinear' ! set if we solve the linear or nonlinear problem
+  Na      = 2           ! number of species (this sets the number of species namelists to be read)
+  mu_x    = 1.0         ! numerical diffusion parameter in the radial direction
+  mu_y    = 1.0         ! numerical diffusion parameter in the binormal direction
+  N_HD    = 4           ! degree of the numerical diffusion
+  mu_z    = 0.0         ! numerical diffusion parameter in the parallel direction (fourth order only)
+  HYP_V   = 'hypcoll'   ! numerical diffusion scheme in the velocity space (experimental)
+  mu_p    = 0.0         ! numerical diffusion parameter in the parallel velocity (experimental)
+  mu_j    = 0.0         ! numerical diffusion parameter in the magnetic moment (experimental)
+  nu      = 0.1         ! collision rate, ~0.5 GENE parameter (better to use it instead of num. diff. in velocities)
+  beta    = 0.0         ! plasma beta (not in percent)
+  ADIAB_E = .f.         ! Use an adiabatic electron model (required if Na = 1)
+  ADIAB_I = .f.         ! Use an adiabatic ion model (experimental)
+/
+&CLOSURE
+  hierarchy_closure='truncation' ! closure scheme
+  dmax = -1             ! set the maximal degree of moment to evolve (-1 evolves all th)
+  nonlinear_closure='anti_laguerre_aliasing' ! set the truncation of the Laguerre convolution (truncation,full_sum,anti_laguerre_aliasing)
+  nmax = -1             ! set the maximal degree in the truncation of the Laguerre convolution (experimental)
+/
+&SPECIES ! Defines the species a (out of Na species)
+ name_ = 'ions' ! name of the species a
+ tau_  = 1.0    ! temperature (Ta/Te)
+ sigma_= 1.0    ! sqrt mass ratio (sqrt(ma/mi))
+ q_    = 1.0    ! charge (qa/e)
+ k_N_  = 2.0    ! density background gradient, grad ln N (omn in GENE)
+ k_T_  = 0.4    ! temperature background gradient, grad ln T (omT in GENE)
+/
+&SPECIES ! defines electrons (if not, use ADIAB_E = true)
+ ! electrons
+ name_ = 'electrons'
+ tau_  = 1.0
+ sigma_= 0.023338
+ q_    =-1.0
+ k_N_  = 2.0
+ k_T_  = 0.4
+/
+&COLLISION
+  collision_model = 'DG'   !LB/DG/SG/PA/LD (Lenhard-Bernstein, Dougherty, Sugama, pitch angle, Landau)
+  GK_CO           = .t.    ! gyrokinetic version of the collision operator (or longwavelength only)
+/
+&INITIAL
+  INIT_OPT         = 'phi' ! initilization of the system ('phi' put noise in electrostatic, 'blob' is like 'ppj' in GENE)
+  init_noiselvl    = 0.005 ! amplitude of the noise init (for 'phi' init only)
+  iseed            = 42    ! random seed
+/
+&TIME_INTEGRATION
+  numerical_scheme = 'RK4' ! numerical scheme for time integration (RK2,3,4 etc.)
+/