Calculates rmajor = 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)
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real, | intent(in), | dimension(-ntgrid:) | :: | rgrid | ||
| real, | intent(in), | dimension(-ntgrid:) | :: | theta | ||
| real, | intent(out), | dimension(-ntgrid:) | :: | rmajor |
subroutine rmajortgrid(rgrid, theta, rmajor)
implicit none
real, dimension(-ntgrid:), intent(in) :: rgrid, theta
real, dimension(-ntgrid:), intent(out) :: rmajor
integer :: i
do i = -nth, nth
rmajor(i) = geom%Rpos(rgrid(i), theta(i))
end do
end subroutine rmajortgrid