invRtgrid Subroutine

private subroutine invRtgrid(rgrid, theta, invR)

Calculates invR = 1/R(theta), where R is the distance to the central axis (rmajor not to be confused with rmaj, r_geo etc) rgrid is the submodule radial grid (r_circ for EFIT & deq, rp otherwise)

Arguments

Type IntentOptional Attributes Name
real, intent(in), dimension(-ntgrid:) :: rgrid
real, intent(in), dimension(-ntgrid:) :: theta
real, intent(out), dimension(-ntgrid:) :: invR

Contents

Source Code

invRtgrid

Source Code

  subroutine invRtgrid(rgrid, theta, invR)
    implicit none
    real, dimension(-ntgrid:), intent(in) :: rgrid, theta
    real, dimension(-ntgrid:), intent(out) :: invR
    integer :: i
    do i = -nth, nth
       invR(i) = geom%invR(rgrid(i), theta(i))
    end do
  end subroutine invRtgrid