dp_dpsi.f90

Go to the documentation of this file.

function dp_dpsi(flux) result(f_dp_dpsi)
!------------------------------------------------------------------
! the normalized profile of the pressure gradient versus flux
! this routine must be initialized by a call to pressure_profile
!------------------------------------------------------------------

  use itm_types
  use mod_profiles
  use helena_spline

  implicit none

  real(r8), intent(in) :: flux
  real(r8) :: f_dp_dpsi

  real(r8), dimension(3) :: ablt   
  real(r8) :: flux_tmp

  flux_tmp = flux

  if (flux >= 1._r8) then
    flux_tmp = 1._r8
  end if

  f_dp_dpsi = spwert(npts, flux_tmp, dp_spline, psivec, ablt, 0)

end function dp_dpsi