Calculate rgrid for the case of eeq and deq: in this case, rgrid is the distance to the magnetic axis as a function of theta.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real, | intent(out), | dimension(-ntgrid:) | :: | rgrid | ||
| real, | intent(in) | :: | rp | |||
| real, | intent(in), | dimension(-ntgrid:) | :: | theta |
subroutine rtgrid(rgrid, rp, theta)
implicit none
real, dimension(-ntgrid:), intent (out) :: rgrid
real, intent (in) :: rp
real, dimension(-ntgrid:), intent (in) :: theta
integer :: i
do i = -nth, nth
rgrid(i) = geom%rfun(rp, theta(i))
end do
end subroutine rtgrid