From 2e3952b0e2382fefa8359a94181609a456d1dedc Mon Sep 17 00:00:00 2001 From: Antoine Hoffmann <antoine.hoffmann@epfl.ch> Date: Wed, 18 Jan 2023 10:54:28 +0100 Subject: [PATCH] Comment all input parameters in fort_example.90 --- fort_00.90 | 85 --------------------------------------------- fort_example.90 | 91 +++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 91 insertions(+), 85 deletions(-) delete mode 100644 fort_00.90 create mode 100644 fort_example.90 diff --git a/fort_00.90 b/fort_00.90 deleted file mode 100644 index ae9b19ba..00000000 --- a/fort_00.90 +++ /dev/null @@ -1,85 +0,0 @@ -&BASIC - nrun = 100000000 - dt = 0.01 - tmax = 500 - maxruntime = 60 ! 1h 14400 !4h -/ -&GRID - pmaxe = 6 - jmaxe = 3 - pmaxi = 6 - jmaxi = 3 - Nx = 200 - Lx = 120 - Ny = 32 - Ly = 60 - Nz = 1 - SG = .f. -/ -&GEOMETRY - geom = 'Z-pinch' - q0 = 0 - shear = 0 - eps = 0 -/ -&OUTPUT_PAR - nsave_0d = 10 - nsave_1d = -1 - nsave_2d = -1 - nsave_3d = 100 - nsave_5d = 10000 - write_doubleprecision = .t. - write_gamma = .t. - write_hf = .t. - write_phi = .t. - write_Na00 = .t. - write_Napj = .t. - write_Sapj = .f. - write_dens = .t. - write_temp = .t. - job2load = -1 -/ -&MODEL_PAR - ! Collisionality - CLOS = 0 - NL_CLOS = 0 - LINEARITY = 'nonlinear' - KIN_E = .t. - mu_x = 1.0 - mu_y = 1.0 - N_HD = 4 - mu_z = 0.0 - mu_p = 0 - mu_j = 0 - nu = 0.01 - tau_e = 1 - tau_i = 1 - sigma_e = 0.023338 - sigma_i = 1 - q_e = -1 - q_i = 1 - K_Ne = 1.7 - K_Ni = 1.7 - K_Te = 0.425 - K_Ti = 0.425 - GradB = 1 - CurvB = 1 - lambdaD = 0 -/ -&COLLISION_PAR - collision_model = 'LD' - gyrokin_CO = .true. - interspecies = .true. - mat_file = 'LDGK_P10_J5_dk_5e-2_km_5_NFLR_30.h5' - collision_kcut = 1.8 -/ -&INITIAL_CON - INIT_OPT = 'phi' - ACT_ON_MODES = 'donothing' - init_background = 0 - init_noiselvl = 0.001 - iseed = 42 -/ -&TIME_INTEGRATION_PAR - numerical_scheme = 'RK4' -/ diff --git a/fort_example.90 b/fort_example.90 new file mode 100644 index 00000000..a681f044 --- /dev/null +++ b/fort_example.90 @@ -0,0 +1,91 @@ +&BASIC + nrun = 100000000 !number of maximal time steps + dt = 0.01 !time step + tmax = 500 !maximal physical time + maxruntime = 60 ! 1h 14400 !4h !maximal run time +/ +&GRID + pmaxe = 6 !maximal degree of Hermite polynomials for e + jmaxe = 3 !maximal degree of Laguerre polynomials for e + pmaxi = 6 !maximal degree of Hermite polynomials for i + jmaxi = 3 !maximal degree of Laguerre polynomials for i + Nx = 200 !resolution in x (=Nkx) + Lx = 120 !box size in x + Ny = 64 !resolution in y (=2(Nky-1)) + Ly = 160 !box size in y + Nz = 24 !resolution in z + Npol = 1 !number of poloidal turns (Lz=2piNpol) + Nexc = 1 !factor to increase Lx in sheared geometry + SG = .f. !staggered grid option (not recommended) +/ +&GEOMETRY + geom = 's-alpha' !magnetic equilibrium geometry (Z-pinch,s-alpha,miller) + q0 = 1.4 !safety factor (s-alpha,miller only) + shear = 0.8 !shear (s-alpha,miller only) + eps = 0.18 !aspect ratio (s-alpha,miller only) + kappa = 1 !elongation (miller only) + delta = 0 !triangularity (miller only) + zeta = 0 !squareness (miller only) + parallel_bc = 'dirichlet' !boundary condition for modes that does not connect due to shear (dirichlet,periodic) +/ +&OUTPUT_PAR + nsave_0d = 50 !period in number of step for time traces + nsave_1d = -1 !unused + nsave_2d = -1 !unused + nsave_3d = 100 !period in number of step for 3D fields (phi,psi,...) + nsave_5d = 1000 !period in number of step for 5D fields (moments) + write_doubleprecision = .t. !for HDF5 output (double precision faster on marconi) + write_gamma = .t. !to write particle flux + write_hf = .t. !to write heat flux + write_phi = .t. !to write ES and EM potentials + write_Na00 = .t. !to write gyrocenter densities + write_Napj = .t. !to write moments + write_Sapj = .f. !to write nonlinear terms + write_dens = .t. !to write particle densities + write_temp = .t. !to write particle temperatures + job2load = -1 !ID of the job to load in a restart (-1 means no restart) +/ +&MODEL_PAR + ! Collisionality + CLOS = 0 !closure model (0: zero-truncation) + NL_CLOS = 0 !NL closure model (-1 is full FLR sum until n=J-j, n>-1 is up to nth term) + LINEARITY = 'nonlinear' !to activate nonlinear term (linear,nonlinear) + KIN_E = .f. !to have a kinetic electron model (adiabatic otherwise) + mu_x = 1.0 !x numerical diffusion coefficient + mu_y = 1.0 !y numerical diffusion coefficient + N_HD = 4 !xy numerical diffusion order + mu_z = 2.0 !z numerical diffusion coefficient (order 4) + mu_p = 0 !p numerical diffusion coefficient (usually not used) + mu_j = 0 !j numerical diffusion coefficient (usually not used) + nu = 0.05 !collision frequency (=0.49*nu_GENE) + tau_e = 1 !electron temperature ratio + tau_i = 1 !ion temperature ratio + sigma_e = 0.023338 !electron mass ratio + sigma_i = 1 !ion mass ratio + q_e = -1 !electron charge + q_i = 1 !ion charge + K_Ne = 2.22 !electron density gradient intensity + K_Te = 6.96 !electron temperature gradient intensity + K_Ni = 2.22 !ion density gradient intensity + K_Ti = 6.96 !ion density temperature intensity + GradB = 1 !magnetic field gradient strength + CurvB = 1 !magnetic curvature strength + lambdaD = 0 !Debye length (not tested when non zero) +/ +&COLLISION_PAR + collision_model = 'DG' !collision model (DG,SG,LD,LR), all need a matrix except DG + gyrokin_CO = .f. !activate gyrokinetic terms in the CO + interspecies = .t. !activate interspecies if CO has some + mat_file = 'LDGK_P10_J5_dk_5e-2_km_5_NFLR_30.h5' !path to find the collision matrix + collision_kcut = 1.8 !maximal wavelength of the CO matrix. For higher kperp, the matrix at collision_kcut will be applied. +/ +&INITIAL_CON + INIT_OPT = 'phi' !initialization option (phi,mom00,allmom,ppj) + ACT_ON_MODES = 'donothing' !to perform numerical experiments (unused) + init_background = 0 !background value for a noise initialization + init_noiselvl = 0.0001 !fluctuation value '' + iseed = 42 !seed of the noise +/ +&TIME_INTEGRATION_PAR + numerical_scheme = 'RK4' !numerical scheme for time-stepping (RK2,RK3,RK4,DOPRI5) +/ -- GitLab