 |
SPEC 3.20
Stepped Pressure Equilibrium Code
|
|
SPEC 3.20
Stepped Pressure Equilibrium Code
|
| Diagnostics to check the code | |
| Free-Boundary Computation | |
| Parallelization | |
| Geometry | |
| Plasma Currents | |
| "global" force | |
| ▼Input namelists and global variables | |
| physicslist | The namelist physicslist controls the geometry, profiles, and numerical resolution |
| numericlist | The namelist numericlist controls internal resolution parameters that the user rarely needs to consider |
| locallist | The namelist locallist controls the construction of the Beltrami fields in each volume |
| globallist | The namelist globallist controls the search for global force-balance |
| diagnosticslist | The namelist diagnosticslist controls post-processor diagnostics, such as Poincaré plot resolution, etc |
| screenlist | The namelist screenlist controls screen output. Every subroutine, e.g. xy00aa.h, has its own write flag, Wxy00aa |
| "local" force | |
| Integrals | |
| Solver/Driver | |
| Build matrices | |
| Metric quantities | |
| Solver for Beltrami (linear) system | |
| Force-driver | |
| "packing" of Beltrami field solution vector | |
| Conjugate-Gradient method | |
| Initialization of the code | |
| Output file(s) | |
| Coordinate axis | |
| Rotational Transform | |
| Plasma volume | |
| Smooth boundary | |
| Enhanced resolution for metric elements | Enhanced resolution is required for the metric elements, \(g_{ij}/\sqrt g\), which is given by mne, ime, and ine. The Fourier resolution here is determined by lMpol=2*Mpol and lNtor=2*Ntor |
| Enhanced resolution for transformation to straight-field line angle | Enhanced resolution is required for the transformation to straight-field line angle on the interfaces, which is given by mns, ims and ins. The Fourier resolution here is determined by iMpol and iNtor |
| ▼Internal Variables | |
| Fourier representation | |
| Interface geometry: iRbc, iZbs etc. | The Fourier harmonics of the interfaces are contained in iRbc(1:mn,0:Mvol) and iZbs(1:mn,0:Mvol), where iRbc(l,j), iZbs(l,j) contains the Fourier harmonics, \(R_j\), \(Z_j\), of the \(l\)-th interface |
| Fourier Transforms | The coordinate geometry and fields are mapped to/from Fourier space and real space using FFTW3. The resolution of the real space grid is given by Nt=Ndiscrete*4*Mpol and Nz=Ndiscrete*4*Ntor |
| Volume-integrated Chebyshev-metrics | These are allocated in dforce(), defined in ma00aa(), and are used in matrix() to construct the matrices |
| Vector potential and the Beltrami linear system | |
| Field matrices: dMA, dMB, dMC, dMD, dME, dMF | |
| Construction of "force" | The force vector is comprised of Bomn and Iomn |
| Covariant field on interfaces: Btemn, Bzemn, Btomn, Bzomn | The covariant field |
| covariant field for Hessian computation: Bloweremn, Bloweromn | |
| Geometrical degrees-of-freedom: LGdof, NGdof | The geometrical degrees-of-freedom |
| Parallel construction of derivative matrix | |
| Derivatives of multiplier and poloidal flux with respect to geometry: dmupfdx | |
| Trigonometric factors | |
| Volume integrals: lBBintegral, lABintegral | |
| Internal global variables | Internal global variables; internal logical variables; default values are provided here; these may be changed according to input values |
| Miscellaneous | The following are miscellaneous flags required for the virtual casing field, external (vacuum) field integration, .. |
1.9.2