jacobian_polar Function

public function jacobian_polar(equi, x, y, phi)

Jacobian of polar coordinate system To compute derivatives on rho, finite difference is used N.B. currently only axisymmetric (2D), proper inclusion of toroidal phi will change this function

Arguments

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

Equilibrium

real(kind=FP), intent(in) :: x

Cartesian coordinate x

real(kind=FP), intent(in) :: y

Cartesian coordinate y

real(kind=FP), intent(in) :: phi

Toroidal angle

Return Value real(kind=fp)