chease_chi_index Function

public elemental function chease_chi_index(nchi, itheta)

Convert from theta-index which is -pi to pi, to chi-index, which is 0 to 2*pi. Assumes nchi is even?

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nchi
integer, intent(in) :: itheta

Return Value integer


Contents

Source Code

chease_chi_index

Source Code

  elemental integer function chease_chi_index(nchi, itheta)
    integer, intent(in) :: nchi, itheta
    ! when ichi = 0, itheta = (ntheta-1)/2 + 1
    ! when itheta = 0, ichi = ntheta / 2
    ! itheta  1 2 3 4 5 6 7 8
    ! ichi    5 6 7 8 1 2 3 4
    if (itheta > nchi/2) then
       chease_chi_index = itheta - nchi/2
    else
       chease_chi_index = itheta + nchi/2
    end if
  end function chease_chi_index