Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(abstract_eqfile_geo_type), | intent(in) | :: | self | |||
real, | intent(in) | :: | pbar |
real function pfun_eqfile (self, pbar) result(pfun)
use splines, only: splint
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
pfun = 0.5 * self%beta_0 * splint(pbar, self%pressure_spl)
end function pfun_eqfile