PENTRC Conversion (WIP)

benchmarks/DIIID_kinetic_example/dcon.toml

@@ -0,0 +1,45 @@
+[DCON_CONTROL]
+bal_flag = false              # Ideal MHD ballooning criterion for short wavelengths
+mat_flag = true               # Construct coefficient matrices for diagnostic purposes
+ode_flag = true              # Integrate ODE's for determining stability of internal long-wavelength mode (must be true for GPEC)
+vac_flag = true               # Compute plasma, vacuum, and total energies for free-boundary modes
+mer_flag = true               # Evaluate the Mercier criterian
+
+set_psilim_via_dmlim = true               # Safety factor (q) limit determined as q_ir+dmlim...
+dmlim = 0.2                   # See sas_flag
+psiedge = 1.00                # If less then psilim, calculates dW(psi)...
+qlow = 1.02                   # Integration initiated at q determined by min(q0, qlow)...
+qhigh = 1e3                   # Integration terminated at q limit determined by min(qa, qhigh)...
+sing_start = 0                # Start integration at the sing_start'th rational from the axis (psilow)
+
+nn = 1                        # Toroidal mode number
+delta_mlow = 8                # Expands lower bound of Fourier harmonics
+delta_mhigh = 8               # Expands upper bound of Fourier harmonics
+delta_mband = 0               # Integration keeps only this wide a band...
+mthvac = 512                  # Number of points used in splines over poloidal angle at plasma-vacuum interface.
+thmax0 = 1                    # Linear multiplier on the automatic choice of theta integration bounds
+
+kin_flag = true              # Kinetic EL equation (default: false)
+con_flag = true              # Continue integration through layers (default: false)
+kinfac1 = 1.0                 # Scale factor for energy contribution (default: 1.0)
+kinfac2 = 1.0                 # Scale factor for torque contribution (default: 1.0)
+kingridtype = 0               # Regular grid method (default: 0)
+passing_flag = true           # Includes passing particle effects (default: false)
+ktanh_flag = true             # Ignore kinetic effects in the core smoothly (default: false)
+ktc = 0.1                     # Parameter for ktanh_flag (default: 0.1)
+ktw = 50.0                    # Parameter for ktanh_flag (default: 50.0)
+ion_flag = true               # Include ion dW_k when kin_flag is true
+electron_flag = false         # Include electron dW_k when kin_flag is true
+
+tol_nr = 1e-6                 # Relative tolerance of dynamic integration steps away from rationals
+tol_r = 1e-7                  # Relative tolerance of dynamic integration steps near rationals
+crossover = 1e-2              # Fractional distance from rational q at which tol switches
+singfac_min = 1e-4            # Fractional distance from rational q at which ideal jump enforced
+ucrit = 1e4                   # Maximum fraction of solutions allowed before re-normalized
+
+use_classic_splines = false   # Use a classical cubic spline instead of tri-diagonal solution
+
+[DCON_OUTPUT]
+write_dcon_out = true         # Write dcon.out which contains equilibrium data and final eigenvalues/vectors
+write_euler_h5 = true         # Write euler.h5 (euler.bin Julia equivalent)
+write_eqdata_h5 = true        # Write eqdata.h5, which contains profile and stability data

benchmarks/DIIID_kinetic_example/equil.toml

@@ -0,0 +1,30 @@
+[EQUIL_CONTROL]
+eq_type = "efit"             # Type of the input 2D equilibrium file
+eq_filename = "g147131.02300_DIIID_KEFIT"  # path to equilibrium file
+
+jac_type = "hamada"          # Coordinate system (hamada, pest, boozer, equal_arc)
+power_bp = 0                 # del.B ~ B_p^power_bp * B^power_b / R^power_r
+power_b = 0                  # del.B ~ B_p^power_bp * B^power_b / R^power_r
+power_r = 0                  # del.B ~ B_p^power_bp * B^power_b / R^power_r
+
+grid_type = "ldp"            # Radial grid packing
+psilow = 1e-4                # Min psi (normalized)
+psihigh = 0.993              # Max psi (normalized)
+mpsi = 128                   # Radial grid intervals
+mtheta = 256                 # Poloidal grid intervals
+
+newq0 = 0                    # Override q(0)
+etol = 1e-7                  # Reconstruction tolerance
+use_classic_splines = false  # Use classical spline (vs. tri-diagonal)
+
+input_only = false           # Quit after input read
+
+[EQUIL_OUTPUT]
+gse_flag = false           # Output G-S equation accuracy diagnostics
+out_eq_1d = false          # ASCII output of 1D eq data
+bin_eq_1d = false          # Binary output of 1D eq data
+out_eq_2d = false          # ASCII output of 2D eq data
+bin_eq_2d = true           # Binary output of 2D eq data (used by GPEC)
+out_2d = false             # ASCII output of processed 2D data
+bin_2d = false             # Binary output of processed 2D data
+dump_flag = false          # Binary dump of equilibrium data