dpfun_eqfile Function

public function dpfun_eqfile(self, pbar) result(dpfun)

Uses

Type Bound

abstract_eqfile_geo_type

Arguments

Type IntentOptional Attributes Name
class(abstract_eqfile_geo_type), intent(in) :: self
real, intent(in) :: pbar

Return Value real


Contents

Source Code


Source Code

  real function dpfun_eqfile (self, pbar) result(dpfun)
    use splines, only: dsplint
    implicit none
    class(abstract_eqfile_geo_type), intent(in) :: self
    real, intent(in) :: pbar
    ! p_N would be B**2/mu_0 => p = beta/2 in our units
    dpfun = 0.5 * self%beta_0 * dsplint(pbar, self%pressure_spl) / (self%psi_a - self%psi_0)
  end function dpfun_eqfile