q_profile.f90

Go to the documentation of this file.

function q_profile(flux) result(f_q_profile)
!------------------------------------------------------------------
! the normalized profile of the current density versus flux
!------------------------------------------------------------------
  
  use itm_types
  use mod_profiles
  
  implicit none
  
  real(r8) :: flux
  real(r8) :: f_q_profile
  
  integer(itm_i4) :: nint
  
  real(r8) :: dps
  
  dps = 1._r8 / dble(npts - 1)
  nint = int((npts - 1) * sqrt(flux)) + 1
  if (flux >= 1._r8) nint = npts - 1
  f_q_profile = q_in(nint) + (sqrt(flux) - dps * (nint - 1)) / dps &
   * (q_in(nint + 1) - q_in(nint))
  
end function q_profile