Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(abstract_eqfile_geo_type), | intent(in) | :: | self | |||
real, | intent(in) | :: | pbar |
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