il_is_trapped Function

public elemental function il_is_trapped(il) result(is_trapped)

Returns true if the passed pitch angle grid point is considered trapped.

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: il

Return Value logical


Contents

Source Code

il_is_trapped

Source Code

  elemental logical function il_is_trapped(il) result(is_trapped)
    implicit none
    integer, intent(in) :: il
    is_trapped = il > ng2 + 1
    ! If we want to improve the consistency of the trapped wfb
    ! treatment we can replace the above with
    ! is_trapped = il > ng2
    ! To be explicit this might look like
    ! if (trapped_wfb) then
    !     is_trapped = il > ng2
    ! else
    !     is_trapped = il > ng2 + 1
    ! end if
  end function il_is_trapped