dgamma_dpsi.f90

Go to the documentation of this file.

function dgamma_dpsi(flux) result(f_dgamma_dpsi)
!------------------------------------------------------------------
! the normalized profile of gradient of gamma versus flux
!------------------------------------------------------------------

  use itm_types
  use mod_profiles
  use helena_spline

  implicit none

  real(r8), intent(in) :: flux
  real(r8) :: f_dgamma_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_dgamma_dpsi = spwert(npts, flux_tmp, gam_spline, psivec, ablt, 0)

end function dgamma_dpsi