safety_factor_m Module

Computes safety factor N.B. This module is currently axisymmetric - can only calculate pitch and safety factor for 2D equilibria PHITEMP parameter needed for equilibrium interface


Uses


Functions

public function safety_factor(equi, rho)

Compute safety factor q(rho)

Arguments

Type IntentOptional Attributes Name
class(equilibrium_t), intent(inout), target :: equi

Equilibrium defining the field line
real(kind=FP), intent(in) :: rho

Flux surface label

Return Value real(kind=fp)


Subroutines

public subroutine solout_none(me, nr, xold, x, y, irtrn, xout)

An empty function needed in order to run the field line tracer

Arguments

Type IntentOptional Attributes Name
class(dop853_pitch_t), intent(inout) :: me

Instance of the type
integer, intent(in) :: nr

Grid point (0,1,...)
real(kind=DP), intent(in) :: xold

The preceding grid point
real(kind=DP), intent(in) :: x

Current grid point
real(kind=DP), intent(in), dimension(:) :: y

State vector [size n]
integer, intent(inout) :: irtrn

Serves to interrupt the integration, if irtrn is set <0
real(kind=DP), intent(out) :: xout

Can be used for efficient intermediate output