salpha_equilibrium_m Module

Module containing the implementation of an salpha like equilibrium type. N.b. this equilibrium is axisymmetric - phi, when required, is not used

Note: This equilibrium is currently implemented in an inconsitent way, i.e. the analytical q-profile does not match the calculated one. To be consistent the q-profile needs to stay as it is while the magnetic field has to be b_x = -rhosin(theta)/(xrhosqrt(1-rho2)) and b_y = rhocos(theta)/(xrhosqrt(1-rho**2)) A consistent equilibrium of this form with a smaller set of parameters is implemented in circular_toroidal_equilibrium_m.f90


Uses

  • module~~salpha_equilibrium_m~~UsesGraph module~salpha_equilibrium_m salpha_equilibrium_m module~comm_handling_m comm_handling_m module~salpha_equilibrium_m->module~comm_handling_m module~constants_m constants_m module~salpha_equilibrium_m->module~constants_m module~descriptors_m descriptors_m module~salpha_equilibrium_m->module~descriptors_m module~equilibrium_m equilibrium_m module~salpha_equilibrium_m->module~equilibrium_m module~error_handling_m error_handling_m module~salpha_equilibrium_m->module~error_handling_m module~params_equi_salpha_m params_equi_salpha_m module~salpha_equilibrium_m->module~params_equi_salpha_m module~precision_m precision_m module~salpha_equilibrium_m->module~precision_m module~screen_io_m screen_io_m module~salpha_equilibrium_m->module~screen_io_m module~status_codes_m status_codes_m module~salpha_equilibrium_m->module~status_codes_m mpi mpi module~comm_handling_m->mpi module~constants_m->module~precision_m module~descriptors_m->module~screen_io_m module~equilibrium_m->module~precision_m module~error_handling_m->module~comm_handling_m module~error_handling_m->module~precision_m module~error_handling_m->module~screen_io_m module~error_handling_m->module~status_codes_m module~error_handling_m->mpi netcdf netcdf module~error_handling_m->netcdf module~params_equi_salpha_m->module~error_handling_m module~params_equi_salpha_m->module~precision_m module~params_equi_salpha_m->module~screen_io_m module~params_equi_salpha_m->module~status_codes_m iso_fortran_env iso_fortran_env module~params_equi_salpha_m->iso_fortran_env iso_c_binding iso_c_binding module~precision_m->iso_c_binding module~precision_m->iso_fortran_env module~precision_m->mpi module~precision_m->netcdf module~screen_io_m->module~precision_m module~screen_io_m->iso_fortran_env module~screen_io_m->netcdf

Used by

  • module~~salpha_equilibrium_m~~UsedByGraph module~salpha_equilibrium_m salpha_equilibrium_m module~equilibrium_factory_m equilibrium_factory_m module~equilibrium_factory_m->module~salpha_equilibrium_m program~benchmark_helmholtz_solvers benchmark_helmholtz_solvers program~benchmark_helmholtz_solvers->module~equilibrium_factory_m program~diagnose_poincare diagnose_poincare program~diagnose_poincare->module~equilibrium_factory_m program~test_diffusion test_diffusion program~test_diffusion->module~equilibrium_factory_m

Derived Types

type, public, extends(equilibrium_t) ::  salpha_t

Type implementing an salpha like equilibrium. The equilibrium has concentric flux surfaces in a toroidal geometry. Shear as well as the q-factor can be specified. To get the correct poloidal flux in Wb = V*s the reference magnetic field strength has to be provided in Tesla and the reference length scale has to be provided in meters.

Components

Type Visibility Attributes Name Initial
logical, public :: initialized = .false.
real(kind=FP), public :: x0

Magnetic axis x = R/R0 (in normalised units)
real(kind=FP), public :: y0

Magnetic axis y = Z/R0 (in normalised units)
real(kind=FP), public :: phi0 = 0.0_FP

Magnetic axis phi
real(kind=FP), public :: xmin

Box limits
real(kind=FP), public :: xmax

Box limits
real(kind=FP), public :: ymin

Box limits
real(kind=FP), public :: ymax

Box limits
real(kind=FP), public :: rhomax

Global limits for rho (rho = normalised psi, n.b. there may also be region-specific limits defined in equi)
real(kind=FP), public :: rhomin

Global limits for rho (rho = normalised psi, n.b. there may also be region-specific limits defined in equi)
real(kind=FP), public :: rho_limiter

Radial position (start) of limiter

Type-Bound Procedures

procedure, public, pass(self) :: absb
procedure, public, pass(self) :: bpol
procedure, public :: init
procedure, public :: display
procedure, public :: debug
procedure, public :: is_axisymmetric
procedure, public :: rho
procedure, public :: psi
procedure, public :: bx
procedure, public :: by
procedure, public :: btor
procedure, public :: jacobian
procedure, public :: epol
procedure, public :: erad
procedure, public :: district
procedure, public :: in_vessel
procedure, public :: mag_axis_loc
procedure, public :: qval
procedure, public, nopass :: theta
procedure, public :: x_point_psi
procedure, public :: o_point_psi