phys_profiles.f90

Go to the documentation of this file.

module phys_profiles

  use itm_types
  use mod_dat, only : verbosity
  use mod_helena_io, only : iu6

  implicit none

  contains

!--------------------------------------------------------------------------

  function j_phi(flux) result(f_j_phi)
!--------------------------------------------------------------------------
! current density versus flux (has to be normalized to 1 on axis)
!--------------------------------------------------------------------------

    use mod_profiles
    use helena_spline

    implicit none

    real(r8) :: flux
    real(r8) :: f_j_phi

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

    if (flux > 1._r8) then
      if (verbosity > 0) &
       write(iu6, *) 'Warning: normalized flux > 1 ! Should not occur.'
      flux = 1._r8
    end if

    f_j_phi = spwert(npts, flux, j_spline, psivec, ablt, 0)

  end function j_phi

!===========================================================================

end module phys_profiles